啊一道好题
感觉写得挺爽的啊这题
这种有一点懵逼然后学了一点东西之后很明朗的感觉真是好!
预处理参考 :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;
}

  

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

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

    410. [NOI2009] 植物大战僵尸 ★★★   输入文件:pvz.in   输出文件:pvz.out   简单对比时间限制:2 s   内存限制:512 MB [问题描述] Plants vs ...

  2. P2805 [NOI2009]植物大战僵尸

    题目地址:P2805 [NOI2009]植物大战僵尸 最大权闭合子图 若有向图 \(G\) 的子图 \(V\) 满足: \(V\) 中顶点的所有出边均指向 \(V\) 内部的顶点,则称 \(V\) 是 ...

  3. COGS410. [NOI2009] 植物大战僵尸

    410. [NOI2009] 植物大战僵尸 ★★★   输入文件:pvz.in   输出文件:pvz.out   简单对比时间限制:2 s   内存限制:512 MB [问题描述] Plants vs ...

  4. BZOJ 1565: [NOI2009]植物大战僵尸

    1565: [NOI2009]植物大战僵尸 Time Limit: 10 Sec  Memory Limit: 64 MBSubmit: 2317  Solved: 1071[Submit][Stat ...

  5. 【刷题】BZOJ 1565 [NOI2009]植物大战僵尸

    Description Plants vs. Zombies(PVZ)是最近十分风靡的一款小游戏.Plants(植物)和Zombies(僵尸)是游戏的主角,其中Plants防守,而Zombies进攻. ...

  6. 【bzoj1565】[NOI2009]植物大战僵尸

    1565: [NOI2009]植物大战僵尸 Time Limit: 10 Sec  Memory Limit: 64 MBSubmit: 2164  Solved: 1001[Submit][Stat ...

  7. 【最大权闭合子图 tarjan】bzoj1565: [NOI2009]植物大战僵尸

    dinic+tarjan板子练手题 Description Plants vs. Zombies(PVZ)是最近十分风靡的一款小游戏.Plants(植物)和Zombies(僵尸)是游戏的主角,其 中P ...

  8. BZOJ1565: [NOI2009]植物大战僵尸

    Description Input Output 仅包含一个整数,表示可以获得的最大能源收入.注意,你也可以选择不进行任何攻击,这样能源收入为0. Sample Input 3 2 10 0 20 0 ...

  9. 【bzoj1565】 NOI2009—植物大战僵尸

    http://www.lydsy.com/JudgeOnline/problem.php?id=1565 (题目链接) 题意 给出$n*m$的棋盘,僵尸攻击每个格子可以获得$v$的分数,每个格子又会保 ...

  10. luogu2805 [NOI2009]植物大战僵尸

    想象一下,要搞掉一个植物,必须先搞掉另一些植物--我们可以发现这是一个最大权闭合子图的问题. 最大权闭合子图的话,太空飞行计划问题是一个入门题,可以一看. 然而我们手玩一下样例就会惊恐地发现,保护关系 ...

随机推荐

  1. 基于Flink秒级计算时CPU监控图表数据中断问题

     基于Flink进行秒级计算时,发现监控图表中CPU有数据中断现象,通过一段时间的跟踪定位,该问题目前已得到有效解决,以下是解决思路:   一.问题现象       以SQL02为例,发现本来10秒一 ...

  2. js 与 ios Android交互

    一.android 交互 1.js调用webview 在android API Level 17及以上的版本中,就会出现js调用不了android的代码,这是版本兼容的问题,需要在调用的方法上面加一个 ...

  3. 【深度学习系列】用PaddlePaddle和Tensorflow实现AlexNet

    上周我们用PaddlePaddle和Tensorflow实现了图像分类,分别用自己手写的一个简单的CNN网络simple_cnn和LeNet-5的CNN网络识别cifar-10数据集.在上周的实验表现 ...

  4. linux磁盘管理系列-软RAID的实现

    1 什么是RAID RAID全称是独立磁盘冗余阵列(Redundant Array of Independent Disks),基本思想是把多个磁盘组合起来,组合一个磁盘阵列组,使得性能大幅提高. R ...

  5. php命令执行脚本

    php -f jiaoben.php &  读入并解释指明的文件.

  6. 小白的Python之路 day1 模块初识

    模块初识 Python的强大之处在于他有非常丰富和强大的标准库和第三方库,几乎你想实现的任何功能都有相应的Python库支持,以后的课程中会深入讲解常用到的各种库,现在,我们先来象征性的学2个简单的. ...

  7. 对datatable进行简单的操作

    筛选出datatable中c_level=1的数据 dataRow[] rows = dt.Select("c_level=0"); 克隆表dt的结构到表dt,并将dt的数据复制到 ...

  8. reverse函数实现指定页面跳转

    需求: 在views中返回的url需要返回到具体的某一篇文章的评论列表 return redirect(reverse('cms_comment_manage',args=(number,))) dj ...

  9. 使用背景图修改radio、checkbox样式

    如果觉得设置样式太麻烦,或者页面上选中的样式太复杂,也可以用背景图去修改样式<div class=""> <label><input type=&qu ...

  10. 我的Android进阶之旅------&gt;Android知识图谱

    Android知识图谱.快来看看哪方面有漏洞? 该图转自:http://blog.csdn.net/xyz_lmn/article/details/41411355