题目

Source

http://www.lydsy.com/JudgeOnline/problem.php?id=1565

3 2
10 0
20 0
-10 0
-5 1 0 0
100 1 2 1
100 0

25

代码

```#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
#include<algorithm>
using namespace std;
#define INF (1<<30)
#define MAXN 666
#define MAXM 666*666*2

struct Edge{
int v,cap,flow,next;
}edge[MAXM];
int vs,vt,NE,NV;

edge[NE].v=v; edge[NE].cap=cap; edge[NE].flow=0;
edge[NE].v=u; edge[NE].cap=0; edge[NE].flow=0;
}

int level[MAXN];
int gap[MAXN];
void bfs(){
memset(level,-1,sizeof(level));
memset(gap,0,sizeof(gap));
level[vt]=0;
gap[level[vt]]++;
queue<int> que;
que.push(vt);
while(!que.empty()){
int u=que.front(); que.pop();
int v=edge[i].v;
if(level[v]!=-1) continue;
level[v]=level[u]+1;
gap[level[v]]++;
que.push(v);
}
}
}

int pre[MAXN];
int cur[MAXN];
int ISAP(){
bfs();
memset(pre,-1,sizeof(pre));
int u=pre[vs]=vs,flow=0,aug=INF;
gap[0]=NV;
while(level[vs]<NV){
bool flag=false;
for(int &i=cur[u]; i!=-1; i=edge[i].next){
int v=edge[i].v;
if(edge[i].cap!=edge[i].flow && level[u]==level[v]+1){
flag=true;
pre[v]=u;
u=v;
//aug=(aug==-1?edge[i].cap:min(aug,edge[i].cap));
aug=min(aug,edge[i].cap-edge[i].flow);
if(v==vt){
flow+=aug;
for(u=pre[v]; v!=vs; v=u,u=pre[u]){
edge[cur[u]].flow+=aug;
edge[cur[u]^1].flow-=aug;
}
//aug=-1;
aug=INF;
}
break;
}
}
if(flag) continue;
int minlevel=NV;
int v=edge[i].v;
if(edge[i].cap!=edge[i].flow && level[v]<minlevel){
minlevel=level[v];
cur[u]=i;
}
}
if(--gap[level[u]]==0) break;
level[u]=minlevel+1;
gap[level[u]]++;
u=pre[u];
}
return flow;
}

vector<int> G[MAXN];
int score[MAXN],deg[MAXN];

int main(){
int n,m,a,b,c;
scanf("%d%d",&n,&m);
for(int i=0; i<n*m; ++i){
scanf("%d",&score[i]);
scanf("%d",&c);
while(c--){
scanf("%d%d",&a,&b);
G[i].push_back(a*m+b);
++deg[a*m+b];
}
if(i%m){
G[i].push_back(i-1);
++deg[i-1];
}
}
vs=n*m; vt=vs+1; NV=vt+1; NE=0;
int tot=0;
queue<int> que;
for(int i=0; i<vs; ++i){
if(deg[i]==0) que.push(i);
}
while(!que.empty()){
int u=que.front(); que.pop();
if(score[u]<0){
}else{
tot+=score[u];
}
for(int i=0; i<G[u].size(); ++i){
int v=G[u][i];
--deg[v];
if(deg[v]==0) que.push(v);
}
}
printf("%d",tot-ISAP());
return 0;
}
```

BZOJ1565 [NOI2009]植物大战僵尸（拓扑排序 + 最大权闭合子图）的更多相关文章

1. b2OJ_1565_[NOI2009]植物大战僵尸_拓扑排序+最大权闭合子图

b2OJ_1565_[NOI2009]植物大战僵尸_拓扑排序+最大权闭合子 题意:n*m个植物,每个植物有分数(可正可负),和能保护植物的位置.只能从右往左吃,并且不能吃正被保护着的,可以一个不吃,求 ...

2. 【BZOJ-1565】植物大战僵尸 拓扑排序 + 最小割

1565: [NOI2009]植物大战僵尸 Time Limit: 10 Sec  Memory Limit: 64 MBSubmit: 1972  Solved: 917[Submit][Statu ...

3. [bzoj1565][NOI2009]植物大战僵尸_网络流_拓扑排序

植物大战僵尸 bzoj1565 题目大意:给你一张网格图,上面种着一些植物.你从网格的最右侧开始进攻.每个植物可以对僵尸提供能量或者消耗僵尸的能量.每个植物可以保护一个特定网格内的植物,如果一个植物被 ...

4. Bzoj 1565: [NOI2009]植物大战僵尸 最大权闭合图,拓扑排序

题目: http://cojs.tk/cogs/problem/problem.php?pid=410 410. [NOI2009] 植物大战僵尸 ★★★   输入文件:pvz.in   输出文件:p ...

5. BZOJ 1565: [NOI2009]植物大战僵尸( 最小割 )

先拓扑排序搞出合法的, 然后就是最大权闭合图模型了.... --------------------------------------------------------------------- ...

6. bzoj 1565 最大权闭合子图

因为每个植物都有保护的点(每排相邻的右面的算是保护左面的),所以连他和保护 的点一条边,然后每个点有自己的权值,要找到最大的权值且满足每个点在访问时他 的前驱一定被访问,那么反向建边,转化为后继必须访 ...

7. HDU 3879 Base Station（最大权闭合子图）

经典例题,好像说可以转化成maxflow(n,n+m),暂时只可以勉强理解maxflow(n+m,n+m)的做法. 题意:输入n个点,m条边的无向图.点权为负,边权为正,点权为代价,边权为获益,输出最 ...

8. [BZOJ 1497][NOI 2006]最大获利（最大权闭合子图）

题目:http://www.lydsy.com:808/JudgeOnline/problem.php?id=1497 分析: 这是在有向图中的问题,且边依赖于点,有向图中存在点.边之间的依赖关系可以 ...

9. HDU4971 A simple brute force problem.（强连通分量缩点 + 最大权闭合子图）

题目 Source http://acm.hdu.edu.cn/showproblem.php?pid=4971 Description There's a company with several ...

随机推荐

1. 解析opencv中Box Filter的实现并提出进一步加速的方案（源码共享）。

说明:本文所有算法的涉及到的优化均指在PC上进行的,对于其他构架是否合适未知,请自行试验. Box Filter,最经典的一种领域操作,在无数的场合中都有着广泛的应用,作为一个很基础的函数,其性能的好 ...

2. Servlet规范简介——web框架是如何注入到Servlet中的

Servlet规范简介--web框架是如何注入到Servlet中的 引言 Web框架一般是通过一个Servlet提供统一的请求入口,将指定的资源映射到这个servlet,在这个servlet中进行框架 ...

题目链接:http://acm.hust.edu.cn/vjudge/contest/126708#problem/E 题意:给做出的题目个数,5个的100分,4个的前n/2的同学95,后n/2的90 ...

5. Midnight.js – 实现奇妙的固定头部切换效果

Midnight.js 是一款 jQuery 插件,在页面滚动的时候实现多个头设计之间的切换,所以你总是有一个头与它下面的内容层叠,看起来效果很不错. Midnight.js 可以让你轻松实现这种切换 ...

6. UWP开发中两种网络图片缓存方法

通常情况下,我们的程序需要从服务器读取图片,但如果需要不止一次读取某一张图片的话,就需要做本地缓存了,这样既为用户省一点流量,又能显得你的APP很快. 假如你已经知道了某一张图片的地址,那么第一件事就 ...

7. CentOS6.4 安装SSDB

1.安装 wget --no-check-certificate https://github.com/ideawu/ssdb/archive/master.zipunzip mastercd ssd ...

8. POJ 3687 Labeling Balls()

Labeling Balls Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 9641 Accepted: 2636 Descri ...

9. UBUNTU查看软件版本

1.查看已安装软件版本aptitude show softwarename 2.查看软件安装目录dpkg -L softwarename

10. 数据绑定表达式(上)：.NET发现之旅(一)

数据绑定表达式(上):.NET发现之旅(一) 2009-06-30 10:29:06 来源:网络转载 作者:佚名 共有评论(0)条 浏览次数:859 作为.NET平台软件开发者,我们频繁与各种各样的数 ...