#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using uint = unsigned int; using ll = long long; using ull = unsigned long long; constexpr ll TEN(int n) { return (n==0) ? 1 : 10*TEN(n-1); } template using V = vector; template using VV = V>; using D = long double; ll rand_int(ll l, ll r) { //[l, r] using Di = uniform_int_distribution; static random_device rd; static mt19937 gen(rd()); return Di(l, r)(gen); } D calc(V v) { D sm = 0; for (ll d: v) sm += sqrt(D(d)); return sm; } const int M = 30; V solve(D _x) { V v_base; { v_base.push_back(1842018); v_base.push_back(1842018); v_base.push_back(1160398); v_base.push_back(1160398); v_base.push_back(731005); v_base.push_back(731005); D x = _x - calc(v_base); for (int i = 6; i < 30; i++) { D y = x/(M-i); v_base.push_back(ll(y*y)); x -= sqrt(v_base.back()); } v_base.back()++; } while (true) { V v = v_base; D x = calc(v) - _x; for (int ph = 0; ph < 2000; ph++) { int i = rand_int(6, M-1); int j = rand_int(6, M-1); if (i == j) continue; if (v[i] > v[j]) swap(v[i], v[j]); D a = sqrt(D(v[i])) + sqrt(D(v[j])); D b = sqrt(D(v[i]-1)) + sqrt(D(v[j]+1)); if (x < a-b) continue; v[i]--; v[j]++; x -= (a-b); } x = calc(v) - _x; for (int z = 2; z >= 0; z--) { while (true) { int i = 2*z; int j = 2*z+1; D a = sqrt(D(v[i])) + sqrt(D(v[j])); D b = sqrt(D(v[i]-1)) + sqrt(D(v[j]+1)); if (x < a-b) break;; v[i]--; v[j]++; x -= (a-b); } } if (abs(calc(v) - _x) < 0.95e-10) { return v; } } assert(false); } int main() { cin.tie(0); ios::sync_with_stdio(false); cout << setprecision(20); // for (int ph = 0; ph < 1000; ph++) { // solve(10000); // } // return 0; int q; cin >> q; for (int ph = 0; ph < q; ph++) { D x; cin >> x; auto v = solve(x); // cout << x << " " << calc(v) << endl; // cout << abs(x - calc(v)) << endl; cout << v.size(); for (auto d: v) { cout << " " << d; } cout << endl; } return 0; }