UVA 10214 - Trees in a Wood.

题目的提示已经很明显了,预处理phi(x)函数,有一部分要直接统计,O(a^2)的复杂度。

细节参见代码:

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxa = 2005;
const int maxb = 2000005;
int n,m,phi[maxa+5];
ll a,b;
void init() {
    memset(phi,0,sizeof(phi));
    phi[1] = 1;
    for(int i=2;i<=maxa;i++) if(!phi[i])
    for(int j=i;j<=maxa;j+=i) {
        if(!phi[j]) phi[j] = j;
        phi[j] = phi[j] / i * (i-1);
    }
}
int gcd(int a,int b) {
    return b == 0 ? a : gcd(b, a%b);
}
int main() {
    init();
    while(~scanf("%lld%lld",&a,&b)) {
        if( !a && !b ) return 0;
        ll ans = 0;
        for(int i=1;i<=a;i++) {
            int v = b / i;
            ans += v*phi[i];
            for(int j=v*i+1;j<=b;j++) 
                if(gcd(i,j)==1) ans++;
        }
        ans = 4*ans + 4;
        ll N = 4*a*b + 2*a + 2*b;
        printf("%.7f\n",ans*1.0/N);
    }
    return 0;
}

你可能感兴趣的:(函数,数论,uva,ACM-ICPC,phi)