結果

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

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
template<class T>bool chmax(T &a, const T &b) { if (a<b) { a=b; return true; } return false; }
template<class T>bool chmin(T &a, const T &b) { if (b<a) { a=b; return true; } return false; }
#define vi vector<int>
#define vl vector<ll>
#define vii vector<pair<int,int>>
#define vll vector<pair<ll,ll>>
#define vvi vector<vector<int>>
#define vvl vector<vector<ll>>
#define vvii vector<vector<pair<int,int>>>
#define vvll vector<vector<pair<ll,ll>>>
#define vst vector<string>
#define pii pair<int,int>
#define pll pair<ll,ll>
#define pb push_back
#define all(x) (x).begin(),(x).end()
#define mkunique(x) sort(all(x));(x).erase(unique(all(x)),(x).end())
#define fi first
#define se second
#define mp make_pair
#define si(x) int(x.size())
const int mod=998244353,MAX=300005,INF=15<<26;

struct UF{
    int n;
    vector<int> par,size,edge;
    
    void init(int n_){
        n=n_;
        par.assign(n,-1);
        size.assign(n,1);
        edge.assign(n,0);
        
        for(int i=0;i<n;i++){
            par[i]=i;
        }
    }
    
    int root(int a){
        if(par[a]==a) return a;
        else return par[a]=root(par[a]);
    }
    
    void unite(int a,int b){
        edge[root(a)]++;
        if(root(a)!=root(b)){
            size[root(a)]+=size[root(b)];
            edge[root(a)]+=edge[root(b)];
            par[root(b)]=root(a);
        }
    }
    
    bool check(int a,int b){
        return root(a)==root(b);
    }
};

map<vi,double> MA;

double solve(vi X){
    if(MA.count(X)) return MA[X];
    int ng=0;
    double sum=1;
    for(int i=0;i<6;i++){
        for(int j=i;j>=0;j--){
            if(X[j]){
                X[j]--;
                sum+=solve(X)/6;
                X[j]++;
                break;
            }
            if(j==0){
                ng++;
            }
        }
    }
    
    if(ng==6) return MA[X]=0;
    return MA[X]=sum*6/(6-ng);
}

int main(){
    
    std::ifstream in("text.txt");
    std::cin.rdbuf(in.rdbuf());
    cin.tie(0);
    ios::sync_with_stdio(false);
    
    int N,M;cin>>N>>M;
    UF uf;uf.init(N);
    vvii S(6);
    for(int i=0;i<M;i++){
        int a,b,c;cin>>a>>b>>c;a--;b--;c--;
        S[c].pb(mp(a,b));
    }
    vi CN(6);
    for(int i=0;i<6;i++){
        for(auto [a,b]:S[i]){
            if(!uf.check(a,b)){
                uf.unite(a,b);
                CN[i]++;
            }
        }
    }
    cout<<fixed<<setprecision(25)<<solve(CN)<<endl;
}
0