結果

問題 No.3318 客に卵をかける
コンテスト
ユーザー みうね
提出日時 2025-10-19 17:36:30
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,790 bytes
コンパイル時間 2,488 ms
コンパイル使用メモリ 277,560 KB
実行使用メモリ 7,716 KB
最終ジャッジ日時 2025-10-26 11:40:29
合計ジャッジ時間 3,664 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other WA * 2
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

using ld = long double;
using ll = long long;

ld f(ll x, ld Q, ld R, int remain, ld p) {
    ld a = (remain - x) * (Q - (1.0L - p) * R);   // もみ消す
    ld b = Q * ((1.0L - powl(p, x)) / (1.0L - p));  // クビになる(ゼロ除算に注意)
    return a + b;
}
ld g(ll x, ld Q, ld R, ll remain, ld p) {
    ld a = -Q + (1 - p) * R;
    ld b = Q * powl(p, x);
    return a + b;
}
ld h(ll t, ll N, ld X, ld Y, ld P, ld Q, ld R) {
    ll skill = X + t * Y;
    ld p = 1.0L - 1.0L / skill;
    ll remain = N - t;
    if (skill == 1) {  // 必ず卵をかけてしまう場合
        if (remain == 1) {
            return P * t + Q;
        } else {
            return P * t + (Q - R) * (remain - 1) + Q;
        }
    }
    ll left = 0, right = remain;
    while (1 < right - left) {
        ll mid = (left + right) / 2;
        if (g(mid, Q, R, remain, p) < 0) {
            right = mid;
        } else {
            left = mid;
        }
    }
    return P * t + f(right, Q, R, remain, p);
}
int main(void) {
    auto solve = []() {
        ll N, X, Y;
        cin >> N >> X >> Y;
        ld P, Q, R;
        cin >> P >> Q >> R;
        ll l = 0, r = N;
        while (r - l > 3) {
            ll m1 = l + (r - l) / 3;
            ll m2 = r - (r - l) / 3;
            ld v1 = h(m1, N, X, Y, P, Q, R);
            ld v2 = h(m2, N, X, Y, P, Q, R);
            if (v1 < v2) {
                l = m1;
            } else {
                r = m2;
            }
        }
        ld ans = 0;
        for (ll t = l; t <= r; ++t) {
            ans = max(ans, h(t, N, X, Y, P, Q, R));
        }
        cout << fixed << setprecision(10) << ans << '\n';
    };
    int T;
    cin >> T;
    for (int i = 0; i < T; i++) {
        solve();
    }
}
0