[NOI2009]植物大战僵尸

这道题跟NOI2006 最大获利其实是很像的

一样都是要搞定一些点才能搞定另一些点,然后有些点正权有些点负权

这种问题,其实是最大权闭合子图

amber的最小割论文有详细的讲解法和证明


闭合子图的定义是,图中每个点所连接的的任何一条边不指向图外,可以有边指向这个图

这实际上就是一个依赖关系,如果我们把a依赖b(在这道题就是b保护a),在图中用一条a指向b的边表示

那么我们求的就是一个闭合子图,就是说不可能出现图中选了a而没有选b的情况,因为这样不满足闭合子图的定义

而最大权闭合子图则可以用网络流来解决

对于原图,我们把每条边的流量设为∞

然后对于每个正权点,从st向点连一条边,边权为点的权值

对于每个负权点,从点向ed连一条边,边权为点的权值的绝对值

那么最大权闭合子图的答案就是所有正权-网络流的最小割

因为是最小割,显然我们不会割到原图的流量为∞的边,这样就满足了闭合子图的定义——并没有破坏依赖关系

然后我们如果割了连接st的边,相当于舍弃掉这个正权点,割了连接ed的边,相当于选择了这个负权点

最后所有正权点的权值和(可能达到的最大值)-最小割(付出的代价),就是最大权闭合子图的答案


对于这道题和NOI2006最大获利,一样的建图即可


但是这道题有一点比最大获利要麻烦,就是这道题可以出现环

如果出现了环,求最大权闭合子图的时候则视情况将这个环里的点一并取走或不取

但这道题则要求无论如何都不能取

因此要先toposort一遍去环再建图才行……


这道题我用了递归版的ISAP(主要是可以少写一个BFS),结果暴慢……虽然也能AC……

GZH神犇写了个非递归的ISAP,比我的Dinic快……

好吧以后再也不写递归ISAP了……老老实实写Dinic……


//Lib
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<ctime>

#include<iostream>
#include<algorithm>
#include<vector>
#include<string>
#include<queue>
#include<stack>
using namespace std;
//Macro
#define	rep(i,a,b)	for(int i=a,tt=b;i<=tt;++i)
#define	drep(i,a,b)	for(int i=a,tt=b;i>=tt;--i)
#define	erep(i,e,x)	for(int i=x;i;i=e[i].next)
#define	irep(i,x)	for(__typedef(x.begin()) i=x.begin();i!=x.end();i++)
#define	read()	(strtol(ipos,&ipos,10))
#define	sqr(x)	((x)*(x))
#define	pb	push_back
#define	PS	system("pause");
typedef	long long	ll;
typedef	pair<int,int>	pii;
const int oo=~0U>>1;
const double inf=1e100;
const double eps=1e-6;
string name="pvz", in=".in", out=".out";
//Var
struct E
{
	int next,node,cap;
}e[1000008],e2[1000008];
queue<int> q;
int n,m,st,ed,tot=1;
int cnt[1008],dis[1008],h[1008],h2[1008],node[28][38],map[28][38],ind[1008],ans,size;
bool vis[1008];
void add(int a,int b,int c)
{
	e[++tot].next=h[a];e[tot].node=b;e[tot].cap=c;h[a]=tot;
	e[++tot].next=h[b];e[tot].node=a;e[tot].cap=0;h[b]=tot;
}
void add2(int a,int b){e2[++tot].next=h2[a];e2[tot].node=b;h2[a]=tot;ind[b]++;}
void Toposort()
{
	stack<int> s;
	int u,v;
	rep(i,1,n)
		rep(j,1,m)
		if(!ind[node[i][j]])s.push(node[i][j]);
	while(!s.empty())
	{
		u=s.top();s.pop();
		vis[u]=true;
		erep(i,e2,h2[u])
		{
			v=e2[i].node;
			if((--ind[v])==0)s.push(v);
		}
	}
}
void Build()
{
	tot=1;int v;
	rep(i,1,n)
		rep(j,1,m)
			if(vis[node[i][j]])
			{
				if(map[i][j]>0)		add(st,node[i][j],map[i][j]),ans+=map[i][j];
				else if(map[i][j]<0)add(node[i][j],ed,-map[i][j]);
				erep(k,e2,h2[node[i][j]])
					if(vis[v=e2[k].node])add(v,node[i][j],oo);
			}
}
void Init()
{
	scanf("%d%d",&n,&m);
	st=n*m+1;ed=n*m+2;int c,t,a,b;
	rep(i,1,n)rep(j,1,m)node[i][j]=++size;
	rep(i,1,n)
	{
		rep(j,1,m)
		{
			scanf("%d%d",&c,&t);
			map[i][j]=c;
			rep(k,1,t)
			{
				scanf("%d%d",&a,&b);a++;b++;
				add2(node[i][j],node[a][b]);
			}
			if(j!=m)add2(node[i][j+1],node[i][j]);
		}
	}
	Toposort();
	Build();
}
bool BFS()
{
	memset(dis,-1,sizeof dis);
	dis[ed]=0;q.push(ed);
	int u,v;bool flag=false;
	while(!q.empty())
	{
		u=q.front();q.pop();
		erep(i,e,h[u])
		{
			if(dis[v=e[i].node]==-1&&e[i^1].cap)
			{
				dis[v]=dis[u]+1;
				q.push(v);
				if(v==st)flag=true;	
			}
		}	
	}
	return flag;
}
int DFS(int u,int low)
{
	if(u==ed)return low;
	int ret=low,tmp,v;
	erep(i,e,h[u])
	{
		if(!e[i].cap||dis[v=e[i].node]!=dis[u]-1)continue;
		tmp=DFS(v,min(low,e[i].cap));
		e[i].cap-=tmp;
		e[i^1].cap+=tmp;
		low-=tmp;
		if(low==0)break;
	}
	if(ret==low)dis[u]=-1;
	return ret-low;
}
void Work()
{
	int flow;
	while(BFS())
		while(flow=DFS(st,oo))
			ans-=flow;
	cout<<ans<<endl;
}
int main()
{
//	freopen((name+in).c_str(),"r",stdin);
//	freopen((name+out).c_str(),"w",stdout);
	Init();
	Work();
	return 0;
}

阅读更多
换一批

没有更多推荐了,返回首页