BZOJ 1412 & 最小割

什么时候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;
}