結果
問題 | No.1683 Robot Guidance |
ユーザー |
|
提出日時 | 2021-09-18 14:25:45 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,800 bytes |
コンパイル時間 | 4,240 ms |
コンパイル使用メモリ | 253,312 KB |
最終ジャッジ日時 | 2025-01-24 15:43:14 |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 26 WA * 12 |
ソースコード
//#define _GLIBCXX_DEBUG #include <bits/stdc++.h> #include <iostream> #include <algorithm> #include <atcoder/all> using namespace std; using namespace atcoder; using mint = modint1000000007; using ll = long long; using pii = pair<int,int>; using pll = pair<ll,ll>; using vi = vector<int>; using vvi = vector<vi>; using vl = vector<ll>; using vvl = vector<vl>; using vb = vector<bool>; using vvb = vector<vb>; using vm = vector<mint>; using vvm = vector<vm>; using vpi = vector<pii>; using vvpi = vector<vpi>; using vpl = vector<pll>; using vvpl = vector<vpl>; const int inf = 1 << 30; const ll INF = 1LL << 60; #define rep(i,m,n) for (int i = m; i < (int)(n); i++) #define rrep(i,m,n) for (int i = m; i > (int)(n); i--) const int MAX = 2000001; //nの最大 vm fac(MAX),finv(MAX); void set_fac(){ fac[0] = fac[1] = 1; finv[0] = finv[1] = 1; rep(i,2,MAX){ fac[i] = fac[i-1]*i; finv[i] = finv[i-1] * mint(i).inv(); } } mint cmb(int n,int k){ if (k < 0 || n < 0 || n < k) return 0; return fac[n]*finv[n-k]*finv[k]; } int main(){ set_fac(); int a,b,x,y; cin >> a >> b >> x >> y; vi s(4,(b+1)/4-1); rep(i,0,(b+1)%4){ s[i]++; } rep(i,0,4){ if (s[i] < 0){ s[i] = 0; } } //cout << s[0] << s[1] << s[2] << s[3] << endl; mint ans = 0; rep(i,0,a+1){ int j = a - i; if ((i+x)%2 == 1 || (j+y)%2 == 1) continue; //cout << i << " " << j << endl; int x0 = (i+x)/2,x2 = (i-x)/2,y1 = (j+y)/2,y3 = (j-y)/2; if (x0 < 0 || x2 < 0 || y1 < 0 || y3 < 0) continue; //cout << x0 << " " << x2 << " " << y1 << " " << y3 << endl; ans += cmb(x0+s[0],s[0])*cmb(x2+s[2],s[2])*cmb(y1+s[1],s[1])*cmb(y3+s[3],s[3]); } cout << ans.val() << endl; }