結果

問題 No.1683 Robot Guidance
ユーザー monnumonnu
提出日時 2021-09-18 11:36:17
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 2,283 bytes
コンパイル時間 3,809 ms
コンパイル使用メモリ 230,668 KB
実行使用メモリ 50,524 KB
最終ジャッジ日時 2024-06-30 08:48:49
合計ジャッジ時間 7,286 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 54 ms
50,304 KB
testcase_01 AC 56 ms
50,304 KB
testcase_02 AC 57 ms
50,432 KB
testcase_03 AC 60 ms
50,304 KB
testcase_04 AC 59 ms
50,304 KB
testcase_05 WA -
testcase_06 AC 66 ms
50,280 KB
testcase_07 WA -
testcase_08 AC 73 ms
50,432 KB
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 AC 71 ms
50,304 KB
testcase_13 AC 2 ms
5,376 KB
testcase_14 WA -
testcase_15 AC 56 ms
50,432 KB
testcase_16 AC 2 ms
5,376 KB
testcase_17 AC 2 ms
5,376 KB
testcase_18 AC 2 ms
5,376 KB
testcase_19 AC 62 ms
50,304 KB
testcase_20 AC 59 ms
50,176 KB
testcase_21 AC 57 ms
50,304 KB
testcase_22 AC 64 ms
50,396 KB
testcase_23 AC 57 ms
50,304 KB
testcase_24 AC 77 ms
50,304 KB
testcase_25 AC 57 ms
50,176 KB
testcase_26 AC 2 ms
5,376 KB
testcase_27 AC 59 ms
50,316 KB
testcase_28 AC 61 ms
50,304 KB
testcase_29 AC 2 ms
5,376 KB
testcase_30 AC 57 ms
50,304 KB
testcase_31 AC 60 ms
50,304 KB
testcase_32 AC 61 ms
50,304 KB
testcase_33 AC 59 ms
50,360 KB
testcase_34 AC 63 ms
50,432 KB
testcase_35 AC 66 ms
50,304 KB
testcase_36 AC 68 ms
50,304 KB
testcase_37 WA -
testcase_38 AC 55 ms
50,292 KB
testcase_39 AC 56 ms
50,304 KB
testcase_40 AC 2 ms
5,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{
      X=-X;
      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{
      Y=-Y;
      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