BZOJ 2190: [SDOI2008]仪仗队 欧拉函数

Description

作为体育委员，C君负责这次运动会仪仗队的训练。仪仗队是由学生组成的N * N的方阵，为了保证队伍在行进中整齐划一，C君会跟在仪仗队的左后方，根据其视线所及的学生人数来判断队伍是否整齐(如下图)。

现在，C君希望你告诉他队伍整齐时能看到的学生人数。

---

Input

　　共一个数N。

---

Output

　　共一个数，即C君应看到的学生人数。

---

Sample Input

　　4

---

Sample Output

　　9

---

HINT

【数据规模和约定】 　　对于 100% 的数据，1 ≤ N ≤ 40000

---

题解：经过观察可以发现对于一个点P（x，y），只要GCD（x，y）不等于1，则说明不可以被看到，所以直接裸欧拉函数就行了

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int MAXN=40005;
const int INF=0x3f3f3f3f;
int phi[MAXN],p[MAXN];
int N,ans;
bool vis[MAXN];
void getphi(){
    phi[1]=1;
    for(register int i=2;i<=N+1;i++){
        if(!vis[i]){
            p[++p[0]]=i;
            phi[i]=i-1;
        }
        for(register int j=1;j<=p[0]&&i*p[j]<=N+1;j++){
            vis[i*p[j]]=true;
            if(i%p[j]==0){
                phi[i*p[j]]=phi[i]*p[j];
                break;
            }else{
                phi[i*p[j]]=phi[i]*(p[j]-1);
            }
        }
    }
}

int main(){
    scanf("%d",&N);
    getphi();
    for(register int i=1;i<=N-1;i++){
        ans+=phi[i];
    }
    printf("%d\n",ans*2+1);
    return 0;
}