Minimal Steiner Tree

最新推荐文章于 2020-06-20 23:04:23 发布

lwbaptx

最新推荐文章于 2020-06-20 23:04:23 发布

阅读量3.7k

点赞数 1

分类专栏：算法研究文章标签： tree ini

本文链接：https://blog.csdn.net/lwbaptx/article/details/7297840

版权

算法研究专栏收录该内容

3 篇文章 0 订阅

订阅专栏

Minimal Steiner Tree 生成树问题。

定义：G(V,E),A是V的一个子集，求至少包含A中所有点的最小子树。

用动态规划解决。

推荐两道题：

（1） poj 3123

（2） Garden

//Minimal Steiner Tree
//G(V,E),A是V的一个子集，求至少包含A中所有点的最小子树
//时间复杂度 O(V^3+V*2^A*(2^A+V))
//INIT: d[][]距离矩阵; id[]置为集合A中点的标号;
//CALL: steiner(int n, int a);

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstdlib>
#include <string>
#include <cmath>
#include <vector>
using namespace std;

const int inf = 2147483647;
const int N = 205; //总顶点数
const int A = 7 ;  //要包含的顶点数
int id[A];         //包含顶点的序号
int d[N][N];       //距离矩阵
int dp[1<<A][N];  //dp[i][j]表示点j到联通的集合i的最短距离

//返回最小值
int Steiner(int n,int a)
{
	int top = 1<<a;
	for(int k=0;k<n;k++)      //vertex 0---n-1
		for(int i=0;i<n;i++)
			for(int j=0;j<n;j++)
				if(d[i][j]>d[i][k]+d[k][j])
					d[i][j]=d[i][k]+d[k][j];

	for(int i=0;i<top;i++)
		for(int j=0;j<n;j++)
			dp[i][j]=inf;

	for(int i=0;i<a;i++)
	{
		for(int j=0;j<n;j++)
			dp[1<<i][j]=d[id[i]][j];
	}
	for(int i=1;i<top;i++)
	{
		if(i&(i-1)==0)
			continue;
		for(int k=0;k<n;k++)   //init
		{
			for(int j=1;j<i;j++)
			{
				if((i|j)==i && dp[i][k]>dp[j][k]+dp[i-j][k])
					dp[i][k]=dp[j][k]+dp[i-j][k];
			}
		}

		for(int j=0;j<n;j++)  //update
		{
			for(int k=0;k<n;k++)
			{
				if(dp[i][j]>dp[i][k]+d[k][j])
					dp[i][j]=dp[i][k]+d[k][j];
			}
		}
	}

	int ans = inf;
	for(int i=0;i<n;i++)
	{
		if(ans > dp[top-1][i])
			ans=dp[top-1][i];
	}
	return ans;
}