什么时候ZJ省选再现一次这么良心的题吧...
题意:
在一个染色的格子画分割线,使其不想连,求最少的线段
SOL:
裸裸的最小割.题目要求两种颜色不想连,我们把他分到两个集合,也就是把所有相连的边切断-----这不就是最小割嘛. 把其中一个颜色与源相连,另一个颜色与汇相连,容量为正无穷,然后中间相连的容量均为1,然后跑下dinic即可.
Code:
/*==========================================================================
# Last modified: 2016-03-11 18:09
# Filename: 1412.cpp
# Description:
==========================================================================*/
#define me AcrossTheSky
#include <cstdio>
#include <cmath>
#include <ctime>
#include <string>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm> #include <set>
#include <map>
#include <stack>
#include <queue>
#include <vector> #define lowbit(x) (x)&(-x)
#define FOR(i,a,b) for((i)=(a);(i)<=(b);(i)++)
#define FORP(i,a,b) for(int i=(a);i<=(b);i++)
#define FORM(i,a,b) for(int i=(a);i>=(b);i--)
#define ls(a,b) (((a)+(b)) << 1)
#define rs(a,b) (((a)+(b)) >> 1)
#define getlc(a) ch[(a)][0]
#define getrc(a) ch[(a)][1] #define maxn 10005
#define maxm 100000
#define pi 3.1415926535898
#define _e 2.718281828459
#define inf 1070000000
using namespace std;
typedef long long ll;
typedef unsigned long long ull; template<class T> inline
void read(T& num) {
bool start=false,neg=false;
char c;
num=0;
while((c=getchar())!=EOF) {
if(c=='-') start=neg=true;
else if(c>='0' && c<='9') {
start=true;
num=num*10+c-'0';
} else if(start) break;
}
if(neg) num=-num;
}
/*==================split line==================*/
using namespace std;
int first[maxn],d[maxn],cur[maxn];
bool vis[maxn];
int cnt=1,n,m;
int xx[4]={0,0,1,-1},yy[4]={1,-1,0,0},mp[105][105];
struct data{int to,next,v;}e[500001];
int T,S;
void ins(int u,int v,int w)
{e[++cnt].to=v;e[cnt].next=first[u];e[cnt].v=w;first[u]=cnt;}
void insert(int u,int v,int w)
{ins(u,v,w);ins(v,u,0);} int bfs(){
queue<int> q;
for(int i=S;i<=T;i++) vis[i]=false;
q.push(0); d[0]=0; vis[0]=true;
while (!q.empty()){
int now=q.front(); q.pop();
for (int i=first[now];i;i=e[i].next)
if (!vis[e[i].to] && e[i].v){
d[e[i].to]=d[now]+1;
vis[e[i].to]=true;
q.push(e[i].to);
}
}
return vis[T];
}
int dfs(int now,int a){
if (now==T || !a) return a;
int f,flow=0;
for (int & i=cur[now];i;i=e[i].next)
if (d[now]+1==d[e[i].to] && (f=dfs(e[i].to,min(a,e[i].v)))>0){
flow+=f; a-=f; e[i].v-=f; e[i^1].v+=f;
if (!a) break;
}
return flow; }
int dinic(){
int ans=0;
while(bfs()){
FORP(i,0,T) cur[i]=first[i];
ans+=dfs(0,inf);
}
return ans;
}
void init()
{
read(n); read(m);
T=n*m+1,S=0;
FORP(i,1,n)
FORP(j,1,m) read(mp[i][j]);
}
void build()
{
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
{
if(mp[i][j]==1)insert(0,(i-1)*m+j,inf);
else if(mp[i][j]==2)insert((i-1)*m+j,T,inf);
for(int k=0;k<4;k++)
{
int nowx=i+xx[k],nowy=j+yy[k];
if(nowx<1||nowx>n||nowy<1||nowy>m||mp[i][j]==2)continue;
if(mp[i][j]!=1||mp[nowx][nowy]!=1)
insert((i-1)*m+j,(nowx-1)*m+nowy,1);
}
}
}
int main()
{
init();
build();
printf("%d",dinic());
return 0;
}