【BZOJ】【2756】【SCOI2012】奇怪的游戏

　　这题……建模部分先略过

　　这道题是会卡时限的T_T俺的Dinic被卡了，在此放几篇很棒的讲网络流算法的文章，至于大家耳熟能详的论文就不放了……

　　http://www.cppblog.com/panzhizhou/articles/172978.html?opt=admin

　　里面的各种超链接也很不错的……

　　好的来重新更新一下……这题因为要二分，需要多次重建跑最大流，所以不能用像lrj大爷的白书上那样用vector存边（太慢），需用前向星= =

　　然后……本蒻由于第一次写前向星，且印象中好像不加【当前弧优化】效率也不会低太多……所以顺利TLE了。事实上当前弧的设计还是和vector存边时差不多的，而且很有必要优化这一下……

round #1 2b版：

  1 /**************************************************************
  2     Problem: 2756
  3     User: ProgrammingApe
  4     Language: C++
  5     Result: Accepted
  6     Time:8000 ms
  7     Memory:2008 kb
  8 ****************************************************************/
  9  
 10 //BZOJ 2756
 11 #include<queue>
 12 #include<cstdio>
 13 #include<vector>
 14 #include<cstring>
 15 #include<cstdlib>
 16 #include<iostream>
 17 #include<algorithm>
 18 #define rep(i,n) for(int i=0;i<n;++i)
 19 #define F(i,j,n) for(int i=j;i<=n;++i)
 20 #define D(i,j,n) for(int i=j;i>=n;--i)
 21 #define pb push_back
 22 using namespace std;
 23 const int N=50,INF=~0u>>2;
 24 const long long infll=~0uLL>>2;
 25 const int fx[]={1,0,-1,0},
 26           fy[]={0,1,0,-1};
 27 typedef long long LL;
 28 //#define debug
 29  
 30 void read(int &v){
 31     v=0; int sig=1;
 32     char ch=getchar();
 33     while(ch<'0'||ch>'9'){ if (ch=='-') sig=-1; ch=getchar();}
 34     while(ch>='0'&&ch<='9'){ v=v*10+ch-'0'; ch=getchar();}
 35     v*=sig;
 36 }
 37  
 38 int n,m,a[N][N],d[N*N],cur[N*N],s,t,cnt;
 39 bool color[N][N];
 40 LL sum[2];
 41 struct edge{
 42     int from,to;
 43     LL cap,flow;
 44     int next;
 45 }E[N*N*10];
 46 int head[2000];
 47 //vector<edge>E;
 48 //vector<int>G[2000];
 49  
 50 void add(int from,int to,LL cap){
 51     E[++cnt]=(edge){from,to,cap,0,head[from]};
 52     head[from]=cnt;
 53     E[++cnt]=(edge){to,from,0,0,head[to]};
 54     head[to]=cnt;
 55 }
 56 //queue<int>Q;//·ÅÔÚmklevelÀïÃæ»áÂýÐí¶à
 57 int Q[N*N];
 58 bool mklevel(){
 59     memset(d,63,sizeof d);
 60     int l=0,r=0;
 61     d[s]=0; Q[r++]=s;
 62     while(l<r){
 63         int x=Q[l++];
 64         for(int i=head[x];i;i=E[i].next){
 65             edge&e=E[i];
 66             if (d[e.to]>200000 && e.cap>e.flow){
 67                 d[e.to]=d[x]+1;
 68                 Q[r++]=e.to;
 69             }
 70         }
 71     }
 72     return d[t]<200000;
 73 }
 74 LL dfs(int x,LL a){
 75     if (x==t||a==0) return a;
 76     LL flow=0;
 77     for(int &i=cur[x];i;i=E[i].next){
 78         edge&e=E[i];
 79         if (d[e.to]!=d[x]+1) continue;
 80         LL f=dfs(e.to,min(a,e.cap-e.flow));
 81         if (f>0){
 82             flow+=f;
 83             e.flow+=f;
 84             E[((i-1)^1)+1].flow-=f;
 85             a-=f;
 86             if (a==0) break;
 87         }
 88     }
 89     return flow;
 90 }
 91 LL dinic(){
 92     LL flow=0;
 93     while(mklevel()){
 94         F(i,s,t) cur[i]=head[i];//µ±Ç°»¡ÓÅ»¯ºÜÖØÒªµÄT_T
 95         flow+=dfs(s,infll);//ÕâÀïÒ²Òª¸Ä³Éinfll
 96     }
 97     return flow;
 98 }
 99 LL Total=0;
100 int maxw=0;
101 void solve1(){
102     memset(E,0,sizeof E);
103     memset(head,0,sizeof head); cnt=0;
104     LL d=sum[0]-sum[1];
105     if (d<maxw) {printf("-1\n"); return;}
106     int x,y;
107     Total=0;
108     F(i,1,n)
109         F(j,1,m){
110 //          if (a[i][j]>d) {printf("-1\n"); return;}
111             if (color[i][j]) add((i-1)*m+j,t,d-a[i][j]);
112             else{
113                 add(s,(i-1)*m+j,d-a[i][j]);
114                 F(k,0,3){
115                     x=i+fx[k],y=j+fy[k];
116                     if (x<1||y<1||x>n||y>m) continue;
117                     add( (i-1)*m+j , (x-1)*m+y , infll);
118                 }
119                 Total+=d-a[i][j];
120             }
121         }
122     LL ans=dinic();
123     if (Total==ans) printf("%lld\n",ans);
124     else printf("-1\n");
125 }
126  
127 bool check(LL x){
128     memset(E,0,sizeof E);
129     memset(head,0,sizeof head);
130     cnt=0;//Çå¿Õ±ß¼¯Êý×éµÄÊ±ºò£¬¼ÇµÃ°Ñ±ß¼¯´óÐ¡cntÒ²Çå¿Õ
131     Total=0;
132     F(i,1,n)
133         F(j,1,m){
134             if (color[i][j]) add( (i-1)*m+j,t,x-a[i][j] );
135             else{
136                 add(s,(i-1)*m+j,x-a[i][j]);
137                 F(k,0,3){
138                     int tx=i+fx[k],ty=j+fy[k];
139                     if (tx<1||ty<1||tx>n||ty>m) continue;
140                     add( (i-1)*m+j , (tx-1)*m+ty , infll);
141                 }
142                 Total+=x-a[i][j];
143             }
144         }
145     LL flow=dinic();
146     return Total==flow;
147 }
148  
149 void solve(){//Èç¹ûÊÇÅ¼Êý£º¶þ·Ö
150     if (sum[0]!=sum[1]){
151         printf("-1\n");
152         return;
153     }
154     LL l=maxw,r=l*2,mid,ans=-1;//¾ÓÈ»ÊÇÕâÀïÍüÁË¸Ä¡
155     while(l<r){
156 //      cout <<l<<" "<<r<<endl;
157 //      printf("l=%lld r=%lld\n",l,r);
158         mid=l+r>>1;
159         if (check(mid)) {ans=Total; r=mid;}
160         else l=mid+1;
161     }
162     if (ans!=-1) printf("%lld\n",ans);
163     else printf("-1\n");
164 }
165  
166 int main(){
167     #ifndef ONLINE_JUDGE
168     freopen("input.txt","r",stdin);
169 //  freopen("output.txt","w",stdout);
170     #endif
171     int T;
172     scanf("%d",&T);
173     while(T--){
174         read(n),read(m);
175         s=0; t=n*m+1;
176         memset(a,0,sizeof a);
177         memset(color,0,sizeof color);
178         sum[0]=sum[1]=maxw=0;
179         F(i,1,n) F(j,1,m) {
180             read(a[i][j]);
181             color[i][j]=(i+j)&1;
182             sum[color[i][j]]+=a[i][j];
183             maxw=max(maxw,a[i][j]);
184         }
185         if ((n&1) && (m&1)) solve1();
186         else solve();
187     }
188     return 0;
189 }

View Code