結果
問題 |
No.108 トリプルカードコンプ
|
ユーザー |
|
提出日時 | 2024-04-05 22:53:08 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 26 ms / 5,000 ms |
コード長 | 2,213 bytes |
コンパイル時間 | 7,318 ms |
コンパイル使用メモリ | 350,240 KB |
実行使用メモリ | 24,576 KB |
最終ジャッジ日時 | 2024-10-01 02:55:27 |
合計ジャッジ時間 | 8,831 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 20 |
ソースコード
#pragma GCC target("avx2") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #include<bits/stdc++.h> #include <atcoder/all> using namespace std;using namespace atcoder;using vst = vector<string>;using ll = long long; using ld = long double;using pll = pair<ll,ll>;using vll = vector<ll>;using vvll = vector<vll>;using vpl = vector<pll>; #define V(a) vector<a> #define rep(i, n) for (ll i = 0; i < n; i++) #define repp(i,k,n) for (ll i = k; i < n; i++) #define per(i,s,e) for(ll i = s; i >= e; i--) #define all(v) v.begin(),v.end() #define yesno(a) a ? cout << "Yes" << '\n' : cout << "No" << '\n' #define YESNO(a) a ? cout << "YES" << '\n' : cout << "NO" << '\n' #define UNOmap unordered_map #define UNOset unordered_set #define join(G,a,b) G[a].push_back(b);G[b].push_back(a); #define PB(a,b) a.emplace_back(b) template<class... T>void in(T&... a){(cin >> ... >> a);} template<class T, class... Ts>void out(const T& a, const Ts&... b){cout << a;((cout << ' ' << b), ...);cout << '\n';} template<class T> void vin(vector<T> &v){for(ll i = 0; i < (ll)v.size(); i++)in(v[i]);} template<class T> void vout(vector<T> &v){for(ll i = 0; i < (ll)v.size(); i++) cout << v[i] << ' ';cout << "\n";} ll INF = 4611686018427387904;ll MOD9 =998244353; ll MOD1 =1000000007; /* INF = 2^62 */ // #include<boost/multiprecision/cpp_int.hpp> // using namespace boost::multiprecision; using mint9 = modint998244353; using mint1 = modint1000000007; using mint = modint; //---------------------------------------------- ll N; vector dp(109,vector(109,vector(109,ld(-1)))); ld dfs(ll x,ll y, ll z){ if(dp[x][y][z] != -1) return dp[x][y][z]; dp[x][y][z] = 0; ld e = 1.0 * N /(N-z); if(x != N) dp[x][y][z] += (dfs(x+1,y,z) + e)*(N-x)/(N-z); if(y != x) dp[x][y][z] += (dfs(x,y+1,z) + e)*(x-y)/(N-z); if(z != y) dp[x][y][z] += (dfs(x,y,z+1) + e)*(y-z)/(N-z); return dp[x][y][z]; } int main(){ ios::sync_with_stdio(false); std::cin.tie(nullptr); cout << fixed << setprecision(15); //============================================== in(N); dp[N][N][N] = 0; vll cnt(4); rep(i,N){ ll a; in(a); rep(j,min(3LL,a)+1){ cnt[j]++; } } out(dfs(cnt[1],cnt[2],cnt[3])); }