bzoj1412-网络流最小割
Posted dancer16
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj1412-网络流最小割相关的知识,希望对你有一定的参考价值。
狼爱上羊啊爱的疯狂,谁让他们真爱了一场;狼爱上羊啊并不荒唐,他们说有爱就有方向......” Orez听到这首歌,心想:狼和羊如此和谐,为什么不尝试羊狼合养呢?说干就干! Orez的羊狼圈可以看作一个n*m个矩阵格子,这个矩阵的边缘已经装上了篱笆。可是Drake很快发现狼再怎么也是狼,它们总是对羊垂涎三尺,那首歌只不过是一个动人的传说而已。所以Orez决定在羊狼圈中再加入一些篱笆,还是要将羊狼分开来养。 通过仔细观察,Orez发现狼和羊都有属于自己领地,若狼和羊们不能呆在自己的领地,那它们就会变得非常暴躁,不利于他们的成长。 Orez想要添加篱笆的尽可能的短。当然这个篱笆首先得保证不能改变狼羊的所属领地,再就是篱笆必须修筑完整,也就是说必须修建在单位格子的边界上并且不能只修建一部分。
Input
文件的第一行包含两个整数n和m。接下来n行每行m个整数,1表示该格子属于狼的领地,2表示属于羊的领地,0表示该格子不是任何一只动物的领地。
Output
文件中仅包含一个整数ans,代表篱笆的最短长度。
Sample Input
2 2
2 2
1 1
2 2
1 1
Sample Output
2
数据范围
10%的数据 n,m≤3
30%的数据 n,m≤20
100%的数据 n,m≤100
数据范围
10%的数据 n,m≤3
30%的数据 n,m≤20
100%的数据 n,m≤100
HINT
Source
本题是经典的网络流最小割模型的应用
T的莫名其妙,明明是dinic模板
#include<iostream> #include<algorithm> #include<cstdio> #include<cstring> #include<cmath> #include<cstdlib> #include<vector> using namespace std; typedef long long ll; typedef long double ld; typedef pair<int,int> pr; const double pi=acos(-1); #define rep(i,a,n) for(int i=a;i<=n;i++) #define per(i,n,a) for(int i=n;i>=a;i--) #define Rep(i,u) for(int i=head[u];i;i=Next[i]) #define clr(a) memset(a,0,sizeof(a)) #define pb push_back #define mp make_pair #define fi first #define sc second #define pq priority_queue #define pqb priority_queue <int, vector<int>, less<int> > #define pqs priority_queue <int, vector<int>, greater<int> > #define vec vector #define T 10001 ld eps=1e-9; ll pp=1000000007; ll mo(ll a,ll pp){if(a>=0 && a<pp)return a;a%=pp;if(a<0)a+=pp;return a;} ll powmod(ll a,ll b,ll pp){ll ans=1;for(;b;b>>=1,a=mo(a*a,pp))if(b&1)ans=mo(ans*a,pp);return ans;} void fre() { freopen("c://test//input.in", "r", stdin); freopen("c://test//output.out", "w", stdout); } //void add(int x,int y,int z){ v[++e]=y; next[e]=head[x]; head[x]=e; cost[e]=z; } int dx[5]={0,-1,1,0,0},dy[5]={0,0,0,-1,1}; ll read(){ ll ans=0; char last=‘ ‘,ch=getchar(); while(ch<‘0‘ || ch>‘9‘)last=ch,ch=getchar(); while(ch>=‘0‘ && ch<=‘9‘)ans=ans*10+ch-‘0‘,ch=getchar(); if(last==‘-‘)ans=-ans; return ans; } #define M 500005 #define N 50005 const int INF=1<<30; int head[N],level[N],q[N],e,n,m; struct E_node{ int v,f,Next; }edge[M]; int Map[1005][1005]; void add(int x,int y,int z){ edge[++e].v=y; edge[e].f=z; edge[e].Next=head[x]; head[x]=e; edge[++e].v=x; edge[e].f=0; edge[e].Next=head[y]; head[y]=e; } int bfs(int s,int t){ memset(level,0,sizeof(level)); level[s]=1; int h=0,t_=1; q[h]=s; while (h<t_){ int x=q[h++]; if (x==t) return 1; for (int i=head[x];i;i=edge[i].Next){ int v=edge[i].v,f=edge[i].f; if (!level[v] && f>0){ level[v]=level[x]+1; q[t_++]=v; } } } return 0; } int dfs(int u,int maxf,int t){ if (u==t) return maxf; int ret=0; for (int i=head[u];i;i=edge[i].Next){ int v=edge[i].v,f=edge[i].f; if (level[v]==level[u]+1 && f>0){ int Min=min(maxf-ret,f); f=dfs(v,Min,t); edge[i].f-=f; edge[i^1].f+=f; ret+=f; if (ret==maxf) return ret; } } return ret; } int Dinic(int s,int t){ int ans=0; while (bfs(s,t)) ans+=dfs(s,INF,t); return ans; } void build() { for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) { if(Map[i][j]==1) add(0,(i-1)*m+j,INF); if(Map[i][j]==2) add((i-1)*m+j,T,INF); for(int k=1;k<=4;k++) { int nowi=i+dx[k],nowj=j+dy[k]; if(nowi<1||nowi>n||nowj>m||nowj<1||Map[i][j]==2) continue; if(Map[i][j]!=1||Map[nowi][nowj]!=1) { cout<<(i-1)*m+j<<‘ ‘<<(nowi-1)*m+nowj<<endl; add((i-1)*m+j,(nowi-1)*m+nowj,1); } } } } int main(){ n=read(),m=read(); for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) scanf("%d",&Map[i][j]); build(); printf("%d\n",Dinic(0,T)); return 0; }
#include<iostream> #include<cstdio> #include<cstring> #define inf 0x7fffffff #define T 10001 using namespace std; int head[10005],q[10005],h[10005]; int cnt=1,ans,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]; void ins(int u,int v,int w) {e[++cnt].to=v;e[cnt].next=head[u];e[cnt].v=w;head[u]=cnt;} void insert(int u,int v,int w) {ins(u,v,w);ins(v,u,0);} bool bfs() { int t=0,w=1,i,now; memset(h,-1,sizeof(h)); q[0]=0;h[0]=0; while(t<w) { now=q[t];t++;i=head[now]; while(i) { if(e[i].v&&h[e[i].to]==-1) { h[e[i].to]=h[now]+1; q[w++]=e[i].to; } i=e[i].next; } } return h[T]==-1? 0:1; } int dfs(int x,int f) { if(x==T)return f; int w,used=0,i; i=head[x]; while(i) { if(e[i].v&&h[e[i].to]==h[x]+1) { w=f-used; w=dfs(e[i].to,min(w,e[i].v)); e[i].v-=w; e[i^1].v+=w; used+=w; if(used==f)return f; } i=e[i].next; } if(!used)h[x]=-1; return used; } void dinic(){while(bfs())ans+=dfs(0,inf);} void ini() { scanf("%d%d",&n,&m); for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) scanf("%d",&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); if(mp[i][j]==2) insert((i-1)*m+j,T,inf); for(int k=0;k<4;k++) { int nowi=i+xx[k],nowj=j+yy[k]; if(nowi<1||nowi>n||nowj>m||nowj<1||mp[i][j]==2) continue; if(mp[i][j]!=1||mp[nowi][nowj]!=1) { cout<<(i-1)*m+j<<‘ ‘<<(nowi-1)*m+nowj<<endl; insert((i-1)*m+j,(nowi-1)*m+nowj,1); } } } } int main() { ini(); build(); dinic(); printf("%d",ans); return 0; }
希望哪位dalao看出错误指教。
以上是关于bzoj1412-网络流最小割的主要内容,如果未能解决你的问题,请参考以下文章
[bzoj2229][Zjoi2011]最小割_网络流_最小割树