結果

問題 No.108 トリプルカードコンプ
ユーザー tsutaj
提出日時 2018-02-19 03:22:32
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 10 ms / 5,000 ms
コード長 1,973 bytes
コンパイル時間 711 ms
コンパイル使用メモリ 100,036 KB
実行使用メモリ 12,800 KB
最終ジャッジ日時 2024-07-05 20:33:13
合計ジャッジ時間 1,564 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

diff # 

// 基本テンプレート
 
#include <iostream>
#include <iomanip>
#include <cstdio>
#include <string>
#include <cstring>
#include <deque>
#include <list>
#include <queue>
#include <stack>
#include <vector>
#include <utility>
#include <algorithm>
#include <map>
#include <set>
#include <complex>
#include <cmath>
#include <limits>
#include <cfloat>
#include <climits>
#include <ctime>
#include <cassert>
#include <numeric>
#include <fstream>
#include <functional>
using namespace std;
 
#define rep(i,a,n) for(int (i)=(a); (i)<(n); (i)++)
#define repq(i,a,n) for(int (i)=(a); (i)<=(n); (i)++)
#define repr(i,a,n) for(int (i)=(a); (i)>=(n); (i)--)
#define int long long int
 
template<typename T> void chmax(T &a, T b) {a = max(a, b);}
template<typename T> void chmin(T &a, T b) {a = min(a, b);}
template<typename T> void chadd(T &a, T b) {a = a + b;}
 
typedef pair<int, int> pii;
typedef long long ll;
 
int dx[] = {0, 0, 1, -1};
int dy[] = {1, -1, 0, 0};
const ll INF = 1001001001001001LL;
const ll MOD = 1000000007LL;

int N;
double dp[110][110][110];
int cnt[4];

double solve(int x, int y, int z) {
    // メモ化しよう
    if(dp[x][y][z] != -1) return dp[x][y][z];
    dp[x][y][z] = 0;
    // 2 枚以下しか持っていないものが出るまで引く回数の期待値
    double exp_num = 1.0 * N / (N - z);
    // (先の期待値) * (それを引く確率) の足しあわせ
    if(x != N) dp[x][y][z] += (solve(x+1, y, z) + exp_num) * (N - x) / (N - z);
    if(y != x) dp[x][y][z] += (solve(x, y+1, z) + exp_num) * (x - y) / (N - z);
    if(z != y) dp[x][y][z] += (solve(x, y, z+1) + exp_num) * (y - z) / (N - z);
    return dp[x][y][z];
}
 
signed main() {
    cin >> N;
    repq(i,0,N) repq(j,0,N) repq(k,0,N) dp[i][j][k] = -1;
    dp[N][N][N] = 0;

    rep(i,0,N) {
        int A; cin >> A;
        cnt[min(A, 3LL)]++;
    }
    repr(i,2,0) {
        cnt[i] += cnt[i+1];
    }

    printf("%.12f\n", solve(cnt[1], cnt[2], cnt[3]));
    return 0;
}
0