結果

問題 No.802 だいたい等差数列
ユーザー 37zigen37zigen
提出日時 2019-03-18 03:01:17
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 1,337 ms / 2,000 ms
コード長 822 bytes
コンパイル時間 2,796 ms
コンパイル使用メモリ 144,444 KB
実行使用メモリ 34,876 KB
最終ジャッジ日時 2023-09-22 17:18:09
合計ジャッジ時間 44,534 ms
ジャッジサーバーID
(参考情報)
judge15 / judge14
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1,160 ms
34,688 KB
testcase_01 AC 1,160 ms
34,592 KB
testcase_02 AC 1,161 ms
34,716 KB
testcase_03 AC 1,161 ms
34,536 KB
testcase_04 AC 1,158 ms
34,544 KB
testcase_05 AC 1,159 ms
34,720 KB
testcase_06 AC 1,160 ms
34,592 KB
testcase_07 AC 1,158 ms
34,528 KB
testcase_08 AC 1,161 ms
34,636 KB
testcase_09 AC 1,157 ms
34,588 KB
testcase_10 AC 1,161 ms
34,528 KB
testcase_11 AC 1,159 ms
34,588 KB
testcase_12 AC 1,158 ms
34,540 KB
testcase_13 AC 1,160 ms
34,592 KB
testcase_14 AC 1,175 ms
34,560 KB
testcase_15 AC 1,205 ms
34,528 KB
testcase_16 AC 1,162 ms
34,540 KB
testcase_17 AC 1,159 ms
34,588 KB
testcase_18 AC 1,159 ms
34,628 KB
testcase_19 AC 1,172 ms
34,876 KB
testcase_20 AC 1,337 ms
34,560 KB
testcase_21 AC 1,292 ms
34,796 KB
testcase_22 AC 1,160 ms
34,624 KB
testcase_23 AC 1,199 ms
34,588 KB
testcase_24 AC 1,158 ms
34,588 KB
testcase_25 AC 1,157 ms
34,640 KB
testcase_26 AC 1,165 ms
34,540 KB
testcase_27 AC 1,158 ms
34,716 KB
testcase_28 AC 1,200 ms
34,680 KB
testcase_29 AC 1,290 ms
34,524 KB
testcase_30 AC 1,160 ms
34,536 KB
testcase_31 AC 1,161 ms
34,628 KB
testcase_32 AC 1,159 ms
34,688 KB
testcase_33 AC 1,161 ms
34,560 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>

long long mo=(long long)1e9+7;

long long fac[2000000];
long long ifac[2000000];

long long modpow(long long a,long long n){
  long long ret=1;
  for(;n>0;n>>=1,a=a*a%mo){
    if(n%2==1){
      ret=ret*a%mo;
    }
  }
  return ret;
}

long long inv(long a){
  return modpow(a,mo-2);
}

long long comb(int n,int k){
  return fac[n]*ifac[k]%mo*ifac[n-k]%mo;
}

void build(){
  fac[0]=1;
  for(int i=1;i<2000000;++i){
    fac[i]=fac[i-1]*i%mo;
  }
  for(int i=0;i<2000000;++i){
    ifac[i]=inv(fac[i]);
  }
}


int main(){
  build();
  long long N,M,D1,D2;
  std::cin>>N>>M>>D1>>D2;
  long long ans=0;
  for(int i=0;i<=N-1;++i){
    long long  res=M-1-D1*(N-1)-i*(D2-D1+1);
    if(res<0)continue;
    ans=(ans+modpow(mo-1,i)*comb(N-1,i)%mo*comb(N+res,res)%mo)%mo;
  }
  std::cout<<ans<<std::endl;
}
0