POJ - 3090 Visible Lattice Points 欧拉函数

题目链接：

https://vjudge.net/problem/POJ-3090

题意：

输入一个n，问从（0,0）点能看到多少个点，坐标范围   0<=i,j<=n;

题解：

通过观察我们能够发现，除了（0,1）（1,1）（1,0）这三个点，其它能看到的点满足：横坐标x和纵坐标y互质，并且经过直线y=x对称。然后我们只要算出y=x直线下的点，乘2就是答案，y=x直线下的点，横坐标为x能看到的点的个数是euler（x）。euler（x）是欧拉函数。可以线性筛求出欧拉函数然后求前缀和就可以了，然后每次查询都是O(1)的复杂度。

代码：

#include<iostream>
#include<cstdio>
using namespace std;
const int maxn=1005;
int v[maxn],prime[maxn],cnt,phi[maxn];
int sum[maxn];
void init(int n){
    for(int i=2;i<maxn;i++){
        if(!v[i]){
            prime[cnt++]=i;
            v[i]=i;
            phi[i]=i-1;
        }
        for(int j=0;j<cnt;j++){
            if(prime[j]>v[i]||prime[j]>n/i) break;
            v[i*prime[j]]=prime[j];
            phi[i*prime[j]]=phi[i]*(i%prime[j]?prime[j]-1:prime[j]);
        }
    }
    sum[1]=3;
    for(int i=2;i<=1000;i++) sum[i]=sum[i-1]+2*phi[i];
}
int main(){
    init(1005);
    int t;
    cin>>t;
    for(int i=1;i<=t;i++){
        int x;
        scanf("%d",&x);
        printf("%d %d %d\n",i,x,sum[x]);
    }
    return 0;
}