[bzoj1001][BeiJing2006]狼抓兔子-题解[平面图最小割转最短路]/[Dinic求最小割]

  1 #include<iostream>
  2 #include<cstdio>
  3 #include<cstdlib>
  4 #include<cstring>
  5 #include<queue>
  6 #define N 1000
  7 #define M 1000
  8 #define mod N*M
  9 #define Max 2147483647
 10 #define ll long long
 11 #define l(A,B) ((A-1)*m+B)
 12 using namespace std;
 13 struct edge
 14 {
 15     int to;
 16     int next;
 17     int dis;
 18 }mp[((N-1)*M+(M-1)*N+(N-1)*(M-1))*2];
 19 int n,m,h,head[N*M+5],que[N*M+5],tip,hd,tl,ans,s,t;
 20 ll tim=1,user[N*M+5];
 21 int read()
 22 {
 23     char ch=getchar();int kin=1,gi=0;
 24     while(ch>‘9‘||ch<‘0‘){if(ch==‘-‘)kin=-1;ch=getchar();}
 25     while(ch>=‘0‘&&ch<=‘9‘){gi=gi*10+ch-48;ch=getchar();}
 26     return kin*gi;
 27 }
 28 void add(int u,int v,int dis)
 29 {
 30     mp[++h].dis=dis;
 31     mp[h].to=v;
 32     mp[h].next=head[u];
 33     head[u]=h;
 34 }
 35 void init()
 36 {
 37     int k;
 38     for(int i=1;i<=n;++i)
 39     for(int j=1;j<m;++j)
 40     {
 41         k=read();
 42         add(l(i,j),l(i,j+1),k);
 43         add(l(i,j+1),l(i,j),k);
 44     }
 45     for(int i=1;i<n;++i)
 46     for(int j=1;j<=m;++j)
 47     {
 48         k=read();
 49         add(l(i,j),l(i+1,j),k);
 50         add(l(i+1,j),l(i,j),k);
 51     }
 52     for(int i=1;i<n;++i)
 53     for(int j=1;j<m;++j)
 54     {
 55         k=read();
 56         add(l(i,j),l(i+1,j+1),k);
 57         add(l(i+1,j+1),l(i,j),k);
 58     }
 59 }
 60 int bfs()
 61 {
 62     int nw;
 63     tl=0;hd=0;
 64     tip=user[s]=tim;que[tl++]=s;tl%=mod;
 65     while(tl!=hd)
 66     {
 67         nw=que[hd++];hd%=mod;
 68         for(int i=head[nw];i;i=mp[i].next)
 69         {
 70             if(user[mp[i].to]<tim&&mp[i].dis>0)
 71             {
 72                 user[mp[i].to]=user[nw]+1;
 73                 if(user[mp[i].to]>tip)tip=user[mp[i].to];
 74                 que[tl++]=mp[i].to;tl%=mod;
 75                 if(mp[i].to==t)return 1;
 76             }
 77         }
 78     }
 79     return 0;
 80 }
 81 int dinic(int x,int y)
 82 {
 83     if(x==t)return y;
 84     int l=y,travel;
 85     for(int i=head[x];i;i=mp[i].next)
 86     {
 87         if(user[mp[i].to]==user[x]+1&&mp[i].dis>0)
 88         {
 89             travel=dinic(mp[i].to,min(mp[i].dis,y));
 90             y-=travel;
 91             mp[i].dis-=travel;
 92             mp[((i-1)^1)+1].dis+=travel;
 93             if(!y)break;
 94         }
 95     }
 96     if(y==l)user[x]=0;
 97     return l-y;
 98 }
 99 void work()
100 {
101     while(bfs())
102     {
103         tim=tip+1;
104         ans+=dinic(s,Max);
105     }
106 }
107 int main()
108 {
109     n=read();m=read();
110     s=1;t=l(n,m);
111     init();
112     work();
113     printf("%d\n",ans);
114     return 0;
115 }

  1 #include<iostream>
  2 #include<cstdio>
  3 #include<cstdlib>
  4 #include<cstring>
  5 #include<queue>
  6 #define Max 2147483640
  7 #define dis first
  8 #define pos second
  9 #define mk(A,B) make_pair(A,B)
 10 #define d(A,B,C) ((A-1)*m*2+m*(C-1)+B)
 11 using namespace std;
 12 typedef pair<int,int> nod;
 13 int read()
 14 {
 15     char ch=getchar();int cn=0,kin=1;
 16     while(ch>‘9‘||ch<‘0‘){if(ch==‘-‘)kin=-1;ch=getchar();}
 17     while(ch>=‘0‘&&ch<=‘9‘){cn=cn*10+ch-48;ch=getchar();}
 18     return cn*kin;
 19 }
 20 struct ni
 21 {
 22     vector<int>s[2];
 23     int ml;
 24 }a[2000005];
 25 int n,m,t;
 26 priority_queue<nod,vector<nod>,greater<nod> >que;
 27 void get(int,int);
 28 void init();
 29 int main()
 30 {
 31     freopen("a.in","r",stdin);
 32     n=read()-1;m=read()-1;
 33     if(n==0||m==0)
 34     {
 35 
 36         if(n<m)swap(n,m);
 37         int la,ans=Max;
 38         for(int i=1;i<=n;++i)
 39         {
 40             la=read();
 41             ans=min(ans,la);
 42         }
 43         cout<<ans<<endl;return 0;
 44     }
 45     t=n*2*m+1;
 46     init();que.push(mk(0,0));
 47     int nw,d;
 48     while(!que.empty())
 49     {
 50         nw=que.top().second;d=que.top().first;que.pop();
 51         for(int i=0;i<a[nw].s[0].size();++i)
 52         {
 53             if(d+a[nw].s[1][i]<a[a[nw].s[0][i]].ml)
 54             {
 55                 a[a[nw].s[0][i]].ml=d+a[nw].s[1][i];
 56                 que.push(mk(a[a[nw].s[0][i]].ml,a[nw].s[0][i]));
 57             }
 58         }
 59     }
 60     cout<<a[t].ml<<endl;
 61     fclose(stdin);
 62     return 0;
 63 }
 64 void get(int x,int y)
 65 {
 66     int t1;
 67     t1=read();
 68     a[x].s[0].push_back(y);
 69     a[x].s[1].push_back(t1);
 70     a[y].s[0].push_back(x);
 71     a[y].s[1].push_back(t1);
 72 }
 73 void init()
 74 {
 75     for(int i=1;i<=n*2*m;++i)
 76     {
 77         a[i].ml=Max;
 78     }a[t].ml=Max;
 79     for(int i=1;i<=n+1;++i)
 80     for(int j=1;j<=m;++j)
 81     {
 82         if(i==1)
 83         get(d(i,j,1),t);
 84         else if(i==n+1)
 85         get(0,d(i-1,j,2));
 86         else
 87         get(d(i,j,1),d(i-1,j,2));
 88     }
 89     for(int i=1;i<=n;++i)
 90     for(int j=1;j<=m+1;++j)
 91     {
 92         if(j==1)
 93         get(d(i,j,2),0);
 94         else if(j==m+1)
 95         get(d(i,j-1,1),t);
 96         else
 97         get(d(i,j,2),d(i,j-1,1));
 98     }
 99     for(int i=1;i<=n;++i)
100     for(int j=1;j<=m;++j)
101     {
102         get(d(i,j,1),d(i,j,2));
103     }
104 }