JZOJ P4313 电话线铺设(最小生成树+倍增求LCA)

Description

Solution

王牌电缆n-2根,李牌电缆1根,使得形成的连通块边权之和最小,因此可以想到最小生成树;
首先对于所有的W边跑一遍Kruskal,得到一棵MST(最小生成树);
然后对于所有的L边:
若开始时得到的MST不完全连通,则依次枚举两个端点分别属于两个连通块的L边,取其最小值;
否则依次枚举所有L边,如果将其加入MST后会形成一个环,因此要删去这条环上最大的W边—对于L边的两个端点,用倍增求它们的最近公共祖先,就能得到两点间线段上最大的边及其序号.
最后输出权值最小值以及MST中处了已经删除的那一条边外其余的W边和一条L边.

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn=100005,maxm=200005,inf=0x3f3f3f3f;
int fa[maxn],head[maxn],dep[maxn],Ans[maxn];
int f[maxn][21],g[maxn][21],h[maxn][21];
int n,W,L,ans=inf,_ans=0,cnt=0,tot=0,maxpos;

struct E{
	int u,v,w,pos;
}e[maxm],o[maxm];

struct Edge{
	int nxt,to,edge,pos;
}ed[maxm];

int read()
{
	int res=0,fh=1;char ch=getchar();
	while(ch>'9'||ch<'0') {if(ch=='-') fh=-1; ch=getchar();}
	while(ch>='0'&&ch<='9') res=(res<<3)+(res<<1)+(ch^'0'),ch=getchar();
	return fh*res;
}
void add(int u,int v,int w,int pos){
	ed[++tot]=(Edge){ head[u], v, w, pos },head[u]=tot;
}
int max(int a,int b) { return a>b?a:b; }
int min(int a,int b) { return a<b?a:b; }
void swap(int &a,int &b) { int t=a;a=b,b=t; }
bool cmp(E a,E b) { return a.w<b.w; }
int find(int x) { return fa[x]!=x ? fa[x]=find(fa[x]) : fa[x] ; }

void dfs(int u,int father)
{
	dep[u]=dep[father]+1;
	for(int i=1;i<=20;++i) {
		f[u][i]=f[f[u][i-1]][i-1];
		if(g[u][i-1]>g[f[u][i-1]][i-1])
			g[u][i]=g[u][i-1],h[u][i]=h[u][i-1];
		else
			g[u][i]=g[f[u][i-1]][i-1],h[u][i]=h[f[u][i-1]][i-1];
	}
	for(int i=head[u];i;i=ed[i].nxt)
	{
		int v=ed[i].to;
		if(v==father) continue;
		f[v][0]=u,g[v][0]=ed[i].edge,h[v][0]=ed[i].pos;
		dfs(v,u);
	}
}

int lca(int x,int y)
{
	int maxedge=0;
	if(dep[x]<dep[y]) swap(x,y);
	for(int i=20;i>=0;--i)
		if(dep[f[x][i]]>=dep[y])
		{
			if(maxedge<g[x][i])
				maxedge=g[x][i],maxpos=h[x][i];
			x=f[x][i];
		}
	if(x==y) return maxedge;
	for(int i=20;i>=0;--i)
		if(f[x][i]!=f[y][i])
		{
			if(g[x][i]>maxedge)
				maxedge=g[x][i],maxpos=h[x][i];
			if(g[y][i]>maxedge)
				maxedge=g[y][i],maxpos=h[y][i];
			x=f[x][i],y=f[y][i];
		}
	if(maxedge<g[x][0]) maxedge=g[x][0],maxpos=h[x][0];
	if(maxedge<g[y][0]) maxedge=g[y][0],maxpos=h[y][0];
	return maxedge;
}

int main()
{
	n=read(),W=read(),L=read();
	for(int i=1;i<=W;++i) {
		int u=read(),v=read(),w=read();
		e[i]=(E){u,v,w,i};
	}
	for(int i=1;i<=L;++i) {
		int u=read(),v=read(),w=read();
		o[i]=(E){u,v,w,i};
	}
	
	sort(e+1,e+W+1,cmp);
	for(int i=1;i<=n;++i) fa[i]=i;
	for(int i=1;i<=W;++i)
	{
		int fx=find(e[i].u),fy=find(e[i].v);
		if(fx==fy) continue;
		fa[fx]=fy,_ans+=e[i].w,Ans[++cnt]=e[i].pos;
		add(e[i].u,e[i].v,e[i].w,e[i].pos);
		add(e[i].v,e[i].u,e[i].w,e[i].pos);
	}
	
	int at=0,del=0;
	if(cnt<n-1)
	{
		for(int i=1;i<=L;++i)
			if(find(o[i].u)!=find(o[i].v))
				if(ans>_ans+o[i].w)
					ans=_ans+o[i].w,at=i;
		printf("%d\n",ans);
		for(int i=1;i<=cnt;++i) printf("%d\n",Ans[i]);
		printf("%d",at);
		return 0;
	}
	dep[1]=1,dfs(1,0);
	for(int i=1;i<=L;++i)
	{
		int u=o[i].u,v=o[i].v,w=o[i].w;
		int maxedge=lca(u,v);
		if(ans>_ans-maxedge+w)
			ans=_ans-maxedge+w,at=i,del=maxpos;
	}
	
	printf("%d\n",ans);
	for(int i=1;i<=cnt;++i)
		if(Ans[i]!=del) printf("%d\n",Ans[i]);
	printf("%d\n",at);
	
	return 0;
}