結果

問題 No.174 カードゲーム（Hard）
ユーザー ei1333333ei1333333
提出日時 2019-02-05 23:06:04
言語 C++17
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 328 ms / 2,000 ms
コード長 1,589 bytes
コンパイル時間 2,262 ms
コンパイル使用メモリ 207,692 KB
実行使用メモリ 19,856 KB
最終ジャッジ日時 2023-08-30 23:58:35
合計ジャッジ時間 5,937 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,376 KB
testcase_01 AC 2 ms
4,384 KB
testcase_02 AC 323 ms
19,644 KB
testcase_03 AC 322 ms
19,664 KB
testcase_04 AC 321 ms
19,760 KB
testcase_05 AC 323 ms
19,756 KB
testcase_06 AC 317 ms
19,856 KB
testcase_07 AC 323 ms
19,596 KB
testcase_08 AC 328 ms
19,764 KB
testcase_09 AC 324 ms
19,668 KB
testcase_10 AC 1 ms
4,376 KB
testcase_11 AC 1 ms
4,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

int main() {
  int N, A[20], B[20];
  double X, Y;

  cin >> N >> X >> Y;
  for(int i = 0; i < N; i++) cin >> A[i];
  for(int i = 0; i < N; i++) cin >> B[i];
  sort(A, A + N);
  sort(B, B + N);

  vector< double > dp1(1 << N), dp2(1 << N);
  vector< vector< double > > dp3(N + 1, vector< double >(N)), dp4(N + 1, vector< double >(N));
  dp1.back() = dp2.back() = 1.0;


  for(int i = (1 << N) - 1; i > 0; i--) {
    vector< int > enable;
    for(int j = 0; j < N; j++) {
      if((i >> j) & 1) enable.emplace_back(j);
    }
    if(enable.size() == 1) {
      dp1[0] += dp1[i];
      dp2[0] += dp2[i];
      dp3[enable.size()][enable[0]] += dp1[i];
      dp4[enable.size()][enable[0]] += dp2[i];
    } else {
      dp1[i & ~(1 << enable[0])] += dp1[i] * X;
      dp2[i & ~(1 << enable[0])] += dp2[i] * Y;
      dp3[enable.size()][enable[0]] += dp1[i] * X;
      dp4[enable.size()][enable[0]] += dp2[i] * Y;
      double Z = 1.0 / ((int) enable.size() - 1);
      for(int j = 1; j < enable.size(); j++) {
        dp1[i & ~(1 << enable[j])] += dp1[i] * Z * (1.0 - X);
        dp2[i & ~(1 << enable[j])] += dp2[i] * Z * (1.0 - Y);
        dp3[enable.size()][enable[j]] += dp1[i] * Z * (1.0 - X);
        dp4[enable.size()][enable[j]] += dp2[i] * Z * (1.0 - Y);
      }
    }
  }

  double ret = 0.0;
  for(int i = 1; i <= N; i++) {
    for(int j = 0; j < N; j++) {
      for(int k = 0; k < N; k++) {
        if(A[j] > B[k]) ret += dp3[i][j] * dp4[i][k] * (A[j] + B[k]);
      }
    }
  }

  cout << fixed << setprecision(10) << ret << endl;
}
0