只看一个象限简化问题，最后答案乘4+4
象限里面枚举x， 在当前这条固定的平行于y轴的直线中
分成长度为x的一段段。符合题目要求的点gcd（x，y） = 1
那么第一段1<= y <= x，个数为phi（x）个，即是x的欧拉函数值
第二段x+1 <= y <= 2x， 因为gcd(x + i, x) = gcd(x, i)
接下来同理
一直到最后一段
kx + 1 <= y <= b要单独一个个算，因为最后一段的长度不一定为x
感觉这道题算每一条直线时分成长度为x的一段段很巧妙，不好想到
#include<cstdio>
#include<cmath>
#define REP(i, a, b) for(int i = (a); i < (b); i++)
using namespace std;typedef long long ll;int phi(int n)
{int ans = n;for(int i = 2; i * i <= n; i++)if(n % i == 0){ans = ans / i * (i - 1);while(n % i == 0) n /= i;} if(n > 1) ans = ans / n * (n - 1);return ans;
}int gcd(int a, int b) { return !b ? a : gcd(b, a % b); }ll f(int a, int b)
{ll ans = 0;REP(x, 1, a + 1){int k = b / x;ans += phi(x) * k;REP(y, k*x + 1, b + 1)if(gcd(x, y) == 1)ans++; } return 4 * ans + 4;
}int main()
{int a, b;while(~scanf("%d%d", &a, &b) && a){ll K = f(a, b);ll N = (ll)(2 * a + 1) * (2 * b + 1) - 1;printf("%.7lf\n", (double)K / N);} return 0;
}