Interconnect UVA - 1390

使用map来记录状态，vector<int> p(31,0)中p[i]代表的是节点的个数为i的子联通图的个数，map建立的是状态到概率的映射。首先对初始的状态进行计算并且记录，然后开始进行深搜，每次判断，对于i个节点的联通图，如果存在，那么手首先这些联通图内部的节点与节点之间进行互联可以得到一个概率，如果i个节点的联通图的个数大于等于2，那么这些联通图可以相互之间进行连接，得到的是节点个数为2*i的联通图，然后再进行递归求解，同时从i+1个节点到最后，如果对应的节点数的联通图存在，那么相互之间也可以进行连接，继续递归计算，最后算出期望并且返回，具体实现见如下代码：

#include<iostream>
#include<vector>
#include<string>
#include<set>
#include<stack>
#include<queue>
#include<map>
#include<algorithm>
#include<cmath>
#include<iomanip>
#include<cstring>
#include<sstream>
#include<cstdio>
#include<deque>
#include<functional>
using namespace std;

const int maxn = 35;
int n, m;
bool visit[maxn];
vector<int> road[maxn];
map<vector<int>, double> record;
double base;

int dfs(int ind){
	visit[ind] = 1;
	int res = 1;
	for (int i = 0; i < road[ind].size(); i++){
		if (!visit[road[ind][i]]){
			res += dfs(road[ind][i]);
		}
	}
	return res;
}

double solve(vector<int> state){
	if (record.find(state) != record.end()) return record[state];
	double dis = 0, con = 0;
	for (int i = 1; i < state.size(); i++){
		if (state[i]){
			dis += 1.0*i*(i-1)/2*state[i] / base;
			if (state[i] >= 2){
				vector<int> temp(state);
				temp[i] -= 2;
				temp[i * 2] += 1;
				con += 1.0*i*i*state[i] * (state[i] - 1) / 2 / base*(solve(temp) + 1);
			}
			for (int j = i + 1; j < state.size(); j++){
				if (state[j]){
					vector<int> temp(state);
					temp[i]--;
					temp[j]--;
					temp[i + j]++;
					con += 1.0*i*j*state[i] * state[j] / base*(solve(temp) + 1);
				}
				
			}
		}
	}
	return record[state] = (dis + con) / (1 - dis);
}

int main(){
	while (cin >> n >> m){
		base = n*(n - 1) / 2.0;
		memset(visit,0,sizeof(visit));
		vector<int> connect(31,0);
		for (int i = 0; i < maxn; i++) road[i].clear();
		record.clear();
		connect[n] = 1;
		record[connect] = 0;
		connect[n] = 0;
		for (int i = 0; i < m; i++){
			int a, b;
			cin >> a >> b;
			road[a].push_back(b);
			road[b].push_back(a);
		}
		for (int i = 1; i <= n; i++){
			if (!visit[i]){
				int t = dfs(i);
				connect[t]++;
			}
		}
		cout <<fixed<< setprecision(6) << solve(connect) << endl;
	}
	return 0;
}