結果

問題 No.174 カードゲーム（Hard）
ユーザー kazumakazuma
提出日時 2018-12-28 13:40:45
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 309 ms / 2,000 ms
コード長 1,674 bytes
コンパイル時間 2,053 ms
コンパイル使用メモリ 211,780 KB
実行使用メモリ 20,188 KB
最終ジャッジ日時 2024-10-01 14:59:06
合計ジャッジ時間 5,093 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 1 ms
6,816 KB
testcase_02 AC 302 ms
19,992 KB
testcase_03 AC 308 ms
19,912 KB
testcase_04 AC 309 ms
20,076 KB
testcase_05 AC 294 ms
19,888 KB
testcase_06 AC 295 ms
20,100 KB
testcase_07 AC 296 ms
20,180 KB
testcase_08 AC 291 ms
20,100 KB
testcase_09 AC 288 ms
20,188 KB
testcase_10 AC 1 ms
5,248 KB
testcase_11 AC 1 ms
5,248 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using ld = long double;

vector<vector<ld>> calc(int N, ld P) {
	vector<ld> dp(1 << N);
	vector<vector<ld>> p(N, vector<ld>(N));
	dp[0] = 1;
	for (int i = 0; i < (1 << N) - 1; i++) {
		int id = 0, cnt = __builtin_popcount(i);
		for (int j = 0; j < N; j++) if (((i >> j) & 1) == 0) {
			id = j;
			break;
		}
		if (cnt + 1 == N) {
			p[cnt][id] += dp[i];
			dp[i | (1 << id)] += dp[i];
			continue;
		}
		p[cnt][id] += P * dp[i];
		dp[i | (1 << id)] += P * dp[i];
		ld pp = (1 - P) / (N - (cnt + 1));
		for (int j = 0; j < N; j++) if (((i >> j) & 1) == 0 && j != id) {
			p[cnt][j] += pp * dp[i];
			dp[i | (1 << j)] += pp * dp[i];
		}
	}
	return p;
}

int main()
{
	cout << fixed << setprecision(10);
	int N;
	ld 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());
	auto p = calc(N, Pb);
	vector<vector<ld>> bsum(N, vector<ld>(N + 1));
	vector<vector<ld>> psum(N, vector<ld>(N + 1));
	for (int i = 0; i < N; i++) {
		for (int j = 0; j < N; j++) {
			bsum[i][j + 1] = bsum[i][j] + p[i][j] * B[j];
			psum[i][j + 1] = psum[i][j] + p[i][j];
		}
	}
	vector<ld> dp(1 << N);
	for (int i = 1; i < (1 << N); i++) {
		int cnt = N - __builtin_popcount(i);
		ld p1 = cnt == N - 1 ? 1 : Pa, p2 = cnt == N - 1 ? 0 : (1 - Pa) / (N - 1 - cnt);
		for (int j = 0, k = 0, c = 0; j < N; j++) if (i & (1 << j)) {
			while (k < N && B[k] < A[j]) ++k;
			dp[i] += (dp[i ^ (1 << j)] + bsum[cnt][k] + A[j] * psum[cnt][k]) * (!c ? p1 : p2);
			++c;
		}
	}
	cout << dp.back() << endl;
	return 0;
}
0