1個數N(N <= 10^9) //實際上應該是n<=10^18

公約數之和

6

15


尤拉函式  
對於 樣例  
1 2 3 2 1 ----6 
有2個1 2個2 和 1個3 
兩個1的情況 1,5與6互質 
兩個2的情況 1(2/2),2(4/2) 與3互質(6/2) 
一個三的情況 1(3/3) 與2(6/3) 互質 
列舉n的所有因子
所以 對於某個因子 i
ans+=i*phi(n/i);

 1 #include <cmath>
 2 #include <cstdio>
 3 #include <cctype>
 4 
 5
 
 typedef long long LL;
 6 
 7 LL n;
 8 
 9 LL Get_Phi(LL x) {
10     LL sum=x;
11     for(LL i=2;i*i<=x;++i) {
12         if(x%i==0) {
13             sum-=sum/i;
14             while(x%i==0) x/=i;
15         }
16     }
17     if(x>1) sum-=sum/x;
18     return sum;
19 }
20 
21 int hh() {
22     scanf("
 
%lld",&n);
23     LL ans=0;
24     LL sq=sqrt(n);
25     for(LL i=1;i<=sq;++i) {
26         if(n%i==0) {
27             LL s=Get_Phi(n/i);
28             ans+=(LL)s*i;
29             LL p=Get_Phi(i);
30             ans+=(LL)p*(n/i);
31         }
32     }
33     if(sq*sq==n) ans-=sq*Get_Phi(sq);
34
 
     printf("%lld\n",ans);
35     return 0;
36 }
37 
38 int sb=hh();
39 int main(int argc,char**argv) {;}