結果
問題 | No.895 MESE |
ユーザー |
![]() |
提出日時 | 2019-09-27 22:23:44 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 24 ms / 2,000 ms |
コード長 | 1,344 bytes |
コンパイル時間 | 1,073 ms |
コンパイル使用メモリ | 93,768 KB |
実行使用メモリ | 10,564 KB |
最終ジャッジ日時 | 2024-09-24 18:26:24 |
合計ジャッジ時間 | 3,168 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 26 |
ソースコード
#include <iostream>#include <algorithm>#include <string>#include <vector>#include <cmath>#include <map>#include <queue>#include <iomanip>#include <set>#include <tuple>#define mkp make_pair#define mkt make_tuple#define rep(i,n) for(int i = 0; i < (n); ++i)using namespace std;typedef long long ll;const ll MOD=1e9+7;#define MAX_N 300030ll inv[MAX_N+10],fac[MAX_N+10],ifac[MAX_N+10];void setComb(){inv[0]=1;inv[1]=1;fac[1]=1;ifac[1]=1;fac[0]=1;ifac[0]=1;for(int i=2;i<MAX_N;i++){inv[i]=(-MOD/i)*inv[MOD%i]%MOD;fac[i]=fac[i-1]*i%MOD;ifac[i]=ifac[i-1]*inv[i]%MOD;inv[i]=(inv[i]+MOD)%MOD;fac[i]=(fac[i]+MOD)%MOD;ifac[i]=(ifac[i]+MOD)%MOD;}return;}ll comb(ll n,ll k){if(n<k||n<0||k<0) return 0;else return ((fac[n]*ifac[k]%MOD*ifac[n-k]%MOD+MOD)%MOD);}ll hcomb(ll n,ll r){if(n==0&&r==0) return 1;else if(n<0||r<0) return 0;else return comb(n+r-1,r-1);}void add(ll &a,ll b){a=(a+b)%MOD;}void mul(ll &a,ll b){a=a*b%MOD;}ll A,B,C;int main(){cin>>A>>B>>C;setComb();ll ans=0;ll two=1;for(ll i=0;i<A+B+C-2;i++){ll p=comb(A+B+C-2,A-1);ll s=A+B+C-2-i;ll q=comb(i,A-1-s);add(p,MOD-q);mul(p,comb(B+C-2,B-1));mul(p,two);add(ans,p);mul(two,2ll);}cout<<ans<<endl;return 0;}