Loj #2321. 「清华集训 2017」无限之环
曾经有一款流行的游戏，叫做 *Infinity Loop***，先来简单的介绍一下这个游戏:
游戏在一个 (n imes m) 的网格状棋盘上进行，其中有些小方格中会有水管，水管可能在方格某些方向的边界的中点有接口，所有水管的粗细都相同，所以如果两个相邻方格的公共边界的中点都有接头，那么可以看作这两个接头互相连接。水管有以下 (15) 种形状：
游戏开始时，棋盘中水管可能存在漏水的地方。
形式化地:如果存在某个接头，没有和其它接头相连接，那么它就是一个漏水的地方。
玩家可以进行一种操作：选定一个含有*非直线型*水管的方格，将其中的水管绕方格中心顺时针或逆时针旋转 (90) 度。
直线型水管是指左图里中间一行的两种水管。
现给出一个初始局面，请问最少进行多少次操作可以使棋盘上不存在漏水的地方。
输入格式
第一行两个正整数 (n,m)，代表网格的大小。
接下来 (n) 行每行 (m) 个数，每个数是 ([0,15]) 中的一个，你可以将其看作一个 (4) 位的二进制数，从低到高每一位分别代表初始局面中这个格子上、右、
下、左方向上是否有 水管接头。
特别地，如果这个数是 (0)，则意味着这个位置没有水管。
比如 (3(0011_{(2)})) 代表上和右有接头，也就是一个 L 型，而 (12(1100_{(2)})) 代表下和左有接头，也就是将 L 型旋转 (180) 度。
输出格式
输出共一行，表示最少操作次数。如果无法达成目标，输出 (-1).
数据范围与提示
(n imes m le 2000)
好神仙啊!
大致思路就是用一个水管的旋转代替所有水管的旋转。
代码：
#include<bits/stdc++.h>
#define ll long long
#define N 10005
using namespace std;
inline int Get() {int x=0,f=1;char ch=getchar();while(ch<'0'||ch>'9') {if(ch=='-') f=-1;ch=getchar();}while('0'<=ch&&ch<='9') {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}return x*f;}
int n,m;
struct road {
int to,next;
int flow,cost;
}s[N*200];
int h[N],cnt=1;
void add(int i,int j,int f,int c) {
s[++cnt]=(road) {j,h[i],f,c};h[i]=cnt;
s[++cnt]=(road) {i,h[j],0,-c};h[j]=cnt;
}
int S,T;
int dis[N];
queue<int>q;
int e[N],fr[N];
bool in[N];
int ans,maxflow;
bool spfa() {
memset(dis,0x3f,sizeof(dis));
dis[S]=0;
q.push(S);
while(!q.empty()) {
int v=q.front();
q.pop();
in[v]=0;
for(int i=h[v];i;i=s[i].next) {
int to=s[i].to;
if(s[i].flow&&dis[to]>dis[v]+s[i].cost) {
dis[to]=dis[v]+s[i].cost;
fr[to]=v;
e[to]=i;
if(!in[to]) in[to]=1,q.push(to);
}
}
}
if(dis[T]>1e9) return 0;
for(int i=T;i!=S;i=fr[i]) {
s[e[i]].flow--;
s[e[i]^1].flow++;
}
maxflow++;
ans+=dis[T];
return 1;
}
vector<int>pipe;
int tot;
int ID[2005][2005][5];
int dx[]={-1,0,1,0},dy[]={0,1,0,-1};
int mp[2005][2005];
int main() {
n=Get(),m=Get();
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
mp[i][j]=Get();
T=n*m+1;
tot=T;
int lpipe=0,rpipe=0;
for(int i=1;i<=n;i++) {
for(int j=1;j<=m;j++) {
int id=(i-1)*m+j;
int a=mp[i][j];
pipe.clear();
for(int k=0;k<4;k++) if(a>>k&1) pipe.push_back(k);
for(int k=0;k<4;k++) ID[i][j][k]=++tot;
if(i+j&1) {
lpipe+=pipe.size();
add(S,id,pipe.size(),0);
for(int k=0;k<pipe.size();k++) add(id,ID[i][j][pipe[k]],1,0);
if(a==10||a==5||a==15) continue ;
if(pipe.size()==1) {
int now=pipe[0];
add(ID[i][j][now],ID[i][j][(now+1)%4],1,1);
add(ID[i][j][now],ID[i][j][(now+2)%4],1,2);
add(ID[i][j][now],ID[i][j][(now+3)%4],1,1);
} else if(pipe.size()==2) {
for(int k=0;k<pipe.size();k++) {
int now=pipe[k];
add(ID[i][j][now],ID[i][j][(now+2)%4],1,1);
}
} else if(pipe.size()==3) {
int now;
for(int k=0;k<4;k++) if(!(a>>k&1)) now=k;
for(int k=0;k<3;k++) {
if((pipe[k]+2)%4==now) add(ID[i][j][pipe[k]],ID[i][j][now],1,2);
else add(ID[i][j][pipe[k]],ID[i][j][now],1,1);
}
}
} else {
rpipe+=pipe.size();
add(id,T,pipe.size(),0);
for(int k=0;k<pipe.size();k++) add(ID[i][j][pipe[k]],id,1,0);
if(a==10||a==5||a==15) continue ;
if(pipe.size()==1) {
int now=pipe[0];
add(ID[i][j][(now+1)%4],ID[i][j][now],1,1);
add(ID[i][j][(now+2)%4],ID[i][j][now],1,2);
add(ID[i][j][(now+3)%4],ID[i][j][now],1,1);
} else if(pipe.size()==2) {
for(int k=0;k<pipe.size();k++) {
int now=pipe[k];
add(ID[i][j][(now+2)%4],ID[i][j][now],1,1);
}
} else if(pipe.size()==3) {
int now;
for(int k=0;k<4;k++) if(!(a>>k&1)) now=k;
for(int k=0;k<3;k++) {
if((pipe[k]+2)%4==now) add(ID[i][j][now],ID[i][j][pipe[k]],1,2);
else add(ID[i][j][now],ID[i][j][pipe[k]],1,1);
}
}
}
}
}
for(int i=1;i<=n;i++) {
for(int j=1;j<=m;j++) {
if(i+j&1) {
for(int d=0;d<4;d++) {
int X=i+dx[d],Y=j+dy[d];
if(1<=X&&X<=n&&1<=Y&&Y<=m) {
add(ID[i][j][d],ID[X][Y][(d+2)%4],1,0);
}
}
}
}
}
if(lpipe!=rpipe) cout<<-1;
else {
while(spfa());
if(maxflow!=lpipe) cout<<-1;
else cout<<ans;
}
return 0;
}