[GXOI/GZOI2019]与或和（位运算，单调栈）

题目链接懒得放了。

题目大意懒得写了。

省选原题哪有找不到的……

---

说实话，其实这题是个大水题，被我十秒钟内口胡出来了。

首先位运算除了拆位还能干啥？以下以与为例，或是差不多的。

我们考虑有多少个子矩阵会对这一位答案产生贡献，其实就是全 $1$ 的子矩阵。

问题转化为计算全 $1$ 子矩阵的个数。

这是一个简单题。考虑枚举右下角，发现包括这个右下角的子矩阵肯定长这样：（画的比较丑，意会就好了）

也就是高度单调递增。

高度可以做到 $O(1)$ 转移（从 $h[i-1][j]$）转移。

至于递增的高度，直接一个单调栈。（设为 $s$）

那么这个点为右下角的矩阵个数为 $(s_1-s_0)h[i][s_1]+(s_2-s_1)h[i][s_2]+\cdots+(s_{top}-s_{top-1})h[i][s_{top}]$。这个也可以入出栈时随便更新一下。

时间复杂度 $O(n^2\log a_i)$。

（然而一开始式子推错了，调了好久，回来再看看发现自己就是个sb……）

#include<bits/stdc++.h>

using namespace std;

const int maxn=,mod=;

#define FOR(i,a,b) for(int i=(a);i<=(b);i++)

#define ROF(i,a,b) for(int i=(a);i>=(b);i--)

#define MEM(x,v) memset(x,v,sizeof(x))

inline int read(){

    char ch=getchar();int x=,f=;

    while(ch<'' || ch>'') f|=ch=='-',ch=getchar();

    while(ch>='' && ch<='') x=x*+ch-'',ch=getchar();

    return f?-x:x;

}

int n,a[maxn][maxn],b[maxn][maxn],ans1,ans2,h[maxn],stk[maxn],tp,sum;

int calc1(){

    int ans=;

    MEM(h,);

    FOR(i,,n){

        MEM(stk,);tp=sum=;

        FOR(j,,n) h[j]=b[i][j]?h[j]+:;

        FOR(j,,n){

            while(tp && h[j]<h[stk[tp]]){

                sum=(sum-1ll*h[stk[tp]]*(stk[tp]-stk[tp-])%mod+mod)%mod;

                tp--;

            }

            stk[++tp]=j;

            sum=(sum+1ll*h[stk[tp]]*(stk[tp]-stk[tp-]))%mod;

            ans=(ans+sum)%mod;

        }

    }

    return ans;

}

int calc2(){

    int ans=;

    MEM(h,);

    FOR(i,,n){

        MEM(stk,);tp=sum=;

        FOR(j,,n) h[j]=b[i][j]?:h[j]+;

        FOR(j,,n){

            while(tp && h[j]<h[stk[tp]]){

                sum=(sum-1ll*h[stk[tp]]*(stk[tp]-stk[tp-])%mod+mod)%mod;

                tp--;

            }

            stk[++tp]=j;

            sum=(sum+1ll*h[stk[tp]]*(stk[tp]-stk[tp-]))%mod;

            ans=(ans+sum)%mod;

        }

    }

    int tot=1ll*n*(n+)*n*(n+)/%mod;

    return (tot-ans+mod)%mod;

}

int main(){

    n=read();

    FOR(i,,n) FOR(j,,n) a[i][j]=read();

    FOR(_,,){

        FOR(i,,n) FOR(j,,n) b[i][j]=(a[i][j]>>_)&;

        ans1=(ans1+1ll*calc1()*(<<_))%mod;

        ans2=(ans2+1ll*calc2()*(<<_))%mod;

    }

    printf("%d %d\n",ans1,ans2);

}