結果

問題 No.3079 Unite Japanese Prefectures
ユーザー GOTKAKO
提出日時 2025-03-28 23:04:07
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,332 ms / 4,000 ms
コード長 1,900 bytes
コンパイル時間 2,592 ms
コンパイル使用メモリ 218,068 KB
実行使用メモリ 49,792 KB
最終ジャッジ日時 2025-03-28 23:04:26
合計ジャッジ時間 12,470 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 27
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

class UnionFind{
    private:
    vector<int> par,siz;
    public:
    UnionFind(int N){
        par.resize(N,-1);
        siz.resize(N,1);
    }
 
    int root(int x){ //連結成分の代表頂点を返す.
        if(par.at(x) == -1) return x;
        else return par.at(x) = root(par.at(x));
    }
    bool unite(int u, int v){ //u,vを連結する 連結してた->false,した->trueを返す.
        u = root(u),v = root(v);
        if(u == v) return false;
 
        if(siz.at(u) < siz.at(v)) swap(u,v); //Union by size.
        par.at(v) = u;
        siz.at(u) += siz.at(v);
        return true;
    }
    bool issame(int u, int v){ //同じ連結成分ならtrue.
        if(root(u) == root(v)) return true;
        else return false;
    }
    int size(int pos){return siz.at(root(pos));} //posの連結成分の大きさを返す.
};

int main(){
    ios_base::sync_with_stdio(false);
    cin.tie(nullptr);

    int N,M; cin >> N >> M;
    vector<vector<pair<int,int>>> UV(6);
    while(M--){
        int u,v,c; cin >> u >> v >> c; u--; v--; c--;
        UV.at(c).push_back({u,v});
    }

    UnionFind Z(N);
    vector<int> C(6);
    for(int i=0; i<6; i++){
        for(auto [u,v] : UV.at(i)) if(Z.issame(u,v) == false) Z.unite(u,v),C.at(i)++;
    } 

    map<vector<int>,double> memo;
    memo[{0,0,0,0,0,0}] = 0;
    auto solve = [&](auto solve,vector<int> &A) -> double {
        if(memo.count(A)) return memo[A];
        int zero = 0,use = -1;
        double ret = 0;
        for(int i=0; i<6; i++){
            if(A.at(i)) use = i;
            if(use == -1){zero++; continue;}

            A.at(use)--; 
            ret += solve(solve,A);
            A.at(use)++; 
        }
        ret /= 6-zero; ret += 6/(6.0-zero);
        memo[A] = ret;
        return ret;
    };
    cout << fixed << setprecision(20) << solve(solve,C) << endl;
}
0