結果

問題 No.1683 Robot Guidance
ユーザー wait_sushiwait_sushi
提出日時 2021-09-17 22:59:10
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 143 ms / 2,000 ms
コード長 7,536 bytes
コンパイル時間 2,205 ms
コンパイル使用メモリ 221,484 KB
実行使用メモリ 30,720 KB
最終ジャッジ日時 2024-06-29 21:46:28
合計ジャッジ時間 9,227 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 128 ms
30,532 KB
testcase_01 AC 132 ms
30,592 KB
testcase_02 AC 140 ms
30,592 KB
testcase_03 AC 129 ms
30,568 KB
testcase_04 AC 125 ms
30,604 KB
testcase_05 AC 126 ms
30,592 KB
testcase_06 AC 130 ms
30,572 KB
testcase_07 AC 126 ms
30,592 KB
testcase_08 AC 139 ms
30,556 KB
testcase_09 AC 129 ms
30,592 KB
testcase_10 AC 133 ms
30,592 KB
testcase_11 AC 134 ms
30,592 KB
testcase_12 AC 139 ms
30,592 KB
testcase_13 AC 126 ms
30,464 KB
testcase_14 AC 138 ms
30,592 KB
testcase_15 AC 129 ms
30,592 KB
testcase_16 AC 130 ms
30,628 KB
testcase_17 AC 129 ms
30,464 KB
testcase_18 AC 128 ms
30,572 KB
testcase_19 AC 127 ms
30,580 KB
testcase_20 AC 129 ms
30,552 KB
testcase_21 AC 127 ms
30,628 KB
testcase_22 AC 128 ms
30,464 KB
testcase_23 AC 130 ms
30,720 KB
testcase_24 AC 143 ms
30,552 KB
testcase_25 AC 128 ms
30,712 KB
testcase_26 AC 127 ms
30,624 KB
testcase_27 AC 129 ms
30,528 KB
testcase_28 AC 127 ms
30,616 KB
testcase_29 AC 127 ms
30,652 KB
testcase_30 AC 127 ms
30,592 KB
testcase_31 AC 132 ms
30,536 KB
testcase_32 AC 128 ms
30,540 KB
testcase_33 AC 129 ms
30,592 KB
testcase_34 AC 134 ms
30,592 KB
testcase_35 AC 129 ms
30,572 KB
testcase_36 AC 129 ms
30,592 KB
testcase_37 AC 129 ms
30,592 KB
testcase_38 AC 126 ms
30,576 KB
testcase_39 AC 129 ms
30,632 KB
testcase_40 AC 125 ms
30,552 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define endl "\n"
using namespace std;
typedef long long ll;
typedef long double ld;
typedef vector<ll> vl;
typedef pair<ll, ll> PP;
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
#pragma GCC target("avx")
#pragma comment(linker, "/stack:200000000")
#define overload4(_1, _2, _3, _4, name, ...) name
#define overload3(_1, _2, _3, name, ...) name
#define rep1(n) for (ll i = 0; i < n; ++i)
#define rep2(i, n) for (ll i = 0; i < n; ++i)
#define rep3(i, a, b) for (ll i = a; i < b; ++i)
#define rep4(i, a, b, c) for (ll i = a; i < b; i += c)
#define rep(...) overload4(__VA_ARGS__, rep4, rep3, rep2, rep1)(__VA_ARGS__)
#define rrep1(n) for (ll i = n; i--;)
#define rrep2(i, n) for (ll i = n; i--;)
#define rrep3(i, a, b) for (ll i = b; i-- > (a);)
#define rrep4(i, a, b, c)                                                      \
  for (ll i = (a) + ((b) - (a)-1) / (c) * (c); i >= (a); i -= c)
#define rrep(...)                                                              \
  overload4(__VA_ARGS__, rrep4, rrep3, rrep2, rrep1)(__VA_ARGS__)
#define elif else if
ll beki(ll i) { return 1ll << i; }
#define skip(M)                                                                \
  if (M)                                                                       \
  continue
#define all(v) v.begin(), v.end()
#define allm(x, y, M) for (auto [(x), (y)] : (M))
#define alls(i, S) for (auto(&i) : (S))
#define pb push_back
#define fi first
#define se second
const ll INF = (1LL << 61) - 1;
const ll mod = 1000000007;
// const ll mod = 998244353;
const ll MAX_N1 = 200010;
const ll MAX_N2 = 500010;
template <class X> void print(X x) { cout << x << endl; }
template <class X> void escape(X x) {
  cout << x << endl;
  return;
}
ll min(int x, ll y) { return min(ll(x), y); }
ll min(ll x, int y) { return min(x, ll(y)); }
ll max(int x, ll y) { return max(ll(x), y); }
ll max(ll x, int y) { return max(x, ll(y)); }
ld DIV(ll x, ll y) { return ld(x) / ld(y); }
ll PLUS(ll x, ll y) {
  if (x < 0 && y < 0)
    return -PLUS(-x, -y);
  if (x > INF - y)
    return INF;
  return x + y;
}
ll MUL(ll x, ll y) {
  if (x == 0 || y == 0)
    return 0;
  if (x < 0 || y < 0)
    return signbit(x) * signbit(y) * (MUL(x, y));
  if (x > INF / y)
    return INF;
  return x * y;
}
void print(vl x) {
  for (ll i : x) {
    cout << i << " ";
  }
  cout << endl;
}
void print(vector<PP> x) {
  for (PP i : x) {
    cout << i.first << " " << i.second << endl;
  }
  cout << endl;
}
template <class X> void printr(X x) { cerr << x << endl; }
void printr(vl x) {
  for (ll i : x) {
    cerr << i << " ";
  }
  cerr << endl;
}
void printr(vector<PP> x) {
  for (PP i : x) {
    cerr << i.first << " " << i.second << endl;
  }
  cout << endl;
}
template <class X> int lbound(vector<X> &Y, X a) {
  return lower_bound(all(Y), a) - Y.begin();
}
template <class X> int ubound(vector<X> &Y, X a) {
  return upper_bound(all(Y), a) - Y.begin();
}
template <class X> X lbound2(vector<X> &Y, X a, int mode = 0) {
  int x = lbound(Y, a) - mode;
  if (x < 0 || x >= Y.size())
    return INF;
  return Y[x];
}
template <class X> X ubound2(vector<X> &Y, X a, int mode = 0) {
  int x = ubound(Y, a) - mode;
  if (x < 0 || x >= Y.size())
    return INF;
  return Y[x];
}
template <class... T> void cl(T &... t) { (..., (t.clear())); }
template <class... T> void in(T &... t) { (..., (cin >> t)); }
template <class... T> void put(vl &V, int n, int mode = 0) {
  ll k;
  if (mode == 0)
    cl(V);
  rep(i, n) {
    cin >> k;
    V.pb(k);
  }
}
int max_index(vl &V) { return max_element(all(V)) - V.begin(); }
int min_index(vl &V) { return min_element(all(V)) - V.begin(); }
ll sum(vl &V) { return accumulate(all(V), 0ll); }
template <typename T> void UNIQUE(vector<T> &v) {
  sort(v.begin(), v.end());
  v.erase(unique(v.begin(), v.end()), v.end());
}
template <typename T> T ADD(T a, T b) {
  T res;
  return __builtin_add_overflow(a, b, &res) ? numeric_limits<T>::max() : res;
}
template <typename T> T MUL(T a, T b) {
  T res;
  return __builtin_mul_overflow(a, b, &res) ? numeric_limits<T>::max() : res;
};
template <typename T> T POW(T a, int n) {
  T res = 1;
  for (; n > 0; n >>= 1, a *= a)
    if (n & 1)
      res *= a;
  return res;
}
template <class X> bool ma(X &a, X b) {
  if (a < b) {
    a = b;
    return true;
  }
  return false;
}
template <class X> bool mi(X &a, X b) {
  if (a > b) {
    a = b;
    return true;
  }
  return false;
}
void vset(vl &A, ll n, ll k) { A = vl(n, k); }
void so(vl &A) { sort(all(A)); }
ll age(ll x, ll y) { return (x + y - 1) / y; }
ll a, b, c, d, aa, bb, cc, dd, h, x, y, z, p, q, n, t, r, k, w, l, ans, m, u, v;
ll codeforces = 1;
string S, T, U;
vl A, B, C, D;
vl g[MAX_N2];
vector<PP> VP, VP2;
set<ll> s1, s2;
map<ll, ll> M;
struct mint {
  using ull = unsigned long long int;
  ull v;
  mint(ll vv = 0) { s(vv % mod + mod); }
  mint &s(ull vv) {
    v = vv < mod ? vv : vv - mod;
    return *this;
  }
  //オーバーロード
  mint operator-() const { return mint() - *this; } // mint型にキャスト
  mint &operator+=(const mint &val) { return s(v + val.v); }
  mint &operator-=(const mint &val) { return s(v + mod - val.v); }
  mint &operator*=(const mint &val) {
    v = ull(v) * val.v % mod;
    return *this;
  }
  mint &operator/=(const mint &val) { return *this *= val.inv(); }
  mint operator+(const mint &val) { return mint(*this) += val; }
  mint operator-(const mint &val) { return mint(*this) -= val; }
  mint operator*(const mint &val) { return mint(*this) *= val; }
  mint operator/(const mint &val) { return mint(*this) /= val; }
  mint pow(ll n) const {
    mint res(1), x(*this);
    while (n) {
      if (n & 1)
        res *= x;
      x *= x;
      n >>= 1ll;
    }
    return res;
  }
  mint inv() const { return pow(mod - 2); }
  //拡張ユークリッドの互除法
  /* mint inv()const{
  int x,y;
  int g=extgcd(v,mod,x,y);
  assert(g==1);
  if(x<0)x+=mod;
  return mint(x);
  }*/
  friend ostream &operator<<(ostream &os, const mint &val) {
    return os << val.v;
  } //出力
  bool operator<(const mint &val) const { return v < val.v; }
  bool operator==(const mint &val) const { return v == val.v; }
  bool operator>(const mint &val) const { return v > val.v; }
};
vector<mint> fact(MAX_N2 * 2 + 1), invfact(MAX_N2 * 2 + 1);
void factcreate(void) {
  fact[0] = 1;
  rep(i, MAX_N2 * 2) { fact[i + 1] = fact[i] * (i + 1); }
  rep(i, MAX_N2 * 2) { invfact[i] = mint(1) / fact[i]; }
}
mint nCk(ll x, ll y) {
  if (x < 0 || y < 0 || x < y)
    return 0;
  return fact[x] * invfact[y] * invfact[x - y];
}
void solve() {
  in(a, b, x, y);
  mint ans = 0;
  t = b + 1;
  q = t / 4;
  rep(4) {
    if (t % 4 != 0) {
      A.pb(q + 1);
      t--;
    } else
      A.pb(q);
  }
  rep(a + 1) {
    cl(B);
    aa = i;
    bb = aa - x;
    if (bb < 0)
      continue;
    p = a - aa - bb;
    if (p < 0)
      continue;
    if ((p + y) % 2 != 0)
      continue;
    cc = (y + p) / 2;
    dd = (-y + p) / 2;
    if (cc < 0 || dd < 0)
      continue;
    mint aaa = 1;
    B.pb(aa);
    B.pb(cc);
    B.pb(bb);
    B.pb(dd);
    // cout << aa << " " << cc << " " << bb << " " << dd << endl;
    rep(j, 4) {
      n = A[j];
      m = B[j];
      if (n == 0 && m == 0)
        continue;
      aaa *= nCk(n + m - 1, n - 1);
    }
    ans += aaa;
  }
  // print(A);
  print(ans);
}
int main() {
  cout << fixed << setprecision(15);
  cin.tie(0);
  ios::sync_with_stdio(false);
  // cin >> codeforces;
  factcreate();
  while (codeforces--) {
    ans = 0;
    solve();
  }
}
0