// [モデル解法 4] // モデル解法 3 (https://yukicoder.me/submissions/915817) に加え、発送部数の計算に D[i] の情報を使うようにする // 具体的には、D[i] の予測値を d とするとき、「(1.05^{人気度} * D[i] * 0.8 / 0.3) の 2 乗」冊の在庫を目指す // ちなみに 0.8 という定数は rand(0.75, 1.25) の中の小さめの値 (確実に 3 割以上売りたいので) // // [D の値をどうやって予測するか?] // 現在のターンまでにおける、(売上 + 0.5) / [売上予測値] の平均 // ここで、売上予測値は 1.05^{週の始まり時点の人気度} * 在庫^{0.5} 冊 // なお、+0.5 が式に含まれている理由は、売上部数の式が小数点以下切り捨てであるため #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]; 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 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++) { // 広告を打つ場合 if (Money >= 1500000 && t >= Prev_Ad + 2) { cout << "2 2" << endl; Prev_Ad = t; } // 広告を打たない場合 else { int NeedMoney = 0; for (int i = 1; i <= N; i++) { int ZaikoNum = pow(pow(1.05, P[i]) * Expected_D[i] * 0.8 / 0.3, 2.0); // どの程度の在庫にすべきか? 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) ÷ 売上予測値 Sum_D[i] += Grade; Cnt_D[i] += 1.0; Expected_D[i] = max(0.5, min(1.5, Sum_D[i] / Cnt_D[i])); // D[i] は 0.5 以上 1.5 以下であるため } } return 0; }