AcWing 220. 最大公约数

题意：

题解：

题目就变成了AcWing 201. 可见的点
当然有微调，因为可见的点里面是从0开始，本题从1开始，所以本题中phi[1]认为是0
AcWing 201. 可见的点的题解

代码：

#include<bits/stdc++.h>
#define debug(a,b) printf("%s = %d\n",a,b);
typedef long long ll;
using namespace std;inline int read(){int s=0,w=1;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}while(ch>='0'&&ch<='9') s=s*10+ch-'0',ch=getchar();//s=(s<<3)+(s<<1)+(ch^48);return s*w;
}
const int maxn=1e7+9;
int prime[maxn],cnt;
bool st[maxn];
int phi[maxn];
ll s[maxn];
void init(int n){phi[1]=0;//在本题中 for(int i=2;i<=n;i++){if(!st[i]){prime[cnt++]=i;phi[i]=i-1;}for(int j=0;prime[j]*i<=n;j++){st[prime[j]*i]=1;if(i%prime[j]==0){phi[i*prime[j]]=phi[i]*prime[j];break;}phi[i*prime[j]]=phi[i]*(prime[j]-1);}}for(int i=1;i<=n;i++)s[i]=s[i-1]+phi[i];
}
int main()
{int n,m;cin>>n;init(n);ll res=0;for(int i=0;i<cnt;i++){int p=prime[i];res+=s[n/p]*2+1;} cout<<res;return 0;
}