bzoj2190: [SDOI2008]仪仗队

最新推荐文章于 2018-03-07 14:01:00 发布

mengbi_er

最新推荐文章于 2018-03-07 14:01:00 发布

阅读量245

点赞数

分类专栏： bzoj ————数学欧拉函数

本文链接：https://blog.csdn.net/mengbi_er/article/details/75136499

版权

bzoj 同时被 3 个专栏收录

107 篇文章 1 订阅

订阅专栏

————数学

17 篇文章 0 订阅

订阅专栏

欧拉函数

4 篇文章 0 订阅

订阅专栏

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int tot,prime[20005],fai[40005];bool check[40005];
int main()
{
    int n;
    scanf("%d",&n);
    unsigned long long ans=0;
    for(int i=2;i<=n;i++)
    {
        if(!check[i])
        {
            prime[tot++]=i;
            fai[i]=i-1;
        }
        for(int j=0;j<tot;j++)
        {
            if(i*prime[j]>n)
            break;
            check[i*prime[j]]=1;
            if(i%prime[j])
            {
                fai[i*prime[j]]=fai[i]*(prime[j]-1);
            }
            else
            {
                fai[i*prime[j]]=fai[i]*prime[j];
                break;
            }
        }
    }
    for(int i=1;i<=n-1;i++)
    ans+=fai[i];
    if(n==1)
    cout<<"0";
    else if(n==2)
    cout<<"3";
    else
    cout<<(ans+2<<1)-1;
}