UVA 1390 Interconnect

https://vjudge.net/problem/UVA-1390

题意：

给出n个点m条边的无向图，

每次随机加一条非自环的边，（加完后可出现重边），

添加每条边的概率是相等的

求使图连通的期望添边次数

只关心图的连通状况，即连通块的个数和大小

所以可以用{a1，a2，a3……an} 表示状态（n个连通块，每个连通块大小为ai）

添加一条边后有两种可能

1、状态不变

2、状态变为 {a1，……ai+aj，……a_n-1}

将状态哈希

dp[i]表示哈希后为i的状态 添边至连通的期望次数

dp[i]= dp[i]*p + dp[k1]*p1 + dp[k2]*p2  + …… + dp[km]*pm

其中p表示状态不变的概率，p1=Σ C（ai，2）/ C（n，2）

pi表示装移到ki这种状态的概率 ，pi= （ai*aj） / C（n，2）

上述式子移项得状态转移方程：dp[i]=（dp[k1]*p1 + dp[k2]*p2  + …… + dp[km]*pm）/（1-p）

#include<cstdio>
#include<cstring>
#include<algorithm>
const int mod=;
int n,m;
int fa[],siz[];
struct sta
{
    int x[];
    bool flag;
    double val;
    void clear() { memset(x,,sizeof(x)); }
    void sort() { std::sort(x,x+); }
    int hashme()
    {
        int v=;
        for(int i=,b=;i&&x[i];i--)
        {
            v+=x[i]*b;
            v%=mod;
            b*=;
            b%=mod;
        }
        return v;
    }
    bool operator == (sta b)
    {
        for(int i=;i<;i++)
            if(x[i]!=b.x[i]) return false;
        return true;
    }
    bool operator != (sta b)
    {
        return *this == b ? false : true;
    }
}st,hash[mod];
int find(int i) { return fa[i]==i ? i : fa[i]=find(fa[i]); }
double gethash(sta st)
{
    int x=st.hashme();
    while(hash[x].flag && hash[x]!=st)
        if(++x==mod) x=;
    return hash[x]==st ? hash[x].val : -;
}
double inhash(sta st)
{
    int x=st.hashme();
    while(hash[x].flag)
        if(++x==mod) x=;
    hash[x]=st; hash[x].flag=true;
}
double dp(sta st)
{
    if(st.hashme()==n) return ;
    double x=gethash(st);
    if(x!=-) return x;
    double tmp=,ans=;
    for(int i=;i<;i++)
        tmp+=st.x[i]*(st.x[i]-)/;
    for(int i=;i<;i++)
        for(int j=i+;j<;j++)
        {
            if(!st.x[i] || !st.x[j]) continue;
            sta tmp=st;
            tmp.x[i]+=tmp.x[j];
            tmp.x[j]=;
            tmp.sort();
            ans+=st.x[i]*st.x[j]*dp(tmp);
        }
    ans/=n*(n-)/;
    ans++;
    ans/=-tmp/(n*(n-)/);
    st.val=ans;
    inhash(st);
    return ans;
}
int main()
{
    int u,v;
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        for(int i=;i<mod;++i) hash[i].flag=false;
        for(int i=;i<=n;i++) fa[i]=i;
        while(m--)
        {
            scanf("%d%d",&u,&v);
            fa[find(u)]=find(v);
        }
        st.clear();
        memset(siz,,sizeof(siz));
        for(int i=;i<=n;i++) siz[find(i)]++;
        int tot=;
        for(int i=;i<=n;i++)
            if(siz[i]) st.x[tot++]=siz[i];
        st.sort();
        printf("%.7lf\n",dp(st));
    }
}