結果
問題 | No.425 ジャンケンの必勝法 |
ユーザー |
|
提出日時 | 2016-12-16 20:14:59 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 110 ms / 2,000 ms |
コード長 | 1,740 bytes |
コンパイル時間 | 1,611 ms |
コンパイル使用メモリ | 167,216 KB |
実行使用メモリ | 14,052 KB |
最終ジャッジ日時 | 2024-11-30 10:38:28 |
合計ジャッジ時間 | 4,755 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 5 |
other | AC * 18 |
ソースコード
#include <bits/stdc++.h>using namespace std;typedef long long ll;typedef vector< int > vi;typedef vector< vi > vvi;typedef vector< ll > vl;typedef vector< vl > vvl;typedef pair< int, int > pii;typedef vector< pii > vp;typedef vector< double > vd;typedef vector< vd > vvd;typedef vector< string > vs;template< class T1, class T2 >int upmin( T1 &x, T2 v ){if( x > v ){x = v;return 1;}return 0;}template< class T1, class T2 >int upmax( T1 &x, T2 v ){if( x < v ){x = v;return 1;}return 0;}const int INF = 0x3f3f3f3f;double P0, Q;void init(){cin >> P0 >> Q;}const int T = 20;double dp[ T ][ 1 << T ];void preprocess(){P0 /= 100.0;Q /= 100.0;dp[ 1 ][ 0 ] = 1.0 / 3.0;for( int i = 1; i + 1 < T; ++i )for( int j = 0; j < 1 << i; ++j ){double skill = P0;for( int k = 1; k < i; ++k ){if( j & 1 << k )skill = min( 1.0, max( 0.0, skill - Q ) );elseskill = min( 1.0, max( 0.0, skill + Q ) );}dp[ i + 1 ][ j | 1 << i ] += dp[ i ][ j ] * skill * 0.5;dp[ i + 1 ][ j ] += dp[ i ][ j ] * ( 1.0 - skill ) / 3.0;}}void solve(){double ans = 1.0 / 3.0;for( int i = 1; i < T; ++i )for( int j = 0; j < 1 << i; ++j ){double skill = P0;for( int k = 1; k < i; ++k ){if( j & 1 << k )skill = min( 1.0, max( 0.0, skill - Q ) );elseskill = min( 1.0, max( 0.0, skill + Q ) );}ans += dp[ i ][ j ] * skill * 0.5;ans += dp[ i ][ j ] * ( 1.0 - skill ) / 3.0;}cout << fixed << setprecision( 9 ) << ans << endl;}signed main(){ios::sync_with_stdio( 0 );init();preprocess();solve();return 0;}