Bzoj1565:[NOI2009]植物大战僵尸:拓扑排序+网络流

13 篇文章 0 订阅
1 篇文章 0 订阅

题目链接:1565:[NOI2009]植物大战僵尸

显然的最大权闭合图模型,但是裸上是错的

每个植物相他所攻击的位置连边,如果形成一个环那么这个换必然无解

进一步,这个环所保护的点都不能取到

所以我们拓扑排序求出可以干掉的植物集合再做网络流即可

#include<queue>
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn=1000010;
const int inf=0x7fffffff/3;
int n,m,tot=1,h[maxn],In[maxn],S,T;
struct edge{int to,next,w;}G[maxn];
int cur[maxn],vis[maxn],val[maxn],ans=0;
bool flag[maxn];

void add(int x,int y,int z){
	In[x]++;
	G[++tot].to=y;G[tot].next=h[x];h[x]=tot;G[tot].w=z;
	G[++tot].to=x;G[tot].next=h[y];h[y]=tot;G[tot].w=0;
}

void toposort(){
	queue<int>q;  
	for (int i=S;i<=T;++i) if (!In[i]) q.push(i);
	while (!q.empty()){
		int u=q.front(); q.pop(); flag[u]=1;
		if (val[u]>0) ans+=val[u];
		for (int i=h[u];i;i=G[i].next){
			int v=G[i].to;
			if (!G[i].w){
				In[v]--;
				if (!In[v]) q.push(v);
			}
			
		}
	}
}

bool bfs(){
    for (int i=S;i<=T;++i) vis[i]=-1;
    queue<int>q; q.push(S); vis[S]=0;
    while (!q.empty()){
        int u=q.front(); q.pop();
        for (int i=h[u];i;i=G[i].next){
            int v=G[i].to;
            if (vis[v]==-1&&G[i].w>0&&flag[v]){
                vis[v]=vis[u]+1;
                q.push(v);
            }
        }
    }return vis[T]!=-1;
}
 
int dfs(int x,int f){
    if (x==T||!f) return f;
    int used=0,w;
    for (int i=cur[x];i;i=G[i].next)
        if (vis[G[i].to]==vis[x]+1){
            w=f-used;
            w=dfs(G[i].to,min(G[i].w,w));
            G[i].w-=w; G[i^1].w+=w;
            used+=w; if (G[i].w) cur[x]=i;
            if (used==f) return used;
        }
    if (!used) vis[x]=-1;
    return used;
}
 
int dinic(){
    int ret=0;
    while (bfs()){
        for (int i=S;i<=T;++i) cur[i]=h[i];
        ret+=dfs(S,inf);
    }return ret;
}

int main(){
	scanf("%d%d",&n,&m);
	S=0; T=n*m+1;
	for (int i=1;i<=n*m;++i){
		scanf("%d",&val[i]);
		if (val[i]>0) add(S,i,val[i]);
		else add(i,T,-val[i]); 
		int a,b,c;
		scanf("%d",&a);
		for (int j=1;j<=a;++j){
			scanf("%d%d",&b,&c);
			add(b*m+c+1,i,inf);
		}
		if (i%m!=0) add(i,i+1,inf);//选了左侧必须先选右侧 
	}
	toposort();
	printf("%d",ans-dinic());
} 



  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值