[bzoj1565][NOI2009]植物大战僵尸_网络流_拓扑排序
植物大战僵尸 bzoj1565
题目大意:给你一张网格图,上面种着一些植物。你从网格的最右侧开始进攻。每个植物可以对僵尸提供能量或者消耗僵尸的能量。每个植物可以保护一个特定网格内的植物,如果一个植物被保护,那么如果僵尸想吃掉该植物就必须先吃掉保护它的植物。问:僵尸最多能获得多少能量。
注释:1<=N(网格的宽)<=20,1<=M(网格的长)<=30,-20,000<=代价和收益<=20,000.
想法:前置题目([NOI2006]最大获利)。这道题和最大获利比较相像,如果不仔细观察甚至可能认为tmCCF出了一道原题?但是,我们发现这个算法是完全错误的,因为如果A植物保护B植物,B植物也保护A植物,那么这两个植物是无敌的。这个无敌的定义就是从进攻开始到进攻结束这两个植物都不会有任何的问题。就是说,如果存在环的话,整个环之内的所有植物都是绝对无敌的。之后.... .... tm咋jb整啊?咳咳,这些无敌的小东西待着也没啥用,删了就得了。怎么删呢?拓扑排序啊!用toposort删掉不可能提供收益或代价的植物,之后裸最大权闭合子图即可。
最后,附上丑陋的代码... ...
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
#define inf 1000000000
#define ll long long
#define p(x,y) (x-1)*m+y
using namespace std;
int n,m,T;
int tot,top,cnt,ans;
int v[605],cur[605],h[605],d[605],q[605];
bool del[605];
struct edge{
int to,next,v;
}e[400005],ed[800005];int last[605],last2[605];
void insert(int u,int v)
{
d[v]++;
e[++cnt].to=v;e[cnt].next=last[u];last[u]=cnt;
}
void insert2(int u,int v,int w)
{
ed[++cnt].to=v;ed[cnt].next=last2[u];last2[u]=cnt;ed[cnt].v=w;
ed[++cnt].to=u;ed[cnt].next=last2[v];last2[v]=cnt;ed[cnt].v=0;
}
void dfs(int x)
{
del[x]=1;
for(int i=last[x];i;i=e[i].next)
if(!del[e[i].to])dfs(e[i].to);
}
void topsort()
{
for(int i=1;i<=n*m;i++)
if(!d[i])q[++top]=i;
else del[i]=1;
while(top)
{
int now=q[top];del[now]=0;
top--;
for(int i=last[now];i;i=e[i].next)
{
d[e[i].to]--;
if(!d[e[i].to])q[++top]=e[i].to;
}
}
for(int i=1;i<=n*m;i++)
if(del[i])dfs(i);
}
void rebuild()
{
cnt=1;T=n*m+1;
for(int x=1;x<=n*m;x++)
if(!del[x])
{
if(v[x]>0)tot+=v[x],insert2(x,T,v[x]);
else insert2(0,x,-v[x]);
for(int i=last[x];i;i=e[i].next)
if(!del[e[i].to])
insert2(x,e[i].to,inf);
}
}
bool bfs()
{
int head=0,tail=1;
for(int i=0;i<=T;i++)h[i]=-1;
h[0]=0;q[0]=0;
while(head!=tail)
{
int now=q[head];head++;
for(int i=last2[now];i;i=ed[i].next)
if(ed[i].v&&h[ed[i].to]==-1)
{
h[ed[i].to]=h[now]+1;
q[tail++]=ed[i].to;
}
}
return h[T]!=-1;
}
int dfs(int x,int f)
{
if(x==T)return f;
int w,used=0;
for(int i=cur[x];i;i=ed[i].next)
if(h[ed[i].to]==h[x]+1)
{
w=dfs(ed[i].to,min(ed[i].v,f-used));
ed[i].v-=w;ed[i^1].v+=w;
if(ed[i].v)cur[x]=i;
used+=w;if(used==f)return f;
}
if(!used)h[x]=-1;
return used;
}
void dinic()
{
while(bfs())
{
for(int i=0;i<=T;i++)cur[i]=last2[i];
ans+=dfs(0,inf);
}
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
{
scanf("%d",&v[p(i,j)]);
int x;
scanf("%d",&x);
while(x--)
{
int a,b;
scanf("%d%d",&a,&b);
a++;b++;
insert(p(i,j),p(a,b));
}
}
for(int i=1;i<=n;i++)
for(int j=m;j>1;j--)
insert(p(i,j),p(i,j-1));
topsort();
rebuild();
dinic();
printf("%d",tot-ans);
return 0;
}
小结:网络流真tm有意思
[bzoj1565][NOI2009]植物大战僵尸_网络流_拓扑排序的更多相关文章
- bzoj1565: [NOI2009]植物大战僵尸 最大权闭合子图,tarjan
bzoj1565: [NOI2009]植物大战僵尸 链接 https://www.lydsy.com/JudgeOnline/problem.php?id=1565 思路 很容易的想到最大权闭合子图 ...
- [BZOJ1565][NOI2009]植物大战僵尸-[网络流-最小割+最大点权闭合子图+拓扑排序]
Description 传送门 Solution em本题知识点是用网络流求最大点权闭合子图. 闭合图定义:图中任何一个点u,若有边u->v,则v必定也在图中. 建图:运用最小割思想,将S向点权 ...
- BZOJ1565 [NOI2009]植物大战僵尸(拓扑排序 + 最大权闭合子图)
题目 Source http://www.lydsy.com/JudgeOnline/problem.php?id=1565 Description Input Output 仅包含一个整数,表示可以 ...
- BZOJ1565[NOI2009]植物大战僵尸——最大权闭合子图+拓扑排序
题目描述 Plants vs. Zombies(PVZ)是最近十分风靡的一款小游戏.Plants(植物)和Zombies(僵尸)是游戏的主角,其中Plants防守,而Zombies进攻.该款游戏包含多 ...
- [luoguP2805] [NOI2009]植物大战僵尸(网络流)
传送门 结论:这是最大权闭合图的模型 因为可能A保护B,B保护A,出现环. 所以由植物A向植物A保护的植物连边,然后拓扑排序,将环去掉. 然后将拓扑排序的边反向连,建立最大权闭合图的模型. 跑最大流( ...
- BZOJ1565 [NOI2009]植物大战僵尸 【最大权闭合子图 + tarjan缩点(或拓扑)】
题目 输入格式 输出格式 仅包含一个整数,表示可以获得的最大能源收入.注意,你也可以选择不进行任何攻击,这样能源收入为0. 输入样例 3 2 10 0 20 0 -10 0 -5 1 0 0 100 ...
- b2OJ_1565_[NOI2009]植物大战僵尸_拓扑排序+最大权闭合子图
b2OJ_1565_[NOI2009]植物大战僵尸_拓扑排序+最大权闭合子 题意:n*m个植物,每个植物有分数(可正可负),和能保护植物的位置.只能从右往左吃,并且不能吃正被保护着的,可以一个不吃,求 ...
- 【BZOJ1565】【NOI2009】植物大战僵尸(网络流)
[BZOJ1565][NOI2009]植物大战僵尸(网络流) 题面 BZOJ 洛谷 题解 做了这么多神仙题,终于有一道能够凭借自己智商能够想出来的题目了.... 好感动. 这就是一个比较裸的最小割模型 ...
- 【BZOJ-1565】植物大战僵尸 拓扑排序 + 最小割
1565: [NOI2009]植物大战僵尸 Time Limit: 10 Sec Memory Limit: 64 MBSubmit: 1972 Solved: 917[Submit][Statu ...
随机推荐
- hi3531结构: VB_CONF_S
定义视频缓存池属性结构体. typedef struct hiVB_CONF_S { HI_U32 u32MaxPoolCnt; /* max count of pools, (0,VB ...
- R语言︱数据分组统计函数族——apply族用法与心得
每每以为攀得众山小,可.每每又切实来到起点,大牛们,缓缓脚步来俺笔记葩分享一下吧,please~ --------------------------- 笔者寄语:apply族功能强大,实用,可以代替 ...
- Java之split()方法
Java之split()方法 1.方法介绍 (1)public String[] split(String regex) 根据给定正则表达式的匹配拆分此字符串 (2)public String[] s ...
- 解决MyEclipse中的Building workspace问题
解决MyEclipse中的Building workspace问题 1.方法一 点击"Project",取消勾选"Build Automatically" 2. ...
- [Err] 1172 - Result consisted of more than one row
1 错误描述 [Err] 1172 - Result consisted of more than one row Procedure execution failed 1172 - Result c ...
- WebService之CXF注解之五(配置文件)
1.web.xml <?xml version="1.0" encoding="UTF-8"?> <web-app version=" ...
- AM335x(TQ335x)学习笔记——WM8960声卡驱动移植
经过一段时间的调试,终于调好了TQ335x的声卡驱动.TQ335x采用的Codec是WM8960,本文来总结下WM8960驱动在AM335x平台上的移植方法.Linux声卡驱动架构有OSS和ALSA两 ...
- .Net Core从命令行读取配置文件
最近在学习博客园腾飞(jesse)的.Net Core视频教程,收益匪浅,在此作推荐 : http://video.jessetalk.cn/ 言归正传,.Net Core应用程序中如何通过命令行读取 ...
- Coins HDU - 2844
Whuacmers use coins.They have coins of value A1,A2,A3...An Silverland dollar. One day Hibix opened p ...
- 爬虫抓包工具Charles设置
1.安装注册 感谢https://www.jianshu.com/p/89111882fa99提供注册码软件去官网下载安装即可. Registered Name:https://zhile.io Li ...