[GXOI/GZOI2019]与或和(单调栈)

想了想决定把这几题也随便水个解题报告...

bzoj  luogu

思路：

首先肯定得拆成二进制30位啊

此后每一位的就是个01矩阵

Q1就是全是1的矩阵个数

Q2就是总矩阵个数减去全是0的矩阵个数

玉蟾宫警告

就是单调栈乱搞对吧

本题完结

事实上有了思路其他的都是多余的对吧所以就不要介意代码了

 #include<cstdio>
 #define mo 1000000007
 ;
 template<typename tp>inline void read(tp &tar)
 {
     tp ret=,f=;char ch=getchar();
     ;ch=getchar();}
     +ch-';ch=getchar();}
     tar=ret*f;
 }
 int n,ai[N][N],a1[N][N],a0[N][N];
 int ans1,ans0;
 int sta1[N],sta0[N];
 int main()
 {
     read(n);
     ;i<=n;i++)
     {
         ;j<=n;j++)
         {
             read(ai[i][j]);
         }
     }
     ,o=;k<;k++,o<<=)
     {
         ;i<=n;i++)
         {
             ;j<=n;j++)
             {
                 if(ai[i][j]&o)
                     a1[i][j]=a1[i-][j]+,a0[i][j]=;
                 else
                     a0[i][j]=a0[i-][j]+,a1[i][j]=;
             }
         }
         ;i<=n;i++)
         {
             ,hop0=,tmp1=,tmp0=;
             ;j<=n;j++)
             {
                 tmp1+=a1[i][j];
                 while(hop1&&a1[i][j]<a1[i][sta1[hop1]])
                 {
                     tmp1+=(sta1[hop1]-sta1[hop1-])*(a1[i][j]-a1[i][sta1[hop1]]);
                     hop1--;
                 }
                 ans1=(ans1+((1ll*tmp1)<<k)%mo)%mo;
                 sta1[++hop1]=j;
                 tmp0+=a0[i][j];
                 while(hop0&&a0[i][j]<a0[i][sta0[hop0]])
                 {
                     tmp0+=(sta0[hop0]-sta0[hop0-])*(a0[i][j]-a0[i][sta0[hop0]]);
                     hop0--;
                 }
                 ans0=(ans0+((1ll*(i*j-tmp0))<<k)%mo)%mo;
                 sta0[++hop0]=j;
             }
         }
     }
     printf("%d %d\n",ans1,ans0);
     ;
 }

我又被卡常了啊啊啊