結果

問題 No.802 だいたい等差数列
コンテスト
ユーザー vjudge1
提出日時 2026-06-19 11:35:46
言語 C++14
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=c++14 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
WA  
実行時間 -
コード長 1,166 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,440 ms
コンパイル使用メモリ 181,712 KB
実行使用メモリ 19,200 KB
最終ジャッジ日時 2026-06-19 11:36:01
合計ジャッジ時間 6,441 ms
ジャッジサーバーID
（参考情報） 		judge2_1 / judge3_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1 WA * 3
other AC * 2 WA * 11 RE * 17
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
using namespace std;
const int mod = 1e9 + 7;
const int N = 1e6 + 10;
long long f[N], g[N];
long long ksm(long long a, long long b) {
    long long r = 1;
    while (b) {
        if (b & 1) r = r * a % mod;
        a = a * a % mod;
        b >>= 1;
    }
    return r;
}
void init(int n) {
    f[0] = 1;
    for (int i = 1; i <= n; i++) f[i] = f[i - 1] * i % mod;
    g[n] = ksm(f[n], mod - 2);
    for (int i = n - 1; i >= 0; i--) g[i] = g[i + 1] * (i + 1) % mod;
}
long long C(int n, int m) {
    if (m < 0 || m > n) return 0;
    return f[n] * g[m] % mod * g[n - m] % mod;
}
int main() {
    int n, m, d1, d2;
    cin >> n >> m >> d1 >> d2;
    if (d1 > d2) {
        cout << 0 << endl;
        return 0;
    }
    if (d1 == 0 && d2 >= m) {
        cout << C(m + n - 1, n) << endl;
        return 0;
    }
    init(m + n);
    long long ans = 0;
    int k = (m - n * d1 + d2 - 1) / d2;
    for (int i = 0; i <= k; i++) {
        int t = m - n * d1 - i * d2;
        if (t < 0) break;
        long long tmp = C(t + n, n);
        if (i & 1) tmp = mod - tmp;
        ans = (ans + tmp) % mod;
    }
    cout << ans << endl;;;
    return 0;
}
0