結果
問題 | No.1683 Robot Guidance |
ユーザー |
|
提出日時 | 2021-09-17 21:56:54 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 40 ms / 2,000 ms |
コード長 | 1,254 bytes |
コンパイル時間 | 607 ms |
コンパイル使用メモリ | 76,220 KB |
実行使用メモリ | 22,816 KB |
最終ジャッジ日時 | 2024-06-29 20:23:15 |
合計ジャッジ時間 | 2,431 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 38 |
コンパイルメッセージ
main.cpp:55:1: warning: ISO C++ forbids declaration of 'main' with no type [-Wreturn-type] 55 | main() | ^~~~
ソースコード
#include<iostream> #include<algorithm> #include<atcoder/modint> using namespace std; #include<vector> template<typename T> struct combination{ vector<T>fac,invfac; combination(size_t N=0):fac(1,1),invfac(1,1) { make_table(N); } void make_table(size_t N) { if(fac.size()>N)return; size_t now=fac.size(); N=max(N,now*2); fac.resize(N+1); invfac.resize(N+1); for(size_t i=now;i<=N;i++)fac[i]=fac[i-1]*i; invfac[N]=1/fac[N]; for(size_t i=N;i-->now;)invfac[i]=invfac[i+1]*(i+1); } T factorial(size_t n) { make_table(n); return fac[n]; } T P(size_t n,size_t k) { if(n<k)return 0; make_table(n); return fac[n]*invfac[n-k]; } T C(size_t n,size_t k) { if(n<k)return 0; make_table(n); return fac[n]*invfac[n-k]*invfac[k]; } T H(size_t n,size_t k) { if(n==0)return k==0?1:0; return C(n-1+k,k); } }; using mint=atcoder::modint1000000007; combination<mint>C; mint f(int p,int n,int X,int T) { if(T<abs(X)||(T-abs(X))%2==1)return 0; int pc=(T+X)/2,nc=(T-X)/2; return C.H(p,pc)*C.H(n,nc); } main() { int A,B,X,Y; cin>>A>>B>>X>>Y; int px=(B+4)/4; int py=(B+3)/4; int nx=(B+2)/4; int ny=(B+1)/4; mint ans=0; for(int ax=0;ax<=A;ax++) { ans+=f(px,nx,X,ax)*f(py,ny,Y,A-ax); } cout<<ans.val()<<endl; }