| 比赛 |
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;
}