結果
問題 |
No.108 トリプルカードコンプ
|
ユーザー |
![]() |
提出日時 | 2016-01-04 15:26:59 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 7 ms / 5,000 ms |
コード長 | 1,502 bytes |
コンパイル時間 | 709 ms |
コンパイル使用メモリ | 86,528 KB |
実行使用メモリ | 11,940 KB |
最終ジャッジ日時 | 2024-09-19 11:10:49 |
合計ジャッジ時間 | 1,474 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 20 |
ソースコード
#include <iostream> #include <iomanip> #include <vector> #include <list> #include <deque> #include <queue> #include <stack> #include <set> #include <map> #include <algorithm> #include <cmath> #include <cstring> #include <limits> using namespace std; #define FOR(x,y) for(int x = 0;x < (y);x++) #define LLI long long int #define FORR(x,arr) for(auto& x:arr) #define ALL(a) (a.begin()),(a.end()) #define _L(x) cout<<(x)<<endl int n; double dp[101][101][101]; double solve(int a, int b, int c) { if(dp[a][b][c] >= 0) { return dp[a][b][c]; } double div = a + b + c; double expect = (double)n / div; if(c > 0) { expect += (double)c / div * solve(a,b,c-1); } if(b > 0) { expect += (double)b / div * solve(a,b-1,c+1); } if(a > 0) { expect += (double)a / div * solve(a-1,b+1,c); } dp[a][b][c] = expect; //cout << a << "," << b << "," << c << "," << expect << endl; return expect; } int main() { cin >> n; vector<int> data(n); int a = 0; int b = 0; int c = 0; FOR(i, n) { cin >> data[i]; if(data[i] == 0) { a++; } else if(data[i] == 1) { b++; } else if(data[i] == 2) { c++; } } FOR(i, 101) FOR(j, 101) FOR(k, 101) { dp[i][j][k] = -1; } dp[0][0][0] = 0; cout << fixed << setprecision(12) << solve(a,b,c) << endl; return 0; }