bzoj-1565 植物大战僵尸

题意：

给出一片n*m的草坪，上面每个点有一个植物；

现在由你来从最右面放出一些僵尸来进攻这些植物；

僵尸到一个植物面前的时候就可以吃掉这个植物，并且得到这个植物的的得分(可正可负)；

每个植物可以攻击几个位置，并且僵尸只能从右面出发，也就是说在吃掉右面的植物之前不可能吃掉左面的；

僵尸是无限的，求最大得分；

n<=20，m<=30；

题解：

原来这个就是最大权闭合子图啊。。似乎以前没刷过的样子= =

总之先建图：

将每个植物向能保护它的植物连流量为无穷的边；

从源点向正得分植物连流量为得分的边；

从负得分植物向汇点连流量为得分相反数的边；

最后用所有正权值的和减去最小割就是了；

这样建图解决最大权闭合子图问题是正确的，感性上讨论一下各种边就差不多理解的了吧；

然而通常的最大权闭合子图都是拓扑图，在这个问题中可以有环存在；

所以有些点是永远也不可能取到的，这些点我们不要，也就是拓扑排序之后排不到的点；

然后建图就是正确的了，跑最大权闭合子图即可；

代码：

#include<queue>
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define N 40
#define E 400000
using namespace std;
int next[E],to[E],flow[E],head[N*N],ce=1;
int dis[N*N],S,T,sum;
int val[N*N];
queue<int>q;
void add(int x,int y,int fl)
{
to[++ce]=y;
flow[ce]=fl;
next[ce]=head[x];
head[x]=ce;
to[++ce]=x;
flow[ce]=0;
next[ce]=head[y];
head[y]=ce;
}
bool BFS()
{
int x,i;
memset(dis,0,sizeof(dis));
dis[S]=1;
q.push(S);
while(!q.empty())
{
x=q.front(),q.pop();
for(i=head[x];i;i=next[i])
{
if(flow[i]&&!dis[to[i]])
{
dis[to[i]]=dis[x]+1;
q.push(to[i]);
}
}
}
return dis[T]!=0;
}
int dfs(int x,int lim)
{
if(x==T) return lim;
int ret=0,i,temp;
for(i=head[x];i;i=next[i])
{
if(flow[i]&&dis[to[i]]==dis[x]+1)
{
temp=dfs(to[i],min(flow[i],lim-ret));
ret+=temp;
flow[i]-=temp,flow[i^1]+=temp;
if(ret==lim)
return ret;
}
}
if(!ret) dis[x]=0;
return ret;
}
namespace Graph
{
int next[E],to[E],head[N*N],ce;
int in[N*N];
bool vis[N*N];
void add(int x,int y)
{
to[++ce]=y;
next[ce]=head[x];
head[x]=ce;
in[y]++;
}
void slove(int n,int m)
{
int i,x,tot=n*m;
for(x=0;x<tot;x++)
{
if(!in[x])
q.push(x);
}
while(!q.empty())
{
x=q.front(),q.pop();
vis[x]=1;
for(i=head[x];i;i=next[i])
{
in[to[i]]--;
if(!in[to[i]])
q.push(to[i]);
}
}
for(x=0;x<tot;x++)
{
if(!vis[x]) continue;
if(val[x]>=0)
::add(S,x,val[x]),sum+=val[x];
else
::add(x,T,-val[x]);
for(i=head[x];i;i=next[i])
{
if(vis[to[i]])
::add(to[i],x,0x3f3f3f3f);
}
}
}
}
int main()
{
int n,m,w,i,j,k,x,y,ans;
scanf("%d%d",&n,&m);
S=n*m,T=S+1;
for(i=0,sum=0;i<n;i++)
{
for(j=0;j<m;j++)
{
if(j) Graph::add(i*m+j,i*m+j-1);
scanf("%d%d",&val[i*m+j],&w);
for(k=1;k<=w;k++)
{
scanf("%d%d",&x,&y);
Graph::add(i*m+j,x*m+y);
}
}
}
Graph::slove(n,m);
ans=0;
while(BFS())
ans+=dfs(S,0x3f3f3f3f);
printf("%d\n",sum-ans);
return 0;
}