LOJ #2318. 「NOIP2017」宝藏

题目描述

参与考古挖掘的小明得到了一份藏宝图，藏宝图上标出了 nnn 个深埋在地下的宝藏屋，也给出了这 nnn 个宝藏屋之间可供开发的 mmm 条道路和它们的长度。

小明决心亲自前往挖掘所有宝藏屋中的宝藏。但是，每个宝藏屋距离地面都很远，也就是说，从地面打通一条到某个宝藏屋的道路是很困难的，而开发宝藏屋之间的道路则相对容易很多。

小明的决心感动了考古挖掘的赞助商， 赞助商决定免费赞助他打通一条从地面到某个宝藏屋的通道，通往哪个宝藏屋则由小明来决定。

在此基础上， 小明还需要考虑如何开凿宝藏屋之间的道路。已经开凿出的道路可以任意通行不消耗代价。每开凿出一条新道路，小明就会与考古队一起挖掘出由该条道路所能到达的宝藏屋的宝藏。另外，小明不想开发无用道路，即两个已经被挖掘过的宝藏屋之间的道路无需再开发。

新开发一条道路的代价是：

这条道路的长度 ×\times× 从赞助商帮你打通的宝藏屋到这条道路起点的宝藏屋所经过的宝藏屋的数量（包括赞助商帮你打通的宝藏屋和这条道路起点的宝藏屋）。

请你编写程序为小明选定由赞助商打通的宝藏屋和之后开凿的道路，使得工程总代价最小，并输出这个最小值。

输入格式

第一行两个用空格分离的正整数 nnn 和 mmm，代表宝藏屋的个数和道路数。

接下来 mmm 行，每行三个用空格分离的正整数，分别是由一条道路连接的两个宝藏屋的编号（编号为 1∼n1\sim n1∼n），和这条道路的长度 vvv。

输出格式

输出共一行，一个正整数，表示最小的总代价。

样例

样例输入1

4 5
1 2 1
1 3 3
1 4 1
2 3 4
3 4 1

样例输出1

4

样例说明1

样例输入2

4 5
1 2 1
1 3 3
1 4 1
2 3 4
3 4 2

样例输出2

5

样例说明2

数据范围与提示

对于 20%20\%20% 的数据：

保证输入是一棵树， 1≤n≤8，v≤5000 且所有的 vvv 都相等。

对于 40%40\%40% 的数据：

1≤n≤8，0≤m≤1000，v≤5000 且所有的 vvv 都相等。

对于 70%70\%70% 的数据：

1≤n≤8，0≤m≤1000，v≤5000。

对于 100%100\%100% 的数据：

1≤n≤12，0≤m≤1000，v≤500000。

 

想都不想就是状压

3^n的写法

p【S】【i】表示点集为S，点i到点集S的最短距离

dp【S】【i】表示点集为S，深度为i的最小值

dp【S1】【i】=dp【S2】【i-1】+sum{ S1-S2,S2   }

//sum{ S1-S2,S2  }表示集合S1-S2每个点到集合S2的最小值的和

注意特判，当m=0时

#include<cstdio>
#include<iostream>
#define ll long long
using namespace std;
inline int read()
{
	int ret=0;
	char ch=getchar();
	while(ch<'0'||ch>'9') ch=getchar();
	while(ch>='0'&&ch<='9')
		ret=(ret<<1)+(ret<<3)+ch-'0',ch=getchar();
	return ret;
}

int n,m;
const int N=13,SS=(1<<12)+5;
int d[N][N],p[SS][N];
ll f[SS][N];

int main()
{
	n=read(),m=read();
	if(m==0) 
	{
		printf("0\n"); return 0;
	}
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)d[i][j]=1e9;
	
	for(int i=1;i<=m;i++)
	{
		int u=read(),v=read(),w=read();
		d[v][u]=d[u][v]=min(d[u][v],w);
	} 
	
	for(int S=1;S<(1<<n)-1;S++)
		for(int i=1;i<=n;i++)
			if(!((1<<(i-1))&S))  
			{
				p[S][i]=1e9;
				for(int j=1;j<=n;j++)
				if((1<<(j-1))&S)p[S][i]=min(p[S][i],d[i][j]);
			}
	
	ll ans=1e18;
	for(int i=1;i<=n;i++)
	{
		for(int S=0;S<(1<<n);S++)
			for(int j=0;j<=n;j++) f[S][j]=1e18;
				
		f[1<<(i-1)][0]=0;
		for(int S1=1;S1<(1<<n);S1++)
			for(int S2=S1&(S1-1);S2;S2=S1&(S2-1))
			{
				ll s=0;
				for(int j=1;j<=n;j++)
					if((1<<(j-1))&(S1-S2)) s+=p[S2][j];
				for(int j=1;j<=n;j++)
				f[S1][j]=min(f[S1][j],f[S2][j-1]+s*j);
			}
		for(int j=1;j<=n;j++)
			ans=min(ans,f[(1<<n)-1][j]);
	}
	printf("%lld\n",ans);
	return 0;
}