JZOJ3400. 【GDOI2014模拟】旅行 (Standard IO)

Description
从前有一位旅者，他想要游遍天下所有的景点。这一天他来到了一个神奇的王国：在这片土地上，有n个城市，从1到n进行编号。王国中有m条道路，第i条道路连接着两个城市ai,bi，由于年代久远，所有的道路都已经不能使用。如果要修复第i条道路，需要wi的时间。为了更好的旅行，旅者想要将某些道路修复，使得1号城市能够到达n号城市，2号城市能够到达n-1号城市…k号城市能够到达n-k+1号城市。为了满足他的要求，请问最少需要多少时间去修复道路。无解请输出-1。
Input
第一行：n,m,k
接下来m行：ai,bi,wi
含义如上所述。
Output
输出共一行：最少需要多少时间修复道路。如果始终无法满足旅者的要求，请输出-1。
Sample Input
5 5 2
1 3 4
3 5 2
2 3 1
3 4 4
2 4 3
Sample Output
9
Data Constraint
20%的数据满足：k <= 2, n<= 10, m <= 20
40%的数据满足：k <= 3, n<=100, m<=1000
70%的数据满足：k<=4, n<=1000, m<=1000
100%的数据满足：k<=4, n<=10000, m<=10000, n >= 2*k, wi<= 1000, 1 <= ai, bi <= n

斯坦纳树入门

先用斯坦纳树处理出F[i][s]代表以i节点联通的树联通状态是S。

令good[s]=min{f[i][s]}，这样我们就知道了将S联通的最小代价

令ans[s]为答案，转移方程即为ans[s]=min{ans[t]+good[s-t]}

注意：若s状态中i与n-i+1不连通，我们删除这个状态。

这样可以保证最后求出的答案ans[n]是由合法状态转移而来

没啦。是不是So easy？

标程：

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
int i,j,k,m,n,o,p,l,s,t,x,y,z,r,ll;
int dp[10001][501],ans[501],good[501];
int f[100001][4],q[100001],h[100001],bz[10001];
void insert(int x,int y,int z)
{ f[++t][1]=y,f[t][2]=q[x],f[t][3]=z,q[x]=t; } 
int pd(int s)
{
	for (int i=1;i<=ll;i++)
	{
		int gg=s&(1<<(i-1)),hh=s&(1<<(ll+i-1));
		if ((!gg&&hh)||(!hh&&gg)) return 0;
	}
	return 1;
}
int main()
{
	freopen("travel.in","r",stdin);
	freopen("travel.out","w",stdout);
	scanf("%d%d%d",&n,&m,&ll);
	for (i=1;i<=m;i++) scanf("%d%d%d",&x,&y,&z),insert(x,y,z),insert(y,x,z);
	int N=(1<<(2*ll))-1;
	memset(dp,60,sizeof(dp));memset(good,60,sizeof(good));memset(ans,60,sizeof(ans));
	for (i=1;i<=ll;i++)
		dp[i][1<<(i-1)]=dp[n-i+1][1<<(ll+i-1)]=0;
	for (s=1;s<=N;s++)
	{
		for (i=1;i<=n;i++)
		{
			for (t=1;t<=s;t++)
				if ((t|s)==s)
				{
					for (k=q[i];k;k=f[k][2])
						dp[i][s]=min(dp[i][s],dp[i][t]+dp[i][s-t]);
				}
		}
		l=0,r=n;memset(bz,0,sizeof(bz));
		for (j=1;j<=n;j++) h[j]=j,bz[j]=1;
		while (l<r)
		{
			l++;
			for (int k=q[h[l]];k;k=f[k][2])
			{
				if (dp[h[l]][s]+f[k][3]<dp[f[k][1]][s])
				{
					dp[f[k][1]][s]=dp[h[l]][s]+f[k][3];
					if (!bz[f[k][1]]) bz[f[k][1]]=1,h[++r]=f[k][1];	
				}	
			}	
			bz[h[l]]=0;
		} 
		for (i=1;i<=n;i++) good[s]=min(good[s],dp[i][s]);
	}
	ans[0]=0;
	for (s=1;s<=N;s++)
	{
		if (s==N)
		{
			int gds=0;
		 } 
		for (t=0;t<s;t++)
			if ((t|s)==s&&pd(s)) ans[s]=min(ans[s],ans[t]+good[s-t]);
	}
	if (ans[N]>1e9) ans[N]=-1;
	printf("%d\n",ans[N]);
}