[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\%\) 的数据，\(m\), \(n\leq 5\) ；
对于 \(60\%\) 的数据，\(m\), \(n\leq 30\) ；
对于 \(100\%\) 的数据，\(1\leq m\), \(n\leq 100\) ，其他数据不小于 \(0\)，不大于 \(10^5\) 。
题解
最暴力的想法：枚举法
先枚举锤子的长与宽（2层循环）；
再枚举锤子敲的位置（2层循环）；
下一步就该判断落锤是否合理合法，具体地，有两种思路。
思路1
假如有一个\(2\times2\)的锤子要砸击下面这个区域（展示的是全图的一部分）。
对于每一次砸击，我们将被砸区域的每个数减\(1\),直到区域中有任何一个数被减至\(0\)，（时间复杂度为：区间内元素的个数\(\times\)砸到有任何一个数变成\(0\)的次数）然后继续枚举下一个要砸的区域。下图中，黄色区域是正准备受砸的区域，红字是上一次砸完后的结果。
将每个能砸的区域都枚举完毕后，再判断一下全图中是否还有元素未被清零。若有，则说明这个型号的锤子无法满足要求，否则，计入最小答案。
思路2
假如有一个\(2\times2\)的锤子要砸击下面这个区域（展示的是全图的一部分）。
一次性砸多次，砸的次数是区域左上角的数大小。（即将左上角的数清零）下图展示按此方法砸前砸后的变化。
但是，如果出现下面这种，左上角的数清零导致区域内其它数变成负数的情况，就可以直接判断出这种型号的锤子无法满足要求，然后直接枚举下一种锤子的型号。
若有一种型号的锤子能锤遍全图，我们再判断一下全图中是否还有元素未被清零。若有，则说明这个型号的锤子无法满足要求，否则，计入最小答案。
显然，由于思路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;
}