Processing math: 100%

●BZOJ 4278 [ONTAK2015]Tasowanie

题链:

http://www.lydsy.com/JudgeOnline/problem.php?id=4278

题解:

把两个串连接起来,用一个大数连接(必须要用大数)。
倍增算法求出后缀排名。
然后两个指针指向A,B串,每次比较两个后缀的大小,
优先取出小的那个后缀的首字母,然后对应指针向后移动一位。

代码:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
#include<cstdio>
#include<cstring>
#include<iostream>
#define MAXN 450000
#define filein(x) freopen(#x".in","r",stdin);
#define fileout(x) freopen(#x".out","w",stdout);
using namespace std;
int sa[MAXN],rak[MAXN],hei[MAXN];
int A[MAXN],B[MAXN],T[MAXN];
int LA,LB;
void build(int N,int M){
    static int cc[MAXN],ta[MAXN],tb[MAXN],*x,*y,h,p;
    x=ta; y=tb; h=0; T[N]=-1;
    for(int i=0;i<M;i++) cc[i]=0;
    for(int i=0;i<N;i++) cc[x[i]=T[i]]++;
    for(int i=1;i<M;i++) cc[i]+=cc[i-1];
    for(int i=N-1;i>=0;i--) sa[--cc[x[i]]]=i;
    for(int k=1;p=0,k<N;k<<=1){
        for(int i=N-k;i<N;i++) y[p++]=i;
        for(int i=0;i<N;i++) if(sa[i]>=k) y[p++]=sa[i]-k;
        for(int i=0;i<M;i++) cc[i]=0;
        for(int i=0;i<N;i++) cc[x[y[i]]]++;
        for(int i=1;i<M;i++) cc[i]+=cc[i-1];
        for(int i=N-1;i>=0;i--) sa[--cc[x[y[i]]]]=y[i];
        swap(x,y); y[N]=-1; x[sa[0]]=0; M=1;
        for(int i=1;i<N;i++)
            x[sa[i]]=y[sa[i]]==y[sa[i-1]]&&y[sa[i]+k]==y[sa[i-1]+k]?M-1:M++;
        if(M>=N) break;
    }
    for(int i=0;i<N;i++) rak[sa[i]]=i;
    for(int i=0,j;i<N;i++){
        if(h) h--;
        if(rak[i]){
            j=sa[rak[i]-1];
            while(T[i+h]==T[j+h]) h++;
        }
        hei[rak[i]]=h;
    }
}
void solve(){
    int i=0,j=0,ii,jj,h;
    while(i<LA&&j<LB){
        ii=i; jj=LA+1+j;
        if(rak[ii]<rak[jj]) printf("%d ",A[i++]);
        else printf("%d ",B[j++]);
    }
    while(i<LA) printf("%d ",A[i++]);
    while(j<LB) printf("%d ",B[j++]);
}
int main()
{
    int N; scanf("%d",&LA);
    for(int i=0;i<LA;i++) scanf("%d",&A[i]),T[i]=A[i];
    T[LA]=1001; scanf("%d",&LB);
    for(int i=0;i<LB;i++) scanf("%d",&B[i]),T[LA+1+i]=B[i];
    N=LA+LB+1;
    build(N,1005);
    solve();
    return 0;
}


 

posted @   *ZJ  阅读(143)  评论(0编辑  收藏  举报
编辑推荐:
· 如何编写易于单元测试的代码
· 10年+ .NET Coder 心语,封装的思维:从隐藏、稳定开始理解其本质意义
· .NET Core 中如何实现缓存的预热?
· 从 HTTP 原因短语缺失研究 HTTP/2 和 HTTP/3 的设计差异
· AI与.NET技术实操系列:向量存储与相似性搜索在 .NET 中的实现
阅读排行:
· 10年+ .NET Coder 心语 ── 封装的思维:从隐藏、稳定开始理解其本质意义
· 地球OL攻略 —— 某应届生求职总结
· 周边上新:园子的第一款马克杯温暖上架
· Open-Sora 2.0 重磅开源!
· 提示词工程——AI应用必不可少的技术
点击右上角即可分享
微信分享提示