#pragma GCC target("avx2") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #include using namespace std; using lint = long long int; using P = pair; using PL = pair; #define FOR(i, begin, end) for(int i=(begin),i##_end_=(end);i=i##_begin_;i--) #define REP(i, n) FOR(i,0,n) #define IREP(i, n) IFOR(i,0,n) #define ALL(a) (a).begin(),(a).end() constexpr int MOD = 1000000007; constexpr lint B1 = 1532834020; constexpr lint M1 = 2147482409; constexpr lint B2 = 1388622299; constexpr lint M2 = 2147478017; constexpr int INF = 2147483647; void yes(bool expr) {cout << (expr ? "Yes" : "No") << "\n";} templatevoid chmax(T &a, const T &b) { if (avoid chmin(T &a, const T &b) { if (b> T >> N; if(!DEBUG) cin >> Money; else Money = 2000000; vector R(N), P(N); vector D(N); if(DEBUG) REP(i, N) cin >> D[i]; vector> Z(T, vector(N)); if(DEBUG) REP(i, T) REP(j, N) cin >> Z[i][j]; vector weight(N); REP(i, N) weight[i] = 1.0 / N; int score = 0; REP(t, T) { /* cerr << t << endl; cerr << Money << endl; REP(i, N) cerr << R[i] << " "; cerr << endl; REP(i, N) cerr << P[i] << " "; cerr << endl; */ int mode = 1; vector L(N); int x = 1; // 行動を決める int MODE_CHANGE_TURN = 40; double THRESHOLD_D = 0.5 + 0.7 * t / T; double THRESHOLD_Z = 0.75; int limL = Money / 500 / N; if(t < MODE_CHANGE_TURN) { REP(i, N) { int l = 3; int r = 1e9; while(r-l > 1) { int m = (l+r)/2; if(sqrt(m) * pow(1.05, P[i]) * THRESHOLD_D * THRESHOLD_Z >= 0.3 * m) l = m; else r = m; } L[i] = max(0, l - R[i]); chmin(L[i], limL); } } else { REP(i, N) if(P[i] < 0) weight[i] = 0; double weight_sum = 0; REP(i, N) weight_sum += weight[i]; REP(i, N) weight[i] /= weight_sum; REP(i, N) L[i] = Money / 500 * weight[i]; } REP(i, N) assert(L[i] >= 0); if(Money >= 5000000) { mode = 2; x = 4; } if(Money >= 9000000) { mode = 2; x = 5; } cout << mode << " "; if(mode == 1) { REP(i, N) cout << L[i] << (i!=N-1 ? " " : ""); REP(i, N) R[i] += L[i]; REP(i, N) Money -= 500 * L[i]; } else { cout << x; REP(i, N) P[i] = min(60, P[i] + x); Money -= 500000 * pow(2, x-1); } cout << "\n"; cout.flush(); assert(Money >= 0); cerr << Money << endl; vector S(N); if(!DEBUG) { cin >> Money; if(Money == -1) { cerr << "Money is -1" << endl; return 0; } REP(i, N) cin >> S[i]; } //結果の受け取り REP(i, N) { int s = S[i]; if(DEBUG) s = min(R[i], (int)floor(sqrt(R[i]) * pow(1.05, P[i]) * D[i] * Z[t][i])); score += s; Money += 1000 * s; //double estimate_D = s / sqrt(R[i]) / pow(1.05, P[i]); //if(estimate_D > 1.0) weight[i] *= 1.0; //else weight[i] /= 1.0; if(R[i] > 0) { if(s >= 0.3 * R[i]) P[i] = min(60, P[i] + 1); else if(s <= 0.1 * R[i]) P[i] = max(-60, P[i] - 1); } R[i] -= s; } if(!DEBUG) { REP(i, N) cin >> P[i]; REP(i, N) cin >> R[i]; } } cerr << score/100 << endl; }