Updated：这只鸽子学会了LaTeX
题意
求
$$sum_{i=1}^{n} sum_{j=1}^{n} phi(gcd(phi(i),phi(j)))$$
在草稿纸上一番乱划推导之后
$$sum_{i=1}^{n} sum_{j=1}^{n} phi(gcd(phi(i),phi(j)))$$
$$=sum_{d=1}^{n} phi(d){ imes}f(d)$$
这里f(d)表示有多少对i，j满足他们的$phi$值的gcd为d
$$f(d)=sum_{i,j} {lbrack} {gcd(phi(i),phi(j))=d} { brack}$$
我们令
$$F(d)=sum_{i,j}{lbrack}{d{vert}gcd(phi(i),phi(j))}{ brack}$$
就有
$$F(n)=sum_{n{vert}d} f(d)$$
反演得到
$$f(n)=sum_{d} F(n{ imes}d){ imes}mu(d)$$
所以原式就等于
$$sum_{i{ imes}j<=n} phi(i){ imes}F(i{ imes}j){ imes}mu(j)$$
只需要预处理出$mu , phi $ 和F即可
code
1 // 2 // main.cpp 3 // 51nod 4 // 5 // Created by gengyf on 2019/8/8. 6 // Copyright © 2019 yifan Geng. All rights reserved. 7 // 8 9 #include <iostream> 10 using namespace std; 11 namespace gengyf{ 12 #define maxn 2000010 13 #define ll long long 14 inline int read(){ 15 int f=1,x=0;char s=getchar(); 16 while(s<'0'||s>'9'){if(s=='-')f=-1;s=getchar();} 17 while(s>='0'&&s<='9'){x=x*10+s-'0';s=getchar();} 18 return f*x; 19 } 20 ll w[maxn>>2],miu[maxn],phi[maxn],c[maxn],vis[maxn]; 21 int n,T,t;ll ans; 22 void init(){ 23 phi[1]=miu[1]=1; 24 for(int i=2;i<=maxn;i++){ 25 if(!vis[i]){//如果是质数 26 miu[++w[++t]=i]=-1; 27 phi[i]=i-1; 28 } 29 for(int j=1,k;(k=i*w[j])<=maxn;j++){ 30 vis[k]=1; 31 if(i%w[j]==0){ 32 phi[k]=phi[i]*w[j]; 33 miu[k]=0;break; 34 } 35 phi[k]=phi[i]*(w[j]-1);miu[k]=-miu[i]; 36 } 37 } 38 } 39 ll F(int x){return c[x]*c[x];} 40 ll calc(int n){ 41 memset(c,0,sizeof(c)); 42 ans=0; 43 for(int i=1;i<=n;i++){ 44 ++c[phi[i]]; 45 } 46 for(int i=1;i<=n;i++) 47 for(int j=i+i;j<=n;j+=i){ 48 c[i]+=c[j]; 49 } 50 for(int i=1;i<=n;i++){ 51 if(miu[i]) 52 for(int j=1;i*j<=n;j++){ 53 ans+=F(i*j)*miu[i]*phi[j]; 54 } 55 } 56 return ans; 57 } 58 int main(){ 59 init(); 60 T=read(); 61 while(T--){ 62 n=read(); 63 printf("%lld ",calc(n)); 64 } 65 return 0; 66 } 67 } 68 int main(){ 69 gengyf::main(); 70 return 0; 71 }
推式子很有趣的，大家一起来推式子啊
手动滑稽.jpg