NOI2009 植物大战僵尸

啊一道好题
感觉写得挺爽的啊这题
这种有一点懵逼然后学了一点东西之后很明朗的感觉真是好!
预处理参考 :http://www.cppblog.com/MatoNo1/archive/2014/11/01/142798.html?opt=admin
建图参考 :http://hihocoder.com/problemset/problem/1398?sid=1010444

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<algorithm>
#include<queue>
#include<vector>
#define N 80000
#define mod
#define inf 300000000
#define rg register
using namespace std;
typedef long long ll;
int n,m,tot=0,cnt=0,tmp=0,s,t;
int dis[N],hd[N],head[N],c[N],in[N],d[(N<<3)],w[N];
bool del[N];
struct hh{
	int to,nxt,w;
}b1[N*100],b[N*100];
int p (int x,int y) { return (x-1)*m+y; }
void link(int x,int y,int w){
	b1[tot].nxt=hd[x];
	b1[tot].to=y;
	b1[tot].w=w;
	hd[x]=tot++;
}
void lnk(int x,int y,int w){
	b[tot].nxt=head[x];
	b[tot].to=y;
	b[tot].w=w;
	head[x]=tot++;
}
void dele(int x){
	del[x]=1;
	for (int i=hd[x];i!=-1;i=b1[i].nxt)
		if (!del[b1[i].to]) dele(b1[i].to);
}
void top(){
	int h=0,t=0;
	for (int i=1;i<=cnt;++i){
		if (!in[i]) d[++t]=i;
		else del[i]=1;
	}
	while (h<=t){
		int x=d[++h];
		for (int i=hd[x];i!=-1;i=b1[i].nxt){
			in[b1[i].to]--,del[b1[i].to]=0;
			if (!in[b1[i].to]) d[++t]=b1[i].to;
		}
	}
	for (int i=1;i<=cnt;++i) if (del[i]) dele(i);
}
bool bfs(){
	queue <int> q;
	memset(dis,-1,sizeof(dis));
	q.push(s);dis[s]=0;
	while (!q.empty()){
		int h=q.front();
		q.pop();
		for (int i=head[h];i!=-1;i=b[i].nxt){
			int v=b[i].to;
			if (dis[v]!=-1||b[i].w<=0) continue;
			dis[v]=dis[h]+1;
			if (v==t) return 1;
			q.push(v);
		}
	}
	return dis[t]==-1?0:1;
}
int dfs(int u,int f){
	if (u==t||!f) return f;
	int res=f;
	for (int &i=c[u];i!=-1;i=b[i].nxt){
		int v=b[i].to;
		if (dis[v]!=dis[u]+1||b[i].w<=0) continue;
		int d=dfs(v,min(res,b[i].w));
		b[i].w-=d,b[i^1].w+=d;
		res-=d;
		if (res<=0) break;
	}
	return f-res;
}
int dinic(){
	int ans=0;
	while (bfs()){
		for (int i=s;i<=t;++i) c[i]=head[i];
		ans+=dfs(s,inf);
	}
	return ans;
}
int main(){
	freopen ("pvz.in","r",stdin);
	freopen ("pvz.out","w",stdout);
	memset(hd,-1,sizeof(hd));
	memset(head,-1,sizeof(head));
	scanf ("%d%d",&n,&m);
	int x,y,z;
	for (int i=1;i<=n;++i)
		for (int j=1;j<=m;++j){
			scanf ("%d%d",&w[p(i,j)],&z);
			for (int k=1;k<=z;++k) scanf ("%d%d",&x,&y),link(p(i,j),p(x+1,y+1),0),in[p(x+1,y+1)]++;
		}
	for (int i=1;i<=n;++i)
		for (int j=2;j<=m;++j) link(p(i,j),p(i,j-1),0),in[p(i,j-1)]++;
	cnt=n*m;
	top();
	s=0,t=cnt+1,tmp=0;
	tot=0;
	for (int i=1;i<=cnt;++i){
		if (del[i]) continue;
		if (w[i]>0) tmp+=w[i],lnk(i,t,w[i]),lnk(t,i,0);
		else lnk(s,i,-w[i]),lnk(i,s,0);
		for (int j=hd[i];j!=-1;j=b1[j].nxt)
			if (!del[b1[j].to]) lnk(i,b1[j].to,inf),lnk(b1[j].to,i,0);
	}
	tmp-=dinic();
	printf ("%d\n",max(0,tmp));
	return 0;
}

  

时间: 2024-10-14 06:27:25

NOI2009 植物大战僵尸的相关文章

NOI2009植物大战僵尸

这题应该分两步来做: 1.拓扑排序,去掉无敌点 2.求最大闭合子图 需要注意几点: 1.拓扑排序时,如果(i,j)可以攻击到(x,y),那么增加(x,y)的入度,而不是(i,j)的入度 因为入度代表着要攻击它需要事先攻击几个点 2.求最大闭合子图时,用所有的正权点-最大流 3.求最大闭合子图时,如果(i,j)可以攻击到(x,y),那么连一条边(x,y)到(i,j),容量为正无穷 因为在最大闭合子图中边(x,y)到(i,j)意味着选(x,y)就必须要选(i,j),这与实际含义相符 4.s到正权点,

COGS410. [NOI2009] 植物大战僵尸

410. [NOI2009] 植物大战僵尸 ★★★   输入文件:pvz.in   输出文件:pvz.out   简单对比时间限制:2 s   内存限制:512 MB [问题描述] Plants vs. Zombies(PVZ)是最近十分风靡的一款小游戏.Plants(植物)和Zombies(僵尸)是游戏的主角,其中Plants防守,而Zombies进攻.该款游戏包含多种不同的挑战系列,比如Protect Your Brain.Bowling等等.其中最为经典的,莫过于玩家通过控制Plants来

图论(网络流):COGS 410. [NOI2009] 植物大战僵尸

410. [NOI2009] 植物大战僵尸 ★★★   输入文件:pvz.in   输出文件:pvz.out   简单对比时间限制:2 s   内存限制:512 MB [问题描述] Plants vs. Zombies(PVZ)是最近十分风靡的一款小游戏.Plants(植物)和Zombies(僵尸)是游戏的主角,其中Plants防守,而Zombies进攻.该款游戏包含多种不同的挑战系列,比如Protect Your Brain.Bowling等等.其中最为经典的,莫过于玩家通过控制Plants来

Bzoj1565 [NOI2009]植物大战僵尸

Time Limit: 10 Sec  Memory Limit: 64 MBSubmit: 2363  Solved: 1092 Description Input Output 仅包含一个整数,表示可以获得的最大能源收入.注意,你也可以选择不进行任何攻击,这样能源收入为0. Sample Input 3 2 10 0 20 0 -10 0 -5 1 0 0 100 1 2 1 100 0 Sample Output 25 HINT 在样例中, 植物P1,1可以攻击位置(0,0), P2, 0

bzoj 1565 [NOI2009]植物大战僵尸

拓扑排序+最小割. 每个植物有它保护的一些植物,等价于选一些点的先决条件是其它一些点,最大权闭合子图问题. 发现图里有环,所以环里的点都不能选,先决条件是环里的点的点也不能选,所以把所有边反向拓扑排序,拓扑不到的点不能选. 无视掉不能选的点后裸最大权闭合子图建图就好了. 一个最大权闭合子图裸题的题解 1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #include<algorithm>

luogu2805 [NOI2009]植物大战僵尸

想象一下,要搞掉一个植物,必须先搞掉另一些植物--我们可以发现这是一个最大权闭合子图的问题. 最大权闭合子图的话,太空飞行计划问题是一个入门题,可以一看. 然而我们手玩一下样例就会惊恐地发现,保护关系出现了环. 比如说: 从 i 到 j 的箭头是说要搞 i 就要先搞 j,这是为了和最大权闭合子图的模型对应. 要搞 3 就要先搞 4,要搞 4 就要先搞 5,要搞 5 就要先搞 3--简直是乱来-- 这就像是一个无冷却食人花前面放一个坚果 所以就要删环了,把边倒过来拓扑排序一下,能被拓扑到的点就是合

[Luogu 2805]NOI2009 植物大战僵尸

<题目链接> 这题是个比较经典的最大权闭合子图,可以建图转化为最小割问题,再根据最大流最小割定理,采用任意一种最大流算法求得. 对于每个点,如果点权w为正,则从源点到这个点连一条边权为w的有向边:否则如果w为负则从这个点到汇点连一条有向边.加边的时候预处理出反图的每个点入度. 其次,每一个被保护的点到保护它的点连一条边权为INF的有向边. 注意同一行的左侧点受到右侧点的间接保护,因此对于每一个不位于当前行最右的点,都要向其右侧的一个点连一条INF有向边(连一条即可,不必连所有). 初始化完成后

98植物大战僵尸OL_僵尸迷阵

最近在玩植物大战僵尸Ol打到僵尸迷阵,打了个700分倒数第一,擦.....俗话说的好,失败是成功的妈妈,于是花了点时间写了一个小玩意,我相信下次我一定能拿第一 代码非常简单: using UnityEngine; using System.Collections; using UnityEngine.UI; public class MyScript : MonoBehaviour { // Use this for initialization void Start () { } // Upd

植物大战僵尸2天空之城安卓版发布

全体注意!一大波飞行僵尸正在接近中——中国独创版<植物大战僵尸2天空之城>安卓版终于发布了.超乎想象的空战体验.全新僵尸植物军团.独特闪电环境效果和史无前例的战舰成长系统,让你和戴夫共同肩负守卫天空之城的荣耀重任.赶紧下载游戏,加入云端激战吧. <植物大战僵尸2天空之城>安桌版上线 [中国区独有版本 打造铿锵空战体验] <植物大战僵尸2>即将迎来两周年生日,为感谢中国玩家一直以来的热情支持,EA/PopCap携手拓维游戏打造了这款中国区独有新版本<植物大战僵尸2天