在一个平面直角坐标系的第一象限内,如果一个点(x,y)与原点(0,0)的连线中没有通过其他任何点,则称该点在原点处是可见的。
例如,点(4,2)就是不可见的,因为它与原点的连线会通过点(2,1)。
部分可见点与原点的连线如下图所示:
编写一个程序,计算给定整数N的情况下,满足0≤x,y≤N0≤x,y≤N的可见点(x,y)的数量(可见点不包括原点)。
题意:从原点发出光线,可以照到的柱子的数量,光线会被第一个柱子挡住。
题解:
每一条直线都是y=kx,,被照到的柱子都是这条直线上的第一个点,
设这个点为x0,y0,那么y=y0/x0x,这个时候,既然是第一个点,那么y0和x0就不能有公约数,不然还可以再化简成更小的点,所以y0和x0是互质的,,,又关于y=x对称轴,两边的点的数量是一样的,所以只算下面的点就行,,,所以答案就是1-n之间的所有数的欧拉函数,即<=i的数与i这个数互质的个数,,,求和最后2+1就行了
这里用了一个欧拉函数线性模板,修该一下欧拉筛就行。
#include <iostream>
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e5+10;
int vis[N+19];
int prime[N+19];
int phi[N+10];
int cnt;
void init(int n) //欧拉函数线性模板
{
phi[1]=1;
vis[0]=vis[1]=1;
cnt=0;
for(int i=2; i<=n; i++)
{
if(vis[i]==0)
{
prime[cnt++]=i;
phi[i]=i-1;
}
for(int j=0; prime[j]*i<=n; j++)
{
vis[prime[j]*i]=1;
if(i%prime[j]==0)
{
phi[i*prime[j]]=phi[i]*prime[j];
break;
}
phi[i*prime[j]]=phi[i]*(prime[j]-1);
}
}
return ;
}
int main()
{
int c;
scanf("%d",&c);
init(2000);
int Case=0;
while(c--)
{ Case++;
int n;
scanf("%d",&n);
int sum=0;
for(int i=1; i<=n; i++)
{
sum+=phi[i];
}
printf("%d %d %d\n",Case,n,sum*2+1);
}
return 0;
}