結果

問題 No.108 トリプルカードコンプ
ユーザー mdj982mdj982
提出日時 2019-03-02 00:29:33
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 28 ms / 5,000 ms
コード長 3,877 bytes
コンパイル時間 1,897 ms
コンパイル使用メモリ 182,196 KB
実行使用メモリ 16,256 KB
最終ジャッジ日時 2024-06-23 12:12:59
合計ジャッジ時間 3,111 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 2 ms
6,940 KB
testcase_02 AC 1 ms
6,940 KB
testcase_03 AC 2 ms
6,940 KB
testcase_04 AC 2 ms
6,940 KB
testcase_05 AC 2 ms
6,940 KB
testcase_06 AC 2 ms
6,944 KB
testcase_07 AC 28 ms
16,256 KB
testcase_08 AC 28 ms
16,256 KB
testcase_09 AC 26 ms
16,128 KB
testcase_10 AC 27 ms
16,128 KB
testcase_11 AC 28 ms
16,256 KB
testcase_12 AC 2 ms
6,944 KB
testcase_13 AC 4 ms
6,944 KB
testcase_14 AC 5 ms
6,940 KB
testcase_15 AC 4 ms
6,940 KB
testcase_16 AC 5 ms
6,940 KB
testcase_17 AC 4 ms
6,944 KB
testcase_18 AC 24 ms
14,848 KB
testcase_19 AC 23 ms
14,208 KB
testcase_20 AC 26 ms
16,000 KB
testcase_21 AC 27 ms
15,872 KB
testcase_22 AC 24 ms
14,720 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

using vi = vector<int>; using vvi = vector<vi>; using vvvi = vector<vvi>;
using ll = long long int;
using vll = vector<ll>; using vvll = vector<vll>; using vvvll = vector<vvll>;
using vd = vector<double>; using vvd = vector<vd>; using vvvd = vector<vvd>;
using P = pair<int, int>;
using Pll = pair<ll, ll>;
using cdouble = complex<double>;

const double INFD = numeric_limits<double>::infinity();
const double PI = 3.14159265358979323846;
#define Loop(i, n) for(int i = 0; i < (int)(n); i++)
#define Loop1(i, n) for(int i = 1; i <= (int)(n); i++)
#define Loopr(i, n) for(int i = (int)(n) - 1; i >= 0; i--)
#define Loopr1(i, n) for(int i = (int)(n); i >= 1; i--)
#define Foreach(buf, container) for(auto buf : container)
#define Loopdiag(i, j, h, w, sum) for(int i = ((sum) >= (h) ? (h) - 1 : (sum)), j = (sum) - i; i >= 0 && j < (w); i--, j++)
#define Loopdiagr(i, j, h, w, sum) for(int j = ((sum) >= (w) ? (w) - 1 : (sum)), i = (sum) - j; j >= 0 && i < (h); j--, i++)
#define Loopdiagsym(i, j, h, w, gap) for (int i = ((gap) >= 0 ? (gap) : 0), j = i - (gap); i < (h) && j < (w); i++, j++)
#define Loopdiagsymr(i, j, h, w, gap) for (int i = ((gap) > (h) - (w) - 1 ? (h) - 1 : (w) - 1 + (gap)), j = i - (gap); i >= 0 && j >= 0; i--, j--)
#define Loopitr(itr, container) for(auto itr = container.begin(); itr != container.end(); itr++)
#define printv(vector) Loop(ex_i, vector.size()) { cout << vector[ex_i] << " "; } cout << endl;
#define printmx(matrix) Loop(ex_i, matrix.size()) { Loop(ex_j, matrix[ex_i].size()) { cout << matrix[ex_i][ex_j] << " "; } cout << endl; }
#define quickio() ios::sync_with_stdio(false); cin.tie(0);
#define bitmanip(m,val) static_cast<bitset<(int)m>>(val)
#define Comp(type_t) bool operator<(const type_t &another) const
#define fst first
#define snd second
bool nearlyeq(double x, double y) { return abs(x - y) < 1e-9; }
bool inrange(int x, int t) { return x >= 0 && x < t; }
bool inrange(vi xs, int t) { Foreach(x, xs) if (!(x >= 0 && x < t)) return false; return true; }
ll rndf(double x) { return (ll)(x + (x >= 0 ? 0.5 : -0.5)); }
ll floorsqrt(ll x) { ll m = (ll)sqrt((double)x); return m + (m * m <= x ? 0 : -1); }
ll ceilsqrt(ll x) { ll m = (ll)sqrt((double)x); return m + (x <= m * m ? 0 : 1); }
ll rnddiv(ll a, ll b) { return (a / b + (a % b * 2 >= b ? 1 : 0)); }
ll ceildiv(ll a, ll b) { return (a / b + (a % b == 0 ? 0 : 1)); }
ll gcd(ll m, ll n) { if (n == 0) return m; else return gcd(n, m % n); }
ll lcm(ll m, ll n) { return m * n / gcd(m, n); }

/*******************************************************/

int n;
vvvd exdp;
vvvi deg;

void solve(int xi, int xj, int xk) {
	if (inrange(xk + 1, n + 1)) {
		exdp[xi][xj][xk + 1] += (exdp[xi][xj][xk] + double(n) / (xi + xj + xk + 1)) * double(xk + 1) / (xi + xj + xk + 1);
		deg[xi][xj][xk + 1]--;
		if (deg[xi][xj][xk + 1] == 0)	solve(xi, xj, xk + 1);
	}
	if (inrange({ xj + 1, xk - 1 }, n + 1)) {
		exdp[xi][xj + 1][xk - 1] += (exdp[xi][xj][xk] + double(n) / (xi + xj + xk)) * double(xj + 1) / (xi + xj + xk);
		deg[xi][xj + 1][xk - 1]--;
		if (deg[xi][xj + 1][xk - 1] == 0)	solve(xi, xj + 1, xk - 1);
	}
	if (inrange({ xi + 1, xj - 1 }, n + 1)) {
		exdp[xi + 1][xj - 1][xk] += (exdp[xi][xj][xk] + double(n) / (xi + xj + xk)) * double(xi + 1) / (xi + xj + xk);
		deg[xi + 1][xj - 1][xk]--;
		if (deg[xi + 1][xj - 1][xk] == 0)	solve(xi + 1, xj - 1, xk);
	}
}

int main() {
	cin >> n;
	exdp = vvvd(n + 1, vvd(n + 1, vd(n + 1)));
	deg = vvvi(n + 1, vvi(n + 1, vi(n + 1)));
	Loop(i, n + 1) {
		Loop(j, n + 1) {
			Loop(k, n + 1) {
				if (i > 0) deg[i][j][k]++;
				if (j > 0) deg[i][j][k]++;
				if (k > 0) deg[i][j][k]++;
			}
		}
	}
	solve(0, 0, 0);
	vi cnt(3);
	Loop(i, n) {
		int x; cin >> x;
		if (x == 0) cnt[0]++;
		if (x == 1) cnt[1]++;
		if (x == 2) cnt[2]++;
	}
	cout << setprecision(12) << exdp[cnt[0]][cnt[1]][cnt[2]] << endl;
}
0