bzoj1500/洛谷P2042 维护数列 splay与无旋treap两种解法

最新推荐文章于 2020-06-06 11:43:00 发布
最新推荐文章于 2020-06-06 11:43:00 发布
阅读量632 收藏 1
点赞数 2
分类专栏： 数据结构 文章标签： splay 伸展树 无旋treap
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/litble/article/details/78804097
版权

splay

终于完成了！这个splay板子非常牛逼…本蒟蒻在做的时候犯了很多错误，调了三天。
数据结构题，当然要面向代码，现在将比较难懂的函数梳理一下。如果您只想看完整代码，请跳到最下面。

要维护的东西

a:原来的序列
b:插入操作中的序列
v:节点的值
siz:以该节点为根的子树大小
son:该节点的左右儿子
f:该节点的父亲
st:存放可以使用的数组下标的一个栈（即将删去了的节点原来的标记存入其中）
same:区间赋值的懒惰标记
rev:旋转的懒惰标记
sum:子树的权值和
mx:子树中的最大子序列和
ls:包含子树所代表的序列左端点的最大子序列（如果左端点是负值，则为０）
rs:包含子树所代表的序列右端点的最大子序列（如果右端点是负值，则为０）
这题不需要特别注意的基本操作都可以看本蒟蒻的另一篇博客：点我起飞

pushup

void up(int x) {
	int l=son[x][0],r=son[x][1];
	sum[x]=sum[l]+sum[r]+v[x],siz[x]=siz[l]+siz[r]+1;
	ls[x]=max(ls[l],sum[l]+v[x]+ls[r]);
	rs[x]=max(rs[r],sum[r]+v[x]+rs[l]);
	mx[x]=max(max(mx[l],mx[r]),rs[l]+ls[r]+v[x]);
}

…好像看代码能懂吧？就不详细讲解了。

pushdown

别的地方都没有什么问题，就是本蒟蒻在打rev标记的时候出过两个错误，和大家分享一下：
错误１：打了rev标记后就直接跑路了，然后在处理到有rev标记的点的时候再将其左右子树反过来。这样会导致打完rev标记后，该节点祖先的ls和rs出错从而导致mx出错。
错误２：在交换完左右子树后，pushup了一下。这是不对的，因为交换完左右子树后只是部分翻转了，而事实上我们应该把这颗子树完全翻转，所以不能pushup而是要交换ls和rs

void pd(int x) {
	int l=son[x][0],r=son[x][1],t=same[x];
	if(same[x]!=inf) {
		if(l) {//注意防止处理０
			same[l]=v[l]=t,sum[l]=siz[l]*t;
			if(t>=0) ls[l]=rs[l]=mx[l]=sum[l];
			else ls[l]=rs[l]=0,mx[l]=t;
		}
		if(r) {
			same[r]=v[r]=t,sum[r]=siz[r]*t;
			if(t>=0) ls[r]=rs[r]=mx[r]=sum[r];
			else ls[r]=rs[r]=0,mx[r]=t;
		}
		same[x]=inf;
	}
	if(rev[x]) {
		if(l) rev[l]^=1,swap(son[l][0],son[l][1]),swap(ls[l],rs[l]);
		if(r) rev[r]^=1,swap(son[r][0],son[r][1]),swap(ls[r],rs[r]);
		rev[x]=0;
	}
}

build和newjd

没错…就是这么暴力的函数名

int newjd() {//寻找可用的新节点编号
	++trs;//splay的大小
	if(top) return st[top--];
	else return trs;
}
int build(int l,int r,int fnum,int fa,int *a) {
	if(l>r) return 0;
	int t=newjd(),mid=(l+r)/2;
	if(l!=r) build(l,mid-1,0,t,a),build(mid+1,r,1,t,a);
	v[t]=a[mid],son[fa][fnum]=t,same[t]=inf,rev[t]=0,f[t]=fa,up(t);
	return t;
}

由于空间限制，我们需要一个st栈，储存所有删除操作中被清理出来的空间，然后优先使用这些空间。newjd就是一个用来寻找这样的空间的函数。
因而我们一定要注意编号和a数组编号的不同为build函数造成的影响。

do系列

即进行各种操作用的函数。
我们需要两个哨兵节点，一个放在最前面，一个放在最后面，防止把这棵树删干净。
然后将需要操作的序列左端点前面的那个节点splay到根，右端点后面的那个节点splay到根的右节点，这样根的右节点的左节点那棵子树就是我们要操作的序列，直接对其操作即可。

完整代码

#include<bits/stdc++.h>
using namespace std;
int read() {
	int q=0,w=1;char ch=' ';
	while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
	if(ch=='-') w=-1,ch=getchar();
	while(ch>='0'&&ch<='9') q=q*10+ch-'0',ch=getchar();
	return w*q;
}
const int N=500005,inf=100000000;
int a[N],b[N],v[N],siz[N],son[N][2],f[N];
int st[N*10],same[N],rev[N],sum[N],mx[N],ls[N],rs[N];
//st:放置可用节点的栈,same:set标记,rev:翻转标记,sum:和
//mx:最大子列,ls:左边最大,rs:右边最大
int n,m,top,rot,trs;
void up(int x) {
	int l=son[x][0],r=son[x][1];
	sum[x]=sum[l]+sum[r]+v[x],siz[x]=siz[l]+siz[r]+1;
	ls[x]=max(ls[l],sum[l]+v[x]+ls[r]);
	rs[x]=max(rs[r],sum[r]+v[x]+rs[l]);
	mx[x]=max(max(mx[l],mx[r]),rs[l]+ls[r]+v[x]);
}
void pd(int x) {
	int l=son[x][0],r=son[x][1],t=same[x];
	if(same[x]!=inf) {
		if(l) {
			same[l]=v[l]=t,sum[l]=siz[l]*t;
			if(t>=0) ls[l]=rs[l]=mx[l]=sum[l];
			else ls[l]=rs[l]=0,mx[l]=t;
		}
		if(r) {
			same[r]=v[r]=t,sum[r]=siz[r]*t;
			if(t>=0) ls[r]=rs[r]=mx[r]=sum[r];
			else ls[r]=rs[r]=0,mx[r]=t;
		}
		same[x]=inf;
	}
	if(rev[x]) {
		if(l) rev[l]^=1,swap(son[l][0],son[l][1]),swap(ls[l],rs[l]);
		if(r) rev[r]^=1,swap(son[r][0],son[r][1]),swap(ls[r],rs[r]);
		rev[x]=0;
	}
}
int is(int x) {return son[f[x]][1]==x;}
void spin(int x,int &mb) {
	int fa=f[x],g=f[fa],t=is(x);
	if(mb==fa) mb=x;
	else son[g][is(fa)]=x;
	f[x]=g,f[fa]=x,f[son[x][t^1]]=fa;
	son[fa][t]=son[x][t^1],son[x][t^1]=fa;
	up(fa),up(x);
}
void splay(int x,int &mb) {
	while(x!=mb) {
		int fa=f[x],g=f[fa];
		if(fa!=mb) {
			if(is(fa)^is(x)) spin(x,mb);
			else spin(fa,mb);
		}
		spin(x,mb);
	}
}
int find(int x,int kth) {
	if(same[x]!=inf||rev[x]) pd(x);
	if(kth==siz[son[x][0]]+1) return x;
	else if(kth<=siz[son[x][0]]) return find(son[x][0],kth);
	else return find(son[x][1],kth-siz[son[x][0]]-1);
}
int newjd() {//寻找可用的新节点编号
	++trs;
	if(top) return st[top--];
	else return trs;
}
int build(int l,int r,int fnum,int fa,int *a) {
	if(l>r) return 0;
	int t=newjd(),mid=(l+r)/2;
	if(l!=r) build(l,mid-1,0,t,a),build(mid+1,r,1,t,a);
	v[t]=a[mid],son[fa][fnum]=t,same[t]=inf,rev[t]=0,f[t]=fa,up(t);
	return t;
}
void do1() {
	int k,num,x,y;
	k=read(),num=read();
	for(int i=1;i<=num;++i) b[i]=read();
	x=find(rot,k+1),y=find(rot,k+2);
	splay(x,rot),splay(y,son[rot][1]);
	son[y][0]=build(1,num,0,y,b);
	up(y),up(rot);
}
void del(int x) {
	if(!x) return;
	v[x]=siz[x]=f[x]=rev[x]=sum[x]=0,same[x]=inf,mx[x]=-inf;
	ls[x]=rs[x]=-inf,st[++top]=x;
	del(son[x][0]),del(son[x][1]);son[x][0]=son[x][1]=0;
}
void do2() {
	int l=read(),r=read()+l-1,x,y;
	x=find(rot,l),y=find(rot,r+2);
	splay(x,rot),splay(y,son[rot][1]);
	trs-=siz[son[y][0]],del(son[y][0]),son[y][0]=0;
	up(y),up(rot);
}
void do3() {
	int l=read(),r=read()+l-1,x,y,kl,tmp=read();
	x=find(rot,l),y=find(rot,r+2);
	splay(x,rot),splay(y,son[rot][1]),kl=son[y][0];
	v[kl]=same[kl]=tmp,sum[kl]=siz[kl]*tmp;
	if(tmp>=0) mx[kl]=ls[kl]=rs[kl]=sum[kl];
	else mx[kl]=tmp,ls[kl]=rs[kl]=0;
	up(y),up(rot);
}
void do4() {
	int l=read(),r=read()+l-1,x,y,kl;
	x=find(rot,l),y=find(rot,r+2);
	splay(x,rot),splay(y,son[rot][1]),kl=son[y][0];
	rev[kl]^=1,swap(son[kl][0],son[kl][1]),swap(ls[kl],rs[kl]);
}
int do5() {
	int l=read(),r=read()+l-1,x,y;
	x=find(rot,l),y=find(rot,r+2);
	splay(x,rot),splay(y,son[rot][1]);
	return sum[son[y][0]];
}
int do6() {
	int x=find(rot,1),y=find(rot,trs);
	splay(x,rot),splay(y,son[rot][1]);
	return mx[son[y][0]];
}
int main()
{
	char ch[20];
	n=read(),m=read();
	mx[0]=-inf,a[1]=-inf,a[n+2]=inf;
	for(int i=2;i<=n+1;++i) a[i]=read();
	rot=build(1,n+2,0,0,a);
	for(int i=1;i<=m;++i) {
		scanf("%s",ch);
		if(ch[0]=='I') do1();
		else if(ch[0]=='D') do2();
		else if(ch[0]=='M'&&ch[2]=='K') do3();
		else if(ch[0]=='R') do4();
		else if(ch[0]=='G') printf("%d\n",do5());
		else if(ch[0]=='M'&&ch[2]=='X') printf("%d\n",do6());
	}
    return 0;
}

无旋treap

和上面那货差不多，按照无旋treap的基本操作来就可以了。
splay写了3天，无旋treap只写了两个小时，看来，treap这东西虽然慢一些，但是还是吼啊！

#include<bits/stdc++.h>
using namespace std;
int read() {
	int q=0,w=1;char ch=' ';
	while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
	if(ch=='-') w=-1,ch=getchar();
	while(ch>='0'&&ch<='9') q=q*10+ch-'0',ch=getchar();
	return q*w;
}
#define mkp make_pair
#define pr pair<int,int>
const int N=500005,inf=0x3f3f3f3f;
int son[N][2],sz[N],pos[N],v[N],a[N],s[N];
int ls[N],rs[N],mx[N],rev[N],sum[N],laz[N];//s!
int n,m,top,rt,sss;
void up(int x) {
	int l=son[x][0],r=son[x][1];
	sz[x]=sz[l]+sz[r]+1,sum[x]=sum[l]+sum[r]+v[x];
	ls[x]=max(ls[l],sum[l]+v[x]+ls[r]);
	rs[x]=max(rs[r],sum[r]+v[x]+rs[l]);
	mx[x]=max(rs[l]+ls[r]+v[x],max(mx[l],mx[r]));
}
void pd(int x) {
	int l=son[x][0],r=son[x][1],t=laz[x];
	if(laz[x]!=inf) {
		if(l) {
			laz[l]=v[l]=t,sum[l]=sz[l]*t;
			if(t>=0) ls[l]=rs[l]=mx[l]=sum[l];
			else ls[l]=rs[l]=0,mx[l]=t;
		}
		if(r) {
			laz[r]=v[r]=t,sum[r]=sz[r]*t;
			if(t>=0) ls[r]=rs[r]=mx[r]=sum[r];
			else ls[r]=rs[r]=0,mx[r]=t;
		}
		laz[x]=inf;
	}
	if(rev[x]) {
		if(l) rev[l]^=1,swap(son[l][0],son[l][1]),swap(ls[l],rs[l]);
		if(r) rev[r]^=1,swap(son[r][0],son[r][1]),swap(ls[r],rs[r]);
		rev[x]=0;
	}
}
pr split(int x,int num) {
	if(!x) return mkp(0,0);
	pd(x);int l=son[x][0],r=son[x][1];
	if(sz[son[x][0]]==num) {son[x][0]=0,up(x);return mkp(l,x);}
	if(sz[son[x][0]]+1==num) {son[x][1]=0,up(x);return mkp(x,r);}
	if(sz[son[x][0]]>num) {
		pr tmp=split(son[x][0],num);
		son[x][0]=tmp.second,up(x);
		return mkp(tmp.first,x);
	}
	else {
		pr tmp=split(son[x][1],num-sz[son[x][0]]-1);
		son[x][1]=tmp.first,up(x);
		return mkp(x,tmp.second);
	}
}
int merge(int a,int b) {
	if(!a) return b;
	if(!b) return a;
	pd(a),pd(b);
	if(pos[a]<pos[b]) {son[a][1]=merge(son[a][1],b),up(a);return a;}
	else {son[b][0]=merge(a,son[b][0]),up(b);return b;}
}
int newjd(int num) {
	int x;
	if(!top) x=++sss;
	else x=s[top--];
	sz[x]=1,pos[x]=rand(),v[x]=sum[x]=mx[x]=num;
	son[x][0]=son[x][1]=0,laz[x]=inf,ls[x]=rs[x]=max(0,num);
	return x;
}
int build(int l,int r) {
	if(l==r) {int x=newjd(a[l]);return x;}
	int mid=(l+r)>>1;
	merge(build(l,mid),build(mid+1,r));
}
void del(int x) {
	if(son[x][0]) del(son[x][0]);
	if(son[x][1]) del(son[x][1]);
	if(top<N-5) s[++top]=x;
}
int main()
{
	char ch[10];int x,tt,num;
	mx[0]=-inf;
	n=read(),m=read(),srand(m);
	for(int i=1;i<=n;++i) a[i]=read();
	rt=build(1,n);
	while(m--) {
		scanf("%s",ch);
		if(ch[0]=='I') {
			x=read(),tt=read();
			pr tmp=split(rt,x);
			for(int i=1;i<=tt;++i) a[i]=read();
			rt=merge(tmp.first,merge(build(1,tt),tmp.second));
		}
		else if(ch[0]=='D') {
			x=read(),tt=read();
			pr t1=split(rt,x+tt-1),t2=split(t1.first,x-1);
			rt=merge(t2.first,t1.second),del(t2.second);
		}
		else if(ch[0]=='M'&&ch[2]=='K') {
			x=read(),tt=read(),num=read();
			pr t1=split(rt,x+tt-1),t2=split(t1.first,x-1);
			int kl=t2.second;
			laz[kl]=v[kl]=num,sum[kl]=num*sz[kl];
			if(num>=0) ls[kl]=rs[kl]=mx[kl]=sum[kl];
			else ls[kl]=rs[kl]=0,mx[kl]=num;
			rt=merge(t2.first,merge(t2.second,t1.second));
		}
		else if(ch[0]=='R') {
			x=read(),tt=read();
			pr t1=split(rt,x+tt-1),t2=split(t1.first,x-1);
			int kl=t2.second;
			rev[kl]^=1,swap(son[kl][0],son[kl][1]),swap(ls[kl],rs[kl]);
			rt=merge(t2.first,merge(t2.second,t1.second));
		}
		else if(ch[0]=='G') {
			x=read(),tt=read();
			pr t1=split(rt,x+tt-1),t2=split(t1.first,x-1);
			printf("%d\n",sum[t2.second]);
			rt=merge(t2.first,merge(t2.second,t1.second));
		}
		else if(ch[0]=='M') printf("%d\n",mx[rt]);
	}
	return 0;
}
litble
关注
专栏目录
非旋treap (fhq treap)
Josh_Shu_的博客
11-21 987
平衡树 非旋treap fhq treap
POJ 3481 Double Queue（无旋treap/Splay）
shadiao.的博客
11-09 300
题目传送门 解题思路 三种操作，插入一个数，取最大值并删除，取最小值并删除。 应该有挺多方法可以做的。 提供无旋treap和Splay这两种平衡树的数组的做法。 代码 无旋treap版本 #include<cstdio> #include<algorithm> using namespace std; #define ll long long #defi...
普通平衡树（Splay，无旋Treap）
Master.Yi的博客
10-21 425
#include <bits/stdc++.h> #define maxn 500005 using namespace std; int n,sz,root,fa[maxn],val[maxn],ch[maxn][2],siz[maxn],cnt[maxn]; inline bool isc(int x){return ch[fa[x]][1]==x;} inline void up...
无旋Treap
蒟蒻QTY
08-03 1333
想要搞带区间的平衡树，要么用splay，要么用无旋treap。（现在只会后面的）     先%一下大佬LadyLex     无旋treap简言之不再有左旋和右旋，而多了拆分子树和合并子树的过程。     slipt为拆树。D 为定义的pair，first是拆下来的树，second是拆下来后剩下的原树。     如果要拆下来的树的大小比左子树小，就拆左子树，并把拆完的树接回来。
非旋Treap总结 ： 快过Splay 好用过传统Treap
weixin_33690367的博客
12-28 457
非旋$Treap$ 其高级名字叫$Fhq\ Treap$，既然叫$Treap$，它一定满足了$Treap$的性质（虽然可能来看这篇的人一定知道$Treap$，但我还是多说几句：$Fhp\ Treap$就是继承了$Treap$的随机系统，在二叉搜索的基础上，每个点加一个随机化处理，这些随机值满足堆的性质……通俗一点讲，就是$Fhp\ Treap$它每个点有至少两个值，一个是val，即存的数值，这些...
BZOJ3585/洛谷P4137 区间mex（主席树）
hipamp的博客
06-06 428
题意 有 nnn 个数 {an}\{a_n\}{an​}，有 mmm 个询问，每个询问给定 l,rl,rl,r，求 al,...,ara_l,...,a_ral​,...,ar​ 的 mexmexmex。（最小的未出现的自然数） 其中 n,m≤2e5,ai≤109n,m\leq 2e5,a_i\leq10^9n,m≤2e5,ai​≤109 分析 令 viv_ivi​ 表示 iii 最后出现的位置。用主席树维护 viv_ivi​。 那么对于一个询问 (l,r)(l,r)(l,r)，我们在 [1,r][1,r]
BZOJ 1500 [NOI2005]维修数列 （无旋treap良心题解）
shadiao.的博客
09-10 277
题目传送门 上面链接沉了幸亏洛谷还有这道题 无旋treap解题的原理： 1.无旋treap可以在维护平衡的同时保持中序遍历不变（treap也可以，应该都可以），因此我们用下标代表插入二叉树的值，那么我们可以按照中序遍历获得原来的数组。 2.无旋treap的分裂操作可以分裂出包含前k个节点的树和剩下元素的树，就相当于可以把数组切开。 你可以切取想要的区间进行操作，操作完再merge回去即可...
BZOJ1500: [NOI2005]维修数列（洛谷P2042）
forezxl的博客
12-09 359
Splay
无旋treap/fhq-treap 模板 （附良心讲解 + 例题：洛谷P3369 普通平衡树）
shadiao.的博客
09-09 730
博文更新： 20190917 更新了一下merge和split函数的递归理解 定义： 无旋treap，就是不通过旋转的维护堆性质的treap，根本的原理还是 tree + heap。 与普通treap的不同之处： 大概就是可以可持久化。 然后区间操作的话，不知道有旋的行不行，用无旋的反正是比较好搞，区间操作举个例子就是让你把一段区间的数反过来等等。 现在听着可能区间操作什么的可能莫名...
树套树之线段树套线段树（BZOJ3110、洛谷P3332）
forezxl的博客
01-03 1165
线段树可以维护一个序列。当需要维护一个矩阵（即二维平面）内的数值，或者有什么奇怪的区间操作时，就需要用到二维线段树，也就是线段树“套上”线段树。
bzoj 1500
diyan5166的博客
03-08 128
再写一次毒瘤题，改变代码风格改用了指针，以为会不太适应，事实上不会撒。。 对拍数据要等20分钟才出来= = 1 //#include<bits/stdc++.h> 2 #include<cstdio> 3 #include<cstring> 4 #include<algorithm> 5 #inc...
bzoj1500
生活不易，多才多艺
03-11 232
我只是…发个博客玩【捂脸】 #include &lt;cstdio&gt; #include &lt;queue&gt; //#include &lt;algorithm&gt; //using namespace std; #define INF 1000000000 #define N 510000 queue&lt;int&gt; q; int n,m,root,tot=0; int ...
BZOJ 1500([NOI2005]维修数列-Splay的数列维护)
nike0good |Oier&ACMer | 熟能生巧
02-13 2024
1500: [NOI2005]维修数列 Time Limit: 10 Sec  Memory Limit: 64 MB Submit: 3087  Solved: 920 [Submit][Status][Discuss] Description Input 输入文件的第1行包含两个数N和M，N表示初始时数列中数的个数，M表示要进行的操作数目。 第2行包含N个数字，描述初始
【BZOJ1500】[NOI2005]维修数列【Splay】
BraketBN
03-16 435
http://www.lydsy.com/JudgeOnline/problem.php?id=1500 这神题终于过了... 主要是： find()那里少写了一个pushdown()。 辣鸡回收写跪了。 /* Footprints In The Blood Soaked Snow */ #include #include using namespace s
[BZOJ]1500 维修数列
Le Petit Prince
07-25 405
1500: [NOI2005]维修数列 Time Limit: 10 Sec  Memory Limit: 64 MB Submit: 14433  Solved: 4701 [Submit][Status][Discuss] Description Input 输入的第1 行包含两个数N 和M(M ≤20 000)，N 表示初始时数列中数的个数，M表示要进行的操作数目
【BZOJ】1500: [NOI2005]维修数列-splay
远行客
12-05 151
splay模板
洛谷 P1640 BZOJ 1854 [SCOI2010]连续攻击游戏
最新发布
05-14
这道题目需要进行博弈论的思考。我们可以使用动态规划来解决这道题目。 首先，我们可以定义 $f_{i,j}$ 表示当前剩下 $i$ 到 $j$ 这些数时，当前玩家最多能获得的分数。我们可以发现，在这个状态中，当前玩家只能从 $i$ 或者 $j$ 中选择一个数进行取走，然后将剩下的部分转化为对手的状态。因此，我们可以得到如下的转移方程： $$ f_{i,j}=\max\{a_i+\min\{f_{i+2,j},f_{i+1,j-1}\},a_j+\min\{f_{i+1,j-1},f_{i,j-2}\}\} $$ 其中 $a_i$ 表示第 $i$ 个数的分值。上述方程的含义是，当前玩家可以选择取走 $i$ 或者 $j$，如果取走 $i$，则转化为对手的状态是 $f_{i+2,j}$ 或者 $f_{i+1,j-1}$，如果取走 $j$，则转化为对手的状态是 $f_{i,j-2}$ 或者 $f_{i+1,j-1}$。当前玩家希望获得的分数是最大的，因此需要取两种情况中的较大值。 最终答案就是 $f_{1,n}$，表示从 $1$ 到 $n$ 这些数中，当前玩家最多能获得的分数。 最后，需要注意的是，我们需要使用逆序的方式来进行动态规划，即从小到大枚举区间长度 $len$，然后枚举区间起点 $i$，最终计算 $f_{1,n}$。这样可以保证在计算 $f_{i,j}$ 时，$f_{i+1,j}$ 和 $f_{i,j-1}$ 已经被计算出来了。
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值