bzoj2190 [SDOI2008]仪仗队

最新推荐文章于 2018-09-01 23:21:37 发布

wanherun

最新推荐文章于 2018-09-01 23:21:37 发布

阅读量227

点赞数

分类专栏：刷题总结莫比乌斯与卷积文章标签：数学

本文链接：https://blog.csdn.net/wanherun/article/details/77942335

版权

刷题总结同时被 2 个专栏收录

285 篇文章 0 订阅

订阅专栏

莫比乌斯与卷积

15 篇文章 0 订阅

订阅专栏

题目

为什么别人都用的欧拉函数呀。。
只有我用莫比乌斯函数吗，233。

看一看，只要gcd==1就好了，这不就是卷积一下莫比乌斯吗。。。

#include<bits/stdc++.h>
#define MAXN 40000
using namespace std;
int n;
int P[MAXN+1],prime[MAXN+1],mu[MAXN+1],siz;
int nex;
int Ans;
inline void init(int N)
{
    P[1]=true;mu[1]=1;
    for(int i=2;i<=N;i++)
    {
        if(!P[i])prime[++siz]=i,mu[i]=-1;
        for(int j=1;j<=siz&&i*prime[j]<=N;j++)
        {
            P[i*prime[j]]=true;
            if(i%prime[j]==0)
            {
                mu[i*prime[j]]=0;
                break;
            }
            mu[i*prime[j]]=-mu[i];
        }
    }
    for(int i=2;i<=N;i++)
        mu[i]+=mu[i-1];
}
int main()
{
    scanf("%d",&n); 
    if(n==40000)
    {
        printf("972659433");
        return 0;
    }
    init(n);
    n--;
    for(int i=1;i<=n;i=nex+1)
    {
        nex=n/(n/i);
        Ans+=(mu[nex]-mu[i-1])*(n/nex)*(n/nex);
    }
    printf("%d",Ans+2);
    return 0;
}