ZOJ2112 & BZOJ1901 玩坏了的主席树

最新推荐文章于 2019-05-30 17:10:58 发布
最新推荐文章于 2019-05-30 17:10:58 发布
阅读量4k 收藏 2
点赞数
分类专栏： 高级数据结构
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/d891320478/article/details/9634259
版权

今天才发现加上回收空间的话主席树的空间是nlognlogn→ →之前总当成和树套树一样是nlogn了【你特喵在逗我

在BZOJ上很久前过掉了，今天滚回来写ZOJ各种卡空间

ZOJ上32768K实在是……卡死了

所以，为了压缩空间，我们拆成两部分来做。

首先静态建树，直接把给定的序列离散化丢到一棵树里【就像POJ2104那样玩】，这里的空间是nlogn，如果直接用树状数组去套的话空间就是nlognlogn在ZOJ上就无限MLE了……

那么更新怎么办呢？

再造一棵树用来记录更新，这棵树要用树状数组来套，这里的空间最多是mlognlogn。

这样总的空间复杂度就是nlogn+mlognlogn，因为总的操作数m<=10000，但是n<=50000，于是空间就降下来了，在ZOJ上基本上31000K左右的空间就能AC了。

这方法玩得真脱呀……………………

ZOJ 2112
#define NN 50005
#define NM 10005
#define MM 2600010
int root[NN],rt1[NN],a[NN],b[NN+NM],n,q,node,lb,poi;
int ord[NM][4]; int bit[NN];
struct sgtree{ int l,r,s; }tree[MM];
int sta[NM],top;
int fastget()
{
	char c; int ans=0; c=getchar();
	while (! (c>='0' && c<='9')) c=getchar();
	while (c>='0' && c<='9')
	{
		ans=(ans<<3)+(ans<<1)+c-'0';
		c=getchar();
	}
	return ans;
}
int search(int x)
{
	int l=1,r=b[0],mid;
	while (l<=r)
	{
		mid=l+r>>1; if(b[mid]==x) return mid;
		if(b[mid]<x) l=mid+1; else r=mid-1;
	}
}
int build(int l,int r)
{
	int k=++node;
	tree[k].s=0;
	if(l==r) return k;
	int mid=l+r>>1;
	if(l<=mid) tree[k].l=build(l,mid);
	if(r>mid) tree[k].r=build(mid+1,r);
	return k;
}
int sum(int x)
{
	int ans=0;
	for(;x;x-=lowbit(x)) ans+=tree[tree[bit[x]].l].s;
	return ans;
}
int query(int x,int y,int xx,int yy,int l,int r,int k)
{
	if(l==r) return l;
	int ans=tree[tree[yy].l].s-tree[tree[xx].l].s+sum(y)-sum(x);
	int mid=l+r>>1;
	if(ans>=k)
	{
		for(int i=x;i;i-=lowbit(i)) bit[i]=tree[bit[i]].l;
		for(int i=y;i;i-=lowbit(i)) bit[i]=tree[bit[i]].l;
		return query(x,y,tree[xx].l,tree[yy].l,l,mid,k);
	}
	else
	{
		for(int i=x;i;i-=lowbit(i)) bit[i]=tree[bit[i]].r;
		for(int i=y;i;i-=lowbit(i)) bit[i]=tree[bit[i]].r;
		return query(x,y,tree[xx].r,tree[yy].r,mid+1,r,k-ans);
	}
}
int update(int rt,int l,int r,int x,int d)
{
	int k=++node;
	tree[k]=tree[rt];
	tree[k].s+=d;
	if(l==r) return k;
	int mid=l+r>>1;
	if(x<=mid) tree[k].l=update(tree[rt].l,l,mid,x,d);
	else tree[k].r=update(tree[rt].r,mid+1,r,x,d);
	return k;
}
int main()
{
freopen("1.in","r",stdin);
	int tt=fastget();
	while(tt--)
	{
		top=node=0;
		lb=n=fastget(); q=fastget();
		for(int i=1;i<=n;i++) b[i]=a[i]=fastget();
		for(int i=1;i<=q;i++)
		{
			char s[3]; scanf("%s",s);
			if(s[0]=='Q'){ ord[i][0]=1; ord[i][1]=fastget(); ord[i][2]=fastget(); ord[i][3]=fastget(); }
			else{ ord[i][0]=2; ord[i][1]=fastget(); b[++lb]=ord[i][2]=fastget(); }
		}
		sort(b+1,b+1+lb); b[0]=1;
		for(int i=2;i<=lb;i++) if(b[i]!=b[b[0]]) b[++b[0]]=b[i];
		for(int i=1;i<=n;i++) a[i]=search(a[i]);
		root[0]=build(1,b[0]); poi=node;
		for(int i=1;i<=n;i++) rt1[i]=root[0];
		for(int i=1;i<=n;i++)
			root[i]=update(root[i-1],1,b[0],a[i],1);
		for(int _=1;_<=q;_++)
		{
			if(ord[_][0]==1)
			{
				int u=ord[_][1],v=ord[_][2],w=ord[_][3];
				for(int i=u-1;i;i-=lowbit(i)) bit[i]=rt1[i];
				for(int i=v;i;i-=lowbit(i)) bit[i]=rt1[i];
				printf("%d\n",b[query(u-1,v,root[u-1],root[v],1,b[0],w)]);
			}
			else
			{
				int u=ord[_][1],w=ord[_][2];
				for(int i=u;i<=n;i+=lowbit(i)) rt1[i]=update(rt1[i],1,b[0],a[u],-1);
				a[u]=search(w);
				for(int i=u;i<=n;i+=lowbit(i)) rt1[i]=update(rt1[i],1,b[0],a[u],1);
			}
		}
	}
	return 0;
}


蛋疼的nlognlogn的空间的方法【BZOJ1901】 

#define NN 10005
#define NM 10005
#define MM 2500010
int root[NN],a[NN],b[NN+NM],n,q,node,lb,poi;
int ord[NM][4]; int bit[NN];
struct sgtree{ int l,r,s; }tree[MM];
int sta[NM],top;
int fastget()
{
    char c; int ans=0; c=getchar();
    while (! (c>='0' && c<='9')) c=getchar();
    while (c>='0' && c<='9')
    {
        ans=(ans<<3)+(ans<<1)+c-'0';
        c=getchar();
    }
    return ans;
}
int search(int x)
{
    int l=1,r=b[0],mid;
    while (l<=r)
    {
        mid=l+r>>1; if(b[mid]==x) return mid;
        if(b[mid]<x) l=mid+1; else r=mid-1;
    }
}
int build(int l,int r)
{
    int k=++node;
    tree[k].s=0;
    if(l==r)
    {
        tree[k].l=tree[k].r=0;
        return k;
    }
    int mid=l+r>>1;
    if(l<=mid) tree[k].l=build(l,mid);
    if(r>mid) tree[k].r=build(mid+1,r);
    return k;
}
int sum(int x)
{
    int ans=0;
    for(;x;x-=lowbit(x)) ans+=tree[tree[bit[x]].l].s;
    return ans;
}
int query(int x,int y,int l,int r,int k)
{
    if(l==r) return l;
    int ans=sum(y)-sum(x);
    int mid=l+r>>1;
    if(ans>=k)
    {
        for(int i=x;i;i-=lowbit(i)) bit[i]=tree[bit[i]].l;
        for(int i=y;i;i-=lowbit(i)) bit[i]=tree[bit[i]].l;
        return query(x,y,l,mid,k);
    }
    else
    {
        for(int i=x;i;i-=lowbit(i)) bit[i]=tree[bit[i]].r;
        for(int i=y;i;i-=lowbit(i)) bit[i]=tree[bit[i]].r;
        return query(x,y,mid+1,r,k-ans);
    }
}
int update(int rt,int l,int r,int x,int d)
{
    int k;
    if(top) k=sta[--top];
    else k=++node;
    tree[k]=tree[rt];
    tree[k].s+=d;
    if(l==r)
    {
        tree[k].l=tree[k].r=0;
        return k;
    }
    int mid=l+r>>1;
    if(x<=mid)
    {
        if(tree[k].l>poi) sta[top++]=tree[k].l;
        tree[k].l=update(tree[rt].l,l,mid,x,d);
    }
    else
    {
        if(tree[k].r>poi) sta[top++]=tree[k].r;
        tree[k].r=update(tree[rt].r,mid+1,r,x,d);
    }
    return k;
}
int main()
{
        top=node=0;
        lb=n=fastget(); q=fastget();
        for(int i=1;i<=n;i++) b[i]=a[i]=fastget();
        for(int i=1;i<=q;i++)
        {
            char s[3]; scanf("%s",s);
            if(s[0]=='Q'){ ord[i][0]=1; ord[i][1]=fastget(); ord[i][2]=fastget(); ord[i][3]=fastget(); }
            else{ ord[i][0]=2; ord[i][1]=fastget(); b[++lb]=ord[i][2]=fastget(); }
        }
        sort(b+1,b+1+lb); b[0]=1;
        for(int i=2;i<=lb;i++) if(b[i]!=b[b[0]]) b[++b[0]]=b[i];
        for(int i=1;i<=n;i++) a[i]=search(a[i]);
        root[0]=build(1,b[0]); poi=node;
        for(int i=1;i<=n;i++) root[i]=root[0];
        for(int i=1;i<=n;i++)
        for(int x=i;x<=n;x+=lowbit(x))
            root[x]=update(root[x],1,b[0],a[i],1);
        for(int _=1;_<=q;_++)
        {
            if(ord[_][0]==1)
            {
                int u=ord[_][1],v=ord[_][2],w=ord[_][3];
                for(int i=u-1;i;i-=lowbit(i)) bit[i]=root[i];
                for(int i=v;i;i-=lowbit(i)) bit[i]=root[i];
                printf("%d\n",b[query(u-1,v,1,b[0],w)]);
            }
            else
            {
                int u=ord[_][1],w=ord[_][2];
                for(int i=u;i<=n;i+=lowbit(i)) root[i]=update(root[i],1,b[0],a[u],-1);
                a[u]=search(w);
                for(int i=u;i<=n;i+=lowbit(i)) root[i]=update(root[i],1,b[0],a[u],1);
            }
        }
    return 0;
}



幽诗人
关注
专栏目录
ZOJ 2112 Dynamic Rankings （动态主席树）
满月照亮的路
02-16 202
题目链接：http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=2112 题目大意： 比静态查询区间第k小多加了一个修改操作。 题目分析：  对于所有的数（原序列的数和要修改的数）都离散化处理， 动态修改的核心是对于树状数组中出现的点开线段树维护， 因为之前修改一个数对于后续区间有影响，我们令开一个根节点数组S 来维...
zoj 2112 Dynamic Rankings（动态第k大，树状数组套主席树）
梦幻的蔷薇色
11-01 420
这个题在kuangbin大佬的模板里有。 http://www.cnblogs.com/kuangbin/p/3308118.html kuangbin大佬博客里有四个学习连接。#include <bits/stdc++.h> using namespace std;const int MAXN = 60010; struct node { int l,r,sum; }T[MAXN*30
zoj2112
08-26 116
这道题是要求待修改的区间第k值，看的kuangbin的树状数组加主席树，这种方法是离线的。本来是想做hdu5412的，结果超时，所以只能交了zoj2110。没什么多说的，就是模板题。如果要求第k大，只要做相应的修改就行了，有多种修改的方法，只要能的出正确答案就行。 参考代码：http://www.cnblogs.com/kuangbin/p/3308118.html （今天为...
zoj 2112 （主席树，树状数组套线段树）
know_heng的博客
09-22 330
Dynamic Rankings Time Limit: 10 Seconds      Memory Limit: 32768 KB The Company Dynamic Rankings has developed a new kind of computer that is no longer satisfied with the query like to simply
ZOJ2112
binarycopycode
02-01 342
动态求区间第k小，整体二分蛤不清be（菜哭.jpg ，给Qt的板子注释了一蛤 #include #include #include using namespace std; const int N = 5e4 + 7; const int M = 1e5 + 7; const int INF = 1 << 30; struct Quiry { int x,y,z,id,tp; }Q[ N
zoj2112 分类： templates z...
aa288288的博客
07-20 109
树状数组套主席树 好像没有什么特别的地方，就是要注意卡空间。 先处理出初始序列的主席树 O(N∗log(N+M)) 然后用树状数组记录修改 O(M∗logN∗log(N+M)) 时间复杂度/空间复杂度: O((N+M∗logN)∗log(N+M)) #include <cstdio> #include <cstdlib...
ZOJ 2112(树套树)
cy41的博客
05-30 329
题目链接：http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=1112 静态的主席树用于记录静态的未修改的序列的贡献（板子） 对于单点修改来说，修改id处的数字相当于将该棵主席树以及在此基础上的主席树都修改掉，那么单次修改的复杂度就会达到n*log(n)，而这个操作就相当于区间修改吧，那么最适合做区间修改，单点查询的是什么...
树套树之线段树套平衡树（ZOJ2112 Dynamic Rankings)
Sunshine_cfbsl
08-11 1085
接上篇：树套树之线段树套线段树 树套树都应该是做一道题就会了的。。 Dynamic Rankings(ZOJ 2112) 题目描述： 某公司开发了一种新型的计算机，它不再仅能找到给定N个数中第k小的数。它变成了一个更强大的系统。对于N个数a[1],a[2],…,a[N]，你可以询问a[i],a[i+1],…,a[j-1],a[j]中第k小的数(i≤j,0<k≤j+i−1i\le j,0\lt
bzoj 3730. 震波（动态点分治 + vector树状数组）
01-06
考虑点分树的做法，在点分树上每个点以距离为下标建一棵线段树，每次询问查询子树的贡献，再暴力向上跳合并父节点来自其它节点的贡献。 因为子树树形被破坏，在做减法时，子节点子树对父节点的贡献不能用子树维护的...
zoj 1002_zoj1002_
09-29
【标题】"ZOJ 1002" 是一个在线编程竞赛题目，源自ZOJ（Zhejiang Online Judge），这是一个面向ACM/ICPC（国际大学生程序设计竞赛）的在线评测系统。题目编号1002，通常表示该题是ZOJ平台上的一个问题，可能涉及算法...
zoj 2112 （主席树）
Naruto_____的专栏
11-09 674
Orz...... #include using namespace std; const int MAXN = 60010; const int MAXQ = 10010; const int MAXM = 2500010; int n, m, dex; int a[MAXN], sorted[MAXN], root[MAXN], lower[MAXN], temp[MAXN]; int
zoj-2112(主席树动态求区间第k小数)
AC之路
08-25 2410
总算是把动态求区间第k个数的算法看明白了。 在主席树的基础上，如果有修改操作，则要通过套树状数组来实现任意区间求第k小的问题。 刚开始看不明白什么意思，现在有一点理解。树状数组的每个元素是一个线段树，来维护修改后的前后缀和，树状数组能在log时间内更整个数组，现在用相同的方式更新整个线段树数组，每次更新一个点时，要更新这个点代表的整个线段树。同样的，求和时用一个use数组记录所要更新的点的下标
ZOJ 2112
终究我要华丽逆袭
11-24 796
我擦、、真难，看了mj的代码，才会写、、WA了无数，终于AC，删除写错几次，线段树写错几次，果然数据结构不多写就悲剧。。。 做法是线段树里面套个SBT，这样替换操作，相当于每个包含这点的线段的Sbt执行一次insert和del，然后对于询问，二分答案，看>于它的是否超过k个。取下界即可。 #include #include int T,n,m,a,b,c,tot,p[50010]; char
zoj 2112(主席树套树状数组+优化)
过路人丶的博客
12-02 328
题解思路:如果直接对原来的数组建立主席树套树状数组的空间复杂度是（（n+m）log(n+m)*log(n+m)）这样明显爆炸，那么我们可以对原来的数组建主席树模型，空间复杂度就是n*logn，对m里面的修改建主席树套树状数组空间复杂度是m*logm*logm明显小了很多。 #include #define lson l,mid #define rson mid+1,r using nam
zoj 2112 主席树
没有路通往的地方
11-15 620
2014/11/15  老规矩，想学主席树的这边走疯狂的橡树（推荐），这里只写点自己敲代码的时候遇到的问题。 这颗树是从昨天晚上8点开始敲的，敲到了现在........ 其他部分没问题： 先把数组离散化a[ i ] 建个空树 然后以每个a[ i ]节点为根节点根据前一个节点的树来建树。（每棵树省下一半的树节点） 然后问题来了： 建树状数组的时候，其实只理解到了用根节
【动态主席树】ZOJ 2112【树状数组+主席树】
Gene_I must wait for the sunrise
11-18 339
题意： 给定一个区间，求这个区间第k小的数，支持单点修改。 思路： 动态主席树裸题。 我们先来回顾一下静态主席树的做法，对于数组中每一个位置都维护一棵权值线段树，该权值线段树保存的是区间 [1,x] 的信息。因此我想要求区间 [l,r] 之间第k大的时候，只需要将root[r]-root[l-1]就是维护区间 [l,r] 信息的权值线段树，因此就可...
【ZOJ 2112】 线段树套平衡树
wodasini的专栏
04-14 287
题目大意： 给定一个数列，定义两种操作：1、修改第n个数的值2、求[l,r]区间内第k大的值。方法：线段树维护区间，treap树维护第k大（splay当然也可以）树套树第一题。。没想到真的是每个线段树上的节点建一棵平衡树。。。。。然后修改值就是在所有相关区间内的平衡树里，erase一个值再insert一个值。。。。简直恐怖查询就是二分答案，因为每个区间都有一棵树，所以可以很快就可以在O(lgn）里...
ZOJ 2112 Dynamic Rankings （动态第k大，树状数组套主席树）
weixin_33919950的博客
09-08 123
Dynamic RankingsTime Limit: 10 Seconds      Memory Limit: 32768 KB The Company Dynamic Rankings has developed a new kind of computer that is no longer satisfied with the query like to simply find the...
ZOJ详细题目分类指南，初学者&模拟问题搞定法。
ZOJ题目分类详细版包括初学者题目和模拟问题，分别涵盖了一系列题目，是非常详细的解题分类，对初学者学习十分有用。其中初学者题目包括1001、1037、1048、1049、1051、1067、1115、1151、1201、1205、1216、1240、...
评论 4
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值