結果

問題 No.125 悪の花弁
ユーザー Iroha_3856
提出日時 2025-02-26 20:10:08
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 155 ms / 5,000 ms
コード長 1,446 bytes
コンパイル時間 3,824 ms
コンパイル使用メモリ 282,732 KB
実行使用メモリ 15,652 KB
最終ジャッジ日時 2025-02-26 20:10:14
合計ジャッジ時間 5,060 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 6
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#include <atcoder/modint>
using mint = atcoder::modint1000000007;

#define rep(i, l, r) for (int i = (int)(l); i<(int)(r); i++)
#define ll long long

struct Combination {
    vector<mint> fac, infac;
    Combination(int n) {
        fac.resize(n+1);
        infac.resize(n+1);
        fac[0] = 1;
        for (int i = 1; i <= n; i++) fac[i] = fac[i-1] * i;
        infac[n] = fac[n].inv();
        for (int i = n; i >= 1; i--) infac[i-1] = infac[i] * i;
    }
    mint operator()(int n, int k) {
        if (k < 0 || k > n) return 0;
        return fac[n] * infac[k] * infac[n-k];
    }
};

void solve() {
    int K; cin >> K;
    vector<int> C(K);
    rep(i, 0, K) cin >> C[i];
    int N = 0;
    rep(i, 0, K) N += C[i];
    Combination comb(N);
    mint ans = 0;
    vector<int> cnt(N+1);
    for (int d = 1; d <= N; d++) {
        int t = gcd(N, d);
        cnt[t]++;
    }
    for (int d = 1; d <= N; d++) {
        if (N%d) continue;
        vector<int> D(K);
        bool f = true;
        rep(i, 0, K) {
            if (C[i]%(N/d) != 0) {
                f = false;
            }
            D[i] = C[i]/(N/d);
        }
        if (f) {
            mint plus = comb.fac[d];
            rep(i, 0, K) plus *= comb.infac[D[i]];
            ans += plus * cnt[d];
        }
    }
    ans /= N;
    cout << ans.val() << endl;
}

int main() {
    int T = 1;
    //cin >> T;
    while(T--) solve();
}
0