結果

問題 No.108 トリプルカードコンプ
ユーザー oxyshoweroxyshower
提出日時 2020-01-20 16:03:04
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 16 ms / 5,000 ms
コード長 1,271 bytes
コンパイル時間 2,127 ms
コンパイル使用メモリ 178,492 KB
実行使用メモリ 11,904 KB
最終ジャッジ日時 2024-07-02 22:32:33
合計ジャッジ時間 2,863 ms
ジャッジサーバーID
(参考情報)
judge4 / judge2
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 2 ms
6,940 KB
testcase_02 AC 2 ms
6,940 KB
testcase_03 AC 2 ms
6,940 KB
testcase_04 AC 2 ms
6,940 KB
testcase_05 AC 2 ms
6,944 KB
testcase_06 AC 2 ms
6,940 KB
testcase_07 AC 16 ms
11,904 KB
testcase_08 AC 10 ms
11,904 KB
testcase_09 AC 11 ms
11,904 KB
testcase_10 AC 11 ms
11,776 KB
testcase_11 AC 11 ms
11,776 KB
testcase_12 AC 2 ms
6,944 KB
testcase_13 AC 4 ms
6,944 KB
testcase_14 AC 3 ms
6,940 KB
testcase_15 AC 4 ms
6,944 KB
testcase_16 AC 3 ms
6,940 KB
testcase_17 AC 3 ms
6,944 KB
testcase_18 AC 16 ms
11,008 KB
testcase_19 AC 11 ms
10,496 KB
testcase_20 AC 11 ms
11,776 KB
testcase_21 AC 13 ms
11,648 KB
testcase_22 AC 10 ms
10,880 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include<bits/stdc++.h>
using namespace std;
#ifdef LOCAL_DEBUG
  #include "LOCAL_DEBUG.hpp"
#endif
#define int long long

template<class T> vector<T> make_vec(size_t a) { return vector<T>(a); }
template<class T, class... Ts> auto make_vec(size_t a, Ts... ts) {
  return vector<decltype(make_vec<T>(ts...))>(a, make_vec<T>(ts...));
}
template<class T, class V>
typename enable_if<is_class<T>::value == 0>::type fill(T &t, const V &v) {
    t = v;
}
template<class T, class V>
typename enable_if<is_class<T>::value != 0>::type fill(T &t, const V &v){
    for (auto &e : t) fill(e, v);
}
// auto v = make_vec<int>(h, w);
// fill(v, 0);

signed main(){

  int n; cin >> n;
  vector<int> a(n), cnt(4, 0);
  for(int i = 0; i < n; i++){
    cin >> a[i];
    cnt[min(3LL, a[i])]++;
  }

  auto dp = make_vec<double>(n+1, n+1, n+1);
  fill(dp, -1);
  dp[0][0][0] = 0;
  function< double(int,int,int) > rec =
  [&](int a, int b, int c){
    if(dp[a][b][c] != -1) return dp[a][b][c];
    double E = 1.0 * n / (a + b + c);
    if(a > 0) E += rec(a-1, b+1, c) * a / (a + b + c);
    if(b > 0) E += rec(a, b-1, c+1) * b / (a + b + c);
    if(c > 0) E += rec(a, b, c-1) * c / (a + b + c);
    return dp[a][b][c] = E;
  };
  printf("%.9f\n",rec(cnt[0], cnt[1], cnt[2]));

  return 0;
}
0