#915992 (C++14) No.5018 Let's Make a Best-seller Book

提出ソース

結果

問題	No.5018 Let's Make a Best-seller Book
ユーザー	e869120
提出日時	2023-09-26 12:11:01
言語	C++14 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	34 ms / 400 ms
コード長	4,965 bytes
コンパイル時間	606 ms
コンパイル使用メモリ	72,888 KB
実行使用メモリ	24,396 KB
スコア	200,660
平均クエリ数	52.00
最終ジャッジ日時	2023-10-01 12:32:53
合計ジャッジ時間	8,114 ms
ジャッジサーバーID （参考情報）	judge12 / judge13
	純コード判定しない問題か言語

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	34 ms 23,520 KB
testcase_01	AC	27 ms 23,364 KB
testcase_02	AC	26 ms 23,376 KB
testcase_03	AC	26 ms 24,264 KB
testcase_04	AC	27 ms 24,372 KB
testcase_05	AC	27 ms 24,384 KB
testcase_06	AC	27 ms 24,384 KB
testcase_07	AC	27 ms 23,508 KB
testcase_08	AC	27 ms 23,520 KB
testcase_09	AC	27 ms 23,844 KB
testcase_10	AC	27 ms 23,508 KB
testcase_11	AC	27 ms 23,532 KB
testcase_12	AC	27 ms 23,832 KB
testcase_13	AC	28 ms 23,652 KB
testcase_14	AC	28 ms 23,664 KB
testcase_15	AC	27 ms 24,324 KB
testcase_16	AC	27 ms 23,388 KB
testcase_17	AC	26 ms 24,012 KB
testcase_18	AC	27 ms 23,508 KB
testcase_19	AC	26 ms 23,616 KB
testcase_20	AC	27 ms 23,412 KB
testcase_21	AC	27 ms 23,364 KB
testcase_22	AC	27 ms 24,012 KB
testcase_23	AC	27 ms 23,400 KB
testcase_24	AC	27 ms 24,384 KB
testcase_25	AC	27 ms 24,360 KB
testcase_26	AC	28 ms 23,508 KB
testcase_27	AC	28 ms 23,508 KB
testcase_28	AC	27 ms 23,508 KB
testcase_29	AC	27 ms 23,628 KB
testcase_30	AC	27 ms 24,324 KB
testcase_31	AC	27 ms 24,024 KB
testcase_32	AC	27 ms 23,364 KB
testcase_33	AC	27 ms 23,520 KB
testcase_34	AC	27 ms 24,312 KB
testcase_35	AC	27 ms 24,372 KB
testcase_36	AC	27 ms 23,712 KB
testcase_37	AC	26 ms 23,616 KB
testcase_38	AC	27 ms 24,384 KB
testcase_39	AC	26 ms 23,376 KB
testcase_40	AC	26 ms 23,640 KB
testcase_41	AC	27 ms 24,324 KB
testcase_42	AC	27 ms 23,520 KB
testcase_43	AC	28 ms 24,252 KB
testcase_44	AC	27 ms 23,376 KB
testcase_45	AC	27 ms 24,336 KB
testcase_46	AC	27 ms 23,664 KB
testcase_47	AC	26 ms 23,628 KB
testcase_48	AC	27 ms 24,324 KB
testcase_49	AC	30 ms 23,628 KB
testcase_50	AC	27 ms 24,384 KB
testcase_51	AC	26 ms 23,376 KB
testcase_52	AC	27 ms 24,024 KB
testcase_53	AC	27 ms 23,832 KB
testcase_54	AC	27 ms 23,640 KB
testcase_55	AC	27 ms 24,264 KB
testcase_56	AC	26 ms 24,024 KB
testcase_57	AC	28 ms 23,532 KB
testcase_58	AC	28 ms 24,324 KB
testcase_59	AC	27 ms 23,520 KB
testcase_60	AC	27 ms 24,036 KB
testcase_61	AC	29 ms 23,640 KB
testcase_62	AC	28 ms 24,396 KB
testcase_63	AC	27 ms 23,628 KB
testcase_64	AC	27 ms 24,036 KB
testcase_65	AC	27 ms 24,060 KB
testcase_66	AC	28 ms 23,628 KB
testcase_67	AC	27 ms 23,520 KB
testcase_68	AC	27 ms 23,364 KB
testcase_69	AC	27 ms 24,012 KB
testcase_70	AC	27 ms 23,520 KB
testcase_71	AC	28 ms 24,372 KB
testcase_72	AC	28 ms 23,520 KB
testcase_73	AC	27 ms 23,640 KB
testcase_74	AC	27 ms 24,264 KB
testcase_75	AC	28 ms 23,376 KB
testcase_76	AC	28 ms 24,012 KB
testcase_77	AC	27 ms 24,264 KB
testcase_78	AC	27 ms 23,832 KB
testcase_79	AC	28 ms 23,628 KB
testcase_80	AC	28 ms 23,412 KB
testcase_81	AC	28 ms 24,060 KB
testcase_82	AC	27 ms 24,384 KB
testcase_83	AC	28 ms 24,012 KB
testcase_84	AC	28 ms 23,520 KB
testcase_85	AC	27 ms 23,412 KB
testcase_86	AC	28 ms 23,532 KB
testcase_87	AC	28 ms 23,508 KB
testcase_88	AC	27 ms 24,024 KB
testcase_89	AC	27 ms 23,376 KB
testcase_90	AC	27 ms 24,264 KB
testcase_91	AC	27 ms 24,372 KB
testcase_92	AC	28 ms 24,348 KB
testcase_93	AC	28 ms 24,264 KB
testcase_94	AC	27 ms 23,820 KB
testcase_95	AC	28 ms 23,652 KB
testcase_96	AC	27 ms 23,376 KB
testcase_97	AC	28 ms 24,384 KB
testcase_98	AC	27 ms 24,036 KB
testcase_99	AC	28 ms 23,520 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

// [モデル解法 10]
// モデル解法 9 (https://yukicoder.me/submissions/915989) に加えてパラメータ調整をする
// そうすると、20 万点を越えることができる
// ベストセラーおめでとうございます！

#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 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;
}