bzoj2241 打地鼠 暴力&&二维差分

       首先这道题目没有二分的性质。。。

       数据不是很强，O(N^4)暴力就能过呢（实际上由于常数很小就算是极限也能过的。。另外O(N^6)暴力都过了真是太强辣！！！虽然说代码并没有简单多少，似乎还更复杂了）。讲一下我的做法O(N^4)。

       首先暴力枚举长和边，然后就是要O(N^2)判断辣！怎么做呢？

       实际上，就是一个二维区间修改+单点查询。。然而你写个暴力总不至于写一个二维线段树吧。。。和一位一样，可以用差分转成前缀和，即(x1,y1,x2,y2)区间+t转化为(x1,y1)+t,(x1,y2+1)-t,(x2+1,y1)-t,(x2+1,y2+1)+t，然后点(x,y)的值转化为Σ(i=1,x)Σ(j=1,y)(i,j)的值。实际上前缀和求一下是很快的。。。

       然后暴力就没了。。正解好像是这样：显然如果(x,y)不行，那么任何(x*p,y*q)显然都不行。然后就可以用线性筛那种东西来搞了。。。（不过似乎并不会快很多？）

AC代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#define N 105
using namespace std;
 
int m,n,tot,a[N][N],s[N][N];
bool ok(int x,int y){
    if (tot%(x*y)) return 0; int i,j;
	memset(s,0,sizeof(s));
	for (i=1; i<=m; i++)
        for (j=1; j<=n; j++){
            s[i][j]+=s[i-1][j]+s[i][j-1]-s[i-1][j-1];
            if (s[i][j]>a[i][j]) return 0;
			int t=a[i][j]-s[i][j];
			if ((i>m-x+1 || j>n-y+1) && t) return 0;
			if (t){ s[i][j]+=t; s[i+x][j]-=t; s[i][j+y]-=t; s[i+x][j+y]+=t; }
        }
    return 1;
}
int main(){
    scanf("%d%d",&m,&n); int i,j;
    for (i=1; i<=m; i++)
        for (j=1; j<=n; j++){
            scanf("%d",&a[i][j]); tot+=a[i][j];
        }
    int ans=0;
    for (i=m; i; i--) for (j=n; j; j--)
        if (i*j>ans && ok(i,j)) ans=i*j;
    printf("%d\n",tot/ans);
    return 0;
}

by lych

2016.2.29