// [モデル解法 10] // モデル解法 9 (https://yukicoder.me/submissions/915989) に加えてパラメータ調整をする // そうすると、20 万点を越えることができる // ベストセラーおめでとうございます! #include #include using namespace std; long long T = 52; long long N = 10; long long Money; long long S[19], P[19], R[19]; long long Prev_Ad = 0; long long Answer[19], Memo[19]; double Expected_D[19]; double Sum_D[19]; double Cnt_D[19]; void NextInput() { cin >> Money; for (int i = 1; i <= N; i++) cin >> S[i]; for (int i = 1; i <= N; i++) cin >> P[i]; for (int i = 1; i <= N; i++) cin >> R[i]; } int GoingUp(int zaiko) { if (zaiko % 10 == 0) return zaiko + 0; if (zaiko % 10 == 1) return zaiko + 2; if (zaiko % 10 == 2) return zaiko + 1; if (zaiko % 10 == 3) return zaiko + 0; if (zaiko % 10 == 4) return zaiko + 2; if (zaiko % 10 == 5) return zaiko + 1; if (zaiko % 10 == 6) return zaiko + 0; if (zaiko % 10 == 7) return zaiko + 3; if (zaiko % 10 == 8) return zaiko + 2; if (zaiko % 10 == 9) return zaiko + 1; return -1; } int main() { // 最初の入力 cin >> T >> N >> Money; for (int i = 1; i <= N; i++) P[i] = 0; for (int i = 1; i <= N; i++) Expected_D[i] = 1.0; // シミュレーション開始 for (int t = 1; t <= T; t++) { // P[i] = 60 が何個あるかを数える int count60 = 0; for (int i = 1; i <= N; i++) { if (P[i] == 60) count60 += 1; } // 広告を打つかどうかを決める bool AdFlag = false; if (t <= 15 && t >= Prev_Ad + 2 && Money >= 1100000 && count60 <= 2) AdFlag = true; if (t <= 25 && t >= Prev_Ad + 2 && Money >= 1200000 && count60 <= 2) AdFlag = true; if (t <= 38 && t >= Prev_Ad + 2 && Money >= 1300000 && count60 <= 2) AdFlag = true; if (t <= 52 && t >= Prev_Ad + 3 && Money >= 1500000 && count60 <= 2) AdFlag = true; // 広告を打つ場合 if (AdFlag == true) { cout << "2 2" << endl; Prev_Ad = t; } // 広告を打たない場合 (終盤) else if (count60 >= 1) { double cl = 0.0, cr = 10000.0, cm; // f'(x) の値で二分探索 for (int i = 0; i < 40; i++) { cm = (cl + cr) / 2.0; int sum = 0; // f'(x) = cm まで在庫を増やすとき、合計何冊発注する必要があるか? for (int j = 1; j <= N; j++) { double keisuu = pow(1.05, P[j]) * Expected_D[j]; double Level3 = pow(pow(1.05, P[j]) * Expected_D[j] * (1.0 - 3.4 * 0.145 / sqrt(t)) * 0.79 / 0.3, 2.0); // 3 割売れの基準 int ZaikoNum = -1; if (P[j] >= 60) ZaikoNum = pow(keisuu / cm, 2.0); if (P[j] >= 50 && P[j] <= 59) ZaikoNum = min(1.15 * Level3, pow(keisuu / cm, 2.0)); if (P[j] >=-60 && P[j] <= 49) ZaikoNum = pow(keisuu / cm, 2.0); Memo[j] = max(0LL, ZaikoNum - R[j]); sum += Memo[j]; } if (500 * sum <= Money) { cr = cm; for (int j = 1; j <= N; j++) Answer[j] = Memo[j]; } else { cl = cm; } } cout << 1; for (int i = 1; i <= N; i++) cout << " " << Answer[i]; cout << endl; } // 広告を打たない場合 (序盤・中盤) else { int NeedMoney = 0; for (int i = 1; i <= N; i++) { double keisuu = 1.0 - 3.4 * 0.145 / sqrt(t); int ZaikoNum = 0; // どの程度の在庫にすべきか? if (P[i] >= 53 && P[i] <= 60) ZaikoNum = pow(pow(1.05, P[i]) * Expected_D[i] * keisuu * 0.79 / 0.1, 2.0) * 1.0; // 1 割売れを目指す if (P[i] >= 45 && P[i] <= 52) ZaikoNum = pow(pow(1.05, P[i]) * Expected_D[i] * keisuu * 0.79 / 0.3, 2.0) * 1.2; // リスキーな 3 割売れを目指す if (P[i] >= 30 && P[i] <= 44) ZaikoNum = pow(pow(1.05, P[i]) * Expected_D[i] * keisuu * 0.79 / 0.3, 2.0) * 1.1; // ややリスキーな 3 割売れを目指す if (P[i] >=-60 && P[i] <= 29) ZaikoNum = pow(pow(1.05, P[i]) * Expected_D[i] * keisuu * 0.79 / 0.3, 2.0) * 1.0; // 安全な 3 割売れを目指す ZaikoNum = GoingUp(ZaikoNum); Answer[i] = max(0LL, ZaikoNum - R[i]); NeedMoney += 500 * Answer[i]; } if (NeedMoney > Money) { for (int i = 1; i <= N; i++) Answer[i] = Answer[i] * Money / NeedMoney; } cout << 1; for (int i = 1; i <= N; i++) cout << " " << Answer[i]; cout << endl; } // 次のターンの入力 + 予測値の更新 NextInput(); for (int i = 1; i <= N; i++) { int Popularity = P[i]; if (10 * S[i] >= 3 * (S[i] + R[i])) Popularity -= 1; // 週の始まりの時点での人気度を計算 if (10 * S[i] < 1 * (S[i] + R[i])) Popularity += 1; // 週の始まりの時点での人気度を計算 double Grade = (0.5 + S[i]) / (pow(1.05, Popularity) * pow(S[i] + R[i], 0.5)); // (S[i] + 0.5) ÷ 売上予測値 double keisuu = 1.0 - pow(1.0 / S[i], 2.0); if (S[i] != 0) { Sum_D[i] += keisuu * Grade; Cnt_D[i] += keisuu; } if (Cnt_D[i] < 0.01) Expected_D[i] = 0.6; else Expected_D[i] = max(0.5, min(1.5, Sum_D[i] / Cnt_D[i])); // D[i] は 0.5 以上 1.5 以下であるため } } return 0; }