[图论]-------斯坦纳树

最新推荐文章于 2023-07-16 21:03:48 发布

T_a_r_j_a_n

最新推荐文章于 2023-07-16 21:03:48 发布

阅读量4.9k

点赞数 2

分类专栏：图论

本文链接：https://blog.csdn.net/Berserker____/article/details/115637081

版权

图论专栏收录该内容

16 篇文章 2 订阅

订阅专栏

斯坦纳树简介

斯坦纳树问题是组合优化问题，与最小生成树相似，是最短网络的一种。最小生成树是在给定的点集和边中寻求最短网络使所有点连通。而最小斯坦纳树允许在给定点外增加额外的点，使生成的最短网络开销最小。-------- 百度百科

斯坦纳树问题的最初形式是这样的：给定平面上的一系列点 $P_1，P_2，P_3，...，P_n$ ，现在要想办法用边将这 $n$ 个点连通。允许在平面上添加一些新的点（这些点称为“斯坦纳点”）作为辅助，添加的边的端点可以是原本存在的点，也可以是斯坦纳点，要求这个网络的总长度最小。

例如，当 $n = 3$ 时，记三个点为 $A,\ B,\ C$ ，有如下结论：如果三角形 $A B C$ 的每个内角都小于 $120\degree$ ，那么增加一个斯坦纳点 $S$ ，当 $SA\ ,SB,\ SC$ 的夹角都为 $120\degree$ 时， $S A + S B + S C$ 最小。如果三角形 $A B C$ 有一个角大于等于 $120\degree$ ，那么斯坦纳点与大于 $120\degree$ 度的角的顶点重合。
在这里插入图片描述

斯坦纳树问题的定义随着历史的发展不断地推广，后来出现了带点权的斯坦纳树问题。但算法竞赛中似乎没有见到过这种几何意义上的斯坦纳树问题，而是斯坦纳树的思想在图论中的一种体现。

图论中的斯坦纳树问题

图论中的斯坦纳树问题失去了原本的几何意义，实在离散的图上定义的。

给定 $N$ 个结点和 $M$ 条边的连通图，选取 $K$ 个关键点，需要保留一些边使得 $K$ 个关键点连通。要求边权和最小。这 $K$ 个点可能不会直接连通，所以需要利用剩下 $N - K$ 个非关键点，这些点就是斯坦纳点。

求解斯坦纳树的方法是，状态压缩动态规划。

状态的定义是 $f (i, S)$ ， $i$ 表示以 $i$ 结点为根。 $S$ 是一个整数表示的状态，表示的是 $K$ 个关键点的连通信息。也就是 $f (i, S)$ 表示以 $i$ 为根，包含集合 $S$ 中所有点的最小边权值和。

状态转移需要经历两个过程，分别是子集转移和当前状态下的边松弛操作。

$f(i, S)=min\{f(i,S)\ ,\ f(i,T)+f(i,\ S-T)\}$

其中 $T$ 是 $S$ 的一个子集， $S - T$ 表示 $T$ 的补集。子集转移需要枚举状态 $S$ 的所有子集，枚举子集可以用如下方法：

for(int sub=S&(S-1);sub;sub=S&(sub-1))
{
	//sub 就是从大到小枚举的 S的子集 
	//S^sub 就是sub的补集
}

可以证明这样做的正确性：首先显然 $S\&(sub-1)$ 一定是 $S$ 的子集，所以只需证明 $s u b$ 和 $\& S$ 之间没有 $S$ 的子集。
减一相当于将 $s u b$ 最右边的 1 变成 0，同时将这个 1 左边的 0 都变成 1。所以 $sub\&(sub-1)$ 一定是小于 $s u b$ 的最大的 $S$ 的子集。

$f(i, S)=min\{f(i,S)\ \ ,\ \ f(j,S)+w(i,j)\}$

可以发现这部分的转移是在同一个状态 $S$ 下进行的，参与转移的只有边权，实际上就转化为了最短路问题，可以直接用 $D i j k s t r a$ 或 $S P F A$ 来转移。

算法核心代码如下：

memset(dp, INF, sizeof(dp));
for(int i=1;i<=K;i++)
{
	scanf("%d", &key[i]);       //第i个关键点
	dp[key[i]][1<<(i-1)] = 0;
}
for(int S=1;S<(1<<K);S++)       //枚举K个关键点的所有状态 
{
	for(int i=1;i<=N;i++)       //枚举根节点 
	{
		for(int sub=S&(S-1);sub;sub=S&(sub-1))     //枚举 S 的所有子集
		    dp[i][S] = min(dp[i][S], dp[i][sub]+dp[i][S^sub]);
		if(dp[i][S]!=INF)
		    q.push((t_Node){dp[i][S], i});         //最短路的优先队列 
	}
	Dijkstra(S);                //对状态 S 进行状态内转移 
}

实际问题中 $K$ 不会太大，我猜这个问题可能是 $N P$ 的，没有关于 $K$ 的多项式时间算法。

斯坦纳森林

有时候问题并不要求 $K$ 个关键点全部连通，而是满足某些要求的点连通即可，也就是说连通后的图为多棵斯坦纳树组成的斯坦纳森林。

求解斯坦纳森林首先要求出正常的斯坦纳树的结果，然后对斯坦纳树再进行一次状态压缩。

令 $f (S)$ 表示关键点连通状态为 $S$ 的最小边权和。枚举 $S$ 的子集进行转移，但前提是 $S$ 需要满足题目所给的要求。

memset(f, INF, sizeof(f));
for(int S=1;S<1<<K;S++)
{
	if(check(S))         //检查状态S是否满足要求 
	{
		for(int i=1;i<=N;i++)
		    f[S] = min(f[S], dp[i][S]);
	}
}
for(int S=1;S<S<1<<K);S++)
{
	for(int sub=S&(S-1);sub;sub=S&(sub-1))
	{
		if(check(sub)&&check(S^sub))
		    f[S] = min(f[S], f[sub]+f[S^sub]);
	}
}
printf("%d\n", f[(1<<K)-1]);