比赛 2026.4.11 评测结果 AAAAAAAAAA
题目名称 与或和 最终得分 100
用户昵称 RpUtl 运行时间 6.097 s
代码语言 C++ 内存使用 13.01 MiB
提交时间 2026-04-11 12:09:09
显示代码纯文本
#include <iostream>
#include <cstdio>
using namespace std;
const int N=1005;
const int mod=1e9+7;
int n,a[N][N],b[N][N];
int nxt[N][N],pre[N][N],f[N][N];
int top,stk[N];
long long ans;
long long calc1(int k){
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            b[i][j]=((a[i][j]>>k)&1);
        }
    }
    for(int i=1;i<=n;i++){
        int lst=0;
        for(int j=1;j<=n;j++){
            if(!b[i][j])lst=j;
            nxt[i][j]=j-lst;
        }
    }
    for(int j=1;j<=n;j++){
        for(int i=n;i>=1;i--){
            while(top&&nxt[stk[top]][j]>nxt[i][j]){
                pre[stk[top]][j]=i;
                top--;
            }
            stk[++top]=i;
        }
        while(top)pre[stk[top]][j]=0,top--;
    }
    long long res=0;
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            f[i][j]=f[pre[i][j]][j];
            f[i][j]+=(i-pre[i][j])*nxt[i][j];
            (res+=f[i][j])%=mod;
        }
    }
    return res;
}
long long calc2(int k){
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            b[i][j]=((a[i][j]>>k)&1);
            b[i][j]^=1;
        }
    }
    for(int i=1;i<=n;i++){
        int lst=0;
        for(int j=1;j<=n;j++){
            if(!b[i][j])lst=j;
            nxt[i][j]=j-lst;
        }
    }
    for(int j=1;j<=n;j++){
        for(int i=n;i>=1;i--){
            while(top&&nxt[stk[top]][j]>nxt[i][j]){
                pre[stk[top]][j]=i;
                top--;
            }
            stk[++top]=i;
        }
        while(top)pre[stk[top]][j]=0,top--;
    }
    long long res=0;
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            f[i][j]=f[pre[i][j]][j];
            f[i][j]+=(i-pre[i][j])*nxt[i][j];
            (res+=i*j-f[i][j])%=mod;
        }
    }
    return res;
}
int main(){ 
    freopen("andorsum.in","r",stdin);
    freopen("andorsum.out","w",stdout);
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            scanf("%d",&a[i][j]);
        }
    }
    for(int i=0;i<=31;i++){
        ans+=1ll*(1ll<<i)*calc1(i)%mod;
        ans%=mod;
    }
    printf("%lld ",ans);
    ans=0;
    for(int i=0;i<=31;i++){
        ans+=1ll*(1ll<<i)*calc2(i)%mod;
        ans%=mod;
    }
    printf("%lld\n",ans);
    return 0;
}