标签: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