洛谷 P2216 [HAOI2007]理想正方形

洛谷

巨说这是一道单调队列好题，但是我并不是用单调队列做的诶。

如果往最暴力的方向去想，肯定是\(n^3\)的\(dp\)了。

\(f[i][j][k]\)代表当前正方形的左上角定点是\((i,j)\)，边长是\(k\)的正方形的最佳答案。

转移方程很简单，但是你一定会妥妥的\(\texttt{TLE}\)。

那么我们怎么做呢？

往倍增的方向去想，设\(f[i][j][k]\)表示左上角为\((i,j)\)，边长为\(2^j\)的正方形的最佳答案。

那么状态就这么转移：

\[mx[i][j]=max(mx[i][j],max(mx[i+(1<<k)][j+(1<<k)],max(mx[i+(1<<k)][j],mx[i][j+(1<<k)]))); \]

\[mn[i][j]=min(mn[i][j],min(mn[i+(1<<k)][j+(1<<k)],min(mn[i+(1<<k)][j],mn[i][j+(1<<k)]))); \]

很简单的转移方程，不过复杂度肯定比不上单调队列的做法了。

显然复杂度\(O(n^2~log(n))\)。

代码：

#include <bits/stdc++.h>
using namespace std;

const int N=1001;
int a,b,n,t[N][N];
int mx[N][N],mn[N][N],ans=2000000000,tmp;

int main()
{
    ios::sync_with_stdio(false);
    cin>>a>>b>>n;
    for (int i=0;i<a;++i)
        for (int j=0;j<b;++j)
            cin>>t[i][j],mn[i][j]=mx[i][j]=t[i][j];
    while ((1<<(tmp+1))<=n) ++tmp;
    for (int k=0;k<tmp;++k)
        for (int i=0;i+(1<<k)<a;++i)
            for (int j=0;j+(1<<k)<b;++j) {
                mx[i][j]=max(mx[i][j],max(mx[i+(1<<k)][j+(1<<k)],max(mx[i+(1<<k)][j],mx[i][j+(1<<k)])));
                mn[i][j]=min(mn[i][j],min(mn[i+(1<<k)][j+(1<<k)],min(mn[i+(1<<k)][j],mn[i][j+(1<<k)])));
            }
    for (int i=0;i<=a-n;++i)
        for (int j=0;j<=b-n;++j) {
            int MX=max(mx[i][j],max(mx[i+n-(1<<tmp)][j+n-(1<<tmp)],max(mx[i+n-(1<<tmp)][j],mx[i][j+n-(1<<tmp)])));
            int MN=min(mn[i][j],min(mn[i+n-(1<<tmp)][j+n-(1<<tmp)],min(mn[i+n-(1<<tmp)][j],mn[i][j+n-(1<<tmp)])));
            ans=min(MX-MN,ans);
        }
    cout<<ans;
    return 0;
}