結果

問題 No.3079 Unite Japanese Prefectures
コンテスト
ユーザー umimel
提出日時 2025-03-28 22:34:02
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,892 ms / 4,000 ms
コード長 2,519 bytes
コンパイル時間 2,194 ms
コンパイル使用メモリ 187,628 KB
実行使用メモリ 56,320 KB
最終ジャッジ日時 2025-03-28 22:34:18
合計ジャッジ時間 15,983 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 27
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘void solve()’:
main.cpp:65:14: warning: structured bindings only available with ‘-std=c++17’ or ‘-std=gnu++17’ [-Wc++17-extensions]
   65 |         auto [c, a, b] = es[i];
      |              ^

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;
using ll = long long;
using pll = pair<ll, ll>;
#define all(a) (a).begin(), (a).end()
#define pb push_back
#define fi first
#define se second
mt19937_64 rng(chrono::system_clock::now().time_since_epoch().count());
const ll MOD1000000007 = 1000000007;
const ll MOD998244353 = 998244353;
const ll MOD[3] = {999727999, 1070777777, 1000000007};
const ll LINF = 1LL << 60LL;
const int IINF = (1 << 30) - 1;


struct union_find{
    vector<int> par;
    vector<int> siz;
    
    union_find(int n) : par(n), siz(n, 1){
        for(int i=0; i<n; i++) par[i] = i;
    }
    
    int root(int x){
        if (par[x] == x) return x;
        return par[x] = root(par[x]);
    }
 
    void unite(int x, int y){
        int rx = root(x);
        int ry = root(y);
        if (rx == ry) return;
        if (siz[rx] < siz[ry]) swap(rx, ry);
        siz[rx] += siz[ry];
        par[ry] = rx;
    }
 
    bool same(int x, int y){
        int rx = root(x);
        int ry = root(y);
        return rx == ry;
    }
 
    int size(int x){
        return siz[root(x)];
    }
};

using ld = long double;

void solve(){
    int n, m; cin >> n >> m;
    vector<tuple<int, int, int>> es(m);
    for (int i = 0; i < m; i++) {
        int a, b, c; cin >> a >> b >> c;
        a--; b--; c--;
        es[i] = {c, a, b};
    }
    sort(all(es));
    union_find uf(n);
    vector<int> init_cnt(6, 0);
    for(int i=0; i<m; i++){
        auto [c, a, b] = es[i];
        if(uf.same(a, b)) continue;
        uf.unite(a, b);
        init_cnt[c]++;
    }

    map<vector<int>, ld> dp;
    dp[{0, 0, 0, 0, 0, 0}] = 0;
    function<void(vector<int>)> dfs = [&](vector<int> cnt){
        if(dp.find(cnt) != dp.end()) return;
        vector<int> pre(6, 0);
        pre[0] = -1;
        if(cnt[0] > 0) pre[0] = 0;
        for(int i=1; i<6; i++){
            pre[i] = pre[i-1];
            if(cnt[i] > 0) pre[i] = i;
        }
        int c = 0;
        dp[cnt] = 1;
        for(int i=0; i<6; i++){
            if(pre[i] == -1){
                c++;
            }else{
                vector<int> nxt = cnt;
                nxt[pre[i]]--;
                dfs(nxt);
                dp[cnt] += dp[nxt]/(ld)6;
                nxt[pre[i]]++;
            }
        }
        dp[cnt] /= (1 - (ld)c/(ld)6);
    }; dfs(init_cnt);

    cout << fixed << setprecision(10) << dp[init_cnt] << '\n';
}

int main(){
    cin.tie(nullptr);
    ios::sync_with_stdio(false);
    
    int T=1;
    //cin >> T;
    while(T--) solve();
}
0