結果

問題 No.174 カードゲーム（Hard）
コンテスト
ユーザー kazuma
提出日時 2018-12-28 13:40:45
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 328 ms / 2,000 ms
コード長 1,674 bytes
コンパイル時間 2,191 ms
コンパイル使用メモリ 204,464 KB
最終ジャッジ日時 2025-01-06 20:05:19
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 12
権限があれば一括ダウンロードができます

ソースコード

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