結果

問題 No.108 トリプルカードコンプ
ユーザー auauaauaua
提出日時 2020-02-08 23:25:19
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 24 ms / 5,000 ms
コード長 1,096 bytes
コンパイル時間 1,701 ms
コンパイル使用メモリ 176,964 KB
実行使用メモリ 24,960 KB
最終ジャッジ日時 2024-10-01 05:30:41
合計ジャッジ時間 2,884 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

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,-1)));
vector<ll>cnt(4,0);
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(y!=x)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)*(y-z)/(n-z);
    return dp[x][y][z];
}
int main(){
    cin>>n;
    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