bzoj 2190: [SDOI2008]仪仗队.

 1 #include<cstdio>
 2 using namespace std;
 3 const int maxn = 40010;
 4 int phi[maxn],prime[maxn],tot,n,ans=0;
 5 bool vis[maxn];
 6 void GetPhi() {
 7     phi[1]=1;
 8     for(int i=2;i<=n;i++) {
 9         if(vis[i]==false) prime[++tot]=i,phi[i]=i-1;
10         for(int j=1;j<=tot&&i*prime[j]<=n;j++) {
11             vis[i*prime[j]]=true;
12             if(i%prime[j]==0) {
13                 phi[i*prime[j]]=phi[i]*prime[j];
14                 break;
15             } else phi[i*prime[j]]=phi[i]*(prime[j]-1);
16         }
17     }
18     return;
19 }
20 int main() {
21     scanf("%d",&n);
22     GetPhi();
23     for(int i=1;i<n;i++) ans+=phi[i];
24     printf("%d\n",(ans<<1)+1);
25     return 0;
26 }