bzoj1565:[NOI2009]植物大战僵尸

传送门

对于一些植物存在保护它的植物，比如能攻击到它的和它右边的植物都能保护它
这就符合最大权闭合子图的定义
按照最大权闭合子图的建模方式，正权连源点，负权连汇点，中间就对于每个植物连向保护它的植物就好了
可惜有些保护是成环的，环内的植物以及被这个环保护的植物你一个也碰不了
所以先反向建图，跑一边拓扑排序，将环处理出来，建模的时候忽略掉这些点就好了
代码：

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<queue>
#include<cstring>
#include<vector>
using namespace std;
void read(int &x) {
    char ch; bool ok;
    for(ok=0,ch=getchar(); !isdigit(ch); ch=getchar()) if(ch=='-') ok=1;
    for(x=0; isdigit(ch); x=x*10+ch-'0',ch=getchar()); if(ok) x=-x;
}
void print(int x){
    if(x<0)putchar('-'),x=-x;
    if(x>9)print(x/10);
    putchar(x%10+'0');
}
#define rg register
const int maxn=1250,inf=1e9;queue<int>q;bool vis[maxn];vector<int>g[maxn];
int dis[maxn],in[maxn],cnt=1,n,m,ans,pre[maxn*maxn],sum,x[30][40],r[30][40][700],c[30][40][700];
int cur[maxn],s,t,a[21][31],nxt[maxn*maxn],h[maxn],v[maxn*maxn];
int id(int x,int y){return (x-1)*m+y;}
void add(int x,int y,int z)
{
    pre[++cnt]=y,nxt[cnt]=h[x],h[x]=cnt,v[cnt]=z;
    pre[++cnt]=x,nxt[cnt]=h[y],h[y]=cnt,v[cnt]=0;
}
bool bfs()
{
    memset(dis,0,sizeof dis);
    q.push(s),dis[s]=1;
    while(!q.empty())
    {
        int x=q.front();q.pop();
        for(rg int i=h[x];i;i=nxt[i])
            if(v[i]&&!dis[pre[i]]&&!in[pre[i]])dis[pre[i]]=dis[x]+1,q.push(pre[i]);
    }
    return dis[t];
}
int dfs(int x,int flow)
{
    if(x==t||!flow)return flow;
    int f=flow;
    for(rg int &i=cur[x];i;i=nxt[i])
        if(dis[pre[i]]>dis[x]&&v[i])
        {
            int y=dfs(pre[i],min(f,v[i]));
            f-=y,v[i]-=y,v[i^1]+=y;
            if(!f)return flow;
        }
    if(f==flow)dis[x]=-1;
    return flow-f;
}
int main()
{
    read(n),read(m);t=n*m+1;
    for(rg int i=1;i<=n;i++)
        for(rg int j=1;j<=m;j++)
        {
            read(a[i][j]),read(x[i][j]);
            for(rg int k=1;k<=x[i][j];k++)read(r[i][j][k]),read(c[i][j][k]),r[i][j][k]++,c[i][j][k]++,in[id(r[i][j][k],c[i][j][k])]++,g[id(i,j)].push_back(id(r[i][j][k],c[i][j][k]));
            if(j!=1)in[id(i,j-1)]++,g[id(i,j)].push_back(id(i,j-1));
        }
    queue<int>qq;
    for(rg int i=1;i<=n*m;i++)if(!in[i])qq.push(i),vis[i]=1;
    while(!qq.empty())
    {
        int x=qq.front(),w=g[x].size();qq.pop();
        for(rg int i=0;i<w;i++)if(!(--in[g[x][i]]))qq.push(g[x][i]),vis[g[x][i]]=1;
    }
    for(rg int i=1;i<=n;i++)
        for(rg int j=1;j<=m;j++)
            if(vis[id(i,j)])
            {
                if(a[i][j]>0)add(s,id(i,j),a[i][j]),sum+=a[i][j];if(j!=m)add(id(i,j),id(i,j+1),inf);
                for(rg int k=1;k<=x[i][j];k++)add(id(r[i][j][k],c[i][j][k]),id(i,j),inf);
                if(a[i][j]<0)add(id(i,j),t,-a[i][j]);
            }
    for(;bfs();ans+=dfs(s,inf))memcpy(cur,h,sizeof h);
    printf("%d\n",sum-ans);
}

转载于:https://www.cnblogs.com/lcxer/p/10533487.html