#include using namespace std; using ll = long long; #define TEST 0 double rnd() { static mt19937 e(1); static double t = pow(0.5, 32.0); return e() * t; } int main() { #if TEST == 0 int t, n, m; cin >> t >> n >> m; { int s[10] = {}, p[10] = {}, r[10] = {}; double s3[10] = {}, r3[10] = {}; int s4[10] = {}; while (t--) { #else int tn = 1000, s0 = 1 << 30, s1 = -1, ss = 0; for (int tt = 0; tt < tn; tt++) { int t = 52, n = 10, m = 2000000; int p[10] = {}, r[10] = {}; double s3[10] = {}, r3[10] = {}; int s4[10] = {}; int sc = 0; double d[10]; for (int i = 0; i < n; i++) { d[i] = 0.5 + rnd(); if (tn == 1) { cout << d[i] << " \n"[i == n - 1]; } } while (t--) { #endif int l[10] = {}; int x = 0; if (t < 10) { double k = 0; for (int i = 0; i < n; i++) { double d1 = s3[i] / r3[i]; k += pow(d1 * pow(1.05, p[i]), 2); } for (int i = 0; i < n; i++) { double d1 = s3[i] / r3[i]; int x = m * pow(d1 * pow(1.05, p[i]), 2) / k / 500; //x = max(x, 5); //l[i] = max(x - r[i], 0); l[i] = x; } } else { for (int i = 0; i < n; i++) { int x = 9; if (r3[i] > 0) { double d1 = s3[i] / r3[i]; x = (int)round(pow(pow(1.05, p[i]) * d1 / 0.5, 2)); //x = (x + 10) / 2; x = max(x, 5); } //cout << x << " \n"[i == n - 1]; l[i] = min(m / (500 * n), max(x - r[i], 0)); } //if (t >= 52 - 1) x = 2; if (m >= 1750000) x = 2; } #if TEST == 0 if (x == 0) { cout << 1; for (int i = 0; i < n; i++) { cout << ' ' << l[i]; } cout << endl; } else { cout << 2 << ' ' << x << endl; } cin >> m; if (m < 0) exit(0); for (int i = 0; i < n; i++) { cin >> s[i]; } int p2[10]; for (int i = 0; i < n; i++) { p2[i] = p[i]; cin >> p[i]; } for (int i = 0; i < n; i++) { cin >> r[i]; s4[i] += s[i]; s3[i] += (s[i] + 0.0) * pow(1.05, -p2[i]); r3[i] += sqrt(r[i]); } } } #else if (x == 0) { for (int i = 0; i < n; i++) { r[i] += l[i]; m -= l[i] * 500; if (m < 0) throw; } } else { m -= 250000 << x; if (m < 0) throw; for (int i = 0; i < n; i++) { p[i] += x; p[i] = min(p[i], 60); } } int s[10] = {}; for (int i = 0; i < n; i++) { s[i] = min(r[i], (int)floor(sqrt(r[i]) * pow(1.05, p[i]) * d[i] * (0.75 + 0.5 * rnd()))); m += s[i] * 1000; sc += s[i]; s4[i] += s[i]; s3[i] += (s[i] + 0.0) * pow(1.05, -p[i]); r3[i] += sqrt(r[i]); } for (int i = 0; i < n; i++) { if (r[i] > 0) { int u = 10 * s[i] / r[i]; p[i] += u >= 3; p[i] -= u < 1; p[i] = clamp(p[i], -60, 60); } } for (int i = 0; i < n; i++) { r[i] -= s[i]; } if (tn == 1) { cout << m << '\n'; for (int i = 0; i < n; i++) { cout << s[i] << " \n"[i == n - 1]; } for (int i = 0; i < n; i++) { cout << p[i] << " \n"[i == n - 1]; } for (int i = 0; i < n; i++) { cout << r[i] << " \n"[i == n - 1]; } } } s0 = min(s0, sc); s1 = max(s1, sc); ss += sc; } cout << s0 << ' ' << ss / tn << ' ' << s1 << endl; #endif return 0; }