https://www.luogu.org/problemnew/show/P2611
题解
\(n\times m\)肯定过不去。。
我们把给定的点看做障碍点，考虑先补集转化为求全空矩阵。
然后我们枚举每一行，令这一行每个点的权值为从这点向上的极大不包含障碍点的连续段。
然后对这个序列建立笛卡尔树，那么答案为：
\[ f[x]=(h[x]-h[fa[x]])*\frac{szie[x]*(size[x]+1)}{2} \]
我们的笛卡尔树上的的每个节点都要维护一个这样的信息。
现在我们还需要扫描每一行，动态维护这颗笛卡尔树。
如果这行没有障碍点，我们整体加个1就好了，这个可以直接打标记。
对于障碍点，相当于这个位置的值变成了0，那么我们把这个点旋转上来就好了，通过手玩我们可以发现\(rotate\)操作不会破坏除了这个点以外的其他点的笛卡尔树结构，于是我们可以一直\(rotate\)把这个点转上去，顺便更新一下答案就好了，因为是随机的数据，所以每次期望操作次数是\(log\)的。
注意如果按照上面的\(\Delta h\)那样算贡献的话如果一个点的父亲改变了的话这个点需要重新\(pushup\)一次。
代码
#include<bits/stdc++.h>
#define N 100009
#define ls tr[x][0]
#define rs tr[x][1]
using namespace std;
typedef long long ll;
vector<int>vec[N];
vector<int>::iterator it;
int tr[N][2],fa[N],size[N],h[N],la[N],n,m,num,rot;
ll dp[N],ans;
inline ll rd(){
ll x=0;char c=getchar();bool f=0;
while(!isdigit(c)){if(c=='-')f=1;c=getchar();}
while(isdigit(c)){x=(x<<1)+(x<<3)+(c^48);c=getchar();}
return f?-x:x;
}
inline ll calc(ll x){return x*(x+1)/2;}
inline bool ge(int x){return tr[fa[x]][1]==x;}
inline void pushup(int x){
size[x]=size[ls]+size[rs]+1;
dp[x]=dp[ls]+dp[rs]+1ll*(h[x]-h[fa[x]])*calc(size[x]);
}
inline void rotate(int x){
int y=fa[x],o=ge(x);
tr[y][o]=tr[x][o^1];fa[tr[y][o]]=y;
if(fa[y])tr[fa[y]][ge(y)]=x;fa[x]=fa[y];
fa[y]=x;tr[x][o^1]=y;
if(tr[y][o])pushup(tr[y][o]);pushup(y);pushup(x);
}
inline void add(int x,int y){
h[x]+=y;la[x]+=y;
pushup(x);
}
inline void pushdown(int x){
if(la[x]){
if(ls)add(ls,la[x]);
if(rs)add(rs,la[x]);
la[x]=0;
}
}
void _pushdown(int x){
if(fa[x])_pushdown(fa[x]);
pushdown(x);
}
int build(int l,int r){
if(l>r)return 0;
int x=(l+r)>>1;
ls=build(l,x-1);rs=build(x+1,r);
if(ls)fa[ls]=x;if(rs)fa[rs]=x;
size[x]=size[ls]+size[rs]+1;
return x;
}
void dfs(int x){
pushdown(x);
if(ls)dfs(ls);
cout<<x<<" "<<ls<<" "<<rs<<" "<<h[ls]<<" "<<h[rs]<<" "<<h[x]<<" "<<dp[x]<<endl;
if(rs)dfs(rs);
}
int main(){
n=rd();m=rd();num=rd();
rot=build(1,m);
for(int i=1;i<=num;++i){
int x,y;
x=rd(),y=rd();
vec[x].push_back(y);
}
for(int i=1;i<=n;++i){
add(rot,1);
for(it=vec[i].begin();it!=vec[i].end();++it){
int x=*it;
_pushdown(x);
while(fa[x])rotate(x);
h[x]=0;
if(ls)pushup(ls);if(rs)pushup(rs);
pushup(x);
rot=x;
}
ans+=dp[rot];
}
printf("%lld",calc(n)*calc(m)-ans);
return 0;
}