#include "rect.h"
#include <bits/stdc++.h>
#define frst first
#define scnd second
using namespace std;
#define MAX 1000000000
typedef vector<int> vi;
typedef pair<int,int> ii;
typedef pair<ii,ii> iii;
typedef vector<ii> vii;
vector<vi> b;
int vis[2510][2510];
int Y[4]={-1,0,1,0};
int X[4]={0,-1,0,1};
int c1,r1,c2,r2;
bool sw1=0;
int flod(int y,int x){
//cout<<y<<" "<<x<<endl;
c1=min(c1,x);
r1=min(r1,y);
c2=max(c2,x);
r2=max(r2,y);
vis[y][x]=1;
int ans=1;
for(int i=0;i<4;i++){
if(y+Y[i]<b.size()-1 and y+Y[i]>=1 and x+X[i]<b[0].size()-1 and x+X[i]>=1){
if(b[y+Y[i]][x+X[i]]==0 and !vis[y+Y[i]][x+X[i]])
ans+=flod(y+Y[i],x+X[i]);
}
else{
if(b[y+Y[i]][x+X[i]]==0){
sw1=1;
}
}
}
return ans;
}
int ma,u,d,l,r,row,i,j,k,o,m;
ii v,x;
bool sw=0;
long long count_rectangles(std::vector<std::vector<int> > a) {
long long c=0;
if((a.size()<=700 and a[0].size()<=700) or (a.size()<=3 and a[0].size()<=2500)){
vector<vii> ud,lr;
ud.resize(a[0].size());
lr.resize(a.size());
for( i=1;i<a[0].size()-1;i++){
for( j=0;j<a.size()-1;j++){
ma=a[j+1][i];
for( k=j+2;k<a.size();k++){
if(ma<a[j][i] and ma<a[k][i]){
ud[i].push_back(ii(j,k));
}
if(ma>a[j][i]) break;
ma=max(ma,a[k][i]);
}
}
}
for(i=1;i<a.size()-1;i++){
for( j=0;j<a[0].size()-1;j++){
ma=a[i][j+1];
for(k=j+2;k<a[0].size();k++){
if(ma<a[i][j] and ma<a[i][k]){
lr[i].push_back(ii(j,k));
}
if(ma>a[i][j]) break;
ma=max(ma,a[i][k]);
}
}
}
for(i=1;i<ud.size();i++){
for(o=0;o<ud[i].size();o++){
v=ud[i][o];
u=v.frst; d=v.scnd;
l=i-1; r=i+1;
for( j=i;ud.size();j++){
sw=0;
if(j!=i){
for(m=0;m<ud[j].size();m++){
x=ud[j][m];
if(u==x.frst and d==x.scnd){
sw=1;
break;
}
if(u<x.first) break;
}
if(sw==0) break;
else r=j+1;
}
row=(d-u)-1;
for( k=u+1;k<lr.size();k++){
sw=0;
for(m=0;m<lr[k].size();m++){
x=lr[k][m];
if(x.first==l and x.second==r){
sw=1;
break;
}
if(l<x.first) break;
}
if(sw==1) row--;
else break;
if(row==0) break;
}
if(row==0) c++;
}
}
}
}
else{
memset(vis,0,sizeof vis);
b=a;
for(int i=1;i<a.size()-1;i++){
for(int j=1;j<a[i].size()-1;j++){
if(a[i][j]==0 and !vis[i][j]){
c1=r1=MAX;
c2=r2=-MAX;
sw1=0;
int ar=flod(i,j);
if((abs(c1-c2)+1)*(abs(r1-r2)+1)==ar and sw1==0) c++;
}
}
}
}
return c;
}
