hdu 3046 Pleasant sheep and big big wolf（sap最大流）

开始在hrbust上做的，数据改为20，用矩阵和普通的最大流算法均可以过！后来查的时候看到了hdu上的，数据比那个大很多，点数最多可达40000，第一次MLE，换成邻接表，第二次TLE，果断换算法，sap轻松秒过！

建图的方法就是远点连2， 终点连1，容量inf， 把每个挨着的格子连起来，容量1，求最小割！ 这里求最小割很好理解，就像从一个2开始走，向1走，然后每次走都会经过一个格子的边界，不管有多少只wolf经过这个边界，只需要一个栅栏即可，也就是每个边界只有一个位置放栅栏，一般有三个或4个边界，这样，就能够建成一个从所有的2开始走向1的图，只要删除最小割上的点，就能阻止从1s和2s的连通！

下面是代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;


const int MAXN=40010;
const int MAXM=2000020;
const int INF = 0x3fffffff;
struct Node
{
	int to,next,cap;
}edge[MAXM];//注意是MAXM
int tol, n, m;
int head[MAXN];
int gap[MAXN],dis[MAXN],pre[MAXN],cur[MAXN];

void addedge(int u,int v,int w,int rw=0)
{
	edge[tol].to=v;edge[tol].cap=w;edge[tol].next=head[u];head[u]=tol++;
	edge[tol].to=u;edge[tol].cap=rw;edge[tol].next=head[v];head[v]=tol++;
}
int sap(int start,int end,int nodenum)
{
	memset(dis,0,sizeof(dis));
	memset(gap,0,sizeof(gap));
	memcpy(cur,head,sizeof(head));
	int u=pre[start]=start,maxflow=0,aug=-1;
	gap[0]=nodenum;
	while(dis[start]<nodenum)
	{
		loop:
		for(int &i=cur[u];i!=-1;i=edge[i].next)
		{
			int v=edge[i].to;
			if(edge[i].cap&&dis[u]==dis[v]+1)
			{
				if(aug==-1||aug>edge[i].cap)
					aug=edge[i].cap;
				pre[v]=u;
				u=v;
				if(v==end)
				{
					maxflow+=aug;
					for(u=pre[u];v!=start;v=u,u=pre[u])
					{
						edge[cur[u]].cap-=aug;
						edge[cur[u]^1].cap+=aug;
					}
					aug=-1;
				}
				goto loop;
			}
		}
		int mindis=nodenum;
		for(int i=head[u];i!=-1;i=edge[i].next)
		{
			int v=edge[i].to;
			if(edge[i].cap&&mindis>dis[v])
			{
				cur[u]=i;
				mindis=dis[v];
			}
		}
		if((--gap[dis[u]])==0)break;
		gap[dis[u]=mindis+1]++;
		u=pre[u];
	}
	return maxflow;
}

int o[205][205];
int main()
{
    int icase = 1;
    while ( scanf("%d%d", &n, &m) != EOF ) {
        int S = 0, T = n*m+1; tol = 0;
        memset( head, -1, sizeof(head));
        for ( int i = 1; i <= n; ++i ) for ( int j = 1; j <= m; ++j ) {
                scanf("%d", &o[i][j]);
                int u = (i-1)*m+j;
                if ( o[i][j] == 1 ) addedge( u, T, INF, 0 );
                else if ( o[i][j] == 2 ) addedge( S, u, INF, 0 );
                if ( i > 1 ) addedge( u, u-m, 1, 0 );
                if ( i < n ) addedge( u, u+m, 1, 0 );
                if ( j > 1 ) addedge( u, u-1, 1, 0 );
                if ( j < m ) addedge( u, u+1, 1, 0 );
        }
        printf("Case %d:\n%d\n", icase++, sap( S, T, T+1));
    }
}