#include #include #include #include #include #include using namespace std; #define getchar getchar_unlocked #define putchar putchar_unlocked bool dot; const int H = 17; const int N = 1 << H; const double L = 0; const double R = log(10); const double LL = -1e100; const double RR = log(R); const double dx = R / (1 << H); const double logDX = log(R) - H * log(2); const double log00001 = -4 * log(10); double lg[N + 1]; bool prime[N + 1]; vector gr[11][11]; int in() { int n, c; while ((c = getchar()) < '0') { assert(c != EOF); } n = c - '0'; while ((c = getchar()) >= '0') { n = n * 10 + c - '0'; } if (c == '.') { dot = true; } return n; } void out(double x) { static char buf[20]; sprintf(buf, "%.9f", x); for (int i = 0; buf[i] != '\0'; i++) { putchar(buf[i]); } putchar('\n'); } void sieve() { lg[0] = log(0); for (int i = 2; i <= N; i++) { if (!prime[i]) { double g = log(i); for (long long j = i; j < N; j *= i) { for (long long k = j; k < N; k += j) { prime[k] = true; lg[k] += g; } } } } for (int i = 0; i <= N; i++) { lg[i] += logDX; } } int binSearch(int a, int b, double logT) { int l = 0; int r = N; while (r - l > 1) { int mid = (l + r) / 2; (a * mid * dx + b * lg[mid] >= logT ? r : l) = mid; } return l; } int main() { int m = in(); sieve(); vector a(m), b(m), c(m), d(m); for (int i = 0; i < m; i++) { a[i] = in(); b[i] = in(); c[i] = in(); d[i] = dot ? in() : 0; c[i] = c[i] * 10000 + d[i]; gr[a[i]][b[i]].push_back(i); } vector ans(m); for (int a = 0; a <= 10; a++) { for (int b = 0; b <= 10; b++) { if (gr[a][b].empty()) { continue; } if (b == 0) { for (int i : gr[a][b]) { const double logT = lg[c[i]] - logDX + log00001; ans[i] = exp(logT / a); } continue; } if (a == 0) { for (int i : gr[a][b]) { const double logT = lg[c[i]] - logDX + log00001; ans[i] = exp(exp(logT / b)); } continue; } if (gr[a][b].size() <= 500) { for (int i : gr[a][b]) { const double logT = lg[c[i]] - logDX + log00001; int l = binSearch(a, b, logT); double s = (lg[l + 1] - lg[l]) / dx; ans[i] = exp((logT + s * b * l * dx - lg[l] * b) / (s * b + a)); } } else { sort(gr[a][b].begin(), gr[a][b].end(), [&](int i, int j) { return c[i] < c[j]; }); int l = 1; for (int i : gr[a][b]) { const double logT = lg[c[i]] - logDX + log00001; while (l < N && a * l * dx + b * lg[l] < logT) { l++; } double s = (lg[l] - lg[l - 1]) / dx; ans[i] = exp((logT + s * b * l * dx - lg[l - 1] * b) / (s * b + a)); } } } } for (int i = 0; i < m; i++) { out(ans[i]); } }