結果

問題 No.108 トリプルカードコンプ
ユーザー auaua
提出日時 2020-02-08 23:19:36
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 1,141 bytes
コンパイル時間 1,417 ms
コンパイル使用メモリ 167,136 KB
実行使用メモリ 25,088 KB
最終ジャッジ日時 2024-10-01 05:30:12
合計ジャッジ時間 2,740 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1 WA * 2
other AC * 4 WA * 16
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
#include<unordered_set>
#include<unordered_map>
using namespace std;
#define REP(i,m,n) for(int i=(m);i<(n);i++)
#define rep(i,n) REP(i,0,n)
#define pb push_back
#define all(a) a.begin(),a.end()
#define rall(c) (c).rbegin(),(c).rend()
#define mp make_pair
#define double long double
typedef long long ll;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
const ll inf=1e9+7;
const ll mod=1e9+7;
ll n;
vector<vector<vector<double> > >dp(110,vector<vector<double> >(110,vector<double>(110)));
vector<ll>cnt(4);
double solve(ll x,ll y,ll z){
    if(dp[x][y][z]!=-1)return dp[x][y][z];
    dp[x][y][z]=0;
    double k=1.0*n/(n-z);
    if(x!=n)dp[x][y][z]+=(solve(x+1,y,z)+k)*(n-x)/(n-z);
    if(x!=y)dp[x][y][z]+=(solve(x,y+1,z)+k)*(x-y)/(n-z);
    if(y!=z)dp[x][y][z]+=(solve(x,y,z+1)+k)*(z-y)/(n-z);
    return dp[x][y][z];
}
int main(){
    cin>>n;
    rep(i,n+1)rep(j,n+1)rep(k,n+1)dp[i][j][k]=-1;
    rep(i,n){
        ll a;cin>>a;
        cnt[min(a,3LL)]++;
    }
    for(int i=2;i>=0;i--){
        cnt[i]+=cnt[i+1];
    }
    dp[n][n][n]=0;
    cout<<fixed<<setprecision(10)<<solve(cnt[1],cnt[2],cnt[3])<<endl;
}
0