hdu 4059 The Boss on Mars（容斥原理）

给你一个数n，求1-n之间所有与n互质的数字的四次方的和。求小于n且与n互质的数字以及求这些数字的和比较简单，但是这里多了个四次方，求四次方和，这就第一次碰到了。不过网上找了下，有求四次方和的公式。

有了这个那就好求了。这题最大的难处是取模，wa了一大片，一定要注意模运算，一不小心就溢出了。

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
const LL mod = 1e9+7;
LL _30rev;
LL mPow(LL a, LL b)
{
    a %= mod;
    LL ret = 1;
    while(b)
    {
        if(b&1) ret = ret*a%mod;
        a = a*a%mod;
        b >>= 1;
    }
    return ret%mod;
}
LL p[65];
LL plen;

LL calc(LL num)
{
    return num*(num+1)%mod*(2*num+1)%mod*(3*num*num%mod+3*num-1)%mod*_30rev%mod;
}

LL solve(LL n)
{
    plen = 0;
    LL temp = n;
    for(LL i = 2; i <= temp/i; ++i)
    {
        if(temp%i == 0)
        {
            p[plen++] = i;
            while(temp%i == 0) temp /= i;
        }
    }
    if(temp > 1) p[plen++] = temp;
    LL res = calc(n);
    for(LL i = 1; i < (1LL<<plen); ++i)
    {
        LL mult = 1;
        LL cnt = 0;
        for(LL j = 0; j < plen; ++j)
        {
            if(i&(1<<j))
            {
                mult *= p[j];
                ++cnt;
            }
        }
        //a*b=n   a^2*b^2=n^2    3+6+9=3*(1+2+3)
        LL num = (mult%mod*mult%mod)%mod*(mult%mod*mult%mod)%mod*calc(n/mult)%mod;
        if(cnt&1) res = ((res - num)%mod+mod)%mod;
        else res = ((res + num)%mod+mod)%mod;
    }
    return res;
}

int main()
{
    int T;
    LL n;
    _30rev = mPow(30,mod-2);
    scanf("%d",&T);
    while(T--)
    {
        scanf("%I64d",&n);
        printf("%I64d\n",solve(n));
    }
    return 0;
}