Description

给定两个数字串A和B,通过将A和B进行二路归并得到一个新的数字串T,请找到字典序最小的T。
 

Input

第一行包含一个正整数n(1<=n<=200000),表示A串的长度。
第二行包含n个正整数,其中第i个数表示A[i](1<=A[i]<=1000)。
第三行包含一个正整数m(1<=m<=200000),表示B串的长度。
第四行包含m个正整数,其中第i个数表示B[i](1<=B[i]<=1000)。
 

Output

输出一行,包含n+m个正整数,即字典序最小的T串。

题解:利用后缀数组比较字典序大小,同 BZOJ 1692

#include <bits/stdc++.h>
#define setIO(s) freopen(s".in", "r", stdin)
#define maxn 4000000
using namespace std;
int n, m, tot;
int arr[maxn], height[maxn], A[maxn];
namespace SA
{
int rk[maxn], tp[maxn], sa[maxn], tax[maxn];
void qsort()
{
for(int i = 0; i <= m ; ++i) tax[i] = 0;
for(int i = 1; i <= n ; ++i) ++tax[rk[i]];
for(int i = 1; i <= m ; ++i) tax[i] += tax[i - 1];
for(int i = n; i >= 1 ; --i) sa[tax[rk[tp[i]]]--] = tp[i];
}
void build()
{
for(int i = 1; i <= n ; ++i) rk[i] = arr[i], tp[i] = i;
qsort();
for(int k = 1; k <= n ; k <<= 1)
{
int p = 0;
for(int i = n - k + 1; i <= n ; ++i) tp[++p] = i;
for(int i = 1; i <= n ; ++i) if(sa[i] > k) tp[++p] = sa[i] - k;
qsort(), swap(rk, tp), rk[sa[1]] = p = 1;
for(int i = 2; i <= n ; ++i)
{
rk[sa[i]] = (tp[sa[i - 1]] == tp[sa[i]] && tp[sa[i - 1] + k] == tp[sa[i] + k]) ? p : ++p;
}
if(n == p) break;
m = p;
}
int k = 0;
for(int i = 1; i <= n ; ++i) rk[sa[i]] = i;
for(int i = 1; i <= n ; ++i)
{
if(k) --k;
int j = sa[rk[i] - 1];
while(arr[i + k] == arr[j + k]) ++k;
height[rk[i]] = k;
}
}
};
int main()
{
// setIO("input");
int a, b;
scanf("%d",&a);
for(int i = 1 ; i <= a; ++i) scanf("%d",&arr[++n]);
arr[++n] = 2000;
scanf("%d",&b);
for(int i = 1 ; i <= b; ++i) scanf("%d",&arr[++n]);
m = 3000;
SA::build();
int l = 1, r = a + 2;
for(int i = 1; i <= a + b ; ++i)
{
if(l > a) printf("%d ", arr[r++]);
else if(r > n) printf("%d ",arr[l++]);
else if(SA::rk[l] < SA::rk[r]) printf("%d ",arr[l++]);
else printf("%d ",arr[r++]);
}
return 0;
}

  

最新文章

  1. linux:vi 替换命令
  2. 深入理解 KVC\KVO 实现机制 — KVO
  3. 简明Vim练级攻略(转载)
  4. udev/raw/asmlib/多路径 配置asm
  5. 多线程程序设计学习(5)balking模式和timed模式
  6. 转载:Restore SQL Server database and overwrite existing database
  7. 在SCVMM2012R2中删除失去联系的VM GateWay
  8. uva 11168 - Airport
  9. Nape 获取碰撞点加特效
  10. Dedecms列表页标签list/pagelist使用方法及pagelist的样式
  11. stm32开发之使用Keil MDK以及标准外设库创建STM32工程
  12. Python库
  13. flutter 解析json
  14. python中的单向链表实现
  15. 初学node.js,安装nodemon,学习debug模式,安装cpu-stat
  16. struts2+mybatis3+log4j2+bulma+mysql8元工程
  17. 关于html2canvas跨域设置
  18. k8s namespace/volume
  19. ijkplayer之.so文件编译过程
  20. 【转帖】Systemd 入门教程:命令篇

热门文章

  1. 使用Visual Studio Code调试Electron主进程
  2. WinCE:在Win7上连接WinCE手持设备
  3. css3中 弹性盒模型布局之box-flex
  4. 为什么是kafka?
  5. @PropertySource&amp;@ImportResource&amp;@Bean
  6. job调度时间格式
  7. jQuery总结02
  8. 一段程序的人生 第10章: server
  9. Ubuntu17.10 下配置caffe 仅CPU i386可以直接apt install caffe-cpu,但是怎么运行mnist代码我懵逼了
  10. B3402 [Usaco2009 Open]Hide and Seek 捉迷藏 最短路