POJ3090 Visible Lattice Points

欧拉函数裸题。

除了（1，1），（0，1），（1，0）这三个点，只有gcd(x,y)=1的点才会被看到。（实际上可以算上（1，1））

对于每一个gcd(x,y)=1的点，对称的gcd(y,x)=1，也是一个可见点，就转换成了求比N小，与N互质的数的个数，×2+3即可。 

#include <cstring>
#include <cmath>
#include <iostream>
#define ms(a,b) memset(a,b,sizeof(a))
#define INF 0x3f3f3f3f
#define LL long long
using namespace std;
const int N=1e6+10;
int phi[N];
LL pre[N];
int prime[N],cnt;//素数筛
bool vis[N];
void eulersieve(int n)
{
	ms(vis,0);
	for(int i=2;i<=n;i++)
	{
		if(!vis[i])
		{
			prime[cnt++]=i;
			phi[i]=i-1;
		}
		for(int j=0;j<cnt&&prime[j]*i<=n;j++)
		{
			vis[prime[j]*i]=1;
			if(i%prime[j]==0)
			{
				phi[i*prime[j]]=prime[j]*phi[i];
				break;
			}
			else
				phi[i*prime[j]]=(prime[j]-1)*phi[i];
		}	
	}
}
int main()
{
    int c,n;
    eulersieve(1000);
    for(int i=2;i<=1000;i++)
        pre[i]=pre[i-1]+phi[i];
    cin>>c;
    for(int it=1;it<=c;it++)
    {
        cin>>n;
        cout<<it<<' '<<n<<' '<<3+2*pre[n]<<endl;
    }
}