2016ECFinal J(无源汇上下界费用流)

题目链接：http://codeforces.com/group/Us3rfLfgWv/contest/101194

题意：给定4个插头，如题图示，每个格子可以填一个插头或者不填，但一些格子一定要填入插头(下称关键点)。问能不能成若干个完整的回路，如果可以，输出最大价值。回路穿过格子时会产生相应的价值。

想了好久。。

一个显然的性质是一个格子如果加插头，只能在纵向上和一个格子相连，在横向上和一个格子相连，那么可以考虑把每个格子拆成2个点i,i'，分别代表横向和纵向，然后去连边，i和i'之间连边，容量为1

对关键点来说，i和i'的流量必须为1，所以有了下界的限制

然后变成了在无向图上跑带上下界的最大循环费用流。。

然而无向并不能跑费用流，得变成跑有向的。。

模拟几遍之后(其实是看了tls的代码后)发现这个回路是具有高度的对称性的，然后可以发现对整个网格黑白染色后，每个回路都由配对的黑白格子构成，那么直接分成黑白情况来考虑就可以建图了。。

对格子i,i'，如果为黑格子，那么i连向i'，否则i'连向i，i和i'在其对应方向上的出入关系即可确定，然后用循环流的套路把关键点处理一下即可。。

然而如果这样直接跑费用流的话只能满足关键点的要求，网格中还剩下一下不经过关键点的回路也能产生贡献，(又看了tls的代码总之后)其实可以把每个点的i和i'之间的连边都定义一个下界，对非关键点来说，直接连一个反向边形成小循环流，这样可以代表不填插头，那么问题就解决了。。

 

 

 

 

/**
 *　　　　　　  　　┏┓　　 　┏┓
 * 　　　　　  　　┏┛┗━━━━━━━┛┗━━━┓
 * 　　　　　  　　┃　　　　　　　┃ 　
 * 　　　　　  　　┃　　　━　　 　┃
 * 　　　　　  　　┃　＞　　　＜　┃
 * 　　　　　  　　┃　　　　　　　┃
 * 　　　　　  　　┃...　⌒　... 　┃
 * 　　　　　　  　┃              ┃
 * 　　　　　　  　┗━┓          ┏━┛
 * 　　　　　　　　　┃          ┃　Code is far away from bug with the animal protecting　　　　　　　　　　
 * 　　　　　　　　　┃          ┃   神兽保佑,代码无bug
 * 　　　　　　　　　┃          ┃　　　　　　　　　　　
 * 　　　　　　　　　┃          ┃  　　　　　　
 * 　　　　　　　　　┃          ┃
 * 　　　　　　　　　┃          ┃　　　　　　　　　　　
 * 　　　　　　　　　┃          ┗━━━┓
 * 　　　　　　　　　┃              ┣┓
 * 　　　　　　　　　┃              ┏┛
 * 　　　　　　　　　┗┓┓┏━━━━━━━━┳┓┏┛
 * 　　　　　　　　　　┃┫┫       ┃┫┫
 * 　　　　　　　　　　┗┻┛       ┗┻┛
 */ 
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
#include<map>
#include<stack>
#include<cmath>
#include<set>
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,l,r) for(int i=l;i>=r;i--)
#define link(x) for(edge *j=h[x];j;j=j->next)
#define mem(a) memset(a,0,sizeof(a))
#define ll long long
#define eps 1e-8
#define succ(x) (1<<x)
#define lowbit(x) (x&(-x))
#define sqr(x) ((x)*(x))
#define mid (x+y>>1)
#define NM 2005
#define nm 20005
#define N 35
#define pi 3.1415926535897931
const int inf=1e9+7;
using namespace std;
ll read(){
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
    return f*x;
}
 
 



struct edge{int t,v,w;edge*next,*rev;}e[nm],*h[NM],*p[NM],*o=e;
void _add(int x,int y,int w,int v){o->t=y;o->v=v;o->w=w;o->next=h[x];h[x]=o++;}
void add(int x,int y,int w,int v){_add(x,y,w,v);_add(y,x,0,-v);h[x]->rev=h[y];h[y]->rev=h[x];}
int n,m,tot,ans,s,d[NM],b[N][N],w[NM],ca,_x,_y;
bool v[NM],a[N][N];
queue<int>q;

int spfa(){
    inc(i,0,n)d[i]=inf,w[i]=0,v[i]=false;
    d[0]=0;w[0]=inf;q.push(0);
    while(!q.empty()){
	int t=q.front();q.pop();v[t]=false;
	link(t)if(j->w&&d[j->t]>d[t]+j->v){
	    d[j->t]=d[t]+j->v;w[j->t]=min(w[t],j->w);p[j->t]=j;
	    if(!v[j->t])v[j->t]++,q.push(j->t);
	}
    }
    return w[n];
}


int main(){
    int _=read();while(_--){
	mem(h);o=e;mem(a);
	n=read();m=read();tot=ans=s=0;
	inc(i,1,n)inc(j,1,m)b[i][j]=++tot;
	inc(i,1,n)inc(j,1,m-1){
	    _x=-read();
	    if((i+j)%2)add(b[i][j],b[i][j+1],1,_x);else add(b[i][j+1],b[i][j],1,_x);
	}
	inc(i,1,n-1)inc(j,1,m){
	    _x=-read();
	    if((i+j)%2==0)add(b[i][j]+tot,b[i+1][j]+tot,1,_x);else add(b[i+1][j]+tot,b[i][j]+tot,1,_x);
	}
	int __=read();while(__--)_x=read(),_y=read(),a[_x][_y]++;
	inc(i,1,n)inc(j,1,m){
	    if((i+j)%2==0){
		add(0,b[i][j]+tot,1,0);
		add(b[i][j],tot<<1|1,1,0);
		if(!a[i][j])add(b[i][j]+tot,b[i][j],1,0);
	    }else{
		add(0,b[i][j],1,0);
		add(b[i][j]+tot,tot<<1|1,1,0);
		if(!a[i][j])add(b[i][j],b[i][j]+tot,1,0);
	    }
	}
	n=tot<<1|1;
	while(spfa()){
	    ans-=d[n]*w[n];s+=w[n];
	    for(int x=n;x;x=p[x]->rev->t)p[x]->w-=w[n],p[x]->rev->w+=w[n];
	}
	//printf("%d %d\n",s,ans);
	if(s==tot)printf("Case #%d: %d\n",++ca,ans);else printf("Case #%d: Impossible\n",++ca);
    }
    return 0;
}