[费用流] ACM-ICPC Asia China-Final Contest J. Mr.Panda and TubeMaster

%%%Manchery

#include <cstdio>#include <iostream>#include <algorithm>#include <cstring>#include <queue>using namespace std;const int N=10010,inf=1<<30;int t,n,m,cnt,S,T,SS,TT,ans,G[N],tag[N];struct edge{  int t,nx,w,f;}E[N<<2];#define x(i,j) (((i-1)*m+j)<<1)#define y(i,j) (((i-1)*m+j)<<1|1)inline void add(int x,int y,int f,int w){  E[++cnt].t=y; E[cnt].nx=G[x]; G[x]=cnt; E[cnt].w=w; E[cnt].f=f;}inline void addedge(int x,int y,int f,int w,int t){  add(x,y,f,w); tag[cnt]=t; add(y,x,0,-w); tag[cnt]=t;}inline void link(int x,int y,int f,int w){  if(w>0){    ans+=f*w; addedge(y,x,f,-w,1); addedge(SS,y,f,0,0); addedge(x,TT,f,0,0);  }  else addedge(x,y,f,w,0);}int vis[35][35];int dis[N],vst[N],frm[N];queue<int> Q;inline bool spfa(){  for(int i=1;i<=TT;i++) dis[i]=-inf,vst[i]=0;  vst[SS]=1; dis[SS]=0; Q.push(SS);  while(!Q.empty()){    int x=Q.front(); vst[x]=0; Q.pop();    for(int i=G[x];i;i=E[i].nx)      if(E[i].f && dis[E[i].t]<dis[x]+E[i].w){    dis[E[i].t]=dis[x]+E[i].w; frm[E[i].t]=i;    if(!vst[E[i].t])      vst[E[i].t]=1,Q.push(E[i].t);      }  }  if(dis[TT]==-inf) return false;  int mf=inf;  for(int i=TT;i!=SS;i=E[frm[i]^1].t) mf=min(mf,E[frm[i]].f);  ans+=dis[TT]*mf;  for(int i=TT;i!=SS;i=E[frm[i]^1].t) E[frm[i]].f-=mf,E[frm[i]^1].f+=mf;  return true;}int main(){  scanf("%d",&t); int Case=0;  while(t--){    memset(G,0,sizeof(G)); cnt=1; ans=0;    scanf("%d%d",&n,&m); S=2*y(n,m)+1; T=2*y(n,m)+2; SS=T+1; TT=SS+1;    for(int i=1;i<=n;i++)      for(int j=1;j<m;j++){    int c; scanf("%d",&c);    if(~(i+j)&1)      link(x(i,j),x(i,j+1),1,c);    else      link(x(i,j+1),x(i,j),1,c);      }    for(int i=1;i<n;i++)      for(int j=1;j<=m;j++){    int c; scanf("%d",&c);    if((i+j)&1)      link(y(i,j),y(i+1,j),1,c);    else      link(y(i+1,j),y(i,j),1,c);      }    for(int i=1;i<=n;i++)for(int j=1;j<=m;j++) vis[i][j]=0;    int e; scanf("%d",&e);    while(e--){      int x,y; scanf("%d%d",&x,&y); vis[x][y]=1;    }    int ct=0;    for(int i=1;i<=n;i++)      for(int j=1;j<=m;j++){    if((i+j)&1){      if(!vis[i][j])        link(x(i,j),y(i,j),1,0);      else        ct++,link(S,y(i,j),1,0),link(x(i,j),T,1,0);    }    else{      if(!vis[i][j])        link(y(i,j),x(i,j),1,0);      else        ct++,link(S,x(i,j),1,0),link(y(i,j),T,1,0);    }      }    while(spfa());    memset(G,0,sizeof(G)); int lst=cnt;    for(int i=2;i<=lst;i++){      if(E[i^1].t==S)    int sss=0;      if(E[i].t==SS || E[i].t==TT || E[i^1].t==SS || E[i^1].t==TT) continue;      if(tag[i])    add(E[i].t,E[i^1].t,E[i^1].f,-E[i].w);      else    add(E[i^1].t,E[i].t,E[i].f,E[i].w);    }    SS=S; TT=T;    while(spfa()) ct--;    if(ct) printf("Case #%d: Impossible\n",++Case);    else printf("Case #%d: %d\n",++Case,ans);  }  return 0;}