BZOJ 2705 [SDOI2012]Longge的问题

来源:互联网 发布:淘宝评价修改链接 编辑:程序博客网 时间:2024/06/11 16:32

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


【题目分析】
因为数据范围过大。不能存下来,只能使用根号n的算法每次都要算。


【代码】

#include <cstdio>#include <cmath>using namespace std;long long m,t,n,ans;long long phi(long long x){    t=x;    for (long long i=2;i<=m;++i)    if (x%i==0){t=t/i*(i-1);while (x%i==0) x/=i;}    if (x>1) t=t/x*(x-1);    return t;}int main(){    scanf("%lld",&n);    m=sqrt(n);    for (long long i=1;i<=m;++i)    if (n%i==0)    {        ans+=i*phi(n/i);        if (i*i<n) ans+=(n/i)*phi(i);    }    printf("%lld",ans);}
0 0