【NOI2009】bzoj1565 植物大战僵尸

Description Input Output

仅包含一个整数，表示可以获得的最大能源收入。注意，你也可以选择不进行任何攻击，这样能源收入为0。

最大权闭合子图。注意如果有环，那么环和环的后继是都不能选的，可以拓扑排序找出来。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int s=1005,t=1006,oo=0x3f3f3f3f;
int f1[1010],n1[2000010],t1[2000010],w1[2000010],que[1010],f[1010],
ok[1010],du[1010],out[1010],f2[1010],n2[2000010],t2[2000010],a[1010],
n,m,tot1,tot2,sum;
void a1(int u,int v,int x)
{
tot1++;
n1[tot1*2]=f1[u];
f1[u]=tot1*2;
t1[tot1*2]=v;
w1[tot1*2]=x;
n1[tot1*2+1]=f1[v];
f1[v]=tot1*2+1;
t1[tot1*2+1]=u;
w1[tot1*2+1]=0;
}
void a2(int u,int v)
{
du[v]++;
n2[++tot2]=f2[u];
f2[u]=tot2;
t2[tot2]=v;
}
void init()
{
int i,j,k,x,y;
scanf("%d%d",&n,&m);
for (i=1;i<=n;i++)
for (j=1;j<=m;j++)
{
scanf("%d%d",&a[i*m+j],&k);
while (k--)
{
scanf("%d%d",&x,&y);
x++;y++;
a1(x*m+y,i*m+j,oo);
a2(i*m+j,x*m+y);
}
}
for (i=1;i<=n;i++)
for (j=1;j<m;j++)
{
a1(i*m+j,i*m+j+1,oo);
a2(i*m+j+1,i*m+j);
}
}
void top()
{
int i,j,u,v,hd=1,tl=0;
/*for (int i=1;i<=t;i++)
for (int j=f2[i];j;j=n2[j])
printf("%d->%d\n",i,t2[j]);*/
for (i=1;i<=n;i++)
for (j=1;j<=m;j++)
if (du[i*m+j]==0)
que[++tl]=i*m+j;
while (hd<=tl)
{
u=que[hd++];
ok[u]=1;
for (i=f2[u];i;i=n2[i])
if (--du[v=t2[i]]==0)
que[++tl]=v;
}
}
void build()
{
int i,j;
for (i=1;i<=n;i++)
for (j=1;j<=m;j++)
if (ok[i*m+j])
{
if (a[i*m+j]>=0)
{
sum+=a[i*m+j];
a1(s,i*m+j,a[i*m+j]);
}
else a1(i*m+j,t,-a[i*m+j]);
}
}
bool bfs()
{
int i,j,u,v,hd,tl;
for (i=1;i<=n;i++)
for (j=1;j<=m;j++)
f[i*m+j]=0;
f[t]=0;
f[que[1]=s]=hd=tl=1;
while (hd<=tl)
{
u=que[hd++];
for (i=f1[u];i;i=n1[i])
if (w1[i]&&!f[v=t1[i]])
{
que[++tl]=v;
f[v]=f[u]+1;
}
}
return f[t];
}
int dfs(int u,int lim)
{
if (u==t) return lim;
int i,v,x,ret=0;
for (i=f1[u];i&&ret<lim;i=n1[i])
if (w1[i]&&f[v=t1[i]]==f[u]+1)
{
x=dfs(v,min(lim-ret,w1[i]));
w1[i]-=x;
w1[i^1]+=x;
ret+=x;
}
if (!ret) f[u]=0;
return ret;
}
int solve()
{
int ans=0,x;
while (bfs())
while (x=dfs(s,oo))
ans+=x;
return sum-ans;
}
int main()
{
init();
top();
build();
/*for (int i=1;i<=t;i++)
for (int j=f1[i];j;j=n1[j])
if (w1[j])
printf("%d->%d:%d\n",i,t1[j],w1[j]);*/
printf("%d\n",solve());
}