結果
問題 | No.425 ジャンケンの必勝法 |
ユーザー |
|
提出日時 | 2016-09-24 19:23:26 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 30 ms / 2,000 ms |
コード長 | 979 bytes |
コンパイル時間 | 1,559 ms |
コンパイル使用メモリ | 167,308 KB |
実行使用メモリ | 6,824 KB |
最終ジャッジ日時 | 2024-11-18 09:39:00 |
合計ジャッジ時間 | 2,909 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 5 |
other | AC * 18 |
ソースコード
#define _CRT_SECURE_NO_WARNINGS#define _USE_MATH_DEFINES#include "bits/stdc++.h"#define REP(i,a,b) for(int i=a;i<b;++i)#define rep(i,n) REP(i,0,n)#define ll long long#define ull unsigned lltypedef long double ld;#define ALL(a) begin(a),end(a)#define ifnot(a) if(not a)#define dump(x) cerr << #x << " = " << (x) << endlusing namespace std;// #define int ll#ifdef _MSC_VERconst bool test = true;#elseconst bool test = false;#endifint dx[] = { 0,1,0,-1 };int dy[] = { 1,0,-1,0 };#define INF (1 << 28)ull mod = (int)1e9 + 7;//.....................#define MAX (int)1e6 + 5double p0, q;double dfs(double a, double p, int cnt) {if (cnt == 22) return 0;double res = 0;res += a * p / 2;res += dfs(a * p / 2, max(p - q, 0.0), cnt + 1);res += a * (1 - p) / 3;res += dfs(a * (1 - p) / 3, min(p + q, 1.0), cnt + 1);return res;}signed main() {cin >> p0 >> q;p0 /= 100; q /= 100;printf("%.20lf\n", 1.0/3 + dfs(1.0/3, p0, 0));return 0;}