题目的提示已经很明显了,预处理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;
}