打地鼠（susliks） 方法记录

[SDOI2011]打地鼠

题目描述

2020.4.29 数据更新。

打地鼠是这样的一个游戏：地面上有一些地鼠洞，地鼠们会不时从洞里探出头来很短时间后又缩回洞中。玩家的目标是在地鼠伸出头时，用锤子砸其头部，砸到的地鼠越多分数也就越高。

游戏中的锤子每次只能打一只地鼠，如果多只地鼠同时探出头，玩家只能通过多次挥舞锤子的方式打掉所有的地鼠。你认为这锤子太没用了，所以你改装了锤子，增加了锤子与地面的接触面积，使其每次可以击打一片区域。如果我们把地面看做 $m\times n$ 的方阵，其每个元素都代表一个地鼠洞，那么锤子可以覆盖 $r\times c$ 区域内的所有地鼠洞。但是改装后的锤子有一个缺点：每次挥舞锤子时，对于这的区域中的所有地洞，锤子会打掉恰好一只地鼠。也就是说锤子覆盖的区域中，每个地洞必须至少有 $1$ 只地鼠，且如果某个地洞中地鼠的个数大于 $1$，那么这个地洞只会有 $1$ 只地鼠被打掉，因此每次挥舞锤子时，恰好有$r\times c$ 只地鼠被打掉。由于锤子的内部结构过于精密，因此在游戏过程中你不能旋转锤子（即不能互换 $r$ 和 $c$）。

你可以任意更改锤子的规格(即你可以任意规定 $r$ 和 $c$ 的大小)，但是改装锤子的工作只能在打地鼠前进行(即你不可以打掉一部分地鼠后，再改变锤子的规格)。你的任务是求出要想打掉所有的地鼠，至少需要挥舞锤子的次数。

Hint：由于你可以把锤子的大小设置为 $1\times 1$，因此本题总是有解的。

输入格式

第一行包含两个正整数 $m$ 和 $n$；

下面 $m$ 行每行 $n$ 个正整数描述地图，每个数字表示相应位置的地洞中地鼠的数量。

输出格式

输出一个整数，表示最少的挥舞次数。

样例 #1

样例输入 #1

3 3
1 2 1
2 4 2
1 2 1

样例输出 #1

4

提示

【样例说明】

使用 $2\times 2$ 的锤子，分别在左上、左下、右上、右下挥舞一次。

【数据规模和约定】

对于 $30\mathrm{\%}$ 的数据，$m$, $n\le 5$ ；

对于 $60\mathrm{\%}$ 的数据，$m$, $n\le 30$ ；

对于 $100\mathrm{\%}$ 的数据，$1\le m$, $n\le 100$ ，其他数据不小于 $0$，不大于 ${10}^5$ 。

题解

最暴力的想法：枚举法

先枚举锤子的长与宽（2层循环）；

再枚举锤子敲的位置（2层循环）；

下一步就该判断落锤是否合理合法，具体地，有两种思路。

思路1

假如有一个$2\times 2$的锤子要砸击下面这个区域（展示的是全图的一部分）。

对于每一次砸击，我们将被砸区域的每个数减$1$,直到区域中有任何一个数被减至$0$，（时间复杂度为：区间内元素的个数$\times$砸到有任何一个数变成$0$的次数）然后继续枚举下一个要砸的区域。下图中，黄色区域是正准备受砸的区域，红字是上一次砸完后的结果。

将每个能砸的区域都枚举完毕后，再判断一下全图中是否还有元素未被清零。若有，则说明这个型号的锤子无法满足要求，否则，计入最小答案。

思路2

假如有一个$2\times 2$的锤子要砸击下面这个区域（展示的是全图的一部分）。

一次性砸多次，砸的次数是区域左上角的数大小。（即将左上角的数清零）下图展示按此方法砸前砸后的变化。

但是，如果出现下面这种，左上角的数清零导致区域内其它数变成负数的情况，就可以直接判断出这种型号的锤子无法满足要求，然后直接枚举下一种锤子的型号。

若有一种型号的锤子能锤遍全图，我们再判断一下全图中是否还有元素未被清零。若有，则说明这个型号的锤子无法满足要求，否则，计入最小答案。

显然，由于思路2需要枚举的砸击次数少，且能及时排除不合要求的锤子，所以思路2更优

然后再考虑其它优化手段。

假设一把 $r\times c$ 的锤子是合格的，它需要砸 $x$ 次，那么它总共砸毙的地鼠数 $sum=r\times c\times x$ ，这个$sum$ 我们可以事先统计出来，即全部地洞里的地鼠数量之和。

换个角度来想，我们已经统计出 $sum$ ，现在正在枚举 $i\times j$ 型号的锤子。若这种锤子满足要求，就应该满足：$sum$ 是$i\times j$ 的整数倍。在此基础上再判断该型号的锤子能否锤尽全图地鼠，若能，则答案 $ans=sum/i/j$ .可以节省不少时间。这里枚举 $i$ 和 $j$ ,就可以从$1$~$n$ 进行，这样 $i$ 和 $j$ 越枚举越大，$ans$ 也就越来越小，省去了比较大小的过程。

除此之外我们还可以进行剪枝。

若当前的 $sum/i/j$ 小于已经算出来的 $ans$ ，则当前情况必不可能为最优情况。

AC代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=105;
int a[N][N],b[N][N],n,m,ans,sum;
bool check(int x,int y)
{
	memcpy(b,a,sizeof(a));
	for(int i=1;i<=n-x+1;i++)
	{
		for(int j=1;j<=m-y+1;j++)
		{
			if(b[i][j])
			{
				int z=b[i][j];
				for(int k=0;k<x;k++)
				{
					for(int l=0;l<y;l++)
					{
						b[i+k][j+l]-=z;
						if(b[i+k][j+l]<0) return 0;
					}
				}
			}
		}
	}
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			if(b[i][j]) return 0;
		}
	}
	return 1;
}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			scanf("%d",&a[i][j]);
			sum+=a[i][j];
		}
	}
	ans=sum;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			if(sum%(i*j)==0&&sum/i/j<ans)
			{
				if(check(i,j)) ans=sum/i/j;
			}
		}
	}
	printf("%d\n",ans);
	return 0;
}