洛谷-2158 [SDOI2008]仪仗队

题目描述
作为体育委员，C君负责这次运动会仪仗队的训练。仪仗队是由学生组成的N * N的方阵，为了保证队伍在行进中整齐划一，C君会跟在仪仗队的左后方，根据其视线所及的学生人数来判断队伍是否整齐(如下图)。 现在，C君希望你告诉他队伍整齐时能看到的学生人数。
输入格式

共一个数N
输出格式
共一个数，即C君应看到的学生人数。

输入输出样例
输入 #1
4

输出 #1
9

说明/提示
【数据规模和约定】
对于 100% 的数据，1 ≤ N ≤ 40000

解释：很简单的结论题，设点(x,y)，那么一条斜线对应这$\frac{y}{x}$,那只要统计$y,x$互质的个数就好了。直接上欧拉函数模版就OK

#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cctype>
#include<cstdio>
#include<vector>
#include<string>
#include<queue>
#include<stack>
#include<cmath>
#include<map>
#include<set>
using namespace std;
const int inf=0x7fffffff;
const double eps=1e-10;
const double pi=acos(-1.0);
inline int read(){
    int x=0,f=1;char ch;ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-') f=0;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+(ch&15);ch=getchar();}
    if(f)return x;else return -x;
}
const int N=40010;
bool vis[N];
long long prim[N],phi[N];
void get_phi(int n){
	phi[1]=1;
  	for(int i=2;i<=n;i++){
    	if(!vis[i]){phi[i]=i-1;prim[++prim[0]]=i;}
    	for(int j=1;j<=prim[0]&&i*prim[j]<=n;j++){
      		vis[i*prim[j]]=1;
      		if(i%prim[j]==0){phi[i*prim[j]]=phi[i]*prim[j];break;}
      		else phi[i*prim[j]]=phi[i]*(prim[j]-1);
    	}
  	}
  	for(int i=1;i<=n;i++) phi[i]+=phi[i-1];
}
int n;
int main(){
	get_phi(40000);
    n=read();
    if(n==1) cout<<0<<endl;
    else cout<<2*phi[n-1]+1<<endl;
    return 0;
}