UVa10214 - Trees in a Wood.(数论)

只看某一个象限 能看到的数 == 一个 象限*4+4

能看到的树既距离原点的距离 gcd(x,y)==1

a 和 b 一大一小 预处理2000以内的phi函数,枚举小的一条边

从1...a 与 a gcd 为 1 的数的个数就是 phi(a)

从 1+a ... 2*a 与 a gcd 为 1 的数的个数 因为 GCD(i,a) = GCD(i+a,a) 所以还是 phi(a)

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<map>
#include<string>
#include<set>
#include<vector>
#include<cmath>
using namespace std;
int phi[2010];
int gcd(int a,int b){
    return b==0?a:gcd(b,a%b);
}
int main()
{
    memset(phi,0,sizeof(phi));
    phi[1]=1;
    for(int i=2;i<=2000;i++) if(phi[i]==0)
    for(int j=i;j<=2000;j+=i){
        if(!phi[j]) phi[j]=j;
        phi[j]=phi[j]/i*(i-1);
    }
    int n,m;
    while(scanf("%d%d",&n,&m)!=EOF&&n&&m){

        double sum=0;
        for(int i=1;i<=n;i++){
            sum+=phi[i]*(m/i);
            int mn=m%i;
            for(int j=1;j<=mn;j++) if(gcd(j,i)==1) sum+=1;
        }
        sum=sum*2+2;
        double sum2=(unsigned long long)n*m*2+n+m;
        double ans=sum/sum2;
        printf("%.7lf\n",ans);
    }
    return 0;
}