結果
問題 | No.162 8020運動 |
ユーザー | koyumeishi |
提出日時 | 2015-03-06 01:07:36 |
言語 | C++11 (gcc 11.4.0) |
結果 |
AC
|
実行時間 | 124 ms / 5,000 ms |
コード長 | 1,692 bytes |
コンパイル時間 | 589 ms |
コンパイル使用メモリ | 78,548 KB |
実行使用メモリ | 4,384 KB |
最終ジャッジ日時 | 2023-09-06 15:27:57 |
合計ジャッジ時間 | 3,956 ms |
ジャッジサーバーID (参考情報) |
judge13 / judge15 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 5 ms
4,376 KB |
testcase_01 | AC | 60 ms
4,380 KB |
testcase_02 | AC | 94 ms
4,380 KB |
testcase_03 | AC | 122 ms
4,380 KB |
testcase_04 | AC | 123 ms
4,380 KB |
testcase_05 | AC | 124 ms
4,376 KB |
testcase_06 | AC | 122 ms
4,380 KB |
testcase_07 | AC | 121 ms
4,376 KB |
testcase_08 | AC | 121 ms
4,380 KB |
testcase_09 | AC | 5 ms
4,376 KB |
testcase_10 | AC | 92 ms
4,376 KB |
testcase_11 | AC | 60 ms
4,384 KB |
testcase_12 | AC | 30 ms
4,376 KB |
testcase_13 | AC | 123 ms
4,376 KB |
testcase_14 | AC | 11 ms
4,380 KB |
testcase_15 | AC | 11 ms
4,376 KB |
testcase_16 | AC | 36 ms
4,380 KB |
testcase_17 | AC | 18 ms
4,380 KB |
testcase_18 | AC | 115 ms
4,376 KB |
testcase_19 | AC | 85 ms
4,380 KB |
testcase_20 | AC | 97 ms
4,380 KB |
testcase_21 | AC | 98 ms
4,380 KB |
testcase_22 | AC | 91 ms
4,376 KB |
testcase_23 | AC | 98 ms
4,376 KB |
testcase_24 | AC | 90 ms
4,376 KB |
testcase_25 | AC | 104 ms
4,376 KB |
testcase_26 | AC | 5 ms
4,376 KB |
testcase_27 | AC | 67 ms
4,380 KB |
testcase_28 | AC | 115 ms
4,376 KB |
ソースコード
#include <iostream> #include <vector> #include <cstdio> #include <sstream> #include <map> #include <string> #include <algorithm> #include <queue> #include <cmath> #include <set> using namespace std; // probability s -> t double prob(vector<int>& P, int s, int t, int len){ s <<= 1; t <<= 1; if( (s|t) > s ) return 0.0; double ret = 1.0; int mask = 0b111; for(int i=1; i<=len; i++){ if( (s&(1<<i)) == 0 ) continue; int my_s = ( s>>(i-1) ) & mask; if( my_s == 0b010 ){ if( (t&(1<<i)) == 0){ ret *= P[0] * 0.01; }else{ ret *= 1 - P[0] * 0.01; } }else if( my_s == 0b111 ){ if( (t&(1<<i)) == 0){ ret *= P[2] * 0.01; }else{ ret *= 1 - P[2] * 0.01; } }else{ if( (t&(1<<i)) == 0){ ret *= P[1] * 0.01; }else{ ret *= 1 - P[1] * 0.01; } } } return ret; } double dfs(vector<int>& P, vector<vector<double>>& memo, int s, int len, int y){ if(memo[y][len] >= 0) return memo[y][len]; if(y==0){ return len; } //部分集合列挙 double ret = 0.0; int sub = s; do{ double p = prob(P, s, sub, len); double E_ = 0.0; int last_len = 0; for(int i=0; i<len; i++){ if((sub >> i) & 1){ last_len++; }else{ E_ += dfs(P, memo, (1<<last_len)-1, last_len, y-1); last_len = 0; } } if(last_len > 0){ E_ += dfs(P, memo, (1<<last_len)-1, last_len, y-1); last_len = 0; } ret += E_ * p; sub = (sub-1) & s; }while( sub != s ); memo[y][len] = ret; return ret; } int main(){ int A; cin >> A; vector<int> P(3); cin >> P[0] >> P[1] >> P[2]; vector<vector<double>> memo(25, vector<double>(15, -1)); double ans = dfs(P,memo, (1<<14)-1, 14, 80-A) * 2; printf("%.10f\n", ans); return 0; }