• [SDOI2008]仪仗队


    传送门:

    分析完题目,不就是要找y=kx有多少条0.0,那就是要找 (n-1)内gcd(i,j)==1的个数吗,最后答案+2就行了,直接走莫比乌斯反演。。。(题解区并没有人用)

    最后特判n==1输出0就行了

    #include <cstdio>
    #include <cstring>
    #include <string>
    #include <cstdlib>
    #include <algorithm>
    #include <iostream>
    using namespace std;
    #define ll long long
    #define re register
    const int N=4e4+10;
    const int M=4e4;
    inline void read(int &a)
    {
        a=0;
        int d=1;
        char ch;
        while(ch=getchar(),ch>'9'||ch<'0')
            if(ch=='-')
                d=-1;
        a=ch^48;
        while(ch=getchar(),ch>='0'&&ch<='9')
            a=(a<<3)+(a<<1)+(ch^48);
        a*=d;
    }
    int mu[N],pri[N],cnt;
    bool vis[N];
    inline void init()
    {
        mu[1]=1;
        for(re int i=2;i<=M;i++)
        {
            if(!vis[i])
                pri[++cnt]=i,mu[i]=-1;
            for(re int j=1;j<=cnt&&i*pri[j]<=M;j++)
            {
                vis[i*pri[j]]=1;
                if(i%pri[j]==0)
                    break;
                mu[i*pri[j]]=-mu[i];
            }
        }
        for(re int i=1;i<=M;i++)
            mu[i]+=mu[i-1];
    }
    inline ll solve(int n)
    {
        ll ans=0;
        for(int l=1,r;l<=n;l=r+1)
        {
            r=n/(n/l);
            ans+=(mu[r]-mu[l-1])*(n/l)*(n/l);
        }
        return ans;
    }
    int main()
    {
        init();
        int n;
        read(n);
        if(n==1)
            return puts("0"),0;
        n--;
        printf("%lld",solve(n)+2);
        return 0;
    }
  • 相关阅读:
    NHibernate中session.update()及session.merge()的区别
    子序列 (All in All,UVa 10340)
    古老的密码 (Ancient Cipher,NEERC 2004,LA 3213)
    例题1 勇者斗恶龙 (The Dragon of Loowater,UVa 11292)
    HDU1869 六度分离
    B. T-primes
    PoJ 1595 PrimeCuts
    poj 3518 Prime Gap
    PKU1988磁铁
    求组合数
  • 原文地址:https://www.cnblogs.com/acm1ruoji/p/10949182.html
Copyright © 2020-2023  润新知