bzoj2705 欧拉函数

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=2705

2705: [SDOI2012]Longge的问题

Time Limit: 3 Sec  Memory Limit: 128 MB
Submit: 3377  Solved: 2099
[Submit][Status][Discuss]

Description

Longge的数学成绩非常好，并且他非常乐于挑战高难度的数学问题。现在问题来了：给定一个整数N，你需要求出∑gcd(i, N)(1<=i <=N)。

Input

一个整数，为N。

Output

一个整数，为所求的答案。

Sample Input

6

Sample Output

15

HINT

【数据范围】

对于60%的数据，0<N<=2^16。

对于100%的数据，0<N<=2^32。

【分析】：有两种方法。

方法1：偏暴力。对与一个数n，我们可以只算与n不互质的数的gcd，最后再加上1的个数即可，n本身先不算在内。

所以统计n的因子的倍数，这些数都是与n不互质的。

比如n=12，因子有6 ，4 ，3 ，2 ，所以与n的gcd==6的数有12/6-1个，gcd==4的数有12/4-1个，

gcd==3的数，理应由12/3-1个，但是还要减掉gcd==6的数的个数！因为gcd（6,12）==6，覆盖了3。

以此类推。

方法2：欧拉函数->http://blog.csdn.net/winter2121/article/details/78411647

用phi（n）代表欧拉函数。

对于gcd（i，n）=k，若k为1，则符合条件的i有phi（i）个

若k不为1，那怎么统计此时i的个数呢？

gcd（i/k，n/k）=1即可。也就是这是的符合条件的i有phi（n/k）个

所以，统计了gcd为k时的i的数量res，乘上k，累加进答案即可。

【代码-方法1】

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
int c[101010];
int d[100010],top;
void divs(int n)
{
    top=0;
    for(int i=2;i*i<=n;i++)if(n%i==0)
    {
        d[top++]=i;
        if(i!=n/i)d[top++]=n/i;
    }
    sort(d,d+top);
}
int main()
{
    int n;
    while(cin>>n)
    {
        divs(n);
        memset(c,0,sizeof(c));
        for(int i=0;i<top;i++)
            c[i]=n/d[i]-1;//第i个因子的个数
        for(int i=top-1;i>=0;i--)
        {
            for(int j=i-1;j>=0;j--)
                if(d[i]%d[j]==0)
                    c[j]-=c[i];//有更大的公因子，所以j的要去掉
        }
        ll res=0;
        ll sum=0;//记录gcd不为1的数量
        for(int i=0;i<top;i++)
        {
            res+=c[i]*d[i];
            sum+=c[i];
        }
        res+=n-sum+n-1;
        printf("%lld\n",res);
    }
}

【代码-方法2】

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
int phi(int n)
{
    int ans=n;
    for(int i=2;i*i<=n;i++)if(n%i==0)
    {
        ans=ans/i*(i-1);
        while(n%i==0)n/=i;
    }
    if(n>1)ans=ans/n*(n-1);
    return ans;
}
int main()
{
    ll n,ans=0;
    cin>>n;
    for(ll k=1;k*k<=n;k++)if(n%k==0)//枚举因子k
    {
        ans+=k*phi(n/k);//k乘上gcd为k的数量
        if(n/k!=k)ans+=n/k*phi(k);
    }
    cout<<ans<<endl;
}