实质上就是求欧拉函数值

书上有个板子挺好，也不难理解。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<climits>
#include<cctype>
#include<queue>
#include<set>using namespace std;typedef long long ll;
const int INF=0x3f3f3f3f;
const int MAXN=5e4+5;
int phi[MAXN];
int n;
ll sum[MAXN];void phi_table()
{phi[1]=1;for(int i=2;i<MAXN;i++){if(!phi[i]){for(int j=i;j<MAXN;j+=i){if(!phi[j]) phi[j]=j;phi[j]=phi[j]/i*(i-1);}}}
}int main()
{phi_table();for(int i=2;i<MAXN;i++){sum[i]=sum[i-1]+phi[i];}while(~scanf("%d",&n) && n){if(n==1){printf("1\n");continue;} printf("%lld\n",sum[n]*2+1);}return 0;
}