uva 11600 期望dp 状压dp

关键词：期望dp，存储已完成某项操作的所有编号——二进制（状压dp）
题意：一共有n个点，由m条双向道路连接，当你从1出发，每天可以走向2-n的任意一个目标点，下一天从目标点出发，再走向除该点之外的任意一点，如果行走的道路不存在，则在行走之后就自动生成。如果每次行走都是随机的，问所有点都连通的期望天数是多少？
经典题目：起点已知，经过所有的点，且顺序无要求—哈密顿路径！因此采取状态压缩dp进行求解
1.将已连通的点缩成一点，形成n1个连通子集。
2.dp[u][st]:现在处于第u个连通子集，st表示当前的连通状态（用二进制第i位为1表示第i个连通子集已加入“行走连通集”（简称未连通的子集））,dp值表示从当前状态开始还需要的期望步数
$dp[u][st]=[s/(n-1)]*(dp[u][st]+1)+sum{[cnt[v]/(n-1)]*（dp[v][st|(2^v)]+1） }$————<1>
注：求期望的公式：p=$\sum_{i=1}^{n}{p_i*x_i}$
下一个状态的选择有已连通的子集、未连通的子集
v表示枚举所有未连通的子集，cnt[v]表示第v个未连通子集的原节点个数
s为已连通子集的元素总数
由期望公式”x=pi*xi”可得到<1>式
此题不同的是求总步数的期望，还要加上转移状态花费的步数1
<1>式可化简为：
dp[u][st]=(n-1)/(n-s)+sum{ [cnt[v]/(n-s)]*dp[v][st|(2^v)] }———–<2>

拓展1：如果每条道路耗费天数不同，如何求解？
简化问题：一开始所有点都不连通，求从1点出发的期望哈密顿路径长度(每次走向不同于当前点的下一点
dp[u][st]=$\sum_{v=1}^{v=n}1/(n-1)*(dp[v][st|(2^v)]+d[u][v])*(v!=u)$

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<algorithm>
#include<map>
#include<vector>
#include<queue>
#include<stack>
#define ll long long
#define INF 0x3f3f3f3f
#define mem(a,b) memset(a,b,sizeof(a))
using namespace std;

const int maxn = 30+5;
map<int,double> d[maxn];
double ans;

int t,n,m;
vector<int> g[maxn];
int num[maxn],pointc;//cnt[i]：第i个连通分量;num[i]:i节点的连通分量编号;tmpc:当前连通分量编号
int cnt[maxn],cas;

void dfs(int u){//缩点子函数
    num[u]=pointc;
    cnt[pointc]=cnt[pointc]+1;
    for(int i=0;i<g[u].size();i++){
        if(num[g[u][i]]==-1) dfs(g[u][i]);
    }
    return ;
}

void con_point(){//将连通分量缩成一点
    pointc=0;
    for(int i=1;i<=n;i++){
        if(num[i]==-1) { dfs(i); pointc++; }
    }
}

void init(){//初始化
    mem(cnt,0);
    for(int i=0;i<maxn;i++) d[i].clear();
    for(int i=1;i<=n;i++) { g[i].clear(); num[i]=-1; }
}

double dp(int u,int st){//dp递归求解
    if(d[u].count(st)) return d[u][st];
    int s=0;
    for(int v=0;v<pointc;v++) if((1<<v)&st) s+=cnt[v];
    //******************
    if(s==n) { return d[u][st]=0; }// 不能忽略!!!!! 保证分母不为0!!!
    //******************
    double tmp=((n-1)*1.0)/((n-s)*1.0);
    for(int v=0;v<pointc;v++){
        if(!((1<<v)&st)){
             tmp+=dp(v,st|(1<<v))*(cnt[v]*1.0)/((n-s)*1.0);
        }
    }
    return d[u][st]=tmp;
}

int main(){
    //freopen("a.txt","r",stdin);
    scanf("%d",&t);
    cas=0;
    while(t--){
        scanf("%d%d",&n,&m);
        init();
        for(int i=1;i<=m;i++){
            int u,v;
            scanf("%d%d",&u,&v);
            g[u].push_back(v),g[v].push_back(u);
        }
        con_point();
        double ans=dp(0,1);
        printf("Case %d: %.6lf\n",++cas,ans);
    }
    return 0;
}