【刷题】BZOJ 4950 [Wf2017]Mission Improbable

Description

那是春日里一个天气晴朗的好日子,你准备去见见你的老朋友Patrick,也是你之前的犯罪同伙。Patrick在编程竞赛上豪赌输掉了一大笔钱,所以他需要再干一票。为此他需要你的帮助,虽然你已经金盆洗手了。你刚开始很不情愿,因为你一点也不想再回到那条老路上了,但是你觉得听一下他的计划也无伤大雅。在附近的一个仓库里有一批货物,包含一些贵重的消费性部件,Patrick企图从中尽可能多地偷些东西出来。这意味着要找一条进去的路,弄晕安保人员,穿过各种各样的激光射线,你懂的,都是常见的抢劫技术。然而,仓库的核心装备了一套Patrick搞不定的安保系统。这也是他需要你帮助他的地方。这批货物被放置在一些巨大的立方体箱里,每个箱子的尺寸都是相同的。这些箱子堆放成许多整齐的堆,每个箱子可以表示成一个三维的网格。安保系统每个小时会用三台相机对这堆货物进行一次拍照,相机分别为:前置相机(front camera),侧置相机(side camera)和顶置相机(top camera)。前置相机的照片显示了每一行最高的那堆箱子的高度,侧置相机显示了每一列最高的那堆箱子的高度,顶置相机显示了每个位置是否存在一堆箱子。如果安保系统发现任何一张照片出现了变化,它会立即拉响警报。一旦 Patrick 进去了,他会确定每堆箱子的高度并且发给你。图1显示了一种网格可能的放置,以及每台相机会得到的视图。

图 1. 网格的高度值与对应的相机视图。

图 2. 洗劫后网格可能的高度值。

Patrick想尽可能多偷走一些箱子。由于他不能弄坏安保系统,他准备重新安排剩余每堆箱子的放置,使得下一次相机取像时会得到相同的照片,从而骗过安保系统。在上面的例子中,他可以偷走九个箱子。图2显示了一种可能的剩余箱子的安置方案能使得安保系统认为与原安置情况相同。Patrick想请你帮他确定在保证能骗过安保系统的情况下他最多能偷走多少个箱子。你会帮他干完这最后一票么?

Input

第一行包含两个整数r(1≤r≤100)和c(1≤n≤100),分别表示网格的行数与列数。

接下来r行,每行包含c个整数,表示对应行上每堆立方体箱的高度(箱子的数量)。

所有的高度在0到10^9之间 (含边界) 。

Output

输出在不被发现的情况下最多能偷走多少箱子。

Sample Input

样例1

5 5

1 4 0 5 2

2 1 2 0 1

0 2 3 4 4

0 3 0 3 1

1 2 2 1 1

样例2

2 3

50 20 3

20 10 3

Sample Output

样例1

9

样例2

30

Solution

贪心选择，每行每列除了最大值之外都可以变成1

然后因为可以移动，如果存在这样一种情况：

40 20 0

20 10 0

那么可以变成这样：

40 1 0

1 20 0

把两个20变成1个20

存在这种情况必定存在点，它所在的行和所在的列的最大值相同，就可以这样操作

但又存在这样的情况：

0 20

20 0

这样每个点都可以操作，但实际上操作一次才是最优的

所以用二分图匹配一下，这样每个点就只能选一次了

#include<bits/stdc++.h>
#define ui unsigned int
#define ll long long
#define db double
#define ld long double
#define ull unsigned long long
const int MAXN=100+10,inf=0x3f3f3f3f;
int e=1,r,c,line[MAXN],row[MAXN],beg[MAXN<<1],nex[MAXN*MAXN*2],to[MAXN*MAXN*2],cap[MAXN*MAXN*2],out[MAXN*MAXN*2],level[MAXN<<1],s,t,G[MAXN][MAXN],snt,cur[MAXN<<1],clk,vis[MAXN<<1];
ll ans;
std::queue<int> q;
template<typename T> inline void read(T &x)
{
	T data=0,w=1;
	char ch=0;
	while(ch!='-'&&(ch<'0'||ch>'9'))ch=getchar();
	if(ch=='-')w=-1,ch=getchar();
	while(ch>='0'&&ch<='9')data=((T)data<<3)+((T)data<<1)+(ch^'0'),ch=getchar();
	x=data*w;
}
template<typename T> inline void write(T x,char ch='\0')
{
	if(x<0)putchar('-'),x=-x;
	if(x>9)write(x/10);
	putchar(x%10+'0');
	if(ch!='\0')putchar(ch);
}
template<typename T> inline void chkmin(T &x,T y){x=(y<x?y:x);}
template<typename T> inline void chkmax(T &x,T y){x=(y>x?y:x);}
template<typename T> inline T min(T x,T y){return x<y?x:y;}
template<typename T> inline T max(T x,T y){return x>y?x:y;}
inline void insert(int x,int y,int z)
{
	to[++e]=y;
	nex[e]=beg[x];
	out[e]=x;
	beg[x]=e;
	cap[e]=z;
	to[++e]=x;
	nex[e]=beg[y];
	out[e]=y;
	beg[y]=e;
	cap[e]=0;
}
inline bool bfs()
{
	memset(level,0,sizeof(level));
	level[s]=1;
	q.push(s);
	while(!q.empty())
	{
		int x=q.front();
		q.pop();
		for(register int i=beg[x];i;i=nex[i])
			if(!level[to[i]]&&cap[i])
			{
				level[to[i]]=level[x]+1;
				q.push(to[i]);
			}
	}
	return level[t];
}
inline int dfs(int x,int maxflow)
{
	if(x==t||!maxflow)return maxflow;
	vis[x]=clk;
	int res=0,f;
	for(register int &i=cur[x];i;i=nex[i])
		if((vis[x]^vis[to[i]])&&cap[i]&&level[to[i]]==level[x]+1)
		{
			f=dfs(to[i],min(maxflow,cap[i]));
			res+=f;
			maxflow-=f;
			cap[i]-=f;
			cap[i^1]+=f;
			if(!maxflow)break;
		}
	vis[x]=0;
	return res;
}
inline int Dinic()
{
	int res=0;
	while(bfs())++clk,memcpy(cur,beg,sizeof(cur)),res+=dfs(s,inf);
	return res;
}
int main()
{
	read(r);read(c);
	for(register int i=1;i<=r;++i)
		for(register int j=1;j<=c;++j)
		{
			read(G[i][j]);
			chkmax(row[i],G[i][j]);
			chkmax(line[j],G[i][j]);
			ans+=max(0,G[i][j]-1);
		}
	for(register int i=1;i<=r;++i)ans-=max(row[i]-1,0);
	for(register int i=1;i<=c;++i)ans-=max(line[i]-1,0);
	for(register int i=1;i<=r;++i)
		for(register int j=1;j<=c;++j)
			if(row[i]==line[j]&&G[i][j])snt++,insert(i,j+r,1);
	s=r+c+1,t=s+1;
	for(register int i=1;i<=r;++i)insert(s,i,1);
	for(register int i=1;i<=c;++i)insert(i+r,t,1);
	Dinic();
	for(register int i=2;i<=(snt<<1);i+=2)
		if(!cap[i])ans+=max(0,row[out[i]]-1);
	write(ans,'\n');
	return 0;
}