bzoj 3110: [Zjoi2013]K大数查询(树套树,整体二分)
发布时间:2021-03-15 22:14:58 所属栏目:大数据 来源:网络整理
导读: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]); } (编辑:晋中站长网) 【声明】本站内容均来自网络,其相关言论仅代表作者个人观点,不代表本站立场。若无意侵犯到您的权利,请及时与联系站长删除相关内容! |