[BZOJ2190]SDOI2008仪仗队|线性筛

答案就是所有满足gcd（x，y）=1的x<n，y<n的点（x，y）。。方法是线性筛欧拉函数，和线性筛素数很像。利用3个结论（p为质数）①phi（p）=p-1②若i%p=0，phi（i*p）=phi(i)*p③若i%p！=0，phi(i*p)=phi(i)*phi(p)。。证明我也不会。。

#include<iostream>
#include<cstdio>
#include<memory.h>
#define N 40005
using namespace std;
int i,n,cnt=0,ans=0,phi[N],u[N],p[N];
void getphi()
{
    int i,j;
    memset(u,0,sizeof(u));
    phi[1]=1;
    for (i=2;i<=n;i++)
    {
        if (!u[i]) p[++cnt]=i,phi[i]=i-1;
        for (j=1;i*p[j]<=n&&j<=cnt;j++)
        {
            u[i*p[j]]=1;
            if (i%p[j]==0) {phi[i*p[j]]=p[j]*phi[i];break;} 
            else phi[i*p[j]]=(p[j]-1)*phi[i];
        }
    }
}
int main()
{
    cin>>n;
    getphi();
    for (i=1;i<n;i++) ans+=phi[i];
    cout<<2*ans+1;
}