結果
問題 | No.584 赤、緑、青の色塗り |
ユーザー | kopricky |
提出日時 | 2017-10-28 18:27:15 |
言語 | C++11 (gcc 11.4.0) |
結果 |
AC
|
実行時間 | 514 ms / 2,000 ms |
コード長 | 2,640 bytes |
コンパイル時間 | 1,588 ms |
コンパイル使用メモリ | 160,336 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-05-01 21:17:10 |
合計ジャッジ時間 | 2,805 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
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,944 KB |
testcase_04 | AC | 1 ms
6,944 KB |
testcase_05 | AC | 2 ms
6,940 KB |
testcase_06 | AC | 2 ms
6,940 KB |
testcase_07 | AC | 1 ms
6,944 KB |
testcase_08 | AC | 2 ms
6,940 KB |
testcase_09 | AC | 2 ms
6,944 KB |
testcase_10 | AC | 2 ms
6,944 KB |
testcase_11 | AC | 1 ms
6,940 KB |
testcase_12 | AC | 2 ms
6,944 KB |
testcase_13 | AC | 2 ms
6,940 KB |
testcase_14 | AC | 4 ms
6,944 KB |
testcase_15 | AC | 2 ms
6,944 KB |
testcase_16 | AC | 2 ms
6,940 KB |
testcase_17 | AC | 3 ms
6,944 KB |
testcase_18 | AC | 2 ms
6,940 KB |
testcase_19 | AC | 514 ms
6,940 KB |
ソースコード
#include <bits/stdc++.h> #define ll long long #define INF 1000000005 #define MOD 1000000007 #define EPS 1e-10 #define rep(i,n) for(int i=0;i<(int)(n);++i) #define rrep(i,n) for(int i=(int)(n)-1;i>=0;--i) #define srep(i,s,t) for(int i=(int)(s);i<(int)(t);++i) #define each(a,b) for(auto (a): (b)) #define all(v) (v).begin(),(v).end() #define len(v) (int)(v).size() #define zip(v) sort(all(v)),v.erase(unique(all(v)),v.end()) #define cmx(x,y) x=max(x,y) #define cmn(x,y) x=min(x,y) #define fi first #define se second #define pb push_back #define show(x) cout<<#x<<" = "<<(x)<<endl #define spair(p) cout<<#p<<": "<<p.fi<<" "<<p.se<<endl #define svec(v) cout<<#v<<":";rep(kbrni,v.size())cout<<" "<<v[kbrni];cout<<endl #define sset(s) cout<<#s<<":";each(kbrni,s)cout<<" "<<kbrni;cout<<endl #define smap(m) cout<<#m<<":";each(kbrni,m)cout<<" {"<<kbrni.first<<":"<<kbrni.second<<"}";cout<<endl using namespace std; typedef pair<int,int> P; typedef vector<int> vi; typedef vector<vi> vvi; typedef vector<ll> vl; typedef vector<double> vd; typedef vector<P> vp; const int MAX_N = 3005; ll dp[MAX_N]; ll inv[MAX_N],fac[MAX_N],finv[MAX_N]; void make() { fac[0] = fac[1] = 1; finv[0] = finv[1] = 1; inv[1] = 1; for(int i=2;i<MAX_N;i++){ inv[i] = MOD - inv[MOD%i] * (MOD/i) % MOD; fac[i] = fac[i-1] * (ll) i % MOD; finv[i] = finv[i-1] * inv[i] % MOD; } } ll comb(int a,int b) { if(a<b){ return 0; } return fac[a] * (finv[b] * finv[a-b] % MOD) % MOD; } ll mod_pow(ll a,ll b) { a %= MOD; ll res = 1; while(b){ if(b & 1){ res = res * a % MOD; } a = a * a %MOD; b >>= 1; } return res; } ll add(ll x,ll y) { return (x + y)%MOD; } ll sub(ll x,ll y) { return (x+MOD-y)%MOD; } ll mul(ll x,ll y) { return x*y%MOD; } int main() { cin.tie(0); ios::sync_with_stdio(false); int n,r,g,b; cin >> n >> r >> g >> b; int a = r + g + b; ll ans = 0; make(); srep(i,(a+1)/2,a+1){ if(i-1 > n-a){ break; } int rem = (n-a)-(i-1); int one = 2*i-a, two = a-i; ll hoge = mul(fac[i],mul(comb(rem+i,i),mod_pow(2,two))); rep(j,min(min(r,g),two)+1){ rep(k,min(min(g-j,b),two-j)+1){ int l = two-j-k; if(l > min(r-j,b-k)){ continue; } ll tm = hoge; tm = mul(mul(tm,finv[r-j-l]),mul(finv[g-j-k],finv[b-k-l])); tm = mul(mul(tm,finv[j]),mul(finv[k],finv[l])); ans = add(ans,tm); } } } cout << ans << "\n"; return 0; }