結果

問題 No.425 ジャンケンの必勝法
ユーザー koyoprokoyopro
提出日時 2020-02-25 00:22:56
言語 C++11
(gcc 11.4.0)
結果
WA  
実行時間 -
コード長 2,172 bytes
コンパイル時間 1,122 ms
コンパイル使用メモリ 157,804 KB
実行使用メモリ 23,168 KB
最終ジャッジ日時 2024-04-20 21:19:18
合計ジャッジ時間 2,686 ms
ジャッジサーバーID
(参考情報)
judge1 / judge5
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 32 ms
22,912 KB
testcase_01 AC 35 ms
23,040 KB
testcase_02 AC 32 ms
23,040 KB
testcase_03 WA -
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 AC 32 ms
23,168 KB
testcase_11 AC 32 ms
22,912 KB
testcase_12 AC 33 ms
22,912 KB
testcase_13 WA -
testcase_14 AC 32 ms
23,168 KB
testcase_15 AC 33 ms
22,912 KB
testcase_16 AC 32 ms
22,912 KB
testcase_17 AC 32 ms
23,040 KB
testcase_18 WA -
testcase_19 WA -
testcase_20 AC 33 ms
22,912 KB
testcase_21 AC 32 ms
23,040 KB
testcase_22 AC 36 ms
23,040 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include "bits/stdc++.h"
using namespace std;
#define int long long
#define FOR(i, a, b) for(int i=(a);i<(b);i++)
#define RFOR(i, a, b) for(int i=(b-1);i>=(a);i--)
#define REP(i, n) for(int i=0; i<(n); i++)
#define RREP(i, n) for(int i=(n-1); i>=0; i--)
#define REP1(i, n) for(int i=1; i<=(n); i++)
#define RREP1(i, n) for(int i=(n); i>=1; i--)
#define ALL(a) (a).begin(),(a).end()
#define UNIQUE_SORT(l) sort(ALL(l)); l.erase(unique(ALL(l)), l.end());
#define CONTAIN(a, b) find(ALL(a), (b)) != (a).end()
#define out(...) printf(__VA_ARGS__)

int dxy[] = {0, 1, 0, -1, 0};

void solve();
signed main()
{
#if DEBUG
    std::ifstream in("input.txt");
    std::cin.rdbuf(in.rdbuf());
#endif
    cin.tie(0);
    ios::sync_with_stdio(false);
    solve();
    return 0;
}

/*================================*/
#if DEBUG
#define SIZE 1234560
#else
#define SIZE 1234560
#endif

double P,Q;
double ans = 0;
double EP[SIZE][2]; // i回目であいこの確率(直前に必勝法を使っていない/使った)
double PP; // Pprev

void solve() {
    cin>>P>>Q;
    P /= 100, Q /= 100;
    
    PP = -1;
    REP(i,SIZE) {
        if (!i) {
            ans += 1.0/3; // win
            EP[i][0] = 1.0/3; // even
            EP[i][1] = 0;
            continue;
        }
        double p0, p1;
        {
            p0 = (PP < 0) ? P : min(1.0, PP + Q);
            // 使わない→使う
            double puse = EP[i-1][0] * p0;
            ans += puse * 0.5; // win
            EP[i][1] += puse * 0.5; // even
            // 使わない→使わない
            double pfree = EP[i-1][0] * (1-p0);
            ans += pfree / 3; // win
            EP[i][0] += pfree / 3; // even
        }
        {
            // 使う→使う
            p1 = max(0.0, PP - Q);
            double puse = EP[i-1][1] * p1;
            ans += puse * 0.5; // win
            EP[i][1] += puse * 0.5; // even
            // 使う→使わない
            double pfree = EP[i-1][1] * (1-p1);
            ans += pfree / 3; // win
            EP[i][0] += pfree / 3; // even
        }
        PP = (EP[i-1][0] * p0 + EP[i-1][1] * p1) / (EP[i-1][0] + EP[i-1][1]);
    }
    out("%.14f\n", ans);
}
0