POJ 3090 Visible Lattice Points

假设t[i]为x=y=i时可以看到的点数

考虑x<y的情况，如果可以看到，那么直线的斜率为y/x，且x,y互素，因为如果两者不互素，令d=gcd(x,y)>1,那么（x/d)/(y/d)这个点一定是之前可以看到的点，(x,y)就被遮住了

 

当x>y时，与x<y一样，是对称的

 

所以t[i]=t[i-1]+2*PHI[i],i=0时，t[i]=3

 

代码：

#include<iostream> #include<memory.h> #include<string> #include<cstdio> #include<algorithm> #include<math.h> using namespace std; int n; int t[1005]; int eular(int x) { int i,j,res=1; for(i=2;i*i<=x;i++) { if(x%i==0) { res*=i-1; x/=i; while(x%i==0) { res*=i; x/=i; } } } if(x>1) res*=x-1; return res; } int main() { int i,j,T,ca=0; cin>>T; t[1]=3; for(i=2;i<=1000;i++) { j=eular(i); t[i]=t[i-1]+j*2; //cout<<i<<" "<<j<<endl; } while(T--) { scanf("%d",&n); printf("%d %d %d/n",++ca,n,t[n]); } return 0; }