码迷,mamicode.com
首页 > 其他好文 > 详细

bzoj2705: [SDOI2012]Longge的问题

时间:2017-10-08 11:21:31      阅读:155      评论:0      收藏:0      [点我收藏+]

标签:它的   mil   string   logs   mem   int   nbsp   linear   bool   

嗯A了道水题(其实在学polya的时候做过类似的)

题意很裸,就是求sigma(gcd(i,n)),那很容易发现很多i和n的gcd是相等的,那我们就枚举gcd,然后将它的phi求出来,那(n/i)*phi(i)就将全部gcd为(n/i)的值给求出来了,同理phi(n/i)也一样,那先预处理一下素数,然后O(sqrt(n))枚举就行了。

#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
typedef long long LL;
bool v[1100000];
int pr,prime[1100000];
void linear_prime()
{
    memset(v,true,sizeof(v));
    for(int i=2;i<=1000000;i++)
    {
        if(v[i]==true)prime[++pr]=i;
        for(int j=1;j<=pr&&i*prime[j]<=1000000;j++)
        {
            v[i*prime[j]]=false;
            if(i%prime[j]==0)break;
        }
    }
}
LL n;
LL phi(LL x)
{
    LL ans=x;
    for(int i=1;prime[i]*prime[i]<=x;i++)
    {
        if(x%prime[i]==0)
        {
            ans=ans-ans/prime[i];
            while(x%prime[i]==0)x/=prime[i];
        }
    }
    if(x!=1)ans=ans-ans/x;
    return ans;
}
int main()
{
    linear_prime();
    scanf("%lld",&n);
    LL ans=0;
    for(int i=1;i*i<=n;i++)//枚举gcd 
    {
        if(n%i==0)
        {
            ans+=phi(i)*(n/i);
            if(i*i!=n)ans+=phi(n/i)*i;
        }
    }
    printf("%lld\n",ans);
    return 0;
}

 

bzoj2705: [SDOI2012]Longge的问题

标签:它的   mil   string   logs   mem   int   nbsp   linear   bool   

原文地址:http://www.cnblogs.com/AKCqhzdy/p/7636681.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!