結果

問題 No.584 赤、緑、青の色塗り
ユーザー ふっぴーふっぴー
提出日時 2017-10-28 12:56:54
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 2,097 bytes
コンパイル時間 1,736 ms
コンパイル使用メモリ 171,112 KB
実行使用メモリ 10,148 KB
最終ジャッジ日時 2024-11-22 03:05:58
合計ジャッジ時間 5,556 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
10,148 KB
testcase_01 AC 2 ms
5,248 KB
testcase_02 AC 2 ms
5,248 KB
testcase_03 AC 2 ms
5,248 KB
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 AC 2 ms
5,248 KB
testcase_08 AC 2 ms
5,248 KB
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 AC 3 ms
5,248 KB
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 TLE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

#define DEBUG(x) cout<<#x<<": "<<x<<endl;
#define DEBUG_VEC(v) cout<<#v<<":";for(int i=0;i<v.size();i++) cout<<" "<<v[i]; cout<<endl

typedef long long ll;
#define vi vector<int>
#define vl vector<ll>
#define vii vector< vector<int> >
#define vll vector< vector<ll> >
#define vs vector<string>
#define pii pair<int,int>
#define pis pair<int,string>
#define psi pair<string,int>
#define pll pair<ll,ll>
const int inf = 1000000001;
const ll INF = 1e18 * 4;
#define MOD 1000000007
#define mod 1000000009
#define pi 3.14159265358979323846
#define Sp(p) cout<<setprecision(15)<<fixed<<p<<endl;
int dx[4] = { 1,0,-1,0 }, dy[4] = { 0,1,0,-1 };
int dx2[8] = { 1,1,0,-1,-1,-1,0,1 }, dy2[8] = { 0,1,1,1,0,-1,-1,-1 };

ll MAXM = 3010;
vl fact(MAXM);

ll mod_pow(ll x, ll p, ll M) {
	ll a = 1;
	while (p) {
		if (p % 2)
			a = a*x%M;
		x = x*x%M;
		p /= 2;
	}
	return a;
}

ll mod_inverse(ll a, ll m) {
	return mod_pow(a, m - 2, m);
}

ll nCr(ll n, ll r) {
	ll ret = fact[n];
	ret = (ret*mod_inverse(fact[r], MOD)) % MOD;
	ret = (ret*mod_inverse(fact[n - r], MOD)) % MOD;
	return ret;
}


int main() {
	int n, r, g, b;
	cin >> n >> r >> g >> b;

	fact[0] = 1;
	for (int i = 1; i <= n; i++) {
		fact[i] = i*fact[i - 1];
	}

	int t1, t2, t3;
	ll ans = 0;
	for (t1 = 0; t1 <= r; t1++) {
		for (t2 = 0; t2 <= min(r, g) && r - t1 - t2 >= 0; t2++) {
			for (t3 = 0; t3 <= min(g, b) && g - t2 - t3 >= 0 && b - r + t1 + t2 - t3 >= 0; t3++) {
				if (n + 1 - 2 * g - 2 * b - r - t1 + t3 > n + 1 || n + 1 - 2 * g - 2 * b - r - t1 + t3 < 0) {
					break;
				}
				int z = n + 1 - 2 * g - 2 * b - r - t1 + t3;
				int x1 = r - t1 - t2, x2 = t2, x3 = t3;
				int y1 = t1, y2 = b - x1 - t3, y3 = g - t2 - t3;
				int m = z + x1 + x2 + x3 + y1 + y2 + y3;
				ans = (ans + (nCr(m, x1) * mod_pow(2, x1, MOD) % MOD * nCr(m - x1, x2) % MOD * mod_pow(2,x2,MOD) % MOD * nCr(m - x1 - x2, x3) % MOD * mod_pow(2,x3,MOD) % MOD * nCr(m - x1 - x2 - x3, y1) % MOD * nCr(m - x1 - x2 - x3 - y1, y2) % MOD*nCr(m - x1 - x2 - x3 - y1 - y2, y3) % MOD)) % MOD;
			}
		}
	}
	cout << ans << endl;
}
0