NOIP2017 DAY2 T1 宝藏

NOIP 2017 DAY2 T1 宝藏

戳这里

题目描述

有n个宝藏，需要找到最小的代价使得所有的宝藏之间联通，当然赞助商会帮助我们免费打通一条道路。

---

思路

很明显搜索当然是可以的，哪些路开通，哪些路不开通。对于免费的，外面再加一层搜索就可以了。当然这样子会超时。
对于数据： n<=12;m<=1000。很显然可以状态压缩一下，状态用哪些选哪些没选。

状态压缩之后，就可以直接 搜索 + 状压dp。f数组来存储选了某一个状态下所最小的代价。G数组来存储宝藏之间是否可以
开通以及开通所需要的代价。

---

#include<cstdio>
#include<cctype>
#include<cstring>
#include<algorithm>
using namespace std;
#define inf 2147483647
#define maxn 2000
inline int read() {
	int x = 0;char c=getchar();
	while(!isdigit(c)) c=getchar();
	while(isdigit(c)) x=(x<<3)+(x<<1)+c-48,c=getchar();return x;
}

int n,m,G[maxn][maxn],dis[maxn];
int f[10000];

void dfs(int s) {
	
	for(int i=1;i<=n;i++) {
		if((1<<(i-1)) & s) {
			for(int j=1;j<=n;j++) {
				if(((1<<(j-1)) & s) == 0 && G[i][j] != inf) {
					if(f[s|(1<<(j-1))] > f[s] + dis[i]*G[i][j]) {
						int tmp = dis[j];
						dis[j] = dis[i] + 1;
						f[s|1<<(j-1)] = f[s] + dis[i]*G[i][j];
						dfs(s|(1<<(j-1)));
						dis[j] = tmp;
					}
				}
			}
		}
	}
}
int main() {
	//freopen("test.in","r",stdin);
	
	n = read(),m = read();
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			G[i][j] = inf;
	for(int i=1;i<=m;i++) {
		int x = read(),y = read(),w = read();
		G[x][y] = min(G[x][y],w),G[y][x] = min(G[x][y],w);
	}
	
	int ans = inf;
	for(int i=1;i<=n;i++) {
		for(int j=1;j<=n;j++) dis[j] = inf;
		for(int j=1;j<=(1<<n)-1;j++) f[j] = inf;
		
		dis[i] = 1;
		f[1<<(i-1)] = 0;
		dfs(1<<(i-1));
		ans = min(ans,f[(1<<n)-1]); 		
	}
	printf("%d",ans);
	return 0;
}
/*
690550
*/