标签:print 矩阵 文件中 put sample 纪念 最小 sizeof log
题解:考前复习最小割模板题
从S向所有1的点连一条∞的边,从所有2的点向T连一条∞的边,从1向相邻的0和2连一条1的边,从0向相邻的所有点连一条1的边
bz400题 1A纪念~
#include <cstdio>
#include <cstring>
#include <iostream>
#include <queue>
#define P(A,B) ((A-1)*m+B)
using namespace std;
int n,m,cnt,ans,S,T;
int dx[]={0,1,0,-1},dy[]={1,0,-1,0};
int head[10010],next[200010],val[200010],to[200010],d[10010];
int map[110][110];
queue<int> q;
void add(int a,int b,int c)
{
to[cnt]=b,val[cnt]=c,next[cnt]=head[a],head[a]=cnt++;
}
int bfs()
{
int i,u;
while(!q.empty()) q.pop();
memset(d,0,sizeof(d));
d[S]=1,q.push(S);
while(!q.empty())
{
u=q.front(),q.pop();
for(i=head[u];i!=-1;i=next[i])
{
if(!d[to[i]]&&val[i])
{
d[to[i]]=d[u]+1;
if(to[i]==T) return 1;
q.push(to[i]);
}
}
}
return 0;
}
int dfs(int x,int mf)
{
if(x==T) return mf;
int i,temp=mf,k;
for(i=head[x];i!=-1;i=next[i])
{
if(d[to[i]]==d[x]+1&&val[i])
{
k=dfs(to[i],min(temp,val[i]));
if(!k) d[to[i]]=0;
val[i]-=k,val[i^1]+=k,temp-=k;
if(!temp) break;
}
}
return mf-temp;
}
int main()
{
scanf("%d%d",&n,&m);
int i,j,k,x,y,a,b,c;
memset(head,-1,sizeof(head));
for(i=1;i<=n;i++) for(j=1;j<=m;j++) scanf("%d",&map[i][j]);
S=0,T=n*m+1;
for(i=1;i<=n;i++)
{
for(j=1;j<=m;j++)
{
if(map[i][j]==1) add(S,P(i,j),1<<30),add(P(i,j),S,0);
if(map[i][j]==2)
{
add(P(i,j),T,1<<30),add(T,P(i,j),0);
continue;
}
for(k=0;k<4;k++)
{
x=i+dx[k],y=j+dy[k];
if(x<=n&&y<=m&&x&&y&&map[x][y]!=1) add(P(i,j),P(x,y),1),add(P(x,y),P(i,j),0);
}
}
}
while(bfs()) ans+=dfs(S,1<<30);
printf("%d",ans);
return 0;
}
【BZOJ1412】[ZJOI2009]狼和羊的故事 最小割
标签:print 矩阵 文件中 put sample 纪念 最小 sizeof log
原文地址:http://www.cnblogs.com/CQzhangyu/p/6729582.html