后缀数组小结--poj2406，poj3693，poj2774，poj1743

最新推荐文章于 2018-09-17 21:07:00 发布
最新推荐文章于 2018-09-17 21:07:00 发布
阅读量2.6k 收藏
点赞数
分类专栏： acm_数据结构 文章标签： ini c
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/ascii991/article/details/7448261
版权

自己搞OI的时候没学（惭愧，都冲进NOI了居然还不会后缀数组！），现在搞ACM了，就学了学人家的论文，做了几题，小结一下。

感觉这些题有些模块比较重要：

首先是求后缀数组：

int cmp(int *r,int a,int b,int l)
{
    return r[a]==r[b] && r[a+l]==r[b+l];
}
void da(char *s,int *sa,int n,int m)
{
     int *x=wa,*y=wb,*t;
     int i,j,p;
     for (i=0;i<m;i++) w[i]=0;
     for (i=0;i<n;i++) w[x[i]=s[i]]++;
     for (i=1;i<m;i++) w[i]+=w[i-1];
     for (i=n-1;i>=0;i--) sa[--w[x[i]]]=i;
     for (j=1,p=1;p<n;j*=2,m=p)
     {
         for (p=0,i=n-j;i<n;i++) y[p++]=i;
         for (i=0;i<n;i++) if (sa[i]>=j) y[p++]=sa[i]-j;
         for (i=0;i<n;i++) v[i]=x[y[i]];
         for (i=0;i<m;i++) w[i]=0;
         for (i=0;i<n;i++) w[v[i]]++;
         for (i=1;i<m;i++) w[i]+=w[i-1];
         for (i=n-1;i>=0;i--) sa[--w[v[i]]]=y[i];
         for (t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1;i<n;i++)
             x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
     }
}


然后是求height数组

void calheight(char *s,int *sa,int n)
{
     int i,j,k=0;
     for (i=1;i<=n;i++) rk[sa[i]]=i;
     for (i=0;i<n;h[rk[i++]]=k)
         for (k?k--:0,j=sa[rk[i]-1];s[i+k]==s[j+k];k++);
}


然后是对height数组搞RMQ

void RMQ()
{
     int i,j;
     er[0]=1;
     for (i=1;i<20;i++)
         er[i]=er[i-1]*2;
     log2[0]=-1;
     for (i=1;i<=n;i++)
         log2[i]=(i&(i-1))?log2[i-1]:log2[i-1]+1;
     for (i=1;i<=n;i++)
         f[i][0]=h[i];
     for (j=1;j<20;j++)
         for (i=1;i+er[j]-1<=n;i++)
             f[i][j]=min(f[i][j-1],f[i+er[j-1]][j-1]);
}


然后是求公共前缀lcp

int lcp(int a,int b)
{
     int x=rk[a],y=rk[b];这一句很重要，我经常忘记
     if (x>y)
     {
             int t;
             t=x;x=y;y=t;
     }
     x++;
     int k=log2[y-x+1];
     return min(f[x][k],f[y-er[k]+1][k]);
}


当然还有根据height数组的分组思想，和二分答案。

poj 3693

给定一个字符串，求重复次数最多的连续重复子串。

枚举长度L，然后求长为L的子串最多能连续出现几次。其实就是求lcp(i,i+L)，不必像论文里写的向前向后匹配。令k=lcp(i,i+L)，t=i-(l-k%l)。如果t>=0 && k%l!=0，那么ans=k/l+2,否则就是k/l+1。

#include <iostream>
#include <cstring>
const int maxn=100003;
char s[maxn];
int wa[maxn],wb[maxn],w[maxn],v[maxn],sa[maxn],rk[maxn],h[maxn],er[20],log2[maxn],a[maxn];
int f[maxn][20];
int n,cc;
using namespace std;
int cmp(int *r,int a,int b,int l)
{
    return r[a]==r[b] && r[a+l]==r[b+l];
}
void da(char *s,int *sa,int n,int m)
{
     int *x=wa,*y=wb,*t;
     int i,j,p;
     for (i=0;i<m;i++) w[i]=0;
     for (i=0;i<n;i++) w[x[i]=s[i]]++;
     for (i=1;i<m;i++) w[i]+=w[i-1];
     for (i=n-1;i>=0;i--) sa[--w[x[i]]]=i;
     for (j=1,p=1;p<n;j*=2,m=p)
     {
         for (p=0,i=n-j;i<n;i++) y[p++]=i;
         for (i=0;i<n;i++) if (sa[i]>=j) y[p++]=sa[i]-j;
         for (i=0;i<n;i++) v[i]=x[y[i]];
         for (i=0;i<m;i++) w[i]=0;
         for (i=0;i<n;i++) w[v[i]]++;
         for (i=1;i<m;i++) w[i]+=w[i-1];
         for (i=n-1;i>=0;i--) sa[--w[v[i]]]=y[i];
         for (t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1;i<n;i++)
             x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
     }
}
void calheight(char *s,int *sa,int n)
{
     int i,j,k=0;
     for (i=1;i<=n;i++) rk[sa[i]]=i;
     for (i=0;i<n;h[rk[i++]]=k)
         for (k?k--:0,j=sa[rk[i]-1];s[i+k]==s[j+k];k++);
}
void RMQ()
{
     int i,j;
     er[0]=1;
     for (i=1;i<20;i++)
         er[i]=er[i-1]*2;
     log2[0]=-1;
     for (i=1;i<=n;i++)
         log2[i]=(i&(i-1))?log2[i-1]:log2[i-1]+1;
     for (i=1;i<=n;i++)
         f[i][0]=h[i];
     for (j=1;j<20;j++)
         for (i=1;i+er[j]-1<=n;i++)
             f[i][j]=min(f[i][j-1],f[i+er[j-1]][j-1]);
}
int lcp(int a,int b)
{
     int x=rk[a],y=rk[b];
     if (x>y)
     {
             int t;
             t=x;x=y;y=t;
     }
     x++;
     int k=log2[y-x+1];
     return min(f[x][k],f[y-er[k]+1][k]);
}
void work()
{
     int i,j,k,t,r,l,mm=0,tot=0;
     for (l=1;l<n;l++)
         for (i=0;i+l<n;i+=l)
         {
             k=lcp(i,i+l);
             r=k/l+1;
             t=i-(l-k%l);
             if (t>=0 && k%l!=0)
                if (lcp(t,t+l)>=k) r++;
             if (r>mm)
             {
                      tot=0;
                      a[++tot]=l;
                      mm=r;
             }
             else if (r==mm) {
                  a[++tot]=l; }        
         }
     int st=0,tl=n;
     for (i=1;i<n;++i)
         for (j=1;j<=tot;++j)
         {
             k=lcp(sa[i],sa[i]+a[j]);
             if (k>=(mm-1)*a[j])
             {
                 st=sa[i];
                 tl=mm*a[j];
                 i=n;
                 break;
             }
         }
     printf("Case %d: ",cc);
     for (i=st;i<st+tl;i++)
         printf("%c",s[i]);
     printf("\n");
}                      
int main()
{
    freopen("pin.txt","r",stdin);
    freopen("pou.txt","w",stdout);
    cc=0;
    while (scanf("%s",s)!=EOF && !(s[0]=='#' && strlen(s)==1))
    {
          cc++;
          n=strlen(s);
          s[n]=0;
          da(s,sa,n+1,128);
          calheight(s,sa,n);
          RMQ();
          work();
    }
    return 0;
}


poj2406

论文解释的很清楚，求height数组中每一个到height[rank[1]]之间的最小值即可。

记住：是rank[1]

#include <iostream>
#include <cstring>
const int maxn=1000003;
char s[maxn];
int wa[maxn],wb[maxn],w[maxn],v[maxn],sa[maxn],rk[maxn],rm[maxn],h[maxn];
int n;
using namespace std;
int cmp(int *r,int a,int b,int l)
{
    return r[a]==r[b] && r[a+l]==r[b+l];
}
void da(char *s,int *sa,int n,int m)
{
     int *x=wa,*y=wb,*t;
     int i,j,p;
     for (i=0;i<m;i++) w[i]=0;
     for (i=0;i<n;i++) w[x[i]=s[i]]++;
     for (i=1;i<m;i++) w[i]+=w[i-1];
     for (i=n-1;i>=0;i--) sa[--w[x[i]]]=i;
     for (j=1,p=1;p<n;j*=2,m=p)
     {
         for (p=0,i=n-j;i<n;i++) y[p++]=i;
         for (i=0;i<n;i++) if (sa[i]>=j) y[p++]=sa[i]-j;
         for (i=0;i<n;i++) v[i]=x[y[i]];
         for (i=0;i<m;i++) w[i]=0;
         for (i=0;i<n;i++) w[v[i]]++;
         for (i=1;i<m;i++) w[i]+=w[i-1];
         for (i=n-1;i>=0;i--) sa[--w[v[i]]]=y[i];
         for (t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1;i<n;i++)
             x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
     }
}
void calheight(char *s,int *sa,int n)
{
     int i,j,k=0;
     for (i=1;i<=n;i++) rk[sa[i]]=i;
     for (i=0;i<n;h[rk[i++]]=k)
         for (k?k--:0,j=sa[rk[i]-1];s[i+k]==s[j+k];k++);
}
void RMQ()
{
     int i,k=rk[0];
     rm[k]=1000000;
     for (i=k-1;i>=0;i--) 
         if (h[i+1]<rm[i+1]) rm[i]=h[i+1];
         else rm[i]=rm[i+1];
     for (i=k+1;i<=n;i++)
         if (h[i]<rm[i-1]) rm[i]=h[i];
         else rm[i]=rm[i-1];         
}
int work()
{
    int i;
    for (i=1;i<=n/2;i++)
    {
        if (n%i) continue;
        if (rm[rk[i]]==n-i) return n/i;
    }
    return 1;
} 
int main()
{
    freopen("pin.txt","r",stdin);
    freopen("pou.txt","w",stdout);
    while (scanf("%s",s)!=EOF && !(s[0]=='.' && strlen(s)==1))
    {
          n=strlen(s);
          s[n]=0;
          da(s,sa,n+1,128);
          calheight(s,sa,n);
          RMQ();
          printf("%d\n",work());
    }
    return 0;
}

         

poj2774

见论文

#include <iostream>
using namespace std;
const int maxn=200004;
int wa[maxn],wb[maxn],w[maxn],v[maxn],a[maxn],sa[maxn],rk[maxn],h[maxn];
char s[maxn],s1[maxn];
int n,n1;
int cmp(int *r,int a,int b,int l)
{
    return r[a]==r[b] && r[a+l]==r[b+l];
}
void da(char *s,int *sa,int n,int m)
{
     int *x=wa,*y=wb,*t,i,j,p;
     for (i=0;i<m;i++) w[i]=0;
     for (i=0;i<n;i++) w[x[i]=s[i]]++;
     for (i=1;i<m;i++) w[i]+=w[i-1];
     for (i=n-1;i>=0;i--) sa[--w[x[i]]]=i;
     for (j=1,p=1;p<n;j*=2,m=p)
     {
         for (p=0,i=n-j;i<n;i++) y[p++]=i;
         for (i=0;i<n;i++) if (sa[i]>=j) y[p++]=sa[i]-j;
         for (i=0;i<n;i++) v[i]=x[y[i]];
         for (i=0;i<m;i++) w[i]=0;
         for (i=0;i<n;i++) w[v[i]]++;
         for (i=1;i<m;i++) w[i]+=w[i-1];
         for (i=n-1;i>=0;i--) sa[--w[v[i]]]=y[i];
         for (t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1;i<n;i++)
             x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
     }
}
void calheight(char *s,int *sa,int n)
{
     int i,j,k=0;
     for (i=1;i<=n;i++) rk[sa[i]]=i;
     for (i=0;i<n;h[rk[i++]]=k)
         for (k?k--:0,j=sa[rk[i]-1];s[i+k]==s[j+k];k++);
}
int work()
{
    int i,ans=0;
    for (i=2;i<=n;i++)
        if (sa[i]<n1) a[i]=0;
        else a[i]=1;
    for (i=3;i<=n;i++)
        if (a[i]!=a[i-1] && h[i]>ans)
           ans=h[i];
    return ans;
}
int main()
{
    freopen("pin.txt","r",stdin);
    freopen("pou.txt","w",stdout);
    scanf("%s",s);
    scanf("%s",s1);
    n1=strlen(s);
    strcat(s,"$");
    strcat(s,s1);
    n=strlen(s);
    s[n]=0;
    da(s,sa,n+1,128);
    calheight(s,sa,n);
    printf("%d\n",work());
}


poj1743

就是应用二分答案和分组思想（我一开始很二叉的用数组记录了分组，表示就是循环的时候判断一下就可以了，详细可以见论文

#include <stdio.h>
using namespace std;
const int maxn=20005;
int wa[maxn],wb[maxn],sa[maxn],rk[maxn],s[maxn],w[maxn],h[maxn],v[maxn];
int n,ans;
void init()
{
    int i; 
    for (i=0;i<n;i++)
        scanf("%d",&s[i]);
    for (i=0;i<n-1;i++)
        s[i]=s[i+1]-s[i]+100;        
    s[--n]=0;
}        
int cmp(int* r,int a,int b,int l)
{
    return r[a]==r[b] && r[a+l]==r[b+l];
}
void da(int* s,int* sa,int n,int m)
{
     int* x=wa;int* y=wb;int* t;
     int i,j,p;
     for (i=0;i<m;i++) w[i]=0;
     for (i=0;i<n;i++) w[x[i]=s[i]]++;
     for (i=1;i<m;i++) w[i]+=w[i-1];
     for (i=n-1;i>=0;i--) sa[--w[x[i]]]=i;
     for (j=1,p=1;p<n;j*=2,m=p)
     {
         for (p=0,i=n-j;i<n;i++) y[p++]=i;
         for (i=0;i<n;i++) if (sa[i]>=j) y[p++]=sa[i]-j;
         for (i=0;i<m;i++) w[i]=0;
         for (i=0;i<n;i++) v[i]=x[y[i]];
         for (i=0;i<n;i++) w[v[i]]++;
         for (i=1;i<m;i++) w[i]+=w[i-1];
         for (i=n-1;i>=0;i--) sa[--w[v[i]]]=y[i];
         for (t=x,x=y,y=t,x[sa[0]]=0,p=1,i=1;i<n;i++)
             x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
     }
}
void calheight(int* s,int* sa,int n)
{
     int i,j,k=0;
     for (i=1;i<=n;i++)
         rk[sa[i]]=i;
     for (i=0;i<n;h[rk[i++]]=k)
         for (k?k--:0,j=sa[rk[i]-1];s[i+k]==s[j+k];k++);
}
int ok(int t)
{
    int smax,smin,i,j;
    smin=smax=sa[1];
    for (i=2;i<=n;i++)
    {
        if (h[i]>=t && i<n)
        {
           if (sa[i]>smax) smax=sa[i];
           if (sa[i]<smin) smin=sa[i];
           continue;
        }
        if (smax-smin>=t) return 1;
        smin=smax=sa[i];
    }
    return 0;
}     
void bin()
{
    int ll=4,rr=n,mid;
    while (ll<=rr)
    {
          mid=(ll+rr)/2;
          if (ok(mid)) ans=mid,ll=mid+1;
          else rr=mid-1;
    }
}  
int main()
{
    freopen("pin.txt","r",stdin);
    freopen("pou.txt","w",stdout);
    while (scanf("%d",&n),n)
    {
          init();
          da(s,sa,n+1,200);
          calheight(s,sa,n);
          bin();
          ans++;
          if (ans<5) ans=0;
          printf("%d\n",ans);
    }
    return 0;
}


ascii991
关注
专栏目录
leetcodeoj调试-My-POJ:我的POJ
07-01
《深入探索LeetCode OJ与My-POJ调试实践》 在编程的世界里，LeetCode OJ（在线判题系统）和My-POJ（个人POJ）是程序员提升技能、锻炼算法的重要平台。这两个工具提供了丰富的编程挑战，帮助开发者在解决实际问题中...
poj2814.rar_5-101_adg_limit_poj 28_poj2814
09-24
Time Limit: 1000ms Memory limit: 65536kB 题目描述 有9个时钟，排成一个3*3的矩阵。 现在需要用最少的移动，将9个时钟的指针都拨到12点的位置。共允许有9种不同的移动。如右表所示，每个移动会将若干个时钟的...
poj Musical Theme(后缀数组)
学习成神之路
09-17 214
前言： jz大佬说初二不会后缀数组就退役吧，然后我看了看高二的自己竟然苟过了三年…哈哈哈哈哈哈 题目大意： 求一个最长重复不相交的字符串，长度必须大于5 思路： 因为原题比较长，也没有上面说的那么裸，所以我们先用每个数剪掉前面一个数，然后就转化成上面那样了， 然后用后缀数据求出来就好了。 具体后缀数组应用： 我们先分析一个判定性问题，如何判断两个字符串是相同且不重叠的。解决这个问题要用到里面的he...
poj 3294(后缀数组)
俯瞰风景
10-18 609
题意：给出n个字符串，求出现次数超过n/2次及以上的最长子串，如果有多组按字典序输出，没有输出?。 题解：多个字符串的处理方式是把所有串连接起来，用没有出现在串中且字典序最大的字符分割，然后二分出一个长度x去判断是否有长度不小于x的公共前缀出现了n/2次以上。需要注意为了判断两个公共前缀不在一个串中，要把每个串的起始位置记录到sum数组，然后连续的height[i]>=x情况下，如果sa[i]所在
poj 2406 Power Strings （后缀数组）
clover_hxy的博客
12-28 736
Power Strings Time Limit: 3000MS   Memory Limit: 65536K Total Submissions: 45998   Accepted: 19230 Description Given two strings a and b we define a*b to be their concate
poj 2758(后缀数组)
俯瞰风景
10-18 1137
题意：给出一个长度为n的字母组成的序列，然后有两种操作，Q a b 询问以第a个字符为开头的后缀字符串和以第b个字符为开头的后缀字符串的lcp(最长公共前缀)的长度，I a b 表示在第b个字符之前插入字符a。 注意询问时的位置a和b是对应初始串，而比对的后缀串是当前串。 题解：思路想到了，但实现不出来 %>_<%，如果不考虑I操作，先用后缀数组处理出sa和rank，然后得到height相邻名次
【后缀数组】poj 3693
Starry_starry_night的博客
08-17 306
题目：给出一个串，求重复次数最多的连续重复子串枚举长度为L，然后看长度为L的字符串最多连续出现几次。既然长度为L的串重复出现，那么str[0],str[l],str[2*l]……中肯定有两个连续的出现在字符串中。 那么就枚举连续的两个，然后从这两个字符前后匹配，看最多能匹配多远。 即以str[i*l],str[i*l+l]前后匹配，这里是通过查询suffix(i*l),suffix(i*l+l)
poj_2682(3).rar_C++ 数组扩充_poj 26_poj 2682_poj26
09-19
"C++ 数组扩充"提示我们问题可能与如何在C++编程语言中处理数组的增长有关，而"poj 26_poj 2682_poj26"似乎是重复提及问题编号，可能是用户在整理文件时的习惯。 描述中提到的“数链思想”可能是指一种处理数组元素...
string-problem(POJ).rar_POJ 19_poj
09-21
标题中的"string-problem(POJ).rar_POJ 19_poj"表明这是一个与ACM编程竞赛相关的压缩包，特别关注的是字符串处理问题。在ACM编程竞赛中，字符串问题是一个常见的类别，通常涉及到字符串的查找、比较、操作、模式匹配...
poj-1028-Web-Navigation.zip_poj
09-21
【标题】"poj-1028-Web-Navigation.zip_poj" 指向的是一个编程竞赛问题，POJ（Programming Online Judge）平台上的第1028题，名为“Web Navigation”。该问题通常涉及到算法设计和实现，可能是为了考察参赛者的编程...
POJ2406之后缀数组
superxtong的博客
08-01 2783
E - Power StringsDescription Given two strings a and b we define a*b to be their concatenation. For example, if a = “abc” and b = “def” then a*b = “abcdef”. If we think of concatenation as multiplicat
POJ 1509 后缀数组
Dan__ge的博客
05-08 1065
点击打开链接 题意：给个串问最小的后缀，但是这个串是个环 思路：是个环也好处理，那么就将串复制一遍加到后面就行了，然后求出后缀数组和高度数组，刚开始我写的是就求出后缀数组sa，然后从排名第1开始找，比如找到了I，排名为I的从第sa[I]开始，那么如果这个值小于输入串的长度，直接输出就行，提交WA，看了看讨论，有这样的数组，aaaa那么应该输出1，而那个程序输出4，然后想到用高度数组，刚开始找到
poj-3693（后缀数组+RMQ）
Static_Ricardo的博客
02-23 805
poj-3693 Maximum repetitiion substring题目： The repetition number of a string is defined as the maximum number R such that the string can be partitioned into R same consecutive substrings. For example,
【POJ3450】后缀数组+二分
lajiyuan的博客
05-17 498
POJ3450 链接 题意为求多个串的最长公共子串，参考上题的做法，我们可以用特殊字符对所有字符串进行拼接，但是问题在于如果找到这个最长公共子串。数据量为n&amp;lt;4000len&amp;lt;200n&amp;lt;4000len&amp;lt;200n2−n2−n2-n遍历height数组，如果有某一段连续的区间heightheightheight数组均大于midmidmid，而且这段区间拥有来自nnn个串的后缀...
POJ - 3415 Common Substrings（后缀数组 + 单调栈）
Top_xiao的博客
08-18 367
题意： 求两个子串的大于等于k 的重复子串有多少个。 { 求两个子串的重复子串我们可以用单调栈来做， 然而，在同一个子串中，我们想要找重复子串的个数，那就是 height数组之和。 } 思路： 两个子串，我们用后缀数组维护，把两个子串合并在一起，然后中间加一个特殊字符。 {  看网上别人的博客。 一个n长度的子串，我们把后边填一个0，就是在位置n上加0， 这个字符串我们从 0 编...
poj 1988 并查集
AChunter的专栏
05-05 2358
【题意】 有几个stack，初始里面有一个cube。支持两种操作：1.move x y： 将x所在的stack移动到y所在stack的顶部。2.count x：数在x所在stack中，在x之下的cube的个数。 【题解】 并查集。与银河英雄传说类似。 维护两个量s[x]和d[x]，表示子树大小和到根的距离。 详情看代码： #include using namespace st
poj 1635 判断树是否同构
AChunter的专栏
04-26 2211
【题意】 给出两棵树，判断树是否同构。 树的表示：从根节点深搜，0表示向下，1表示向上回溯，得到的01串。 【题解】 将树最小表示，排序 如何得到子树：子树的01串中0和1的数量相等。 所以只需递归求解。 【代码】 #include #include #include #include using namespace std; const int maxn=3005;
poj 1625 AC自动机上的DP
AChunter的专栏
04-15 1642
好题！搞了N久终于AC！ 【题意】 给定N个字符集合和P个禁用的字符串，求长度为M的合法字符串个数。 【题解】 显然是AC自动机上的DP。 令dp[i][j]表示长度为i，状态为j的字符串个数。 按自动机的图转移即可，但要注意细节。 贡献几组数据。 50 50 10 qwertyuiop[]\asdfghjkl;'zxcvbnm,./ QWERTYUIOP{}|A
POJ题目标签与算法分类总结
- 递归和分治法：通过将问题分解为更小的子问题来求解，如POJ中的某些题目。 - 递推：利用已知条件推导出未知项，常见于数学问题的解决中。 - 构造法：直接构建满足条件的解，如POJ3295。 - 模拟法：按照实际...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值