ZOJ 3955 Saddle Point

#include<bits/stdc++.h>
using namespace std;

int n,m;
long long mod = 1e9+7;
int a[1200][1200];
int num1[1200][1200];
int num2[1200][1200];

long long b[1200];

struct X
{
int id,val;
}s[1200];

bool cmp1(X a,X b)
{
return a.val<b.val;
}

bool cmp2(X a,X b)
{
return a.val>b.val;
}

int main()
{
b[0]=1;
for(int i=1;i<=1000;i++) b[i]=(b[i-1]*2)%mod;

int T; scanf("%d",&T);

while(T--)
{
scanf("%d%d",&n,&m);

memset(num1,0,sizeof num1);
memset(num2,0,sizeof num2);

for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) scanf("%d",&a[i][j]);

for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++) s[j].id=j, s[j].val=a[i][j];
sort(s+1,s+1+m,cmp2);
for (int j=1,k;j<=m;j=k)
{
for (k=j;k<=m&&s[k].val==s[j].val;k++)
{
num1[i][s[k].id]=j-1;
}
}
}

for(int j=1;j<=m;j++)
{
for(int i=1;i<=n;i++) s[i].id=i, s[i].val=a[i][j];
sort(s+1,s+1+n,cmp1);
for (int i=1,k;i<=n;i=k)
{
for (k=i;k<=n&&s[k].val==s[i].val;k++)
{
num2[s[k].id][j]=i-1;
}
}
}

long long ans=0;

for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
long long A=b[num1[i][j]],B=b[num2[i][j]];

long long p=A*B%mod;
ans=(ans+p)%mod;

}
}

printf("%lld\n",ans);

}
return 0;
}