题目大意

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

题解

我们从归并排序的角度去想,每次把两者之一较小的取出来

遇到相等的元素的时候,排序时取那个都是一样的.

但是在这道题中我们是求字典序最小,这就不一样了.

我们不能任意取,这一步的决策其实是由后面的数字的大小决定的

怎么说呢...

方法蕴藏在这三张图里了,自己去悟吧。

黑色箭头表示选取,黑色圆圈表示确定选取序列的理由.

迭代第一次



迭代第二次



迭代第三次



剩下的略...

这个东西我们可以用后缀数组(lcp)来维护

Code

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
inline void read(int &x){
x=0;char ch;bool flag = false;
while(ch=getchar(),ch<'!');if(ch=='-')ch=getchar(),flag = true;
while(x=10*x+ch-'0',ch=getchar(),ch>'!');if(flag) x=-x;
}
const int maxn = 410010;
int a[maxn],b[maxn],c[maxn];
int wa[maxn],wb[maxn],ws[maxn];
int sa[maxn],rank[maxn],height[maxn];
inline bool cmp(int *r,int i,int j,int l){
return r[i] == r[j] && r[i+l] == r[j+l];
}
void da(int *r,int n,int m){
int i,j,p,*x = wa,*y = wb;
for(i=0;i<m;++i) ws[i] = 0;
for(i=0;i<n;++i) ws[x[i] = r[i]]++;
for(i=1;i<m;++i) ws[i] += ws[i-1];
for(i=n-1;i>=0;--i) sa[--ws[x[i]]] = i;
for(j=1,p=1;p<n;j<<=1,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) ws[i] = 0;
for(i=0;i<n;++i) ws[x[y[i]]]++;
for(i=1;i<m;++i) ws[i] += ws[i-1];
for(i=n-1;i>=0;--i) sa[--ws[x[y[i]]]] = y[i];
for(swap(x,y),p=1,i=1,x[sa[0]] = 0;i<n;++i)
x[sa[i]] = cmp(y,sa[i-1],sa[i],j) ? p-1 : p++;
}
}
void get_h(int *r,int n){
int i,j,k=0;for(i=1;i<=n;++i) rank[sa[i]] = i;
for(i=0;i<n;height[rank[i++]] = k)
for(k ? --k : 0,j = sa[rank[i]-1];r[i+k] == r[j+k];++k);
}
int loger[maxn],minn[maxn][22];
void pre(int *r,int n){
loger[1] = 0;
for(int i=2;i<=n;++i){
loger[i] = loger[i-1];
if( (1 << loger[i]+1 ) == i) ++loger[i];
}
for(int i=n;i>=1;--i){
minn[i][0] = r[i];
for(int j=1;(i + (1<<j) - 1) <= n;++j){
minn[i][j] = min(minn[i][j-1],minn[i+(1<<j-1)][j-1]);
}
}
}
int lcp(int s,int t){
s = rank[s];t = rank[t];
if(s+1 > t) swap(s,t);++s;
int k = loger[t-s+1];
return min(minn[s][k],minn[t-(1<<k)+1][k]);
}
int n,m,len;
inline bool order(int i,int j){
int k = lcp(i,j);
int x = (i+k) == n ? 0x7f7f7f7f : c[i+k];
int y = (j+k) == len ? 0x7f7f7f7f : c[j+k];
return x <= y;
}
int q[maxn];
int main(){
read(n);
for(int i=0;i<n;++i){
read(a[i]);
c[i] = a[i];
}c[n] = 1005;
read(m);
for(int i=0;i<m;++i){
read(b[i]);
c[n+i+1] = b[i];
}len = n+m+1;c[len] = 0;
da(c,len+1,1010);
get_h(c,len);
pre(height,len);
int i=0,j=0,k=0;
while(i < n && j < m){
int x = order(i,j+n+1);
if(x == 1) q[k++] = a[i++];
else q[k++] = b[j++];
}
while(i < n) q[k++] = a[i++];
while(j < m) q[k++] = b[j++];
for(i=0;i<k;++i) printf("%d ",q[i]);
getchar();getchar();
return 0;
}

bzoj 4278 Tasowanie 后缀数组+贪心的更多相关文章

  1. bzoj 4278 Tasowanie —— 后缀数组

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4278 每次取两个后缀中字典序较小的那个的首字符: 注意超出去的部分是 inf 而不是 0,因 ...

  2. BZOJ 4278: [ONTAK2015]Tasowanie 后缀数组 + 贪心 + 细节

    Code: #include <bits/stdc++.h> #define setIO(s) freopen(s".in", "r", stdin ...

  3. BZOJ4278 [ONTAK2015]Tasowanie[后缀数组+贪心]

    题目 求两数组归并后的数组最小字典序排列. 嘛,可能本人在贪心这块还是太弱了(或者说什么都弱),如果不知道是字符串题估计也想不起来用sa. 显然看得出归并时字典序小的那个数组先往里面加,这就是要比较两 ...

  4. BZOJ 4278 [ONTAK2015]Tasowanie (后缀数组)

    题目链接: https://www.lydsy.com/JudgeOnline/problem.php?id=4278 题解: 居然把后缀数组写成n^2的..我真厉害.. 想了无数种方法,最后发现就是 ...

  5. BZOJ 4278: [ONTAK2015]Tasowanie (后缀数组 / 二分+hash)

    直接归并,然后如果哪边的后缀字典序比较小就去哪边,然后就可以后缀数组 博客传送门- 但是本蒟蒻不会后缀数组 Upd:Upd:Upd:现在会了233.一道差不多的题:BZOJ 1692: [Usaco2 ...

  6. [BZOJ 1692] [Usaco2007 Dec] 队列变换 【后缀数组 + 贪心】

    ---恢复内容开始--- 题目链接:BZOJ - 1692 题目分析 首先,有个比较简单的贪心思路:如果当前剩余字符串的两端字母不同,就选取小的字母,这样显然是正确的. 然而若两端字母相同,我们怎么选 ...

  7. BZOJ 1692: [Usaco2007 Dec]队列变换 [后缀数组 贪心]

    1692: [Usaco2007 Dec]队列变换 Time Limit: 5 Sec  Memory Limit: 64 MBSubmit: 1383  Solved: 582[Submit][St ...

  8. 【BZOJ】1692 & 1640: [Usaco2007 Dec]队列变换(后缀数组+贪心)

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

  9. 【BZOJ-4278】Tasowanie 后缀数组 + 归并

    4278: [ONTAK2015]Tasowanie Time Limit: 10 Sec  Memory Limit: 256 MBSubmit: 164  Solved: 80[Submit][S ...

随机推荐

  1. Java编码辅助工具:Mapstruct—— Java对象转换框架

    项目开发中,业务分层会涉及不同类型的Bean之间需要相互转换,如PO与DTO之间,PO与VO之间等.手动编码setter/getter各个对应属性,会显得臃肿繁琐.通过Mapstruct框架可简单方便 ...

  2. centOS解决乱码问题

    问题描述:输入javac出现乱码,部分字符不能显示解决方法 echo 'export LANG=en_US.UTF-8' >> ~/.bashrc

  3. linux下proc里关于磁盘性能的参数(转)

    我们在磁盘写操作持续繁忙的服务器上曾经碰到一个特殊的性能问题.每隔 30 秒,服务器就会遇到磁盘写活动高峰,导致请求处理延迟非常大(超过3秒).后来上网查了一下资料,通过调整内核参数,将写活动的高峰分 ...

  4. 【Java并发编程实战】—–“J.U.C”:ReentrantLock之二lock方法分析

    前一篇博客简介了ReentrantLock的定义和与synchronized的差别,以下尾随LZ的笔记来扒扒ReentrantLock的lock方法.我们知道ReentrantLock有公平锁.非公平 ...

  5. 【BZOJ4173】数学 欧拉函数神题

    [BZOJ4173]数学 Description Input 输入文件的第一行输入两个正整数 . Output 如题 Sample Input 5 6 Sample Output 240 HINT N ...

  6. LRM-00109: could not open parameter file

    SQL>startup                                                                                       ...

  7. detached HEAD state

    1 detached HEAD state指的是什么 正常情况下,HEAD指向一个branch,而branch又指向一个commit. detached HEAD state指的是HEAD指针没有指向 ...

  8. php总结1 ——php简介、工作原理、运行环境、文件构成、语法结构、注释

    1.1 PHP  超文本预处理程序.实际就是制作网站的脚本程序 1.2 运行环境: wamp——windowns+apache+mySQL+php    常用于开发.学习和研究 lamp ——linu ...

  9. SPFA的两个优化:SLF与LLL

    先举出个例题:洛谷P3371 [模板]单源最短路径 一眼扫去:最短路径. spfa不接受反驳... 附上代码: #include<iostream> #include<algorit ...

  10. 列举你了解的Python较其他语言的优势

    1.简单易学 2.开发速度快 3.拥有最成熟的程序包资源库(第三方库)