今天zky学长讲数论，上午水，舒爽的不行。。后来下午直接while（true）{懵逼；}死循全程懵逼。。。。（可怕）Thinking Bear。
2190: [SDOI2008]仪仗队
Time Limit: 10 Sec Memory Limit: 259 MB
Submit: 2092 Solved: 1325
[Submit][Status][Discuss]
Description
作为体育委员，C君负责这次运动会仪仗队的训练。仪仗队是由学生组成的N * N的方阵，为了保证队伍在行进中整齐划一，C君会跟在仪仗队的左后方，根据其视线所及的学生人数来判断队伍是否整齐(如下图)。
现在，C君希望你告诉他队伍整齐时能看到的学生人数。
Input
共一个数N。
Output
共一个数，即C君应看到的学生人数。
Sample Input
4
Sample Output
9
HINT
【数据规模和约定】
对于 100% 的数据，1 ≤ N ≤ 40000
Source
数论
这道题啊，算是我懵逼过程中最水的题了。。于是下课 顺利一遍水过（哦，第一次复制错了。。。CE。。不算！）
这道题，先小数据捯饬捯饬。。发现了规律：
1.除了坐标带0的点，其他的点都满足gcd（x，y）=1，所以累加欧拉函数φ。。然后加上坐标带0的点即可（2个，一成不变）。。
2.然后发现具有对称性。。于是只需要搞上方的三角即可。但是对角线那个点重复了一次所以-1
于是 ans=Σ（φ【1..n-1】）*2+1（+1其实是+2-1）
ps（程序里，欧拉函数的处理没有处理φ【1】所以有个初始值。。思路一样。。）
那么就是筛出欧拉函数即可。。（根据性质：）
性质① m是素数时，有φ(m)=m-1
性质② 当m、n互素时，φ(m*n)=φ(m)*φ(n)
性质③ 对一切正整数n，有φ(p^n)=[p^(n-1)]*(p-1)
奇怪的代码：
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define maxn 40010
long long prime[maxn],phi[maxn];
bool flag[maxn]={0};
int n;
long long ans;
void get_phi()
{
memset(flag,0,sizeof(flag));
flag[1]=1;
int cnt=0;
for (int i=2; i<=maxn; i++)
{
if (!flag[i])
prime[cnt++]=i,phi[i]=i-1;
for (int j=0; j<cnt && i*prime[j]<maxn; j++)
{
flag[i*prime[j]]=1;
if (i%prime[j]==0)
{
phi[i*prime[j]]=phi[i]*prime[j];
break;
}
else
phi[i*prime[j]]=phi[i]*(prime[j]-1);
}
}
}
int main()
{
scanf("%d",&n);
get_phi();
ans=2;
for (int i=2; i<n; i++)
ans+=phi[i];
printf("%lld",ans*2-1);
return 0;
}