2016.7.9 的数学题！(@﹏@)~

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

想的是自己写来着 实在是 改不出来了。

分析：
枚举公约数k，求gcd（m,n)为 k的 m个数。
即gcd（m,n)=k;
那么(m/k,n/k)=1对不对！【【对~
并且m/k,n/k互质，这个时候我们想到了什么！没错就是欧拉函数【【要灵活运用嘞，队长说这玩意 noip之后挺有用 的？

所以求一下ψ(x)喽
这个等于 x*(1-1/p1)(1-1/p2)…..
p1…是 x的素因数！！
具体怎么求 看代码喽。另外 %%%%%黄学长

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define ll long long
using namespace std;
ll n;
ll ans;
ll phi(ll x)
{
    ll sum;
    sum=x;
    //int tp=x;
    int t=sqrt(x);
    for(int i=2;i<=t;i++)
    {
        if(x%i == 0)   
        {
            sum=(ll)sum*(i-1)/i;
            while(x%i == 0) x=x/i;    //举个例子 12  就可以把2*2 都删掉 就保证 了剩下的 可以整除的 
        }
    }
    if(x>1) sum=(ll)sum*(x-1)/x; 
    return sum;
}
int main()
{
    scanf("%lld",&n);
    int t=sqrt(n);
    for(int i=1;i<=t;i++)
    {
        if(n % i == 0)
        {
            ans+=(ll)(i*phi(n/i));
            if(i*i<n)ans+=(ll)(n/i)*phi(i);    //i*(n/i)=i对不对！  因为 公约数吗，所以sqrt（n)左右两边可能都有。 
        }
    }
    printf("%lld\n",ans);
}

先回家喽，等会再写题！