結果

問題 No.3079 Unite Japanese Prefectures
コンテスト
ユーザー Iroha_3856
提出日時 2025-03-19 17:19:42
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 2,488 bytes
コンパイル時間 3,643 ms
コンパイル使用メモリ 298,472 KB
実行使用メモリ 524,752 KB
最終ジャッジ日時 2025-03-19 17:19:53
合計ジャッジ時間 10,988 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 11 WA * 1 MLE * 1 -- * 14
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#include <atcoder/dsu>

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

ll CurrentHash;

struct PolynomialHash {
    size_t operator()(const array<int, 6>& seq) const {
        return CurrentHash;
    }
};
struct ArrayEqual {
    bool operator()(const array<int, 6>& a, const array<int, 6>& b) const {
        return a == b;
    }
};

int main() {
    int N, M; cin >> N >> M;

    vector<tuple<int, int, int>> Edges(M);
    rep(i, 0, M) {
        int u, v, c; cin >> u >> v >> c;
        u--; v--; c--;
        Edges[i] = {c, u, v};
    }

    array<int, 6> cnt = {0, 0, 0, 0, 0, 0};
    atcoder::dsu d(N);

    sort(Edges.begin(), Edges.end());
    for (int i = 0; i < M; i++) {
        auto[c, u, v] = Edges[i];
        if (!d.same(u, v)) {
            d.merge(u, v);
            cnt[c]++;
        }
    }

    //累積
    for (int i = 4; i >= 0; i--) cnt[i] += cnt[i+1];

    int MOD = 998244353;
    ll power[6];
    power[0] = 1;
    rep(i, 1, 6) power[i] = power[i-1] * 53 % MOD;

    //dp[t] : 今まででたのがtとなるようにするまでにかかるコストの期待値
    unordered_map<array<int, 6>, double, PolynomialHash, ArrayEqual> dp;
    //6, 5, 4, 4, 3, 2 -> 7, 6, 4, 4, 3, 2 : 辞書順で後にいく
    array<int, 6> A;
    for (A[0] = cnt[0]; A[0] >= 0; A[0]--) for (A[1] = min(A[0], cnt[1]); A[1] >= 0; A[1]--) for (A[2] = min(A[1], cnt[2]); A[2] >= 0; A[2]--) {
        for (A[3] = min(A[2], cnt[3]); A[3] >= 0; A[3]--) for (A[4] = min(A[3], cnt[4]); A[4] >= 0; A[4]--) for (A[5] = min(A[4], cnt[5]); A[5] >= 0; A[5]--) {

            double t = 0;
            int loop = 0;
            array<int, 6> B = A;
            CurrentHash = 0;
            rep(i, 0, 6) {
                CurrentHash += power[i]*A[i];
                CurrentHash %= MOD;
            }
            ll AHash = CurrentHash;

            bool f = false;
            rep(i, 0, 6) {
            	if (A[i] != cnt[i]) {
                    f = true;
                    CurrentHash += power[i];
            	    B[i]++;
                }
            	if (f) t += dp[B];
                else loop++;
            }
            CurrentHash = AHash;
            if (loop == 6) {
                dp[A] = 0;
                continue;
            }
            t += 6;
            t /= (6 - loop);
            dp[A] = t;
        }
    }
    cout << fixed << setprecision(15) << dp[{0, 0, 0, 0, 0, 0}] << endl;
}
0