BZOJ 4657 tower(最小割棋盘覆盖模型)

题目

考场上打了一发上下界最小费用流，55pts尚可。
这个题是最小割的模型。

因为炮塔不能瞄准炮塔，所以两个炮塔的攻击区域只会垂直相交。
将炮塔横向覆盖的格子作为X部，纵向覆盖的格子作为Y部。
横向覆盖的格子，相邻之间连无穷大的边，方向指向炮塔。
纵向覆盖的格子，相邻之间连无穷大的边，方向指向炮塔的反方向。
源点向横向覆盖的格子连边，容量为该格子的前缀最大值与它前一个格子的前缀最大值的增量。前是指向炮塔的方向。
纵向覆盖的格子向汇点连边，容量为该格子的前缀最大值与它前一个格子的前缀最大值的增量。前仍然是指向炮塔的方向。
跑一遍最大流。
答案是每个炮塔能消灭的beta的最大值之和-最小割。

我们这里的最小割割的是最大值的增量。
然后就是有一些增量不能同时选，同时我们对于每行每列应该选前缀增量和，也就是一个选了另一个也要选。
就这样连边求前缀增量和 - 最小割就行了。
只需要分两层点。

AC Code:

#include<bits/stdc++.h>
#define maxn 5505
#define inf 100000000
#define maxm 100005
using namespace std;

int n,m;
int buf[maxn],info[maxn],Prev[maxm],to[maxm],cap[maxm],cnt_e=1;
void Node(int u,int v,int c){ Prev[++cnt_e]=info[u],info[u]=cnt_e,to[cnt_e]=v,cap[cnt_e]=c; }
void Line(int u,int v,int c){ Node(u,v,c),Node(v,u,0); }
int mp[100][100];
int S,T,h[maxn];

bool BFS(){
	static int q[maxn],L,R;
	memset(h,-1,sizeof h);
	L = R = 0 , q[R++] = T;
	h[T] = 0;
	for(int now;L<R;){
		now = q[L++];
		for(int i=info[now];i;i=Prev[i])
			if(cap[i^1] && h[to[i]] == -1)
				h[q[R++]=to[i]]=h[now]+1;
	}
	return h[S] != h[0];
}

int aug(int now,int Max){
	if(now == T) return Max;
	int inc , st = Max;
	for(int &i=info[now];i;i=Prev[i])
		if(cap[i] && h[to[i]]+1 == h[now]){
			inc = aug(to[i],min(cap[i],st));
			st-=inc,cap[i]-=inc,cap[i^1]+=inc;
			if(!st) break;
		}
	return Max - st;
}

int id1(int x,int y){ return (x-1)*m+y; }
int id2(int x,int y){ return (x-1)*m+y+n*m; }

int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			scanf("%d",&mp[i][j]);
	S = n * m * 2 + 1, T = S+1;
	int ans = 0;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			if(mp[i][j] < 0)
			{
				if(mp[i][j] == -1){
					int pr=0,mx=0;
					for(int x=i;x>=1;x--){
						mx = max(mx , mp[x][j]);
						if(mx-pr) Line(S,id1(x,j),mx-pr);
						if(x<i) Line(id1(x,j),id1(x+1,j),inf);
						pr = mx;
					}
					ans += mx;
				}
				if(mp[i][j] == -2){
					int pr=0,mx=0;
					for(int x=i;x<=n;x++){
						mx = max(mx , mp[x][j]);
						if(mx-pr) Line(S,id1(x,j),mx-pr);
						if(x>i) Line(id1(x,j),id1(x-1,j),inf);
						pr = mx;
					}
					ans += mx;
				}
				if(mp[i][j] == -3){
					int pr=0,mx=0;
					for(int y=j;y>=1;y--){
						mx = max(mx , mp[i][y]);
						if(mx-pr) Line(id2(i,y),T,mx-pr);
						if(y<j) Line(id2(i,y+1),id2(i,y),inf);
						pr = mx;
					}
					ans += mx;
				}
				if(mp[i][j] == -4){
					int pr=0,mx=0;
					for(int y=j;y<=m;y++){
						mx = max(mx,mp[i][y]);
						if(mx-pr) Line(id2(i,y),T,mx-pr);
						if(y>j) Line(id2(i,y-1),id2(i,y),inf);
						pr = mx;
					}
					ans += mx;
				}
			}
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			Line(id1(i,j),id2(i,j),inf);
	memcpy(buf,info,sizeof info);
	for(;BFS();)
		ans -= aug(S,inf),
		memcpy(info,buf,sizeof buf);
	printf("%d\n",ans);
}