首先看到与或,很显然想到按照位拆分运算。然后就变成了0/1矩阵,要使矩阵在当前位与为1,则矩阵全为1,如果是或为1,则是矩阵不全为0,然后求全为0/1的矩阵个数即可。记录c[i][j]表示以a[i][j]在该位向上0/1的长度。然后对于每一行,单调栈求解即可。
#include<bits/stdc++.h> using namespace std; const int N=1003,mod=1e9+7; int n,ans1,ans2,top,a[N][N],b[N][N],c[N][N],st[N],sum[N]; int calc() { for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) c[i][j]=b[i][j]?c[i-1][j]+1:0; int ret=0; for(int i=1;i<=n;i++) { st[0]=top=0; for(int j=1;j<=n;j++) if(!c[i][j])st[0]=j,top=0; else{ while(top&&c[i][j]<=c[i][st[top]])top--; st[++top]=j,sum[top]=(sum[top-1]+1ll*(j-st[top-1])*c[i][j])%mod; ret=(ret+sum[top])%mod; } } return ret; } int main() { scanf("%d",&n); for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) scanf("%d",&a[i][j]); for(int t=0;t<=30;t++) { for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) b[i][j]=(a[i][j]>>t)&1; ans1=(ans1+(1ll<<t)*calc())%mod; for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) b[i][j]^=1; for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) ans2=(ans2+1ll*(n-i+1)*(n-j+1)%mod*(1ll<<t))%mod; ans2=(ans2-(1ll<<t)*calc()%mod+mod)%mod; } printf("%d %d",ans1,ans2); }