题目链接 https://cn.vjudge.net/problem/UVA-10214
【题意】
平面坐标系,在满足 |x|<=a,|y|<=b (a<2000,b<=2000000)|x|<=a,|y|<=b (a<2000,b<=2000000) 的网格中,除了原点之外的整点各种一棵树,树的半径忽略不计,但是可以互相遮挡,求从原点能看到多少颗树,输出百分比
【思路】
四个坐标轴只能看到一棵树,所以只考虑第一象限,最后的答案乘以4再加4,第一象限所有 x,yx,y 为正整数,能看到 (x,y)(x,y) 当且仅当 gcd(x,y)=1gcd(x,y)=1,由于 xx 范围小,考虑枚举 xx,第 xx 列能看到的树的个数为 0<y<=b0<y<=b 满足 gcd(x,y)=1gcd(x,y)=1 的个数,可以分块计算
1<=y<=x1<=y<=x 有 phi(i)phi(i) 个
x+1<=y<=2xx+1<=y<=2x 也有 phi(i)phi(i) 个,因为 gcd(x+i,x)=gcd(x,i)gcd(x+i,x)=gcd(x,i)
2x+1<=y<=3x2x+1<=y<=3x 也有 phi(i)phi(i) 个,因为 gcd(2x+i,x)=gcd(x,i)gcd(2x+i,x)=gcd(x,i)
……
kx+1<=y<=bkx+1<=y<=b 直接暴力枚举
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = 2005;
int gcd(int a,int b) {
if(0==b)
return a;
else
return gcd(b,a%b);
}
int phi[maxn];
void phi_table(int n) {
for (int i=2; i<=n; ++i)
phi[i]=0;
phi[1]=1;
for(int i=2; i<=n; ++i) {
if(0==phi[i]) {
for (int j=i; j<=n; j+=i) {
if(0==phi[j])
phi[j]=j;
phi[j]=phi[j]/i*(i-1);
}
}
}
}
ll a,b;
int main() {
phi_table(maxn-1);
while(scanf("%lld%lld",&a,&b)==2) {
if(0==a && 0==b)
break;
ll ans=0;
ll all=(2*a+1)*(2*b+1)-1;
for(int x=1; x<=a; ++x) {
int k=b/x;
int num=k*phi[x];
for(int y=k*x+1; y<=b; ++y) {
if(gcd(y,x)==1)
++num;
}
ans+=num;
}
ans=ans*4+4;
printf("%.7lf\n",(double)ans/all);
}
return 0;
}