博客园 首页 私信博主 显示目录 隐藏目录 管理 动画

BZOJ.2521.[SHOI2010]最小生成树(最小割ISAP/Dinic)

题目链接

一条边不变其它边减少可以看做一条边增加其它边不变。
假设要加的边lab为(A->B,v),那么肯定是要使除这条边外,A->B的每条路径上的最小权值都\(>v\),这样在连通A,B时(即Kruskal中Union())才一定会选择这条边。
要求路径上最小边的权值\(>v\),即要求在路径上有任意一边权值\(\leq v\)时不连通。于是求最小割(使它不连通),割掉一条边的代价即\(v[lab]-v[i]+1\)

无向图建双向边。

status里的怎么都那么快?复制了一份20的也60+

ISAP:

//892kb	64ms
//ISAP lev[]的上限是n不是des。。
#include <cstdio>
#include <cctype>
#include <algorithm>
#define gc() getchar()
const int N=505,M=3215;

int n,m,lab,s[805],t[805],v[805];
int src,des,Enum,cur[N],H[N],nxt[M],fr[M],to[M],cap[M],pre[N],lev[N],que[N],num[N];

inline int read()
{
	int now=0;register char c=gc();
	for(;!isdigit(c);c=gc());
	for(;isdigit(c);now=now*10+c-'0',c=gc());
	return now;
}
inline void AddEdge(int u,int v,int w)
{
	to[++Enum]=v, fr[Enum]=u, nxt[Enum]=H[u], H[u]=Enum, cap[Enum]=w;
	to[++Enum]=u, fr[Enum]=v, nxt[Enum]=H[v], H[v]=Enum, cap[Enum]=0;
}
bool BFS()
{
	for(int i=1; i<=n; ++i) lev[i]=n+1;
	lev[des]=0, que[0]=des; int h=0,t=1;
	while(h<t)
	{
		int x=que[h++];
		for(int i=H[x]; i; i=nxt[i])
			if(lev[to[i]]==n+1 && cap[i^1])
				lev[to[i]]=lev[x]+1, que[t++]=to[i];
	}
	return lev[src]<=n;
}
int Augment()
{
	int mn=1e9;
	for(int i=des; i!=src; i=fr[pre[i]])
		mn=std::min(mn,cap[pre[i]]);
	for(int i=des; i!=src; i=fr[pre[i]])
		cap[pre[i]]-=mn, cap[pre[i]^1]+=mn;
	return mn;
}
int ISAP()
{
	if(!BFS()) return 0;
	for(int i=1; i<=n; ++i) ++num[lev[i]],cur[i]=H[i];
	int x=src,res=0;
	while(lev[src]<=n)
	{
		if(x==des) x=src,res+=Augment();
		bool can=0;
		for(int i=cur[x]; i; i=nxt[i])
			if(lev[to[i]]==lev[x]-1 && cap[i])
			{
				can=1, cur[x]=i, pre[x=to[i]]=i;
				break;
			}
		if(!can)
		{
			int mn=n;
			for(int i=H[x]; i; i=nxt[i])
				if(cap[i]) mn=std::min(mn,lev[to[i]]);
			if(!--num[lev[x]]) break;
			++num[lev[x]=mn+1], cur[x]=H[x];
			if(x!=src) x=fr[pre[x]];
		}
	}
	return res;
}

int main()
{
	n=read(),m=read(),lab=read(),Enum=1;
	for(int i=1; i<=m; ++i) s[i]=read(),t[i]=read(),v[i]=read();
	src=s[lab], des=t[lab];
	for(int i=1; i<=m; ++i)
		if(i!=lab&&v[i]<=v[lab]) AddEdge(s[i],t[i],v[lab]-v[i]+1),AddEdge(t[i],s[i],v[lab]-v[i]+1);
	printf("%d",ISAP());

	return 0;
}

Dinic:

//876kb	60ms 果然差不多
#include <cstdio>
#include <cctype>
#include <algorithm>
#define gc() getchar()
const int N=505,M=3215;

int n,m,lab,s[805],t[805],v[805];
int src,des,Enum,cur[N],H[N],nxt[M],to[M],cap[M],lev[N],que[N];

inline int read()
{
	int now=0;register char c=gc();
	for(;!isdigit(c);c=gc());
	for(;isdigit(c);now=now*10+c-'0',c=gc());
	return now;
}
inline void AddEdge(int u,int v,int w)
{
	to[++Enum]=v, nxt[Enum]=H[u], H[u]=Enum, cap[Enum]=w;
	to[++Enum]=u, nxt[Enum]=H[v], H[v]=Enum, cap[Enum]=0;
}
bool BFS()
{
	for(int i=1; i<=n; ++i) lev[i]=0,cur[i]=H[i];
	lev[src]=1, que[0]=src; int h=0,t=1;
	while(h<t)
	{
		int x=que[h++];
		for(int i=H[x]; i; i=nxt[i])
			if(!lev[to[i]] && cap[i])
			{
				lev[to[i]]=lev[x]+1, que[t++]=to[i];
				if(to[i]==des) return 1;
			}
	}
	return 0;
}
int Dinic(int x,int flow)
{
	if(x==des) return flow;
	int used(0);
	for(int &i=cur[x]; i; i=nxt[i])
		if(lev[to[i]]==lev[x]+1&&cap[i])
		{
			int delta=Dinic(to[i],std::min(flow-used,cap[i]));
			if(delta)
			{
				cap[i]-=delta, cap[i^1]+=delta;
				if((used+=delta)==flow) return flow;
			}
		}
	lev[x]=0;
	return used;
}

int main()
{
	n=read(),m=read(),lab=read(),Enum=1;
	for(int i=1; i<=m; ++i) s[i]=read(),t[i]=read(),v[i]=read();
	src=s[lab], des=t[lab];
	for(int i=1; i<=m; ++i)
		if(i!=lab&&v[i]<=v[lab]) AddEdge(s[i],t[i],v[lab]-v[i]+1),AddEdge(t[i],s[i],v[lab]-v[i]+1);
	int res=0;
	while(BFS()) res+=Dinic(src,1e9);
	printf("%d",res);

	return 0;
}
posted @ 2018-03-30 09:57  SovietPower  阅读(175)  评论(0编辑  收藏  举报