結果

問題 No.5018 Let's Make a Best-seller Book
コンテスト
ユーザー e869120
提出日時 2023-09-26 11:50:44
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 35 ms / 400 ms
コード長 4,123 bytes
コンパイル時間 610 ms
コンパイル使用メモリ 71,564 KB
実行使用メモリ 24,396 KB
スコア 185,576
平均クエリ数 52.00
最終ジャッジ日時 2023-10-01 12:32:49
合計ジャッジ時間 8,442 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge15
純コード判定しない問題か言語
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 100
権限があれば一括ダウンロードができます

ソースコード

diff # 

// [モデル解法 8]
// 序盤に関しては、広告をできるだけ速く打つのが得
// しかし中盤になると、広告の次のターンでお金が足りなくなる場合がある
// そこで、以下のようなアルゴリズムにしたがって広告を打つことにする
// 
// 15 ターン目まで: 1100000 円以上で前週広告を打っていなければ打つ
// 25 ターン目まで: 1200000 円以上で前週広告を打っていなければ打つ
// 38 ターン目まで: 1300000 円以上で前週広告を打っていなければ打つ
// 44 ターン目まで: 1500000 円以上で前週広告を打っていなければ打つ

#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++) {
		bool AdFlag = false;
		if (t <= 15 && t >= Prev_Ad + 2 && Money >= 1100000) AdFlag = true;
		if (t <= 25 && t >= Prev_Ad + 2 && Money >= 1200000) AdFlag = true;
		if (t <= 38 && t >= Prev_Ad + 2 && Money >= 1300000) AdFlag = true;
		if (t <= 44 && t >= Prev_Ad + 2 && Money >= 1500000) AdFlag = true;

		// 広告を打つ場合
		if (AdFlag == true) {
			cout << "2 2" << endl;
			Prev_Ad = t;
		}

		// 広告を打たない場合 (終盤)
		else if (t >= 44) {
			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++) {
				double keisuu = 1.0 - 3.0 * 0.145 / sqrt(t);
				int ZaikoNum = pow(pow(1.05, P[i]) * Expected_D[i] * keisuu * 0.8 / 0.3, 2.0); // どの程度の在庫にすべきか？
				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) ÷ 売上予測値
			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