这道题很简单,直接AC就可以了。
#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;
}