BZOJ 4261: 建设游乐场
4261: 建设游乐场
Time Limit: 50 Sec Memory Limit: 256 MB
Submit: 38 Solved: 16
[Submit][Status][Discuss]
Description

Input
Output
Sample Input
1 1 1
1 0 0
1 0 0
48 94 1
78 78 81
1 12 60
Sample Output
HINT
N<=150,M<=30,Vi,j<=100
Source
分析:
好题~~~
考虑构成若干个环的必要条件:每个格子的管道都要链接相邻的两个格子...然后把一个格子拆成两个管子,如果两个管子方向不同,则可以得到价值,如果相同就无法得到...
既然是网格图,那么我们考虑黑白染色,对于每个必须放的格子,我们把这个格子拆成三个点,一个点叫做控制点,一个代表竖向管道,一个代表横向管道...
然后对于每个黑点,从$S$向黑点的控制点连边$<S,id[i][j][0],2,0>$,代表连接两个格子,$<id[i][j][0],id[i][j][1],1,val>$,$<id[i][j][0],id[i][j][1],id[i][j][2],1,val>$,$<id[i][j][0],id[i][j][1],1,0>$,$<id[i][j][0],id[i][j][2],1,0>$...
对于每个白点,从白点的控制点向$T$连边$<id[i][j][0],T,2,0>$,然后$1$、$2$的连边和黑点相反...
然后对于每个黑点,我们把竖向点向上下的相邻格子个竖点连边$<1,0>$,横向点向左右的横向点连边$<1,0>$...
这样,如果满流并且黑白点相同的话就一定是合法解,因为代表每个格子都连接了相邻的两个格子...然后如果一个各自选择了相同方向的管道,那么会得到$val$的收益,如果选择不同的,会得到$2*val$的收益,所以求出最大费用之后减去$\sum val$就好了...
代码:
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<queue>
//by NeighThorn
#define inf 0x3f3f3f3f
#define min(a,b) a>b?b:a
using namespace std; const int N=150+5,maxn=50000+5,maxm=1000000+5; int n,m,S,T,cnt,tat,tot,sum,mp[N][N],id[N][N][3],val[N][N];
int w[maxm],fl[maxm],hd[maxn],to[maxm],nxt[maxm],Min[maxn],dis[maxn],vis[maxn],from[maxn];
int mv[4][2]={1,0,-1,0,0,1,0,-1}; inline int read(void){
char ch=getchar();int x=0;
while(!(ch>='0'&&ch<='9')) ch=getchar();
while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
return x;
} inline void add(int x,int y,int l,int s){
w[cnt]= s;fl[cnt]=l;to[cnt]=y;nxt[cnt]=hd[x];hd[x]=cnt++;
w[cnt]=-s;fl[cnt]=0;to[cnt]=x;nxt[cnt]=hd[y];hd[y]=cnt++;
} inline bool spfa(void){
memset(Min,inf,sizeof(Min));
memset(dis,inf,sizeof(dis));
queue<int> q;q.push(S);dis[S]=0;vis[S]=1;
while(!q.empty()){
int top=q.front();q.pop();vis[top]=0;
for(int i=hd[top];i!=-1;i=nxt[i])
if(dis[to[i]]>dis[top]+w[i]&&fl[i]){
from[to[i]]=i;
dis[to[i]]=dis[top]+w[i];
Min[to[i]]=min(Min[top],fl[i]);
if(!vis[to[i]])
vis[to[i]]=1,q.push(to[i]);
}
}
return dis[T]!=inf;
} inline int find(void){
for(int i=T;i!=S;i=to[from[i]^1])
fl[from[i]]-=Min[T],fl[from[i]^1]+=Min[T];
return Min[T]*dis[T];
} inline int mcmf(void){
int res=0,flow=0;
while(spfa())
res+=find(),flow+=Min[T];
if(flow!=tot) return -1;
return -res-sum;
} signed main(void){
n=read();m=read();
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
mp[i][j]=read();
for(int i=0;i<=n+1;i++) mp[i][0]=mp[i][m+1]=1;
for(int i=0;i<=m+1;i++) mp[0][i]=mp[n+1][i]=1;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++){
val[i][j]=read();
if(!mp[i][j])
sum+=val[i][j];
}
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
if(!mp[i][j])
for(int k=0;k<3;k++)
id[i][j][k]=++tot;
S=0,T=++tot;tot=0;
memset(hd,-1,sizeof(hd));
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
if(!mp[i][j]){
if((i+j)&1)
tat+=2,
add(S,id[i][j][0],2,0),
add(id[i][j][0],id[i][j][1],1,0),
add(id[i][j][0],id[i][j][2],1,0),
add(id[i][j][0],id[i][j][1],1,-val[i][j]),
add(id[i][j][0],id[i][j][2],1,-val[i][j]);
else
tot+=2,
add(id[i][j][0],T,2,0),
add(id[i][j][1],id[i][j][0],1,0),
add(id[i][j][2],id[i][j][0],1,0),
add(id[i][j][1],id[i][j][0],1,-val[i][j]),
add(id[i][j][2],id[i][j][0],1,-val[i][j]);
}
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
if(!mp[i][j]&&((i+j)&1)){
for(int k=0,x,y;k<4;k++){
x=i+mv[k][0],y=j+mv[k][1];
if(!mp[x][y]){
if(k<2)
add(id[i][j][1],id[x][y][1],1,0);
else
add(id[i][j][2],id[x][y][2],1,0);
}
}
}
if(tat!=tot) return puts("-1"),0;
printf("%d\n",mcmf());
return 0;
}
By NeighThorn
BZOJ 4261: 建设游乐场的更多相关文章
- bzoj 4261: 建设游乐场 费用流
题目 现在有一大块土地,可以看成N*M的方格.在这块土地上,有些格子内是崎岖的山地,无法建造任何东西:其他格子都是平原.现在打算在这块土地上建设一个游乐园.游乐园由若干条闭合的过山车轨道组成,每个平原 ...
- 【BZOJ-4261】建设游乐场 最大费用最大流
4261: 建设游乐场 Time Limit: 50 Sec Memory Limit: 256 MBSubmit: 21 Solved: 8[Submit][Status][Discuss] D ...
- BZOJ4261 : 建设游乐场
将图黑白染色,每个点拆成两个点,分别表示水平和竖直方向,再增加一个点以控制流量,那么每个格子都需要找两个方向去连接. $S$到每个黑点的控制点连边,流量$2$,费用$0$: 控制点向两个方向的点各连两 ...
- bzoj AC倒序
Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...
- BZOJ 1096 【ZJOI2007】 仓库建设
Description L公司有N个工厂,由高到底分布在一座山上.如图所示,工厂1在山顶,工厂N在山脚.由于这座山处于高原内陆地区(干燥少雨),L公司一般把产品直接堆放在露天,以节省费用.突然有一天, ...
- 【BZOJ 1096】 [ZJOI2007]仓库建设 (斜率优化)
1096: [ZJOI2007]仓库建设 Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 3940 Solved: 1736 Description ...
- bzoj 1096 仓库建设 -斜率优化
L公司有N个工厂,由高到底分布在一座山上.如图所示,工厂1在山顶,工厂N在山脚.由于这座山处于高原内陆地区(干燥少雨),L公司一般把产品直接堆放在露天,以节省费用.突然有一天,L公司的总裁L先生接到气 ...
- 【BZOJ】1096: [ZJOI2007]仓库建设(dp+斜率优化)
http://www.lydsy.com/JudgeOnline/problem.php?id=1096 首先得到dp方程(我竟然自己都每推出了QAQ)$$d[i]=min\{d[j]+cost(j+ ...
- BZOJ 2001: [Hnoi2010]City 城市建设
2001: [Hnoi2010]City 城市建设 Time Limit: 20 Sec Memory Limit: 162 MBSubmit: 1132 Solved: 555[Submit][ ...
随机推荐
- 【转载】java 客户端链接不上redis解决方案 (jedis)
本文出自:http://blog.csdn.net/lulidaitian/article/details/51946169 出现问题描述: 1.Could not get a resource fr ...
- Django自带后台管理配置
Django自带后台管理的配置 创建项目和应用 修改配置文件 数据库配置 DATABASES = { 'default': { 'ENGINE': 'django.db.backends.mysql' ...
- C语言函数篇(五)静态库和动态库的创建和使用
使用库函数是源码的一种保护??? <我猜的.> 库函数其实不是新鲜的东西,我们一直都在用,比如C库. 我们执行pringf() 这个函数的时候,就是调用C库的函数. 下面记录静态库和动态库 ...
- 第三章习题 C++ Primer 第六版
1.使用一个整数输入自己的身高(单位为cm),并将此身高转化为米和厘米共同表示的形式,使用下划线字符来指示输入的位置,使用一个const符号常量来表示转换因子. #include<iostrea ...
- Aizu:0189-Convenient Location
Convenient Location Time limit 1000 ms Memory limit 131072 kB Problem Description 明年毕业的A为就业而搬家.就职的公司 ...
- Nginx模块详解
Nginx模块介绍 核心模块:core module 标准模块: HTTP modules: Standard HTTP modules Optional HTTP modules Mail modu ...
- 大话卷积神经网络(CNN)
这几年深度学习快速发展,在图像识别.语音识别.物体识别等各种场景上取得了巨大的成功,例如AlphaGo击败世界围棋冠军,iPhone X内置了人脸识别解锁功能等等,很多AI产品在世界上引起了很大的 ...
- ArrayList & Vector的源码实现
#ArrayList & Vector #####前言: 本来按照计划,ArrayList和Vector是分开讲的,但是当我阅读了ArrayList和Vector的源码以后,我就改变了注意,把 ...
- java实时监听日志写入kafka(多目录)
目的 实时监听多个目录下的日志文件,如有新文件切换到新文件,并同步写入kafka,同时记录日志文件的行位置,以应对进程异常退出,能从上次的文件位置开始读取(考虑到效率,这里是每100条记一次,可调整) ...
- 3,Flask 中的模板语言 Jinja2 及 render_template 的深度用法
Flask中默认的模板语言是Jinja2 现在我们来一步一步的学习一下 Jinja2 捎带手把 render_template 中留下的疑问解决一下 首先我们要在后端定义几个字符串,用于传递到前端 S ...