結果
問題 | No.174 カードゲーム(Hard) |
ユーザー | sugim48 |
提出日時 | 2015-03-27 00:06:13 |
言語 | C++11 (gcc 11.4.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 2,139 bytes |
コンパイル時間 | 2,546 ms |
コンパイル使用メモリ | 98,312 KB |
実行使用メモリ | 4,380 KB |
最終ジャッジ日時 | 2023-09-21 02:21:48 |
合計ジャッジ時間 | 2,634 ms |
ジャッジサーバーID (参考情報) |
judge14 / judge11 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 1 ms
4,380 KB |
testcase_01 | AC | 2 ms
4,380 KB |
testcase_02 | AC | 2 ms
4,380 KB |
testcase_03 | AC | 1 ms
4,380 KB |
testcase_04 | AC | 2 ms
4,380 KB |
testcase_05 | AC | 1 ms
4,376 KB |
testcase_06 | AC | 2 ms
4,376 KB |
testcase_07 | AC | 2 ms
4,376 KB |
testcase_08 | AC | 2 ms
4,380 KB |
testcase_09 | AC | 2 ms
4,376 KB |
testcase_10 | AC | 1 ms
4,376 KB |
testcase_11 | AC | 1 ms
4,380 KB |
ソースコード
#define _USE_MATH_DEFINES #include <algorithm> #include <cstdio> #include <functional> #include <iostream> #include <cfloat> #include <climits> #include <cstring> #include <cmath> #include <map> #include <queue> #include <random> #include <set> #include <sstream> #include <stack> #include <string> #include <time.h> #include <vector> using namespace std; typedef long long ll; typedef unsigned long long ull; typedef pair<int, int> i_i; typedef pair<ll, int> ll_i; typedef pair<double, int> d_i; typedef pair<ll, ll> ll_ll; typedef pair<double, double> d_d; struct edge { int u, v; ll w; }; ll MOD = 1000000007; ll _MOD = 1000000009; double EPS = 1e-10; int INF = INT_MAX / 2; int main() { int N; double PA, PB; cin >> N >> PA >> PB; vector<int> A(N), B(N); for (int i = 0; i < N; i++) cin >> A[i]; for (int i = 0; i < N; i++) cin >> B[i]; sort(A.begin(), A.end()); sort(B.begin(), B.end()); vector<vector<double> > pa(N, vector<double>(N)), pb = pa; for (int i = 0; i < N; i++) { vector<double> dp(N); dp[i] = 1; for (int t = 0; t < N; t++) { if (t == N - 1) { pa[i][t] += dp[0]; continue; } vector<double> _dp(N); pa[i][t] += dp[0] * PA; _dp[0] += dp[0] * (1 - PA); for (int j = 1; j < N - t; j++) { pa[i][t] += dp[j] * (1 - PA) / (N - t - 1); _dp[j - 1] += dp[j] * ((1 - PA) / (N - t - 1) * (j - 1) + PA); _dp[j] += dp[j] * (1 - PA) / (N - t - 1) * (N - t - j - 1); } dp = _dp; } } for (int i = 0; i < N; i++) { vector<double> dp(N); dp[i] = 1; for (int t = 0; t < N; t++) { if (t == N - 1) { pb[i][t] += dp[0]; continue; } vector<double> _dp(N); pb[i][t] += dp[0] * PB; _dp[0] += dp[0] * (1 - PB); for (int j = 1; j < N - t; j++) { pb[i][t] += dp[j] * (1 - PB) / (N - t - 1); _dp[j - 1] += dp[j] * ((1 - PB) / (N - t - 1) * (j - 1) + PB); _dp[j] += dp[j] * (1 - PB) / (N - t - 1) * (N - t - j - 1); } dp = _dp; } } double e = 0; for (int i = 0; i < N; i++) for (int j = 0; j < N; j++) if (A[i] > B[j]) for (int t = 0; t < N; t++) e += (A[i] + B[j]) * pa[i][t] * pb[j][t]; printf("%.10f\n", e); }