洛谷.5300.[GXOI/GZOI2019]与或和(单调栈)

LOJ
BZOJ
洛谷

想了一个奇葩的单调栈，算的时候要在中间取\(\min\)，感觉不靠谱不写了=-=
调了十分钟发现输出没取模=v=

BZOJ好逗逼啊 题面连pdf都不挂了 哈哈哈哈

---

枚举每一位。在这一位上与之后得到\(1\)的就是全\(1\)子矩形个数。或之后得到\(1\)的就是总举行个数减去全\(0\)子矩形个数。
单调栈算一下就好啦。
维护一个单调递增的栈。如果在右下角统计贡献，每次遇到往上延伸长度\(l\leq sk[top]\)时，会把\(sk[top]\)的一部分截断。用个变量\(sum\)维护前面所有柱子当前的延伸距离即可。
算惹我不说惹 看代码叭= =

---

//9208kb    7720ms
#include 
#include 
#include 
#define mod 1000000007
//#define gc() getchar()
#define MAXIN 500000
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
typedef long long LL;
const int N=1005;

int A[N][N];
char IN[MAXIN],*SS=IN,*TT=IN;

inline int read()
{
    int now=0;register char c=gc();
    for(;!isdigit(c);c=gc());
    for(;isdigit(c);now=now*10+c-48,c=gc());
    return now;
}
int Calc(const int n,const int bit,const int equ)
{
    static int sk[N],sz[N],mx[N][N];
    LL ans=0; sk[0]=-1;
    for(int i=1; i<=n; ++i)
    {
        for(int j=1; j<=n; ++j) mx[i][j]=((A[i][j]>>bit&1)==equ)?mx[i-1][j]+1:0;
        LL sum=0;
        for(int j=1,top=0; j<=n; ++j)
        {
            int tot=0;
            while(sk[top]>=mx[i][j]) sum-=sk[top]*sz[top], tot+=sz[top--];
            sk[++top]=mx[i][j], sz[top]=tot+1;
            sum+=sk[top]*sz[top], ans+=sum;
        }
    }
    return ans%mod;
}

int main()
{
    const int n=read(); int mx=0;
    for(int i=1; i<=n; ++i)
        for(int j=1; j<=n; ++j) mx=std::max(mx,A[i][j]=read());
    LL ans1=0,ans2=0,tmp=(1ll*n*(n+1)*n*(n+1)>>2)%mod;
    for(int i=0; 1ll<