結果

問題 No.1683 Robot Guidance
ユーザー monnumonnu
提出日時 2021-09-18 11:55:17
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 73 ms / 2,000 ms
コード長 2,263 bytes
コンパイル時間 3,612 ms
コンパイル使用メモリ 227,064 KB
実行使用メモリ 50,520 KB
最終ジャッジ日時 2023-09-12 21:39:41
合計ジャッジ時間 7,223 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 56 ms
50,248 KB
testcase_01 AC 54 ms
50,220 KB
testcase_02 AC 54 ms
50,264 KB
testcase_03 AC 57 ms
50,292 KB
testcase_04 AC 54 ms
50,332 KB
testcase_05 AC 55 ms
50,240 KB
testcase_06 AC 57 ms
50,304 KB
testcase_07 AC 57 ms
50,304 KB
testcase_08 AC 72 ms
50,240 KB
testcase_09 AC 52 ms
50,332 KB
testcase_10 AC 65 ms
50,220 KB
testcase_11 AC 62 ms
50,328 KB
testcase_12 AC 68 ms
50,240 KB
testcase_13 AC 2 ms
4,380 KB
testcase_14 AC 69 ms
50,520 KB
testcase_15 AC 57 ms
50,264 KB
testcase_16 AC 1 ms
4,380 KB
testcase_17 AC 1 ms
4,380 KB
testcase_18 AC 2 ms
4,376 KB
testcase_19 AC 56 ms
50,356 KB
testcase_20 AC 58 ms
50,224 KB
testcase_21 AC 56 ms
50,280 KB
testcase_22 AC 59 ms
50,272 KB
testcase_23 AC 56 ms
50,360 KB
testcase_24 AC 73 ms
50,220 KB
testcase_25 AC 55 ms
50,288 KB
testcase_26 AC 1 ms
4,376 KB
testcase_27 AC 58 ms
50,220 KB
testcase_28 AC 57 ms
50,244 KB
testcase_29 AC 1 ms
4,380 KB
testcase_30 AC 58 ms
50,240 KB
testcase_31 AC 54 ms
50,432 KB
testcase_32 AC 55 ms
50,216 KB
testcase_33 AC 59 ms
50,236 KB
testcase_34 AC 55 ms
50,304 KB
testcase_35 AC 57 ms
50,372 KB
testcase_36 AC 60 ms
50,280 KB
testcase_37 AC 54 ms
50,220 KB
testcase_38 AC 54 ms
50,236 KB
testcase_39 AC 55 ms
50,420 KB
testcase_40 AC 1 ms
4,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#include <atcoder/all>
using namespace atcoder;
using ll=long long;
using Graph=vector<vector<int>>;
#define MAX 2000005
#define MOD 1000000007
//#define MOD 998244353
#define INF 1000000000
//#define INF 1000000000000000000

ll fac[MAX],finv[MAX],inv[MAX];

void COMinit(){
  fac[0]=fac[1]=1;
  finv[0]=finv[1]=1;
  inv[1]=1;
  for(int i=2;i<MAX;i++){
    fac[i]=(fac[i-1]*i)%MOD;
    inv[i]=MOD-inv[MOD%i]*(MOD/i)%MOD;
    finv[i]=finv[i-1]*inv[i]%MOD;
  }
}

ll COM(int n,int k){
  if(n<k){
    return 0;
  }
  if(n<0||k<0){
    return 0;
  }
  return fac[n]*(finv[k]*finv[n-k]%MOD)%MOD;
}


int main(){
  int A,B,X,Y;
  cin>>A>>B>>X>>Y;
  if(abs(X)+abs(Y)>A){
    cout<<0<<'\n';
    return 0;
  }else if((abs(X)+abs(Y))%2!=A%2){
    cout<<0<<'\n';
    return 0;
  }

  COMinit();

  ll ans=0;
  for(int x=abs(X);x<=A-abs(Y);x+=2){
    int y=A-x;
    int x_plus=1;
    int y_plus=0;
    int x_minus=0;
    int y_minus=0;
    if(B%4==0){
      x_plus+=B/4;
      y_plus+=B/4;
      x_minus+=B/4;
      y_minus+=B/4;
    }else if(B%4==1){
      x_plus+=B/4;
      y_plus+=B/4+1;
      x_minus+=B/4;
      y_minus+=B/4;
    }else if(B%4==2){
      x_plus+=B/4;
      y_plus+=B/4+1;
      x_minus+=B/4+1;
      y_minus+=B/4;
    }else{
      x_plus+=B/4;
      y_plus+=B/4+1;
      x_minus+=B/4+1;
      y_minus+=B/4+1;
    }

    ll ret=1;
    if(X>0){
      if(X+(x-X)/2>0){
        ret*=COM(X+(x-X)/2+x_plus-1,x_plus-1);
        ret%=MOD;
      }
      if((x-X)/2>0){
        ret*=COM((x-X)/2+x_minus-1,x_minus-1);
        ret%=MOD;
      }
    }else{
      if(-X+(x+X)/2>0){
        ret*=COM(-X+(x+X)/2+x_minus-1,x_minus-1);
        ret%=MOD;
      }
      if((x+X)/2>0){
        ret*=COM((x+X)/2+x_plus-1,x_plus-1);
        ret%=MOD;
      }
    }
    if(Y>0){
      if(Y+(y-Y)/2>0){
        ret*=COM(Y+(y-Y)/2+y_plus-1,y_plus-1);
        ret%=MOD;
      }
      if((y-Y)/2>0){
        ret*=COM((y-Y)/2+y_minus-1,y_minus-1);
        ret%=MOD;
      }
    }else{
      if(-Y+(y+Y)/2>0){
        ret*=COM(-Y+(y+Y)/2+y_minus-1,y_minus-1);
        ret%=MOD;
      }
      if((y+Y)/2>0){
        ret*=COM((y+Y)/2+y_plus-1,y_plus-1);
        ret%=MOD;
      }
    }
    ans+=ret;
    ans%=MOD;
  }
  cout<<ans<<'\n';
}
0