HDU 3046 Pleasant sheep and big big wolf 最小割

题意：

给定n*m个点的矩阵

0为空点。1为羊。2为狼

相邻点之间有一条路。

问要使得狼与羊不连通最少要去掉几条边

最小割

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<queue>
#include<vector>
using namespace std;

#define ll int 

#define N 40010
#define M 105000
#define inf 1073741824
#define inf64 1152921504606846976
struct Edge{  
	ll from, to, cap, nex, max;  
}edge[M*4];//注意这个一定要够大 不然会re 还有反向弧  

ll head[N], edgenum;  
void add(ll u, ll v, ll cap){  
	Edge E = { u, v, cap, head[u],cap};  
	edge[ edgenum ] = E;  
	head[u] = edgenum ++;  

	Edge E2= { v, u, 0,   head[v],cap};  
	edge[ edgenum ] = E2;  
	head[v] = edgenum ++;  
}  
ll sign[N];  
bool BFS(ll from, ll to){  
	memset(sign, -1, sizeof(sign));  
	sign[from] = 0;  

	queue<ll>q;  
	q.push(from);  
	while( !q.empty() ){  
		int u = q.front(); q.pop();  
		for(ll i = head[u]; i!=-1; i = edge[i].nex)  
		{  
			ll v = edge[i].to;  
			if(sign[v]==-1 && edge[i].cap)  
			{  
				sign[v] = sign[u] + 1, q.push(v);  
				if(sign[to] != -1)return true;  
			}  
		}  
	}  
	return false;  
}  
ll Stack[N], top, cur[N];  
ll dinic(ll from, ll to){  

	ll ans = 0;  
	while( BFS(from, to) )  
	{  
		memcpy(cur, head, sizeof(head));  
		ll u = from;      top = 0;  
		while(1)  
		{  
			if(u == to)  
			{  
				ll flow = inf, loc;//loc 表示 Stack 中 cap 最小的边  
				for(ll i = 0; i < top; i++)  
					if(flow > edge[ Stack[i] ].cap)  
					{  
						flow = edge[Stack[i]].cap;  
						loc = i;  
					}  

					for(ll i = 0; i < top; i++)  
					{  
						edge[ Stack[i] ].cap -= flow;  
						edge[Stack[i]^1].cap += flow;  
					}  
					ans += flow;  
					top = loc;  
					u = edge[Stack[top]].from;  
			}  
			for(ll i = cur[u]; i!=-1; cur[u] = i = edge[i].nex)//cur[u] 表示u所在能增广的边的下标  
				if(edge[i].cap && (sign[u] + 1 == sign[ edge[i].to ]))break;  
			if(cur[u] != -1)  
			{  
				Stack[top++] = cur[u];  
				u = edge[ cur[u] ].to;  
			}  
			else  
			{  
				if( top == 0 )break;  
				sign[u] = -1;  
				u = edge[ Stack[--top] ].from;  
			}  
		}  
	}  
	return ans;  
}  

void init(){memset(head,-1,sizeof head);edgenum = 0;}

int mp[205][205];
int n, m;
int Hash(int x,int y){return (x-1)*m+y;}
bool inmap(int x,int y){return 1<=x&&x<=n&&1<=y&&y<=m;}
int w[N], s[N], step[4][2] = {1,0,-1,0,0,1,0,-1};
int main()
{
    int u, v, d, i, j, Cas = 1;
	while(~scanf("%d %d",&n,&m)){
		init();
		int wt = 0, st = 0;
		for(i = 1; i <= n; i++)
		{
			for(j = 1; j <= m; j++)
			{
				scanf("%d",&mp[i][j]);
				int now = Hash(i,j);
				if(mp[i][j]==1)s[st++] = now;
				if(mp[i][j]==2)w[wt++] = now;
				for(int k = 0; k < 4; k++){
					if(!inmap(i+step[k][0], j+step[k][1]))continue;
					int go = Hash(i+step[k][0], j+step[k][1]);
					add(now, go, 1);
				}
			}
		}
		int from = 0, to = Hash(n,m)+1;
		for(i = 0; i < wt; i++)add(from, w[i], inf);
		for(i = 0; i < st; i++)add(s[i], to, inf);
		printf("Case %d:\n%d\n",Cas++,dinic(from, to));
	}
    return 0;
}
/*
4 6
1 0 0 1 0 0
0 1 1 0 0 0
2 0 0 0 0 0
0 2 0 1 1 0

*/