結果

問題 No.108 トリプルカードコンプ
コンテスト
ユーザー moti
提出日時 2015-07-25 12:32:24
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 9 ms / 5,000 ms
コード長 1,050 bytes
コンパイル時間 599 ms
コンパイル使用メモリ 85,044 KB
実行使用メモリ 14,448 KB
最終ジャッジ日時 2024-07-16 04:52:40
合計ジャッジ時間 1,419 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <iostream>
#include <algorithm>
#include <cmath>
#include <vector>
#include <complex>
#include <queue>
#include <set>
#include <iomanip>
#include <cstring>

using namespace std;

#define REP(i,a,b) for(int i=a;i<(int)b;i++)
#define rep(i,n) REP(i,0,n)

typedef long long ll;

int N;
double dp[111][111][111];

double solve(int i, int j, int k) {

  if(i < 0 || j < 0 || k < 0) { return 0.0; }

  double& ret = dp[i][j][k];
  if(i == 0 && j == 0 && k == 0) { return ret = 1.0; }

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

  ret =  (double)N / (i + j + k); // 自己ループ

  ret += solve(i-1, j+1, k)  * i / (i + j + k);
  ret += solve(i, j-1, k+1)  * j / (i + j + k);
  ret += solve(i, j, k-1)    * k / (i + j + k);

  return ret;
}

int main() {

  memset(dp, -1, sizeof dp);

  cin >> N;
  int cnt0 = 0, cnt1 = 0, cnt2 = 0;
  rep(i, N) {
    int a; cin >> a;
    if(a == 0) { cnt0 ++; }
    if(a == 1) { cnt1 ++; }
    if(a == 2) { cnt2 ++; }
  }

  cout << setprecision(7) << solve(cnt0, cnt1, cnt2) - 1.0 << endl;

  return 0;
}
0