BZOJ 5120: [2017国家集训队测试]无限之环(费用流)
解题思路
神仙题。调了一个晚上+半个上午。。这道咋看咋都不像图论的题竟然用费用流做,将行+列为奇数的点和偶数的点分开,也就是匹配问题,然后把一个点复制四份,分别代表这个点的上下左右接头,如果有这个接头就加一个费用为\(0\),流量为\(1\)的边,如果没有要分情况讨论,因为从源点到这个点的流量是固定的,当只有一个接头时,可以让这个点向自己其余三个点连费用为\(1\),流量为\(1\)的边,当有两个接头并且两个接头相邻时,让这个点的两个接头分别与对应的方向连边,当有三个接头时,让那个没有的接头向相邻的连费用为\(1\)的边,向相对的连费用为\(2\)的边。然后一边费用流就行了。代码比较丑。建图写挫了。。
代码
#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<cstdlib>
using namespace std;
const int MAXN = 10005;
const int MAXM = 100005;
const int inf = 0x3f3f3f3f;
const int zz[4] = {1,2,4,8}; // 0 1 2 3
inline int rd(){
int x=0,f=1;char ch=getchar();
while(!isdigit(ch)) f=ch=='0'?0:1,ch=getchar();
while(isdigit(ch)) x=(x<<1)+(x<<3)+ch-'0',ch=getchar();
return f?x:-x;
}
int n,m,head[MAXN],to[MAXM<<1],nxt[MAXM<<1],val[MAXM<<1],cost[MAXM<<1],pre[MAXN];
int S,T,cnt=1,maxflow,ans,dis[MAXN],ch[MAXN][5],num,tot,sum,pos[2005][2005];
int incf[MAXN],all;
bool vis[MAXN];
queue<int> Q;
inline void add(int bg,int ed,int z,int w){
to[++cnt]=ed,nxt[cnt]=head[bg],val[cnt]=z,cost[cnt]=w,head[bg]=cnt;
}
inline void solve1(int p,int x,int op){
for(int i=0;i<=3;i++)
if((x&zz[i])) {
if(op&1) add(p,ch[p][i],1,0),add(ch[p][i],p,0,0);
else add(ch[p][i],p,1,0),add(p,ch[p][i],0,0);
if(i<2) {
if(op&1) add(p,ch[p][i+2],1,2),add(ch[p][i+2],p,0,-2);
else add(ch[p][i+2],p,1,2),add(p,ch[p][i+2],0,-2);
}
else {
if(op&1) add(p,ch[p][i-2],1,2),add(ch[p][i-2],p,0,-2);
else add(ch[p][i-2],p,1,2),add(p,ch[p][i-2],0,-2);
}
if(op&1){
add(p,ch[p][(i+1)%4],1,1),add(ch[p][(i+1)%4],p,0,-1);
add(p,ch[p][(i+3)%4],1,1),add(ch[p][(i+3)%4],p,0,-1);
}
else{
add(ch[p][(i+1)%4],p,1,1),add(p,ch[p][(i+1)%4],0,-1);
add(ch[p][(i+3)%4],p,1,1),add(p,ch[p][(i+3)%4],0,-1);
}
}
}
inline void solve2(int p,int x,int op){
if(x==5 || x==10){
for(int i=0;i<=3;i++)
if(x&zz[i]){
if(op&1) add(p,ch[p][i],1,0),add(ch[p][i],p,0,0);
else add(ch[p][i],p,1,0),add(p,ch[p][i],0,0);
}
return ;
}
for(int i=0;i<=3;i++){
if(x&zz[i]){
if(op&1) add(p,ch[p][i],1,0),add(ch[p][i],p,0,0);
else add(ch[p][i],p,1,0),add(p,ch[p][i],0,0);
}
else{
if(i<2) {
if((op&1))
add(ch[p][i+2],ch[p][i],1,1),add(ch[p][i],ch[p][i+2],0,-1);
else
add(ch[p][i],ch[p][i+2],1,1),add(ch[p][i+2],ch[p][i],0,-1);
}
else {
if((op&1))
add(ch[p][i-2],ch[p][i],1,1),add(ch[p][i],ch[p][i-2],0,-1);
else
add(ch[p][i],ch[p][i-2],1,1),add(ch[p][i-2],ch[p][i],0,-1);
}
}
}
}
inline void solve3(int p,int x,int op){
for(int i=0;i<=3;i++){
if(x&zz[i]){
if(op&1) add(p,ch[p][i],1,0),add(ch[p][i],p,0,0);
else add(ch[p][i],p,1,0),add(p,ch[p][i],0,0);
}
else {
if((op&1)){
add(ch[p][(i+1)%4],ch[p][i],1,1);
add(ch[p][i],ch[p][(i+1)%4],0,-1);
add(ch[p][(i+3)%4],ch[p][i],1,1);
add(ch[p][i],ch[p][(i+3)%4],0,-1);
if(i<2) {
add(ch[p][i+2],ch[p][i],1,2);
add(ch[p][i],ch[p][i+2],0,-2);
}
else {
add(ch[p][i-2],ch[p][i],1,2);
add(ch[p][i],ch[p][i-2],0,-2);
}
}
else{
add(ch[p][i],ch[p][(i+1)%4],1,1);
add(ch[p][(i+1)%4],ch[p][i],0,-1);
add(ch[p][i],ch[p][(i+3)%4],1,1);
add(ch[p][(i+3)%4],ch[p][i],0,-1);
if(i<2) {
add(ch[p][i],ch[p][i+2],1,2);
add(ch[p][i+2],ch[p][i],0,-2);
}
else {
add(ch[p][i],ch[p][i-2],1,2);
add(ch[p][i-2],ch[p][i],0,-2);
}
}
}
}
}
inline void solve4(int p,int op){
for(int i=0;i<=3;i++){
if(op&1) add(p,ch[p][i],1,0),add(ch[p][i],p,0,0);
else add(ch[p][i],p,1,0),add(p,ch[p][i],0,0);
}
}
bool spfa(){
while(Q.size()) Q.pop();
memset(dis,0x3f,sizeof(dis));
memset(vis,false,sizeof(vis));
vis[S]=1;dis[S]=0;incf[S]=inf;Q.push(S);
while(Q.size()){
int x=Q.front();Q.pop();vis[x]=0;
for(int i=head[x];i;i=nxt[i]){
int u=to[i];
if(dis[u]>cost[i]+dis[x] && val[i]) {
dis[u]=cost[i]+dis[x];
incf[u]=min(incf[x],val[i]);
pre[u]=i;
if(!vis[u]) vis[u]=1,Q.push(u);
}
}
}
return (dis[T]==inf)?false:true;
}
inline void update(){
int x=T,i;
while(x!=S){
i=pre[x];
val[i]-=incf[T];
val[i^1]+=incf[T];
x=to[i^1];
}
maxflow+=incf[T];
ans+=incf[T]*dis[T];
}
int main(){
n=rd(),m=rd();S=++num;T=++num;int x,p;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++){
x=rd();p=++num;pos[i][j]=p;
for(int k=0;k<=3;k++) ch[p][k]=++num;
tot=__builtin_popcount(x);
if((i+j)&1) add(S,p,tot,0),add(p,S,0,0);
else add(p,T,tot,0),add(T,p,0,0);
if(tot==1) solve1(p,x,i+j);
else if(tot==2) solve2(p,x,i+j);
else if(tot==3) solve3(p,x,i+j);
else solve4(p,i+j);
sum+=tot;
}
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++){
if(!((i+j)&1)) continue;
if(i!=1){
add(ch[pos[i][j]][0],ch[pos[i-1][j]][2],1,0);
add(ch[pos[i-1][j]][2],ch[pos[i][j]][0],0,0);
}
if(j!=1){
add(ch[pos[i][j]][3],ch[pos[i][j-1]][1],1,0);
add(ch[pos[i][j-1]][1],ch[pos[i][j]][3],0,0);
}
if(i!=n){
add(ch[pos[i][j]][2],ch[pos[i+1][j]][0],1,0);
add(ch[pos[i+1][j]][0],ch[pos[i][j]][2],0,0);
}
if(j!=m){
add(ch[pos[i][j]][1],ch[pos[i][j+1]][3],1,0);
add(ch[pos[i][j+1]][3],ch[pos[i][j]][1],0,0);
}
}
if(sum&1) {puts("-1");return 0;}
sum>>=1;
while(spfa()) update();
if(maxflow!=sum) puts("-1");
else printf("%d\n",ans);
return 0;
}
BZOJ 5120: [2017国家集训队测试]无限之环(费用流)的更多相关文章
- bzoj 5120 [2017国家集训队测试]无限之环——网络流
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=5120 旋转的话相当于去掉一个插头.新增一个插头,所以在这两个插头之间连边并带上费用即可. 网 ...
- bzoj 5120: [2017国家集训队测试]无限之环【最小费用最大流】
玄妙的建图-- 这种平衡度数的题按套路是先黑白染色然后分别连ST点,相邻格子连黑向白连费用1流量0的边,然后考虑费用怎么表示 把一个点拆成五个,上下左右中,中间点黑白染色连ST, 对于连S的点,中点连 ...
- BZOJ5120 [2017国家集训队测试]无限之环 费用流
欢迎访问~原文出处——博客园-zhouzhendong 去博客园看该题解 题目传送门 - BZOJ5120 题意概括 原题挺简略的. 题解 本题好难. 听了任轩笛大佬<国家队神犇>的讲课才 ...
- [BZOJ5120] [2017国家集训队测试]无限之环
Description 曾经有一款流行的游戏,叫做InfinityLoop,先来简单的介绍一下这个游戏: 游戏在一个n×m的网格状棋盘上进行,其中有些小方格中会有水管,水管可能在方格某些方向的边界的中 ...
- BZOJ.5120.[清华集训2017]无限之环(费用流zkw 黑白染色)
题目链接 LOJ 洛谷 容易想到最小费用最大流分配度数. 因为水管形态固定,每个点还是要拆成4个点,分别当前格子表示向上右下左方向. 然后能比较容易地得到每种状态向其它状态转移的费用(比如原向上的可以 ...
- BZOJ_2622_[2012国家集训队测试]深入虎穴_最短路
BZOJ_2622_[2012国家集训队测试]深入虎穴_最短路 Description 虎是中国传统文化中一个独特的意象.我们既会把老虎的形象用到喜庆的节日装饰画上,也可能把它视作一种邪恶的可怕的动物 ...
- 【BZOJ2622】[2012国家集训队测试]深入虎穴 次短路
[BZOJ2622][2012国家集训队测试]深入虎穴 Description 虎是中国传统文化中一个独特的意象.我们既会把老虎的形象用到喜庆的节日装饰画上,也可能把它视作一种邪恶的可怕的动物,例如“ ...
- 2017国家集训队作业Atcoder题目试做
2017国家集训队作业Atcoder题目试做 虽然远没有达到这个水平,但是据说Atcoder思维难度大,代码难度小,适合我这种不会打字的选手,所以试着做一做 不知道能做几题啊 在完全自己做出来的题前面 ...
- 2017国家集训队作业[agc016b]Color Hats
2017国家集训队作业[agc016b]Color Hats 题意: 有\(N\)个人,每个人有一顶帽子.帽子有不同的颜色.现在,每个人都告诉你,他看到的所有其它人的帽子共有多少种颜色,问有没有符合所 ...
随机推荐
- 饿了么监控系统 EMonitor 与美团点评 CAT 的对比
背景介绍 饿了么监控系统EMonitor:是一款服务于饿了么所有技术部门的一站式监控系统,覆盖了系统监控.容器监控.网络监控.中间件监控.业务监控.接入层监控以及前端监控的数据存储与查询.每日处理总数 ...
- Arthas 开源一周年,GitHub Star 16 K ,我们一直在坚持什么?
缘起 最近看到一个很流行的标题,<开源XX年,star XXX,我是如何坚持的>.看到这样的标题,忽然发觉 Arthas 从 2018 年 9 月开源以来,刚好一年了,正好在这个秋高气爽的 ...
- Mybatis基于SqlSession实现CRUD
之前我们讲的基于XML还是接口注解配置都是使用接口实现CRUD,本文我们将要讲解通过splsession来实现CRUD,这种方法比较灵活. 基本配置 <!-- spring和MyBatis完美整 ...
- I2C用户态驱动设计
一.用户态驱动模型 1.1 I2C通用驱动代码 i2c_dev_init: static int __init i2c_dev_init(void) { int res; printk(KERN_IN ...
- Android学习--apk打包过程
1. 使用aapt工具,给所有的res目录下的资源文件生成对应的id,id会被放进R.java文件中 2. JavaC编译器,将所有Java文件转换为Class文件,其中,内部类会分别生成.class ...
- [MySQL] innobackupex在线备份及恢复(全量和增量)
安装percona-xtrabackup 方法1: percona-xtrabackup-2.1.9-744-Linux-x86_64.tar.gz(D:\share\src\linux-mysql) ...
- MySQL 下载,安装,配置windows 服务
本次使用的是压缩包的方式是可以纯手动自己折腾各种配置... ok,闲话少叙,我们准备发车... 一.先要去mysql官网去下载压缩包咯 ①下载地址:https://dev.mysql.com/down ...
- SecureCRT key登录linux ssh设置
一.首先用secureCrt创建密钥 1.使用SecureCRT创建私钥和公钥. SecureCRT quick Connect-> Authentiation -> Public Key ...
- 力扣算法题—145BinartTreePostorderTraversal
Given a binary tree, return the postorder traversal of its nodes' values. Example: Input: [1,null,2, ...
- 如何将英文版的Firefox添加中文版语言包
http://ftp.mozilla.org/pub/firefox/releases/ xpi中下载zh_CN.xpi 文件 , 把文件拖拽进火狐浏览器 在地址栏输入”about:config”,回 ...