結果

問題 No.5018 Let's Make a Best-seller Book
ユーザー e869120e869120
提出日時 2023-09-25 00:15:07
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 34 ms / 400 ms
コード長 3,703 bytes
コンパイル時間 646 ms
コンパイル使用メモリ 70,700 KB
実行使用メモリ 24,480 KB
スコア 99,435
平均クエリ数 52.00
最終ジャッジ日時 2023-10-01 12:32:34
合計ジャッジ時間 8,395 ms
ジャッジサーバーID
(参考情報)
judge13 / judge11
純コード判定しない問題か言語
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 34 ms
24,264 KB
testcase_01 AC 28 ms
24,036 KB
testcase_02 AC 28 ms
24,336 KB
testcase_03 AC 28 ms
24,348 KB
testcase_04 AC 29 ms
24,360 KB
testcase_05 AC 29 ms
24,276 KB
testcase_06 AC 28 ms
24,408 KB
testcase_07 AC 29 ms
24,036 KB
testcase_08 AC 28 ms
23,628 KB
testcase_09 AC 28 ms
23,376 KB
testcase_10 AC 29 ms
23,376 KB
testcase_11 AC 29 ms
23,436 KB
testcase_12 AC 31 ms
23,640 KB
testcase_13 AC 29 ms
24,264 KB
testcase_14 AC 28 ms
23,652 KB
testcase_15 AC 30 ms
24,000 KB
testcase_16 AC 28 ms
24,012 KB
testcase_17 AC 30 ms
24,024 KB
testcase_18 AC 30 ms
23,652 KB
testcase_19 AC 29 ms
23,616 KB
testcase_20 AC 28 ms
24,012 KB
testcase_21 AC 29 ms
23,508 KB
testcase_22 AC 29 ms
23,832 KB
testcase_23 AC 28 ms
23,424 KB
testcase_24 AC 30 ms
24,384 KB
testcase_25 AC 29 ms
23,412 KB
testcase_26 AC 28 ms
23,832 KB
testcase_27 AC 28 ms
24,012 KB
testcase_28 AC 28 ms
24,036 KB
testcase_29 AC 29 ms
24,324 KB
testcase_30 AC 29 ms
23,628 KB
testcase_31 AC 29 ms
23,364 KB
testcase_32 AC 28 ms
24,036 KB
testcase_33 AC 28 ms
23,664 KB
testcase_34 AC 29 ms
23,532 KB
testcase_35 AC 29 ms
23,892 KB
testcase_36 AC 29 ms
23,376 KB
testcase_37 AC 28 ms
23,532 KB
testcase_38 AC 28 ms
23,376 KB
testcase_39 AC 28 ms
23,628 KB
testcase_40 AC 28 ms
24,312 KB
testcase_41 AC 27 ms
23,604 KB
testcase_42 AC 28 ms
23,520 KB
testcase_43 AC 28 ms
23,508 KB
testcase_44 AC 29 ms
24,060 KB
testcase_45 AC 28 ms
23,400 KB
testcase_46 AC 29 ms
24,012 KB
testcase_47 AC 28 ms
24,480 KB
testcase_48 AC 29 ms
23,376 KB
testcase_49 AC 28 ms
24,024 KB
testcase_50 AC 28 ms
24,336 KB
testcase_51 AC 28 ms
23,832 KB
testcase_52 AC 27 ms
23,508 KB
testcase_53 AC 28 ms
23,520 KB
testcase_54 AC 29 ms
24,324 KB
testcase_55 AC 28 ms
23,364 KB
testcase_56 AC 28 ms
24,480 KB
testcase_57 AC 28 ms
23,844 KB
testcase_58 AC 29 ms
23,628 KB
testcase_59 AC 29 ms
23,376 KB
testcase_60 AC 28 ms
24,012 KB
testcase_61 AC 28 ms
23,640 KB
testcase_62 AC 30 ms
24,048 KB
testcase_63 AC 28 ms
23,628 KB
testcase_64 AC 28 ms
23,376 KB
testcase_65 AC 28 ms
24,012 KB
testcase_66 AC 28 ms
23,424 KB
testcase_67 AC 28 ms
24,024 KB
testcase_68 AC 27 ms
24,024 KB
testcase_69 AC 28 ms
24,360 KB
testcase_70 AC 29 ms
24,372 KB
testcase_71 AC 28 ms
24,024 KB
testcase_72 AC 27 ms
23,424 KB
testcase_73 AC 28 ms
24,048 KB
testcase_74 AC 30 ms
24,336 KB
testcase_75 AC 28 ms
24,000 KB
testcase_76 AC 29 ms
23,628 KB
testcase_77 AC 29 ms
23,520 KB
testcase_78 AC 29 ms
24,276 KB
testcase_79 AC 28 ms
23,664 KB
testcase_80 AC 29 ms
24,036 KB
testcase_81 AC 28 ms
23,376 KB
testcase_82 AC 30 ms
23,244 KB
testcase_83 AC 28 ms
24,336 KB
testcase_84 AC 30 ms
23,640 KB
testcase_85 AC 28 ms
24,012 KB
testcase_86 AC 28 ms
23,412 KB
testcase_87 AC 28 ms
23,376 KB
testcase_88 AC 29 ms
24,336 KB
testcase_89 AC 28 ms
23,376 KB
testcase_90 AC 29 ms
23,640 KB
testcase_91 AC 29 ms
23,664 KB
testcase_92 AC 29 ms
24,348 KB
testcase_93 AC 29 ms
24,360 KB
testcase_94 AC 28 ms
24,012 KB
testcase_95 AC 30 ms
24,384 KB
testcase_96 AC 29 ms
24,000 KB
testcase_97 AC 29 ms
23,640 KB
testcase_98 AC 29 ms
23,376 KB
testcase_99 AC 29 ms
23,412 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

// [モデル解法 5]
// モデル解法 4 (https://yukicoder.me/submissions/915818) では、常に 3 割以上の売れ行きを目指していた
// しかし終盤は、人気度アップ (3 割以上売れ) にこだわらず、売上を貪欲に最大化するのが良い
// それでは、限られた資金で「次のターンの」売上部数の予測値を最大化するにはどうすれば良いか?
// (資金が限られているので、本の発送部数も制限される)
// 
// 現在の在庫を x とし、次週の予測売上を f(x) とする
// このとき、微分をすると f'(x) = x^{-0.5} * 1.05^{人気度} * D[i] となる
// そこで売上部数を最大化するには、どの書店も f'(x) の値が同じになるようにすれば良い
// このときの発送部数は、f'(x) の値で二分探索することで計算できる
// 
// なお、終盤の貪欲最大化は、45 ターン目以降に適用することにする
// また、47 ターン目以降には広告も打たないようにする

#include <iostream>
#include <cmath>
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 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 && t <= 46) {
			cout << "2 2" << endl;
			Prev_Ad = t;
		}

		// 広告を打たない場合 (終盤)
		else if (t >= 45) {
			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];
					int 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++) {
				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;
}
0