結果

問題 No.108 トリプルカードコンプ
ユーザー たこしたこし
提出日時 2019-10-01 22:09:53
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 1,793 bytes
コンパイル時間 2,084 ms
コンパイル使用メモリ 195,696 KB
実行使用メモリ 13,004 KB
最終ジャッジ日時 2024-04-14 08:40:39
合計ジャッジ時間 3,324 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 5 ms
12,820 KB
testcase_01 AC 4 ms
12,892 KB
testcase_02 AC 6 ms
12,740 KB
testcase_03 AC 6 ms
12,592 KB
testcase_04 AC 6 ms
12,928 KB
testcase_05 AC 5 ms
12,604 KB
testcase_06 AC 6 ms
12,892 KB
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 AC 5 ms
12,724 KB
testcase_11 AC 5 ms
12,832 KB
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 WA -
testcase_20 WA -
testcase_21 WA -
testcase_22 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

#define INF 100000000
#define YJ 1145141919
#define INF_INT_MAX 2147483647
#define INF_LL 9223372036854775
#define INF_LL_MAX 9223372036854775807
#define EPS 1e-10
#define MOD 1000000007
#define MOD9 998244353
#define Pi acos(-1)
#define LL long long
#define ULL unsigned long long
#define LD long double

#define int long long

using II = pair<int, int>;

int gcd(int a, int b) { return b != 0 ? gcd(b, a % b) : a; }
int lcm(int a, int b) { return a * b / gcd(a, b); }
int extgcd(int a, int b, int &x, int &y) { int g = a; x = 1; y = 0; if (b != 0) g = extgcd(b, a % b, y, x), y -= (a / b) * x; return g; }

#define FOR(i,a,b) for(int i=(a);i<(b);++i)
#define REP(i,n) FOR(i,0,n)
#define ALL(a)  begin((a)), end((a))
#define RALL(a) (a).rbegin(), (a).rend()
#define PB push_back
#define MP make_pair
#define SZ(a) int((a).size())

const int MAX_N = 105;
int N;
int A[MAX_N];
int sumA[3];

double dp[MAX_N][MAX_N][MAX_N];

double solve(int i, int j, int k) {
  if(i == sumA[0] && j == sumA[1] & k == sumA[2]) {
    return 0.0;
  }

  if(i == 0 && j == 0 && k == 0) {
    return (solve(i,j,k+1) + 1) / N;
  }

  if(dp[i][j][k] >= 0.0) {
    return dp[i][j][k];
  }

  double ret = 0.0;

  if(j-1 >= 0) {
    ret += solve(i+1,j-1,k) * (i+1) + (i+1);
  }
  if(k-1 >= 0) {
    ret += solve(i,j+1,k-1) * (j+1) + (j+1);
  }
  if(i+j+k+1 <= sumA[0] + sumA[1] + sumA[2]) {
    ret += solve(i,j,k+1) * (k+1) + (k+1);
  }

  return dp[i][j][k] = ret / (i+j+k);
}

signed main()
{
  cin >> N;
  REP(n,N) {
    cin >> A[n];
    if(A[n] >= 3) {
      continue;
    } else {
      sumA[A[n]]++;
    }
  }

  REP(i,MAX_N) {
    REP(j,MAX_N) {
      REP(k,MAX_N) {
        dp[i][j][k] = -1.0;
      }
    }
  }

  printf("%.15f\n", solve(0,0,0));

  return 0;
}
0