bzoj4006(斯坦纳树+状压DP)

最新推荐文章于 2020-08-02 15:29:00 发布

qkoqhh

最新推荐文章于 2020-08-02 15:29:00 发布

阅读量328

点赞数

分类专栏：斯坦纳树 DP

本文链接：https://blog.csdn.net/qkoqhh/article/details/86671290

版权

DP 同时被 2 个专栏收录

116 篇文章 1 订阅

订阅专栏

斯坦纳树

1 篇文章 0 订阅

订阅专栏

参考链接：https://www.cnblogs.com/ECJTUACM-873284962/p/7643445.html#autoid-0-1-0

资料比较少，主要看上面学的（博客挺漂亮的

斯坦纳树的定义是给定m个关键点，使得这些关键点相互联通的最小权值的树（自己写的

MST是斯坦纳树的特殊情况，只要把关键点定义成所有点就可以了。。

求解斯坦纳树是NP问题，所以要关键点的个数要少(一般m<10)，然后才可以状压

令d[i][j]为覆盖了i点集，j为根的树的最小权值

那么状态转移有两重：

第一重：把两个子状态合并成一个状态

$d[i][j]=\max\limits_{t\subseteq i} (d[t][j]+d[i-t][j])$

复杂度 O(V3^m)

枚举子状态的时候可以参考挑战的位运算枚举方法

第二重：在同一状态中进行转移（就是SPFA

d[state][j]=max(d[state]][j],d[state][i]+e[i][j])

复杂度 O(E2^m)

然后就把斯坦纳树求解出来了。。

对这题，他要同一频率的点联通，那么再对频率做一个状压DP就可以了，每个状态的初始值从斯坦纳树里面抽取出来，然后做个合并就可以了。。

/**
 *　　　　　　  　　┏┓　　 　┏┓
 * 　　　　　  　　┏┛┗━━━━━━━┛┗━━━┓
 * 　　　　　  　　┃　　　　　　　┃ 　
 * 　　　　　  　　┃　　　━　　 　┃
 * 　　　　　  　　┃　＞　　　＜　┃
 * 　　　　　  　　┃　　　　　　　┃
 * 　　　　　  　　┃...　⌒　... 　┃
 * 　　　　　　  　┃              ┃
 * 　　　　　　  　┗━┓          ┏━┛
 * 　　　　　　　　　┃          ┃　Code is far away from bug with the animal protecting　　　　　　　　　　
 * 　　　　　　　　　┃          ┃   神兽保佑,代码无bug
 * 　　　　　　　　　┃          ┃　　　　　　　　　　　
 * 　　　　　　　　　┃          ┃  　　　　　　
 * 　　　　　　　　　┃          ┃
 * 　　　　　　　　　┃          ┃　　　　　　　　　　　
 * 　　　　　　　　　┃          ┗━━━┓
 * 　　　　　　　　　┃              ┣┓
 * 　　　　　　　　　┃              ┏┛
 * 　　　　　　　　　┗┓┓┏━━━━━━━━┳┓┏┛
 * 　　　　　　　　　　┃┫┫       ┃┫┫
 * 　　　　　　　　　　┗┻┛       ┗┻┛
 */ 
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
#include<map>
#include<stack>
#include<cmath>
#include<set>
#include<bitset>
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,l,r) for(int i=l;i>=r;i--)
#define link(x) for(edge *j=h[x];j;j=j->next)
#define mem(a) memset(a,0,sizeof(a))
#define ll long long
#define eps 1e-8
#define succ(x) (1<<x)
#define lowbit(x) (x&(-x))
#define sqr(x) ((x)*(x))
#define mid (x+y>>1)
#define NM 1025
#define nm 6005
#define pi 3.1415926535897931
const int inf=998244353;
using namespace std;
ll read(){
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
    return f*x;
}
 
 



struct edge{int t,v;edge*next;}e[nm],*h[NM],*o=e;
void add(int x,int y,int v){o->t=y;o->v=v;o->next=h[x];h[x]=o++;}
int n,m,p,tot,d[NM][NM],_x,_y,_t,g[NM],b[NM];
pair<int,int>c[NM];
queue<int>q;
bool v[NM];


void spfa(int k){
    while(!q.empty()){
	int t=q.front();q.pop();v[t]=false;
	link(t)if(d[k][j->t]>d[k][t]+j->v){
	    d[k][j->t]=d[k][t]+j->v;
	    if(!v[j->t])v[j->t]++,q.push(j->t);
	}
    }
}

int main(){
    freopen("data.in","r",stdin);
    n=read();m=read();p=read();tot=succ(p)-1;
    inc(i,1,m){_x=read();_y=read();_t=read();add(_x,_y,_t);add(_y,_x,_t);}
    inc(i,1,p)c[i].first=read(),c[i].second=read();
    sort(c+1,c+1+p);
    inc(i,1,p)if(c[i].first-c[i-1].first>1){
	int t=c[i].first,j=i;
	while(c[j].first==t)c[j++].first=c[i-1].first+1;
	i=j-1;
    }
    inc(i,1,tot)inc(j,1,n)d[i][j]=inf;
    inc(i,1,p)d[succ(i-1)][c[i].second]=0;
    inc(i,1,tot){
	mem(v);
	inc(j,1,n){
	    for(int k=i&(i-1);k;k=(k-1)&i)
		d[i][j]=min(d[i][j],d[k][j]+d[k^i][j]);
	    if(d[i][j]<inf)q.push(j),v[j]++;
	}
	spfa(i);
    }
    inc(i,1,p){
	int j=i;
	while(c[j].first==c[i+1].first)i++;
	b[c[j].first]=succ(i)-succ(j-1);
    }
    p=c[p].first;tot=succ(p)-1;
    inc(i,1,tot){
	int t=0;
	g[i]=inf;
	inc(j,1,p)if(succ(j-1)&i)t^=b[j];
	inc(j,1,n)g[i]=min(g[i],d[t][j]);
    }
    inc(i,1,tot)
	for(int j=i&(i-1);j;j=(j-1)&i)
	    g[i]=min(g[i],g[j]+g[j^i]);
    return 0*printf("%d\n",g[tot]);
}