結果

問題 No.122 傾向と対策:門松列(その3)
ユーザー はまやんはまやんはまやんはまやん
提出日時 2017-04-18 16:41:17
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 54 ms / 5,000 ms
コード長 2,537 bytes
コンパイル時間 1,551 ms
コンパイル使用メモリ 173,204 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-07-19 07:38:07
合計ジャッジ時間 2,390 ms
ジャッジサーバーID
(参考情報)
judge1 / judge3
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 36 ms
5,248 KB
testcase_01 AC 44 ms
5,248 KB
testcase_02 AC 53 ms
5,376 KB
testcase_03 AC 45 ms
5,376 KB
testcase_04 AC 45 ms
5,376 KB
testcase_05 AC 54 ms
5,376 KB
testcase_06 AC 37 ms
5,376 KB
testcase_07 AC 36 ms
5,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include<bits/stdc++.h>
using namespace std;
#define rep(i,a,b) for(int i=a;i<b;i++)


int mod = 1000000007;
int add(int x, int y) { return (x += y) >= mod ? x - mod : x; }
template<class... T> int add(int x, T... y) { return add(x, add(y...)); }
int mul(int x, int y) { return 1LL * x * y % mod; }
template<class... T> int mul(int x, T... y) { return mul(x, mul(y...)); }
int sub(int x, int y) { return add(x, mod - y); }
int modpow(int a, long long b) { int ret = 1; while (b > 0) { if (b & 1) ret = 1LL * ret * a % mod; a = 1LL * a * a % mod; b >>= 1; } return ret; }
int modinv(int a) { return modpow(a, mod - 2); }
int fac[201010], ifac[201010];
void initfac() { fac[0] = ifac[0] = 1; rep(i, 1, 201010) ifac[i] = modinv(fac[i] = mul(i, fac[i - 1])); }
int aCb(int a, int b) { if (b < 0 || a < b) return 0; return mul(fac[a], ifac[a - b], ifac[b]); }
//-----------------------------------------------------------------------------------
int MIN[7], MAX[7];
//-----------------------------------------------------------------------------------
int dp[8][20101];
int count(vector<int> v) {
    rep(i, 0, 8) rep(j, 0, 20101) dp[i][j] = 0;

    dp[0][0] = 1;
    rep(i, 1, 8) {
        int vi = v[i - 1];
        int sm = 0;
        rep(j, 0, MIN[vi]) sm = add(sm, dp[i - 1][j]);
        rep(j, MIN[vi], MAX[vi] + 1) {
            dp[i][j] = sm, sm = add(sm, dp[i - 1][j]);
        }
    }

    int res = 0;
    rep(j, 0, 20101) res = add(res, dp[7][j]);
    return res;
}
//-----------------------------------------------------------------------------------
int main() {
    rep(i, 0, 7) scanf("%d%d", &MIN[i], &MAX[i]);

    int ans = 0;

    // max(a,c,e,g) < min(b,d,f)
    {
        vector<int> a({ 0,2,4,6 }), b({ 1,3,5 });
        do {
            do {
                vector<int> v;
                copy(a.begin(), a.end(), back_inserter(v));
                copy(b.begin(), b.end(), back_inserter(v));
                ans = add(ans, count(v));
            } while (next_permutation(b.begin(), b.end()));
        } while (next_permutation(a.begin(), a.end()));
    }

    // max(b,d,f) < min(a,c,e,g)
    {
        vector<int> a({ 0,2,4,6 }), b({ 1,3,5 });
        do {
            do {
                vector<int> v;
                copy(b.begin(), b.end(), back_inserter(v));
                copy(a.begin(), a.end(), back_inserter(v));
                ans = add(ans, count(v));
            } while (next_permutation(b.begin(), b.end()));
        } while (next_permutation(a.begin(), a.end()));
    }

    cout << ans << endl;
}
0