結果

問題 No.802 だいたい等差数列
ユーザー E31415926E31415926
提出日時 2019-03-19 23:56:30
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 1,996 bytes
コンパイル時間 3,173 ms
コンパイル使用メモリ 156,672 KB
実行使用メモリ 27,104 KB
最終ジャッジ日時 2023-10-12 15:28:57
合計ジャッジ時間 5,974 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 33 ms
26,804 KB
testcase_01 AC 33 ms
26,784 KB
testcase_02 WA -
testcase_03 AC 33 ms
26,752 KB
testcase_04 WA -
testcase_05 AC 33 ms
26,912 KB
testcase_06 AC 34 ms
26,796 KB
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 AC 36 ms
26,744 KB
testcase_15 WA -
testcase_16 AC 36 ms
27,028 KB
testcase_17 AC 35 ms
26,832 KB
testcase_18 AC 35 ms
26,796 KB
testcase_19 AC 35 ms
26,836 KB
testcase_20 AC 36 ms
26,796 KB
testcase_21 WA -
testcase_22 AC 34 ms
26,852 KB
testcase_23 AC 36 ms
26,800 KB
testcase_24 AC 35 ms
26,796 KB
testcase_25 AC 33 ms
26,796 KB
testcase_26 AC 35 ms
26,800 KB
testcase_27 AC 35 ms
26,740 KB
testcase_28 AC 36 ms
26,740 KB
testcase_29 AC 37 ms
26,740 KB
testcase_30 AC 34 ms
26,888 KB
testcase_31 AC 33 ms
26,780 KB
testcase_32 AC 34 ms
26,796 KB
testcase_33 AC 34 ms
26,784 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<vector>
#include<algorithm>
#include<queue>
#include<stack>
#include<cstring>
#include<utility>
#include<cmath>
#include<assert.h>
#include<set>
#include<map>
#include<unordered_set>
#include<unordered_map>
#include<complex>

#define int long long
using namespace std;
#define rep(i, n) for(int i=0;i<(n);++i)
typedef pair<int, int> pii;
const int INF = 1l << 60;
#define u_b upper_bound
#define l_b lower_bound

//const int mod = //998244353;
const int mod = 1000000007;
const int FacN = 1500000;//使わないなら0

constexpr int power(int a, int b) {
    int res = 1;
    while (b) {
        if (b & 1)res *= a;
        a *= a;
        res %= mod;
        a %= mod;
        b /= 2;
    }
    return res;
}

constexpr int inv(int x) {
    return power(x, mod - 2);
}

const struct Fac {
    int val[FacN];

    constexpr Fac() : val() {
        if (FacN == 0)return;
        val[0] = 1;
        for (int i = 1; i < FacN; ++i) val[i] = val[i - 1] * i % mod;
    }

    const int &operator[](size_t i) const {
        return val[i];
    }
} fac;

const struct FacInv {
    int val[FacN];

    constexpr FacInv() : val() {
        if (FacN == 0)return;
        val[FacN - 1] = 1;
        for (int i = 1; i < FacN; ++i)val[FacN - 1] = val[FacN - 1] * i % mod;
        val[FacN - 1] = inv(val[FacN - 1]);
        for (int i = FacN - 1; i >= 1; --i)val[i - 1] = val[i] * i % mod;
    }

    const int &operator[](size_t i) const {
        return val[i];
    }
} fac_inv;

int Perm(int n, int r) {
    return fac[n] * fac_inv[n - r] % mod;
}

int Comb(int n, int r) {
    if (n < r)return 0;
    return Perm(n, r) * fac_inv[r] % mod;
}

int N, M, D1, D2;

signed main() {
    cin >> N >> M >> D1 >> D2;
    int ans = 0;
    rep(i, N) {
        int c = Comb(N - 1, i);
        if (i % 2 == 1)c *= -1;
        int k = D1 * (N - 1 - i) + (D2 + 1) * i + 1;
        c *= Comb(M - k + N, N);
        ans += c;
        ans %= mod;
    }
    cout << ans << endl;
    return 0;
}
0