蒟蒻求助与或和

#include <stdio.h>
#include <string.h>
int n;
int ans;
int height[1005];
int stack[1005];
bool matrix[32][1005][1005];
int main()
{
    scanf("%d",&n);
    for(int i=1; i<=n; i++)
    {
        for(int j=1; j<=n; j++)
        {
            int k;
            scanf("%d",&k);
            for(int bit=0; bit<32; bit++)
            {
                matrix[bit][i][j]=k&1;
                k>>=1;
            }
        }
    }

    for(int bit=0;bit<32;bit++){
        memset(height,0,sizeof height);
        int temp=0;
        for(int i=1;i<=n;i++){
            int top=0;
            for(int j=1;j<=n;j++){
                height[j]=matrix[bit][i][j]?height[j]+1:0;
                temp+=matrix[bit][i][j];
                while(top && height[stack[top]]>=height[j]){
                    temp-=(stack[top]-stack[top-1])*(height[j]-height[stack[top]]);
                    top--;
                }
                stack[++top]=j;
            }
        }
        ans+=temp*(1<<bit);
//        printf("%d ",ans);
    }
    printf("%d ",ans);
    return 0;
}