hdu 4059

来源：互联网发布：淘宝怎么修改差评编辑：程序博客网时间：2024/04/29 12:46

http://acm.hdu.edu.cn/showproblem.php?pid=4059

题意很好理解，熔池原理，还有公式 http://math2.org/math/expansion/power.htm

 #include<cstdio> #include<cstring> #include <iostream>using namespace std;const __int64 iv30 = 233333335;const __int64 mod=1000000007;__int64 n;__int64 cal(__int64 x) {     __int64 y=n/x;     return (x*x)%mod*x%mod*x%mod*((6*y%mod*y%mod*y%mod*y%mod*y%mod+(15*y%mod*y%mod*y%mod*y)%mod)%mod+10*y%mod*y%mod*y%mod-y)%mod*iv30%mod;}int main(){     int ca,p[100];     scanf("%d",&ca);     while(ca--)     {         scanf("%I64d",&n);         int cnt=0;         __int64 tmp=n;         for(int i=2;i*i<=tmp;i++)          if(tmp%i==0)          {              p[cnt++]=i;              while(tmp%i==0)   tmp/=i;          }         if(tmp!=1) p[cnt++]=tmp;         __int64 ans=0;         for(__int64 i=1;i<(1<<cnt);i++)         {              int num=0;                  __int64 mul=1;              for(__int64 j=0;j<cnt;j++)                         if(i&(1<<j)) num++,mul=mul*p[j]%mod;              if(num&1)                  ans=(ans+cal(mul)+mod)%mod;              else                ans=(ans-cal(mul)+mod)%mod;         }         ans=(cal(1)-ans+mod)%mod;         printf("%I64d\n",ans);     }      return 0; }