先令编号从\(1\)开始。我们要求\([1,i]\)这些数字能否构成一个矩形。
考虑能否用线段树维护，让每个叶子节点\(i\)表示前\(i\)个数能否构成矩形。
一种方法是维护前\(i\)个点最左上点和最右下点的坐标，直接判断这两个点构成的矩形面积是否是\(i\)。
发现修改的时候这个最值不好维护，每次修改可能是\(O(n)\)的。
考虑合法矩形的特征。把前\(i\)个点标记为黑点，其余点是白点。那么前\(i\)个点构成了一个矩形当且仅当：
正确性比较显然...？（雾）不说了。
记左边上边都是白点的黑点数量为\(t1\)，上下左右有两个及以上黑点是白点数量为\(t2\)。注意到\(t1>0\)，\(t2\)非负，那么\(i\)合法当且仅当\(t1+t2=1\)，所以只在叶节点处维护前\(i\)个点为黑点时，\(t1+t2\)的值就好了。非叶节点就维护区间最小值及最小值的个数。
考虑修改时如何维护。
记\(l\)为点\(i\)周围点（上下左右）编号的次小值，点\(i\)作为白点时，会对\([l,i-1]\)这些位置有贡献。
记\(r\)为点\(i\)左边、上边的点的编号的最小值，那么点\(i\)作为黑点时，会对\([i,r-1]\)这些位置有贡献。
每次交换两个点\(x,y\)，最多只会影响\(10\)个点的\(l,r\)，所以把这些点取出来，减掉在线段树上的贡献，交换\(x,y\)之后再把它们的贡献加上即可。
注意要对这些点判重。
//27496ms 1688K
#include "seats.h"
#include <cstdio>
#include <cctype>
#include <algorithm>
#define F(p,i) (p+Way[i])
#define ID(x,y) ((x-1)*m+y)
#define Check(x,y) (x>=1&&x<=n&&y>=1&&y<=m)
#define gc() getchar()
const int N=1e6+5,Way[]={1,0,-1,0,1};//down left up right
int n,m,tot,A[N],X[N],Y[N],val[N];
struct Segment_Tree
{
#define ls rt<<1
#define rs rt<<1|1
#define lson l,m,ls
#define rson m+1,r,rs
#define S N<<2
int mn[S],cnt[S],tag[S];
#undef S
#define Upd(rt,v) tag[rt]+=v, mn[rt]+=v
#define Update(rt) mn[rt]=std::min(mn[ls],mn[rs]), cnt[rt]=(mn[rt]==mn[ls]?cnt[ls]:0)+(mn[rt]==mn[rs]?cnt[rs]:0)
inline void PushDown(int rt)
{
Upd(ls,tag[rt]), Upd(rs,tag[rt]), tag[rt]=0;
}
void Build(int l,int r,int rt)
{
if(l==r)
{
mn[rt]=val[l], cnt[rt]=1;
return;
}
int m=l+r>>1;
Build(lson), Build(rson), Update(rt);
}
void Modify(int l,int r,int rt,int L,int R,int v)
{
if(L<=l && r<=R) {Upd(rt,v); return;}
if(tag[rt]) PushDown(rt);
int m=l+r>>1;
if(L<=m) Modify(lson,L,R,v);
if(m<R) Modify(rson,L,R,v);
Update(rt);
}
}T;
inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-48,c=gc());
return now;
}
inline int CalcW(int p)
{
int mn1=tot+1,mn2=mn1,x=X[p],y=Y[p];
for(int i=0,xn,yn; i<4; ++i)
if(xn=F(x,i),yn=F(y,i+1),Check(xn,yn))
{
int w=A[ID(xn,yn)];
if(w<mn1) mn2=mn1, mn1=w;
else if(w<mn2) mn2=w;
}
return mn2;
}
inline int CalcB(int p)
{
int x=X[p],y=Y[p],xn1=F(x,1),yn1=F(y,2),xn2=F(x,2),yn2=F(y,3);
return std::min(Check(xn1,yn1)?A[ID(xn1,yn1)]:tot+1,Check(xn2,yn2)?A[ID(xn2,yn2)]:tot+1);
}
#define S 1,tot,1
void give_initial_chart(int n,int m,std::vector<int> R,std::vector<int> C)
{
int tot=n*m; ::n=n, ::m=m, ::tot=tot;
for(int i=1; i<=tot; ++i) X[i]=R[i-1]+1,Y[i]=C[i-1]+1,A[ID(X[i],Y[i])]=i;
// for(int i=1; i<=tot; ++i) X[i]=read()+1,Y[i]=read()+1,A[ID(X[i],Y[i])]=i;
for(int i=1; i<=tot; ++i)
{
val[i]=val[i-1];
if(CalcW(i)<i) --val[i];
if(CalcB(i)>i) ++val[i];
for(int j=0,x=X[i],y=Y[i],xn,yn,w; j<4; ++j)
if(xn=F(x,j),yn=F(y,j+1),Check(xn,yn))
{
if((w=A[ID(xn,yn)])<i && CalcB(w)==i) --val[i];
else if(w>i && CalcW(w)==i) ++val[i];
}
}
T.Build(S);
}
int swap_seats(int a,int b)
{
static int B[12];
++a, ++b;
int x=X[a],y=Y[a],t=2; B[1]=a, B[2]=b;
for(int i=0,xn,yn; i<4; ++i)
if(xn=F(x,i),yn=F(y,i+1),Check(xn,yn)) B[++t]=A[ID(xn,yn)];
x=X[b],y=Y[b];
for(int i=0,xn,yn; i<4; ++i)
if(xn=F(x,i),yn=F(y,i+1),Check(xn,yn)) B[++t]=A[ID(xn,yn)];
std::sort(B+1,B+1+t);
for(int i=1; i<=t; ++i)
if(B[i]!=B[i-1])
{
int p=B[i],l=CalcW(p),r=CalcB(p);
if(l<p) T.Modify(S,l,p-1,-1);
if(r>p) T.Modify(S,p,r-1,-1);
}
std::swap(A[ID(X[a],Y[a])],A[ID(X[b],Y[b])]);
std::swap(X[a],X[b]), std::swap(Y[a],Y[b]);
for(int i=1; i<=t; ++i)
if(B[i]!=B[i-1])
{
int p=B[i],l=CalcW(p),r=CalcB(p);
if(l<p) T.Modify(S,l,p-1,1);
if(r>p) T.Modify(S,p,r-1,1);
}
return T.cnt[1];
}
//int main()
//{
// int n=read(),m=read(),Q=read(); ::n=n, ::m=m;
// give_initial_chart(n,m);
// while(Q--) printf("%d\n",swap_seats(read(),read()));
//
// return 0;
//}原文：https://www.cnblogs.com/SovietPower/p/10679775.html