HDU5479 Colmerauer 单调栈+暴力优化

http://acm.hdu.edu.cn/showproblem.php?pid=5749

思路：

bestcoder 84

贡献：所有可能的子矩阵的面积和

//len1：子矩阵所有长的和
for(int i=1;i<=L;i++){
    for(int j=1;j<=R;j++){
        len1+=i+j-1;//-1是因为1*1单元格也是鞍点
    }
}

// #pragma comment(linker, "/STACK:102c000000,102c000000")
#include <iostream>
#include <cstdio>
#include <cstring>
#include <sstream>
#include <string>
#include <algorithm>
#include <list>
#include <map>
#include <vector>
#include <queue>
#include <stack>
#include <cmath>
#include <cstdlib>
#include <conio.h>
using namespace std;
#define clc(a,b) memset(a,b,sizeof(a))
#define inf 0x3f3f3f3f
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
const int N = 1010;
const int M = 1e6+10;
const long long MOD = 1LL<<32;
#define LL long long
#define LB long double
#define mi() (l+r)>>1
double const pi = acos(-1);
const double eps = 1e-8;
void fre() {
    freopen("in.txt","r",stdin);
}

int a[N][N];
int l[N][N],r[N][N],u[N][N],d[N][N];
stack<int>st;
LL ans;
int n,m;
void calc(){
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++){
            int L=j-l[i][j],R=r[i][j]-j,U=i-u[i][j],D=d[i][j]-i;
            LL tem=(LL)L*R*U*D*(L+R)*(U+D)/4;
            ans=(ans+(LL)a[i][j]*tem)%MOD;
        }
    printf("%I64d
",ans);
}

int main(){
    int T;
    scanf("%d",&T);
    while(T--){
        scanf("%d%d",&n,&m);
        for(int i=1;i<=n;i++){
            for(int j=1;j<=m;j++){
                scanf("%d",&a[i][j]);
            }
        }
        while(!st.empty()) st.pop();
        for(int i=1;i<=n;i++){
            a[i][0]=a[i][m+1]=-inf;
            while(!st.empty()) st.pop();st.push(0);
            for(int j=1;j<=m;j++){
                while(!st.empty()&&a[i][st.top()]>a[i][j]) st.pop();
                l[i][j]=st.top();
                st.push(j);
            }
            while(!st.empty()) st.pop();st.push(m+1);
            for(int j=m;j>=1;j--){
               while(!st.empty()&&a[i][st.top()]>a[i][j]) st.pop();
               r[i][j]=st.top();
               st.push(j);
            }
        }
        for(int j=1;j<=m;j++){
            a[0][j]=a[n+1][j]=inf;
            while(!st.empty()) st.pop();st.push(0);
            for(int i=1;i<=n;i++){
                while(!st.empty()&&a[st.top()][j]<a[i][j]) st.pop();
                u[i][j]=st.top();
                st.push(i);
            }
            while(!st.empty()) st.pop();st.push(n+1);
            for(int i=n;i>=1;i--){
                while(!st.empty()&&a[st.top()][j]<a[i][j]) st.pop();
                d[i][j]=st.top();
                st.push(i);
            }
        }
        ans=0;
        calc();
    }
    return 0;
}