結果

問題 No.1683 Robot Guidance
ユーザー monnumonnu
提出日時 2021-09-18 11:38:28
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 2,315 bytes
コンパイル時間 4,142 ms
コンパイル使用メモリ 230,608 KB
実行使用メモリ 50,532 KB
最終ジャッジ日時 2024-06-30 08:51:19
合計ジャッジ時間 6,758 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 AC 47 ms
50,276 KB
testcase_02 AC 47 ms
50,432 KB
testcase_03 AC 52 ms
50,340 KB
testcase_04 AC 46 ms
50,328 KB
testcase_05 WA -
testcase_06 AC 50 ms
50,388 KB
testcase_07 WA -
testcase_08 AC 63 ms
50,432 KB
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 AC 61 ms
50,304 KB
testcase_13 AC 2 ms
6,940 KB
testcase_14 WA -
testcase_15 AC 51 ms
50,312 KB
testcase_16 AC 2 ms
6,944 KB
testcase_17 AC 2 ms
6,940 KB
testcase_18 AC 2 ms
6,944 KB
testcase_19 WA -
testcase_20 AC 50 ms
50,316 KB
testcase_21 WA -
testcase_22 AC 47 ms
50,272 KB
testcase_23 WA -
testcase_24 AC 68 ms
50,292 KB
testcase_25 WA -
testcase_26 AC 2 ms
6,940 KB
testcase_27 AC 45 ms
50,308 KB
testcase_28 AC 56 ms
50,380 KB
testcase_29 AC 2 ms
6,944 KB
testcase_30 AC 52 ms
50,376 KB
testcase_31 AC 53 ms
50,376 KB
testcase_32 WA -
testcase_33 AC 55 ms
50,320 KB
testcase_34 WA -
testcase_35 AC 54 ms
50,368 KB
testcase_36 WA -
testcase_37 AC 58 ms
50,384 KB
testcase_38 AC 86 ms
50,364 KB
testcase_39 AC 50 ms
50,424 KB
testcase_40 AC 2 ms
6,944 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