luogu P6192 【模板】最小斯坦纳树

代码:

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>

using namespace std;

const int N=109,INF=1<<30;
int n,m,k,head[N],cnt,f[N][N*10],Max,p[N],vis[N];
priority_queue <pair<int,int> > q;
struct Edge
{
	int nxt,to,w;
}g[N*10];

void add(int from,int to,int w)
{
	g[++cnt].nxt=head[from];
	g[cnt].to=to;
	g[cnt].w=w;
	head[from]=cnt;
}

void init()
{
	scanf("%d %d %d",&n,&m,&k);
	for (int x,y,z,i=1;i<=m;i++)
	{
		scanf("%d %d %d",&x,&y,&z);
		add(x,y,z),add(y,x,z);
	}
	Max=1<<k;
	for (int i=1;i<=n;i++)
		for (int j=0;j<Max;j++)
			f[i][j]=INF;
	for (int i=1;i<=k;i++)
	{
		scanf("%d",&p[i]);
		f[p[i]][1<<i-1]=0;
	}
}

void Dijkstra(int s)
{
	memset(vis,0,sizeof(vis));
	while(!q.empty())
	{
		int x=q.top().second;q.pop();
		if(vis[x]) continue;
		vis[x]=1;
		for (int i=head[x];i;i=g[i].nxt)
		{
			int v=g[i].to;
			if(f[v][s]>f[x][s]+g[i].w)
			{
				f[v][s]=f[x][s]+g[i].w;
				q.push(make_pair(-f[v][s],v));
			}
		}
	}
}

void work()
{
	for (int s=1;s<Max;s++)
	{
		for (int i=1;i<=n;i++)
		{
			for (int sub=s&(s-1);sub;sub=(sub-1)&s)
				f[i][s]=min(f[i][s],f[i][sub]+f[i][s^sub]);
			if(f[i][s]!=INF) q.push(make_pair(-f[i][s],i));
		}
		Dijkstra(s);
	}
	printf("%d\n",f[p[1]][Max-1]);
}

int main()
{
	init();
	work();
	return 0;
}
posted @ 2020-07-12 21:59  With_penguin  阅读(165)  评论(0编辑  收藏  举报