#include #include #include #include #include using namespace std; const double pi = acos(-1); mt19937 mt(123456); void solve() { int t; double p, w, v, gx, gy; cin >> t >> p >> w >> v >> gx >> gy; double gr = hypot(gx, gy); double mn = 1e100; for (int ii = 0; ii < 500; ii++) { double x = 1; double y = 0; vector vx(t); vector vy(t); double det = hypot(1 + v, w); double theta = uniform_real_distribution(0, 2*pi)(mt); double c = det * exp(-(double)ii / 500 * 10); vx[0] = cos(theta) - 1; vy[0] = sin(theta); for (int i = 0; i < t; i++) { double x0 = x; double y0 = y; x = x0 - y0*w + x0*v + vx[i]; y = y0 + x0*w + y0*v + vy[i]; double r = hypot(x, y); double tr; tr = 1; if (i == t - 2) tr = gr * c; double d = tr - r; vx[i] += d*x / r; vy[i] += d*y / r; x += d*x / r; y += d*y / r; } vx[t - 1] += gx - x; vy[t - 1] += gy - y; double s = 0; for (int i = 0; i < t; i++) { s += vx[i]*vx[i] + vy[i]*vy[i]; } mn = min(mn, s); if (s < p - 0.001) { for (int i = 0; i < t; i++) { printf("%.16f %.16f\n", vx[i], vy[i]); } fprintf(stderr, "Sum: %.16f\n", s); fprintf(stderr, "p: %.16f\n", p); return; } } cout << "FAIL" << endl; cout << endl; cout << mn << endl; } int main() { int T; cin >> T; while (T--) solve(); }