bzoj4165: 矩阵（堆+hash）

　　求第k大用堆维护最值并出堆的时候扩展的经典题...

　　因为只有正数，所以一个矩阵的权值肯定比它的任意子矩阵的权值大，那么一开始把所有满足条件的最小矩阵加进堆里，弹出的时候上下左右扩展一行加进堆，用hash判重即可

#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<queue>
#define uint unsigned int
#define ull unsigned long long
using namespace std;
const int maxn=, mod=;
struct poi{int x1, yy1, x2, y2; uint sum;};
struct tjm{ull too; int pre;}e[];
priority_queue<poi>q;
bool operator<(poi a, poi b){return a.sum>b.sum;}
int n, m, mna, mnb, k, tot;
int last[mod];
uint mp[maxn][maxn];
void read(uint &k)
{
    int f=; k=; char c=getchar();
    while(c<'' || c>'') c=='-' && (f=-), c=getchar();
    while(c<='' && c>='') k=k*+c-'', c=getchar();
    k*=f;
}
inline void add(int x, ull y){e[++tot].too=y; e[tot].pre=last[x]; last[x]=tot;}
inline int makehs(ull x)
{
    int tmp=x%mod;
    for(int i=last[tmp];i;i=e[i].pre)
    if(e[i].too==x) return i;
    add(tmp, x);
    return tot;
}
inline int hs(ull x)
{
    int tmp=x%mod;
    for(int i=last[tmp];i;i=e[i].pre)
    if(e[i].too==x) return i;
    return -;
}
inline void update(int x1, int yy1, int x2, int y2)
{
    ull tmp=;
    tmp=tmp*+x1; tmp=tmp*+yy1; tmp=tmp*+x2; tmp=tmp*+y2;
    if(hs(tmp)!=-) return;
    makehs(tmp);
    q.push((poi){x1, yy1, x2, y2, mp[x2][y2]-mp[x2][yy1-]-mp[x1-][y2]+mp[x1-][yy1-]});
}
int main()
{
    scanf("%d%d%d%d%d", &n, &m, &mna, &mnb, &k);
    for(int i=;i<=n;i++)
        for(int j=;j<=m;j++)
            read(mp[i][j]), mp[i][j]+=mp[i-][j]+mp[i][j-]-mp[i-][j-];
    for(int i=mna;i<=n;i++)
        for(int j=mnb;j<=m;j++)
        {
            ull tmp=;
            int x1=i-mna+, yy1=j-mnb+, x2=i, y2=j;
            tmp=tmp*+x1; tmp=tmp*+yy1; tmp=tmp*+x2; tmp=tmp*+y2;
            makehs(tmp);
            q.push((poi){x1, yy1, x2, y2, mp[x2][y2]-mp[x2][yy1-]-mp[x1-][y2]+mp[x1-][yy1-]});
        }
    for(int i=;i<=k;i++)
    {
        if(q.empty()) return puts("-1"), ;
        poi t=q.top(); q.pop();
        int x1=t.x1, yy1=t.yy1, x2=t.x2, y2=t.y2;
        if(x1>) update(x1-, yy1, x2, y2);
        if(yy1>) update(x1, yy1-, x2, y2);
        if(x2<n) update(x1, yy1, x2+, y2);
        if(y2<m) update(x1, yy1, x2, y2+);
        if(i==k) printf("%u\n", t.sum);
    }
}