胡策DAY2

施肥
【问题描述】
学校里有一块田地，所有庄稼排成了 N 行 M 列。恰好化学老师研制出了一种高效肥料——
金坷垃。初始时，每棵庄稼都有一个自己的高度 hi,j。每次可以使用 1mol 的金克拉使一棵庄稼的
高度增加 1。现在有 Q 个询问，化学老师每次想知道最少需要使用多少 mol 的金克拉，才能使田
里出现一块庄稼高度一致，大小为 ai*bi 的田地。
其中每个询问是独立的，也就是说这次询问使用的金克拉不会影响到下一次询问。
【输入格式】
输入第一行包含两个正整数 N 和 M。
接下来 N 行每行包括 M 个正整数，表示每棵庄稼最初的高度 hi,j。
接下来 Q 行每行包括两个正整数 ai 和 bi。
【输出格式】
输出共包括 Q 行，对于每个询问输出使用金克拉的最少 mol 数。
【输入样例】
3 4
1 8 3 4
5 2 3 1
3 6 2 2
4
1 1
2 2
2 3
3 2
【输出样例】
0
4
15
9
【输入输出样例说明】
对于第四个询问，最优方案之一为选择以下庄稼：
3 4
3 1
2 2
【数据说明】
对于 30%的数据 1≤n,m≤30
对于 50%的数据 1≤n,m≤300
对于 100%的数据 1≤n,m,hi,j≤1000 1≤ai≤n 1≤bi≤m 1≤Q≤30
题解：这个题吧怎么做很好想，可是怎么降低时间复杂度是个问题。找一遍每个符合啊ai*bi的田地求最少需要多少化肥。那么需要先知道这个块田地的最大值是什么，所有的和是什么。数据范围1000，需要n2的算法，所以最大值得O(n)求，可以维护一个二维前缀和。O（n）怎么求呢。二维？No No No 不用，可以先处理行在处理列，两边一维解决，接下来就是怎么O（n）求。我们可以维护一个单调递减队列，正确性（lc说显然）。然后问题就解决了O(∩_∩)O~~。
（O（n）求最值可以去做poj滑动窗口 or codves窗口）

#include <iostream>
#include <cstdio>
#include <algorithm>
using namespace std;
void read(int &a)
{
    a=0;
    char c=getchar();
    while(c<'0'||c>'9')
    c=getchar();
    while(c>='0'&&c<='9')
    {
        a*=10;
        a+=c-'0';
        c=getchar();
    }
}
const int lc=1006;
int num[lc][lc],S[lc][lc],sum[lc][lc],mx[lc][lc];
int a,b,n,m,H,T,Q;
int q[lc];
void qzh()
{
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
        {
            read(num[i][j]);
            S[i][j]=S[i-1][j]+S[i][j-1]-S[i-1][j-1]+num[i][j];
        }
    return ;
}
void xw()
{
    while(Q--)
    {
        read(a);read(b);
        if(a==1&&b==1)
        {
            puts("0");
            continue;
        }
        int ans=1e9+7;
        for(int i=a;i<=n;i++)
            for(int j=b;j<=m;j++)
            sum[i][j]=S[i][j]-S[i-a][j]-S[i][j-b]+S[i-a][j-b];
        for(int i=1;i<=m;i++)
        {
            H=1,T=0;
            for(int j=1;j<=n;j++)
            {
                while (H<=T&&num[j][i]>num[q[T]][i]) T--;
                q[++T]=j;
                if(j>a)
                while (H<=T&&q[H]<=j-a) H++;
                mx[j][i]=num[q[H]][i];
            }
        }
        for(int i=a;i<=n;i++)
        {
            H=1,T=0;
            int Max=0;
            for(int j=1;j<=m;j++)
            {
                while (H<=T&&mx[i][j]>mx[i][q[T]]) T--;
                q[++T]=j;
                if(j>b)
                while (H<=T&&q[H]<=j-b) H++;
                Max=mx[i][q[H]];
                if(j>=b)
                ans=min(ans,Max*a*b-sum[i][j]);
            }
        }
        printf("%d\n",ans);
    }
    return;
}
void work()
{
    read(n);read(m);
    qzh();
    read(Q);
    xw();
    return ;
}
int main()
{
    //freopen("fertile.in","r",stdin);
    //freopen("fertile.out","w",stdout);
    work();
    return 0;
}

下面是30暴力。

#include <iostream>
#include <cstdio>
using namespace std;
int mp[1050][1050];
int n,m;
int q,a,b;
int zda(int x,int y)
{
    int maxx=-1;
    for(int i=x;i<=x+a-1;i++)
        for(int j=y;j<=y+b-1;j++)
        maxx=max(maxx,mp[i][j]);
    int ans=0;
    for(int i=x;i<=x+a-1;i++)
        for(int j=y;j<=y+b-1;j++)
        ans+=maxx-mp[i][j];
    return ans;
}
int main()
{
    //freopen("fertile.in","r",stdin);
    //freopen("fertile.out","w",stdout);
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)   
        scanf("%d",&mp[i][j]);
    scanf("%d",&q);
    while (q!=0)
    {
        q--;
        int dan=1e9+7;
        scanf("%d%d",&a,&b);
        if(a==1&&b==1)
        printf("0\n");
        else 
        {
            for(int i=1;i<=n;i++)
                for(int j=1;j<=m;j++)
                {
                    if(i+a-1<=n&&j+b-1<=m)
                    dan=min(dan,zda(i,j));
                }
            printf("%d\n",dan);
        }
    }
    //fclose(stdin);
    //fclose(stdout);
    return 0;