結果

問題 No.2158 X日後に全完するhibit君
コンテスト
ユーザー KowerKoint2010
提出日時 2022-11-21 21:52:18
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 157 ms / 2,000 ms
コード長 2,376 bytes
コンパイル時間 4,177 ms
コンパイル使用メモリ 200,752 KB
最終ジャッジ日時 2025-02-08 23:08:39
ジャッジサーバーID
（参考情報） 		judge5 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 25
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

int main() {
    int n, k; cin >> n >> k;
    vector<int> p(n), s(n), t(n);
    for(int i = 0; i < n; i++) cin >> p[i] >> s[i] >> t[i];
    vector<int> prodp(n+1);
    prodp[0] = 1;
    for(int i = 0; i < n; i++) prodp[i+1] = prodp[i] * (p[i] + 1);
    auto idx_encode = [&](vector<int> data) {
        int idx = 0;
        for(int i = 0; i < n; i++) {
            idx *= p[i] + 1;
            idx += data[i];
        }
        return idx;
    };
    auto idx_decode = [&](int idx) {
        vector<int> data(n);
        for(int i = n-1; i >= 0; i--) {
            data[i] = idx % (p[i] + 1);
            idx /= p[i] + 1;
        }
        return data;
    };
    bool winable = false;
    vector<int> win(1<<n);
    for(int i = 0; i < 1 << n; i++) {
        int tmp = 0;
        for(int j = 0; j < n; j++) {
            if(i >> j & 1) tmp += t[j];
            else tmp += s[j];
        }
        if(tmp <= 60) {
            win[i] = 1;
            winable = true;
        }
    }
    if(!winable) {
        cout << "-1\n";
        return 0;
    }
    vector<int> mx(prodp[n]), mn(prodp[n], 1e9);
    vector<f128> ex(prodp[n]);
    for(int i = (prodp[n])-1; i >= 0; i--) {
        vector<int> data = idx_decode(i);
        int numcases = 0;
        int numselfcases = 0;
        f128 sumsteps = 0;
        for(int j = 0; j < 1<<n; j++) {
            int cases = 1;
            vector<int> ndata = data;
            for(int l = 0; l < n; l++) {
                if(j >> l & 1) cases *= max(0, data[l] - k);
                else {
                    cases *= p[l] - data[l];
                    ndata[l]++;
                }
            }
            if(cases == 0) continue;
            int ni = idx_encode(ndata);
            numcases += cases;
            sumsteps += cases;
            if(win[j]) {
                mn[i] = min(mn[i], 1);
                mx[i] = max(mx[i], 1);
            } else if(ni == i) {
                numselfcases += cases;
            } else {
                sumsteps += cases * ex[ni];
                mn[i] = min(mn[i], mn[ni] + 1);
                mx[i] = max(mx[i], mx[ni] + 1);
            }
        }
        ex[i] = sumsteps / (numcases - numselfcases);
    }
    cout << fixed << setprecision(20);
    cout << mn[0] << " " << mx[0] << " " << (double)ex[0] << '\n';
}
0