結果

問題 No.802 だいたい等差数列
ユーザー 37zigen
提出日時 2019-03-18 03:01:17
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 1,185 ms / 2,000 ms
コード長 822 bytes
コンパイル時間 1,182 ms
コンパイル使用メモリ 159,740 KB
実行使用メモリ 34,688 KB
最終ジャッジ日時 2024-07-08 08:39:30
合計ジャッジ時間 39,156 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

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