Bzoj:1324:Exca王者之剑:网络流，最小割

题目链接：1324:Exca王者之剑

发现我们必然是偶数秒到一颗宝石上才能拿到它，同时他周围的宝石都没了，也就是一颗宝石和他周围的宝石并不能同时拿到

因为一颗宝石和他周围的宝石不能同时拿到，所以我们将原图黑白染色，相邻的黑白格子不能同时到达

那么黑格子向S连边，白格子向T连边，相邻的格子连边后跑一个最小割，代表最少要舍弃多少，相减即可

尝试变了下码风，神清气爽QwQ

#include<queue>
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define rep(i,x,y) for (int i=x;i<=y;++i)
#define rep2(i,x,y) for (int i=x;i<=y;i+=2)
#define nxt(i,x) for (int i=x;i;i=G[i].next)
using namespace std;
const int inf=0x7fffffff/2-1;
const int maxs=11110;
int n,m,h[maxs],tot=1,vis[maxs];
struct edge{int to,next,w;}G[maxs*10];
int a[110][110],b[110][110];
int S,T,cur[maxs],cnt=0;
int dx[]={-1,0,1,0};
int dy[]={0,-1,0,1};

bool bfs(){
rep(i,0,T) vis[i]=-1;
queue<int>q;q.push(S);vis[S]=0;
while (!q.empty()){
int u=q.front(); q.pop();
nxt(i,h[u]){
int v=G[i].to;
if (vis[v]==-1&&G[i].w){
vis[v]=vis[u]+1;
q.push(v);
}
}
}return vis[T]!=-1;
}

int dfs(int x,int f){
if(x==T||f==0)return f; int w,used=0;
nxt(i,cur[x])
if(G[i].w&&vis[G[i].to]==vis[x]+1){
w=f-used; w=dfs(G[i].to,min(w,G[i].w));
G[i].w-=w; G[i^1].w+=w; used+=w;
if (G[i].w) cur[x]=i;
if (used==f) return used;
}if (!used) vis[x]=-1;
return used;
}

void dinic(){
int ret=0;
while (bfs()){
rep(i,0,T) cur[i]=h[i];
ret+=dfs(S,inf);
}printf("%d",cnt-ret);
}

int calc(int x,int y){return (x-1)*m+y;}

void add(int x,int y,int z){
G[++tot].to=y;G[tot].next=h[x];h[x]=tot;G[tot].w=z;
G[++tot].to=x;G[tot].next=h[y];h[y]=tot;G[tot].w=0;
}

int main(){
scanf("%d%d",&n,&m); S=0; T=n*m+1;
rep(i,1,n) rep(j,1,m) scanf("%d",&a[i][j]),cnt+=a[i][j];
rep(i,1,n) rep2(j,i%2+1,m) b[i][j]=1;
rep(i,1,n) rep(j,1,m)
if (b[i][j]) add(S,calc(i,j),a[i][j]);
else add(calc(i,j),T,a[i][j]);
rep(i,1,n) rep(j,1,m)if(b[i][j])rep(k,0,3){
int x=i+dx[k],y=j+dy[k];
if (!x||!y||x>n||y>m) continue;
add(calc(i,j),calc(x,y),inf);
}
dinic();
}