結果

問題 No.802 だいたい等差数列
ユーザー goodbatongoodbaton
提出日時 2019-03-18 08:58:24
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 70 ms / 2,000 ms
コード長 2,996 bytes
コンパイル時間 924 ms
コンパイル使用メモリ 98,964 KB
実行使用メモリ 23,524 KB
最終ジャッジ日時 2023-09-24 20:48:44
合計ジャッジ時間 3,248 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,380 KB
testcase_01 AC 2 ms
4,376 KB
testcase_02 AC 2 ms
4,376 KB
testcase_03 AC 2 ms
4,376 KB
testcase_04 AC 1 ms
4,376 KB
testcase_05 AC 1 ms
4,380 KB
testcase_06 AC 1 ms
4,380 KB
testcase_07 AC 1 ms
4,380 KB
testcase_08 AC 2 ms
4,380 KB
testcase_09 AC 2 ms
4,376 KB
testcase_10 AC 41 ms
18,724 KB
testcase_11 AC 51 ms
23,068 KB
testcase_12 AC 47 ms
21,128 KB
testcase_13 AC 50 ms
22,544 KB
testcase_14 AC 52 ms
21,884 KB
testcase_15 AC 51 ms
20,392 KB
testcase_16 AC 52 ms
22,328 KB
testcase_17 AC 50 ms
22,172 KB
testcase_18 AC 50 ms
22,356 KB
testcase_19 AC 55 ms
23,340 KB
testcase_20 AC 69 ms
23,376 KB
testcase_21 AC 70 ms
23,344 KB
testcase_22 AC 52 ms
23,316 KB
testcase_23 AC 58 ms
23,524 KB
testcase_24 AC 52 ms
23,316 KB
testcase_25 AC 42 ms
18,516 KB
testcase_26 AC 16 ms
8,008 KB
testcase_27 AC 51 ms
22,848 KB
testcase_28 AC 31 ms
12,248 KB
testcase_29 AC 69 ms
23,336 KB
testcase_30 AC 2 ms
4,376 KB
testcase_31 AC 1 ms
4,376 KB
testcase_32 AC 2 ms
4,376 KB
testcase_33 AC 27 ms
12,396 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>

#include <iostream>
#include <complex>
#include <string>
#include <algorithm>
#include <numeric>
#include <vector>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <unordered_map>
#include <unordered_set>

#include <functional>
#include <cassert>

typedef long long ll;
using namespace std;

#ifndef LOCAL
#define debug(x) ;
#else
#define debug(x) cerr << __LINE__ << " : " << #x << " = " << (x) << endl;

template <typename T1, typename T2>
ostream &operator<<(ostream &out, const pair<T1, T2> &p) {
  out << "{" << p.first << ", " << p.second << "}";
  return out;
}

template <typename T>
ostream &operator<<(ostream &out, const vector<T> &v) {
  out << '{';
  for (const T &item : v) out << item << ", ";
  out << "\b\b}";
  return out;
}
#endif

#define mod 1000000007 //1e9+7(prime number)
#define INF 1000000000 //1e9
#define LLINF 2000000000000000000LL //2e18
#define SIZE 200010

vector<ll> factmemo, factmemoInv;
ll factmemoMod = -1;

ll factorial(int n, int M){
  if(factmemoMod == M) return factmemo[n];
  if(n <= 1) return 1;

  ll res = 1;
  for(int i=1; i<=n; i++) res = res * i % M;
  return res;
}

ll power(int k, int n, int M){
  if(n == 0) return 1;
  if(n == 1) return (ll)k;

  ll res = power(k, n/2, M);

  res = res * res % M;
  return n%2 == 1 ? res * k % M : res;
}

void initFactorial(int n, int M){
  factmemo.assign(n+1, 0);
  factmemoInv.assign(n+1, 0);
  factmemoMod = M;
  factmemo[0] = 1;
  for(int i=1;i<=n;i++) factmemo[i] = factmemo[i-1] * i % M;
  factmemoInv[n] = power(factmemo[n], M-2, M);
  for(int i=n;i>0;i--) factmemoInv[i-1] = factmemoInv[i] * i % M;
}

//nCm nPm nHm (mod M)

/*Combination*/
ll C(int n, int m, int M){
  if(n < m) return 0;
  if(m == 0 || n == m) return 1;

  if(factmemoMod == M)
    return factmemo[n] * factmemoInv[m] % M * factmemoInv[n-m] % M;

  ll numer = factorial(n, M);
  ll denom = factorial(m, M) * factorial(n-m, M) % M;

  denom = power((int)denom, M-2, M);

  return numer * denom % M;
}
/*Permutation*/
ll P(int n, int m, int M){
  if(n < m) return 0;
  if(m == 0) return 1;

  if(factmemoMod == M)
    return factmemo[n] * factmemoInv[n-m] % M;

  ll numer = factorial(n, M);
  ll denom = factorial(n-m, M);

  denom = power((int)denom, M-2, M);

  return numer * denom % M;
}
/*Combination with Repetitions*/
ll H(int n, int m, int M){
  if(n == 0 && m == 0) return 1;
  return C(n+m-1, m, M);
}


int main(){
  ll N, M, D1, D2;

  cin >> N >> M >> D1 >> D2;

  initFactorial(N + M, mod);

  if(M < D1 * (N - 1)) {
    cout << 0 << endl;
    return 0;
  }
  M --;
  M -= D1 * (N - 1);
  D2 -= D1;

  ll ans = 0;
  int sign = 1;

  debug(M); debug(D2);

  for(int i=0;i<N;i++){
    ll m = M - (D2 + 1) * i;
    if(m < 0) break;
    ll p = C(m + N, N, mod);
    ll q = C(N-1, i, mod);
    debug(p); debug(q);
    ans += (mod + sign * p * q % mod) % mod;
    sign *= -1;
  }

  cout << ans % mod << endl;

  return 0;
}
0