bzoj 3110: [Zjoi2013]K大数查询(树套树,整体二分)
发布时间:2021-03-19 01:54:07 所属栏目:大数据 来源:网络整理
导读:副标题#e# 3110: [Zjoi2013]K大数查询 Time Limit:?20 Sec?? Memory Limit:?512 MB Submit:?4020?? Solved:?1547 [ Submit][ Status][ Discuss] Description 有N个位置,M个操作。操作有两种,每次操作如果是1 a b c的形式表示在第a个位置到第b个位置,每个
|
注意因为每到达一个新的答案区间,线段树都要清零,但是o(n)的时间复杂度太高,所以我们对于线段树中的第一个节点(就是区间[1,n])打上lazy标记并清零,然后每使用到一个节点如果他有lazy标记,就清他的左右儿子,并下放lazy标记。这样每次只需要清需要用的节点,节省时间。 #include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 500003
#define ul unsigned int
using namespace std;
int n,ans[N];
ul tr[N*4],delta[N*4],pd[N*4];
struct data
{
int l,id,num,pd,x;
}a[N];
void clear(int now)
{
tr[now]=delta[now]=0;
pd[now]=1;
}
int cmp(data a,data b)
{
return a.num<b.num;
}
void pushdown(int now,int r)
{
if (pd[now])
{
//tr[now]=delta[now]=0;
clear(now<<1); clear(now<<1|1);
pd[now]=0;
}
int mid=(l+r)/2;
if (delta[now])
{
delta[now<<1]+=delta[now];
delta[now<<1|1]+=delta[now];
tr[now<<1]+=delta[now]*(mid-l+1);
tr[now<<1|1]+=delta[now]*(r-mid);
delta[now]=0;
}
}
void update(int now)
{
tr[now]=tr[now<<1]+tr[now<<1|1];
}
void qjchange(int now,int rr)
{
if (ll<=l&&r<=rr)
{
tr[now]+=(r-l+1);
delta[now]++;
return;
}
int mid=(l+r)/2;
pushdown(now,r);
if(ll<=mid) qjchange(now<<1,rr);
if (rr>mid) qjchange(now<<1|1,rr);
update(now);
}
ul qjsum(int now,int rr)
{
if (ll<=l&&r<=rr)
return tr[now];
int mid=(l+r)/2;
pushdown(now,r);
ul ans=0;
if (ll<=mid) ans+=qjsum(now<<1,rr);
if (rr>mid) ans+=qjsum(now<<1|1,rr);
return ans;
}
void divide (int l,int x,int y)
{
if (l==r)
{
for (int i=x;i<=y;i++)
if (a[i].pd==2)
ans[a[i].id]=l;
return;
}
int mid,pl,pr;
mid=(l+r)/2;
pl=0; pr=y-x+1;
clear(1);
for (int i=x;i<=y;i++)
if (a[i].pd==1)
{
if (a[i].x<=mid) a[i].num=++pl;
else
{
qjchange(1,a[i].l,a[i].r);
// cout<<"!"<<" "<<a[i].l<<" "<<a[i].r<<endl;
a[i].num=++pr;
}
}
else
{
ul t=qjsum(1,a[i].r);
//cout<<t<<endl;
if (t>=a[i].x) a[i].num=++pr;
else
{
a[i].x=a[i].x-t;
a[i].num=++pl;
}
}
sort(a+x,a+y+1,cmp);
divide(l,x,x+pl-1);
divide(mid+1,x+pl,y);
}
int main()
{
//freopen("a.in","r",stdin);
scanf("%d%d",&m);
for (int i=1;i<=m;i++)
scanf("%d%d%d%d",&a[i].pd,&a[i].l,&a[i].r,&a[i].x),a[i].id=i;
divide(0,m);
for (int i=1;i<=m;i++)
if (ans[i]) printf("%dn",ans[i]);
}
(编辑:佛山站长网) 【声明】本站内容均来自网络,其相关言论仅代表作者个人观点,不代表本站立场。若无意侵犯到您的权利,请及时与联系站长删除相关内容! |



