小蓝的好友
Description
终于到达了这次选拔赛的最后一题,想必你已经厌倦了小蓝和小白的故事,为了回馈各位比赛选手,此题的主角是贯穿这次比赛的关键人物——小蓝的好友。
在帮小蓝确定了旅游路线后,小蓝的好友也不会浪费这个难得的暑假。与小蓝不同,小蓝的好友并不想将时间花在旅游上,而是盯上了最近发行的即时战略游戏——SangoCraft。但在前往通关之路的道路上,一个小游戏挡住了小蓝的好友的步伐。
“国家的战争其本质是抢夺资源的战争”是整款游戏的核心理念,这个小游戏也不例外。简单来说,用户需要在给定的长方形土地上选出一块子矩形,而系统随机生成了N个资源点,位于用户所选的长方形土地上的资源点越多,给予用户的奖励也越多。悲剧的是,小蓝的好友虽然拥有着极其优秀的能力,但同时也有着极差的RP,小蓝的好友所选的区域总是没有一个资源点。
终于有一天,小蓝的好友决定投诉这款游戏的制造厂商,为了搜集证据,小蓝的好友想算出至少包含一个资源点的区域的数量。作为小蓝的好友,这自然是你分内之事。
Input
每个输入文件中仅包含一个测试数据。
第一行包含两个由空格隔开的正整数R,C,N,表示游戏在一块[1,R]X[1,C]的地图上生成了N个资源点。
接下来有N行,每行包含两个整数 x,y,表示这个资源点的坐标(1<=x<=R,1<=Y<=c)。
Output
输出文件应仅包含一个整数,表示至少包含一个资源点的区域的数量。具体的说,设N个资源点的坐标为(i=1..n),你需要计算有多少个四元组(LB,DB,RB,UB)满足1<=LB<=RB<=R,1<=DB<=UB<=C,且存在一个i使得LB<=Xi<=RB,DB<=Yi<=UB均成立。
Sample Input
5 5 4
1 2
2 3
3 5
4 1
Sample Output
139
HINT
【数据范围】
对于100%的数据,R,C<=40000,N<=100000,资源点的位置两两不同,且位置为随机生成。
一道被搬进校内互测的神题,现场无人AC……
当时搬题人给了个70pts的
O(n2)
部分分,其他dalao都用了单调栈,就咱这种数据结构学傻选手强行用主席树加了一个
log
……
思路:
考虑补集转化,也就是计算不合法的矩形数量。
考虑逐行扫描,每次计算以这一行为底边的空矩形个数。
可以发现这跟每一个位置向上第一个遇到的资源点位置有关。
那么可以发现每一时刻咱们要维护的都是这样一个东西:
于是把这个奇怪的东西找个数据结构维护一下就好了~
咱校内互测的时候就用了主席树去维护~
然而主席树会T。
考虑换成平衡树,以高度为关键字建出出一棵笛卡尔树。
此时记每个节点高度为
h[i]
,子树大小为
siz[i]
,那么单个节点的贡献为:
这样只需要支持整棵树 h[i]+1 ,和将某个节点的 h[i] 改为 0 两个操作辣~
这样做貌似复杂度依旧不对。
但是,由于数据随机,这其实也是棵Treap……
于是复杂度得到了保证。
代码比咱的主席树还好写qwq
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=40009;
const int Q=100009;
inline int read()
{
int x=0;char ch=getchar();
while(ch<'0' || '9'<ch)ch=getchar();
while('0'<=ch && ch<='9')x=x*10+(ch^48),ch=getchar();
return x;
}
int rt;
namespace treap
{
const int M=N;
int ch[M][2],siz[M],fa[M];
ll val[M],sum[M],tag[M];
inline void calc(int x)
{
sum[x]=sum[ch[x][0]]+sum[ch[x][1]]+(ll)(siz[x]+1ll)*siz[x]/2*(val[x]-val[fa[x]]);
}
inline void mark(int x,int v){tag[x]+=v;val[x]+=v;calc(x);}
inline void update(int x)
{
if(ch[x][0])calc(ch[x][0]);
if(ch[x][1])calc(ch[x][1]);
siz[x]=siz[ch[x][0]]+siz[ch[x][1]]+1;
calc(x);
}
inline void push(int x)
{
if(tag[x])
{
if(ch[x][0])mark(ch[x][0],tag[x]);
if(ch[x][1])mark(ch[x][1],tag[x]);
tag[x]=0;
}
}
inline void rotate(int x)
{
int y=fa[x],z=fa[y];
int l=(ch[y][1]==x);
if(z)ch[z][ch[z][1]==y]=x;
fa[x]=z;fa[y]=x;fa[ch[x][l^1]]=y;
ch[y][l]=ch[x][l^1];ch[x][l^1]=y;
update(y);update(x);
if(y==rt)rt=x;
}
inline int build(int l,int r)
{
if(l>r)return 0;
int mid=l+r>>1,lc,rc;
if(ch[mid][0]=build(l,mid-1))
fa[ch[mid][0]]=mid;
if(ch[mid][1]=build(mid+1,r))
fa[ch[mid][1]]=mid;
siz[mid]=r-l+1;
return mid;
}
inline void modify(int x)
{
static int stk[M],top;
stk[top=1]=x;
for(int i=x;i!=rt;i=fa[i])stk[++top]=fa[i];
while(top)push(stk[top--]);
val[x]=0;
while(x!=rt)rotate(x);
update(x);
}
}
int n,m,q;
vector<int> g[N];
int main()
{
n=read();m=read();q=read();
for(int i=1;i<=q;i++)
{
int x=read(),y=read();
g[x].push_back(y);
}
for(int i=1;i<=n;i++)
if(g[i].size())
sort(g[i].begin(),g[i].end());
rt=treap::build(1,m);
ll ans=0;
for(int i=1;i<=n;i++)
{
treap::mark(rt,1);
for(int j=0,ej=g[i].size();j<ej;j++)
treap::modify(g[i][j]);
ans+=treap::sum[rt];
}
printf("%lld\n",(n*(n+1ll)/2ll)*(m*(m+1ll)/2ll)-ans);
return 0;
}