結果
問題 | No.802 だいたい等差数列 |
ユーザー | risujiroh |
提出日時 | 2019-03-18 09:27:35 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
MLE
|
実行時間 | - |
コード長 | 3,025 bytes |
コンパイル時間 | 1,624 ms |
コンパイル使用メモリ | 171,888 KB |
実行使用メモリ | 814,440 KB |
最終ジャッジ日時 | 2024-07-17 22:52:24 |
合計ジャッジ時間 | 3,996 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
5,248 KB |
testcase_01 | AC | 2 ms
5,376 KB |
testcase_02 | AC | 2 ms
5,376 KB |
testcase_03 | AC | 2 ms
5,376 KB |
testcase_04 | AC | 2 ms
5,376 KB |
testcase_05 | AC | 1 ms
5,376 KB |
testcase_06 | AC | 2 ms
5,376 KB |
testcase_07 | AC | 2 ms
5,376 KB |
testcase_08 | AC | 2 ms
5,376 KB |
testcase_09 | AC | 2 ms
5,376 KB |
testcase_10 | AC | 16 ms
11,008 KB |
testcase_11 | AC | 18 ms
13,056 KB |
testcase_12 | AC | 10 ms
7,680 KB |
testcase_13 | AC | 18 ms
12,824 KB |
testcase_14 | AC | 15 ms
10,112 KB |
testcase_15 | AC | 10 ms
6,912 KB |
testcase_16 | AC | 16 ms
10,504 KB |
testcase_17 | AC | 2 ms
5,376 KB |
testcase_18 | AC | 2 ms
5,376 KB |
testcase_19 | AC | 12 ms
8,576 KB |
testcase_20 | AC | 17 ms
11,008 KB |
testcase_21 | AC | 22 ms
13,312 KB |
testcase_22 | AC | 19 ms
13,324 KB |
testcase_23 | AC | 9 ms
6,400 KB |
testcase_24 | MLE | - |
testcase_25 | -- | - |
testcase_26 | -- | - |
testcase_27 | -- | - |
testcase_28 | -- | - |
testcase_29 | -- | - |
testcase_30 | -- | - |
testcase_31 | -- | - |
testcase_32 | -- | - |
testcase_33 | -- | - |
ソースコード
#include <bits/stdc++.h> using namespace std; using lint = long long; template<class T = int> using V = vector<T>; template<class T = int> using VV = V< V<T> >; template<unsigned P> struct ModInt { using M = ModInt; unsigned v; ModInt() : v(0) {} template<class Z> ModInt(Z x) : v(x >= 0 ? x % P : -x % P ? P - -x % P : 0) {} constexpr ModInt(unsigned v, int) : v(v) {} static constexpr unsigned p() { return P; } M operator+() const { return *this; } M operator-() const { return {v ? P - v : 0, 0}; } explicit operator bool() const noexcept { return v; } bool operator!() const noexcept { return !(bool) *this; } M operator*(M r) const { return M(*this) *= r; } M operator/(M r) const { return M(*this) /= r; } M operator+(M r) const { return M(*this) += r; } M operator-(M r) const { return M(*this) -= r; } bool operator==(M r) const { return v == r.v; } bool operator!=(M r) const { return !(*this == r); } M& operator*=(M r) { v = (uint64_t) v * r.v % P; return *this; } M& operator/=(M r) { return *this *= r.inv(); } M& operator+=(M r) { v = r.v < P - v ? v + r.v : v - (P - r.v); return *this; } M& operator-=(M r) { v = r.v <= v ? v - r.v : v + (P - r.v); return *this; } M inv() const { int a = v, b = P, x = 1, u = 0; while (b) { int q = a / b; swap(a -= q * b, b); swap(x -= q * u, u); } assert(a == 1); return x; } template<class Z> M pow(Z n) const { n = n >= 0 ? n % (P - 1) : P - 1 - -n % (P - 1); M res = 1; for (M a = *this; n; a *= a, n >>= 1) if (n & 1) res *= a; return res; } template<class Z> friend M operator*(Z l, M r) { return M(l) *= r; } template<class Z> friend M operator/(Z l, M r) { return M(l) /= r; } template<class Z> friend M operator+(Z l, M r) { return M(l) += r; } template<class Z> friend M operator-(Z l, M r) { return M(l) -= r; } friend ostream& operator<<(ostream& os, M r) { return os << r.v; } friend istream& operator>>(istream& is, M& r) { lint x; is >> x; r = x; return is; } template<class Z> friend bool operator==(Z l, M r) { return M(l) == r; } template<class Z> friend bool operator!=(Z l, M r) { return !(l == r); } }; using Mint = ModInt<(unsigned) 1e9 + 7>; V<Mint> fact, ifact; void init(int n) { fact.resize(n + 1); fact[0] = 1; for (int i = 1; i <= n; ++i) { fact[i] = i * fact[i - 1]; } ifact.resize(n + 1); ifact[n] = fact[n].inv(); for (int i = n; i > 0; --i) { ifact[i - 1] = i * ifact[i]; } } Mint comb(int n, int r) { if (r < 0 or r > n) return 0; return fact[n] * ifact[r] * ifact[n - r]; } int main() { cin.tie(nullptr); ios::sync_with_stdio(false); int n, m, d0, d1; cin >> n >> m >> d0 >> d1, --m; if ((m -= d0 * (n - 1)) < 0) return cout << 0 << '\n', 0; init(n + m); Mint res; lint d = d1 - d0 + 1; for (int x = 0; x < n; ++x) { if (m - d * x < 0) break; Mint curr = comb(n - 1, x); curr *= comb(m - d * x + n, n); res += x & 1 ? -curr : curr; } cout << res << '\n'; }