后缀数组 + Hash + 二分 or Hash + 二分 + 双指针 求 LCP ---- 2017icpc 青岛 J Suffix (假题!!)_j.suffix 二分-程序员宅基地

技术标签: 算法  c语言  # 后缀数组  哈希算法  

题目链接


题目大意:

就是给你n个串每个串取一个后缀,要求把串拼起来要求字典序最小!!
s u m _ l e n g t h _ o f _ n ≤ 5 e 5 sum\_length\_of\_n\leq 5e5 sum_length_of_n5e5


MY Slove :

  1. 首先我们知道对于最后一个串肯定是取最小后缀的

  2. 那么我们可以把最后一个串的结果接到倒数第二个串上面在求一次最小后缀就可以得倒数第二个串的结果,依次以此类推就好了。

  3. 但是这样我们算算复杂度:假如我们每次都把答案拼到下一个串上面求个后缀数组那么我们假设最坏就是 a + z z z z z . . . . . . a+zzzzz...... a+zzzzz......取均摊的复杂度就是 n \sqrt n n 那么就是每个串长度是 n \sqrt n n ,有 n \sqrt n n 个串,答案最坏的情况是全部取,那么你每次都拼接的话,最后一个串要动 n \sqrt n n 次,倒数第二个要动 n − 1 \sqrt n-1 n 1次…那么最坏就是 O ( n ( 1 + n ) / 2 + n ( ∑ i = 1 n i ∗ l o g ( i ∗ n ) ) ) O(\sqrt n(1+\sqrt n)/2+ \sqrt n(\sum_{i=1}^{\sqrt n}i*log(i*\sqrt n))) O(n (1+n )/2+n (i=1n ilog(in ))) 铁定T飞了(但是实践证明 O ( n 2 ) O(n^2) O(n2)都可以过的假题)

  4. 我们考虑优化:

  5. 我们想想后面那一大坨的问题就是:后缀数组求解拼接后字符串的时间复杂度是很大的 n ( ∑ i = 1 n i ∗ l o g ( i ∗ n ) ) \sqrt n(\sum_{i=1}^{\sqrt n}i*log(i*\sqrt n)) n (i=1n ilog(in )) 那么我们可不可以考虑只对原来的字符求后缀数组,去虚假拼接求呢?

  6. 我们看看后缀排序后的后缀长什么样子:
    eg:aba
    a
    aba
    b

  7. 通过观察我们知道要拼接的答案肯定是以最小后缀为前缀的后缀的其中一个这个贪心看肯定是这样的
    我们看看这个数据就知道了

3
ababac
b
drl
ans:ababacbdrl

那么我们可以顺着这个串的所有后缀比下去,但是这个串所有后缀字符个数是 O ( n 2 ) O(n^2) O(n2)的?肯定不是这样比呀!!
在这里插入图片描述
因为我们知道前面是一样的,实际上我们每次只比较后个串多出来的那部分,那部分的长度是 O ( n ) O(n) O(n)

  1. 但是如果前面都一样怎么办?我们就要看红色重叠的部分了,但是如果暴力比较红色部分复杂度还是将不下来!!因为最坏就是红色的长度了那么有 n n n个后缀就是 O ( n 2 ) O(n^2) O(n2)
  2. 因为后面是答案串我们可以对答案串进行hash,然后二分那个第一个hash值第一个hash值不一样的位置进行判断字典序!
  3. 我们在动态维护答案串的hash,不能正序hash,因为字符是头插入的,那么我们逆着hash

AC code

细节巨多,局难写,写了一晚上

#include <bits/stdc++.h>
#define mid ((l + r) >> 1) 
#define Lson rt << 1, l , mid
#define Rson rt << 1|1, mid + 1, r
#define ms(a,al) memset(a,al,sizeof(a))
#define log2(a) log(a)/log(2)
#define _for(i,a,b) for( int i = (a); i < (b); ++i)
#define _rep(i,a,b) for( int i = (a); i <= (b); ++i)
#define for_(i,a,b) for( int i = (a); i >= (b); -- i)
#define rep_(i,a,b) for( int i = (a); i > (b); -- i)
#define lowbit(x) ((-x) & x)
#define IOS std::ios::sync_with_stdio(0); cin.tie(0); cout.tie(0)
#define INF 0x3f3f3f3f
#define LLF 0x3f3f3f3f3f3f3f3f
#define hash Hash
#define next Next
#define pb push_back
#define f first
#define s second
using namespace std;
const int N = 1e6 + 10, mod = 1e9 + 9;
const int maxn = 5e5+10;
// const long double eps = 1e-5;
const int base = 233;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> PII;
typedef pair<ll,ll> PLL;
typedef pair<double,double> PDD;
template<typename T> void read(T &x) {
    
    x = 0;char ch = getchar();ll f = 1;
    while(!isdigit(ch)){
    if(ch == '-')f*=-1;ch=getchar();}
    while(isdigit(ch)){
    x = x*10+ch-48;ch=getchar();}x*=f;
}
template<typename T, typename... Args> void read(T &first, Args& ... args)  {
    
    read(first);
    read(args...);
}
//......................这里是后缀数组
// sa[l]排序是lth的后缀的开始位置
//ra[l]是起点是l的后缀排名多少
//lcp(suf(i),suf(j)) = min(H(ra[i] + 1),....H(ra[j]));
//区间最小值倍增求
// H(i)是rk[i]和rk[i-1]的lcp
//这个板子下标是从0开始,rank从1开始
//传进来的是int数组但是不能有0 
struct SA {
     
    int sa[maxn], ra[maxn], height[maxn];
    int t1[maxn], t2[maxn], c[maxn];
    int shift[maxn];
    inline void init(const string &s) {
    
        for(int i = 0; i < s.size(); ++ i) shift[i] = (int)(s[i]-'a'+1);
        build(shift,s.size(),30);
    }
    void build(int *str, int n, int m) {
    //字符数组,长度,字符种类 
        str[n] = 0;
        n++;
        int i, j, p, *x = t1, *y = t2;
        for (i = 0; i < m; i++) c[i] = 0;//排序的桶
        for (i = 0; i < n; i++) c[x[i] = str[i]]++;
        for (i = 1; i < m; i++) c[i] += c[i - 1];
        for (i = n - 1; i >= 0; i--) sa[--c[x[i]]] = i;
        for (j = 1; j <= n; j <<= 1) {
    //所有长度为1,2,4,8....的子串的排序
       //长度够长的话就是所有的后缀排序
            p = 0;
            for (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++) c[i] = 0;
            for (i = 0; i < n; i++) c[x[y[i]]]++;
            for (i = 1; i < m; i++) c[i] += c[i - 1];
            for (i = n - 1; i >= 0; i--) sa[--c[x[y[i]]]] = y[i];
            swap(x, y);
            p = 1;
            x[sa[0]] = 0;
            for (i = 1; i < n; i++)
                x[sa[i]] = (y[sa[i - 1]] == y[sa[i]] && y[sa[i - 1] + j] == y[sa[i] + j]) ? p - 1 : p++;
            if (p >= n) break;
            m = p;
        }
        n--;
        for (int i = 0; i <= n; i++) ra[sa[i]] = i;
        for (int i = 0, j = 0, k = 0; i <= n; i++) {
    
            if (k) k--;
            j = sa[ra[i] - 1];
            while (str[i + k] == str[j + k]) k++;
            height[ra[i]] = k;
        }
        st_init(height, n);
    }
    int lg[maxn], table[23][maxn];
    void st_init(int *arr, int n) {
    
        if (!lg[0]) {
    
            lg[0] = -1;
            for (int i = 1; i < maxn; i++)
                lg[i] = lg[i / 2] + 1;
        }
        for (int i = 1; i <= n; ++i)
            table[0][i] = arr[i];
        for (int i = 1; i <= lg[n]; ++i)
            for (int j = 1; j <= n; ++j)
                if (j + (1 << i) - 1 <= n)
                    table[i][j] = min(table[i - 1][j], table[i - 1][j + (1 << (i - 1))]);
                    //从j开始长度是2^i次方得st表
    }
    int lcp(int l, int r) {
    
        l = ra[l], r = ra[r];
        if (l > r) swap(l, r);
        ++l;
        int t = lg[r - l + 1];//区间长度
        return min(table[t][l], table[t][r - (1 << t) + 1]);
    }
} sa;
//............................... 这里是Hash
string s[maxn], ans;
ull pw[maxn];
vector<ull> Hash;

inline void insert(char a) {
    
    ull last = Hash.empty() ? 0ull : Hash.back();
    last = last * base + a - 'a' + 1;
    Hash.push_back(last);
}

inline ull gethash(int l,int r) {
    
    if(r < l) return 0; 
    return (ull)Hash[r]-(l - 1 < 0 ? 0ull : Hash[l-1])*pw[r-l+1];
}
//........................
inline void update(int i, int anspos) {
     // 更新答案
    for(int m = s[i].size()-ans.size()-1; m >= anspos; -- m) insert(s[i][m]), ans += s[i][m];
}

int main() {
    
    IOS;
    pw[0] = 1;
    for(int i = 1; i < maxn; ++ i) pw[i] = pw[i-1] * base;
    int _;
    cin >> _;
    while (_--) {
    
        Hash.clear();
        int n;
        cin >> n;
        for(int i = 1; i <= n; ++ i) cin >> s[i];
        sa.init(s[n]);
        ans = s[n].substr(sa.sa[1]); // 先找到最后一个串的最小后缀
        
        for(int i = ans.size()-1; i >= 0; -- i) insert(ans[i]); // 逆着Hash
        reverse(ans.begin(),ans.end()); // 将答案串反过来
        
        for(int i = n - 1; i >= 1; -- i) {
     // 枚举第i个串
            if(s[i].size()==1) {
     insert(s[i][0]); ans += s[i]; continue;}// 特判长度为1的串
            // 对这个串建立后缀数组
            sa.init(s[i]);
            // 把答案串拼到新串后面
            for(int j = (int)ans.size()-1; j >= 0; -- j) s[i] += ans[j];
            
            int anspos = sa.sa[1]; //anspos = 答案后缀的开头位置
            for(int k = 2; k <= s[i].size()-ans.size(); ++ k) {
     // 枚举这个串的每个后缀
                if(sa.lcp(sa.sa[k],sa.sa[k-1]) != s[i].size()-ans.size()-anspos) {
     // 如果不是以前面后缀为前缀那就结束了后面答案不会更优
                    update(i,anspos);
                    break; 
                } 
                int eps = s[i].size() - ans.size() - sa.sa[k-1];
                int is = 0;
                // 只比较突出的部分
                int j;
                for(j = 0; j+anspos+eps < s[i].size() && j+sa.sa[k]+eps<s[i].size()-ans.size(); ++ j) 
                   if(s[i][j+anspos+eps] < s[i][j+sa.sa[k]+eps]) {
    // 如果在黑
                      is = 1;
                      update(i,anspos);
                      break;
                   } else if(s[i][j+anspos+eps] > s[i][j+sa.sa[k]+eps]) {
    
                      anspos = sa.sa[k];
                      if(k == s[i].size()-ans.size()) update(i,anspos);//如果到了最后一个串了可以直接更新了
                      is = 2;
                      break;
                   } else if(j+anspos+eps == s[i].size()-1) {
    
                      is = 1;
                      update(i,anspos);
                      break;                       
                   }
                
                if(is == 1) break; // is == 1 表示可以确定答案了提前退出就好了
                if(is == 2) continue;
                //...................................二分红色部分
                int l = 0, r = s[i].size()-(j+anspos+eps)-1;
                //注意我是把答案串拼到了新串后面
                int tmp = j+anspos+eps-(s[i].size()-ans.size());
                while(l < r) {
    //注意Hash是逆序存的,要从末尾开始看
                    if(gethash(Hash.size()-1-tmp-mid,Hash.size()-1-tmp)==gethash(Hash.size()-1-mid,Hash.size()-1)) l = mid+1;
                    else r = mid;
                }
                if(s[i][j+anspos+eps+l]>s[i][j+sa.sa[k]+eps+l]) anspos = sa.sa[k];
                if(k == s[i].size()-ans.size()) update(i,anspos);//如果到了最后一个串了可以直接更新了
            } 
        }
        for(int i = ans.size()-1;~i;i--) cout << ans[i];
        cout << endl;        
    }
    return 0;
}


Solution 2

贪心从后往前看,最后一个串一定选择字典序最小的后缀,然后把这个后缀拼接到第n - 1个串,重复这个步骤就行了。

具体实现:

从后往前遍历,每次找当前串的最小后缀,这个可以对于当前下标和当前最小后缀下标二分+hash找到lcp,看lcp后一个位置的大小,然后更新最小后缀的下标这里是关键)。然后把最小后缀拼接到前一个串即可。

理论时间复杂度跟后缀数组差不多,但后缀数组常数和编码量更大。


AC code

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef  unsigned long long ull;
const ull base=131;
#define mod
#define maxn 605000
ull p[maxn];ull h[maxn];
ull get_hash(ull l,ull mid)
{
    
    return h[l]-h[l-mid]*p[mid];
}
string str[maxn];char ans[maxn];
int len[maxn];
ull check(ull x,ull y)
{
    
    ull l=1,r=y;
    while(l<=r)
    {
    
 
        ull mid=(l+r)/2;
        //cout<<"l="<<l<<" r="<<r<<" mid="<<mid<<'\n';
        if(get_hash(x,mid)==get_hash(y,mid))
        {
    
            l=mid+1;
        }
        else r=mid-1;
    }
    if(l==y+1) return 0;
    return ans[x-r]<ans[y-r];
}
int main()
{
    
    ull T;
    p[0]=1;
    for(ull i=1;i<=maxn-100;i++)
    {
    
        p[i]=p[i-1]*base;
    }
    cin>>T;
    while(T--)
    {
    
        ull n;
        cin>>n;
        for(ull i=1;i<=n;i++)
        {
    
            cin>>str[i];
            len[i]=str[i].size();
        }
        ull id=0;
        for(ull i=n;i>=1;i--)
        {
    
            //cout<<"i="<<i<<" id="<<id<<'\n';
           ans[++id]=str[i][len[i]-1];
           h[id]=h[id-1]*base+str[i][len[i]-1];
           ull to=1;
           ull t=id;
           for(int j=(int)(len[i])-2;j>=0;j--)
           {
    
               //cout<<"i="<<i<<" j="<<j<<'\n';
               ans[id+to]=str[i][j];
               h[id+to]=h[id+to-1]*base+str[i][j];
               if(check(id+to,t))
               {
    
                   t=id+to;
               }
               to++;
           }
           id=t;
        }
        for(ull i=id;i>=1;i--) cout<<ans[i];
        cout<<'\n';
    }
}
/*
3
3
bbb
aaa
ccc
*/
版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/weixin_45630722/article/details/121007642

智能推荐

彻底扒光 通过智能路由器拆解看其本质-程序员宅基地

文章浏览阅读1.7k次。可以看到很多联发科的MT芯片摘自:https://net.zol.com.cn/531/5312999.html彻底扒光 通过智能路由器拆解看其本质2015-07-23 00:40:00[中关村在线 原创] 作者:陈赫|责编:白宁收藏文章 分享到 评论(24)关注智能路由器拆解的朋友们注意啦!我们已经将这五款产品彻底扒开,将主板的真容展现在了大家的眼前。网友们可以看见这些智能路由器主板的做工和用料,我们还为网友们展示了主要的电子元器件,供大家品评观赏。..._路由器拆解

Java--深入JDK和hotspot底层源码剖析Thread的run()、start()方法执行过程_jdk的源码hotspot跟jdk是分开的-程序员宅基地

文章浏览阅读2.1k次,点赞101次,收藏78次。【学习背景】今天主要是来了解Java线程Thread中的run()、start()两个方法的执行有哪些区别,会给出一个简单的测试代码样例,快速理解两者的区别,再从源码层面去追溯start()底层是如何最终调用Thread#run()方法的,个人觉得这样的学习不论对面试,还是实际编程来说都是比较有帮助的。进入正文~学习目录一、代码测试二、源码分析2.1 run()方法2.2 start()方法三、使用总结一、代码测试执行Thread的run()、start()方法的测试代码如下:public_jdk的源码hotspot跟jdk是分开的

透视俄乌网络战之一:数据擦除软件_俄乌网络战观察(一)-程序员宅基地

文章浏览阅读4.4k次,点赞90次,收藏85次。俄乌冲突中,各方势力通过数据擦除恶意软件破坏关键信息基础设施计算机的数据,达到深度致瘫的效果,同时窃取重要敏感信息。_俄乌网络战观察(一)

Maven私服仓库配置-Nexus详解_nexus maven-程序员宅基地

文章浏览阅读1.7w次,点赞23次,收藏139次。Maven 私服是一种特殊的Maven远程仓库,它是架设在局域网内的仓库服务,用来代理位于外部的远程仓库(中央仓库、其他远程公共仓库)。当然也并不是说私服只能建立在局域网,也有很多公司会直接把私服部署到公网,具体还是得看公司业务的性质是否是保密的等等,因为局域网的话只能在公司用,部署到公网的话员工在家里也可以办公使用。_nexus maven

基于AI的计算机视觉识别在Java项目中的使用 (四) —— 准备训练数据_java ocr ai识别训练-程序员宅基地

文章浏览阅读934次。我先用所有的样本数据对模型做几轮初步训练,让深度神经模型基本拟合(数万条记录的训练集,识别率到99%左右),具备初步的识别能力,这时的模型就是“直男”。相较于训练很多轮、拟合程度很高的“油腻男”,它的拟合程度较低,还是“直男愣头青”。..............._java ocr ai识别训练

hibernate 数据库类型 date没有时分秒解决_hibernate解析時間只有年月日沒有時分秒-程序员宅基地

文章浏览阅读688次。一、问题现象:  在数据库表中日期字段中存的日期光有年月日,没有时分秒。二、产生原因:三 解决办法   检查表的相应映射xml文件。 <property name="operateDate" type="Date">如果同上面所写,那问题出在 type类型上了正确写法 :<property name="operateDate" type="java.util..._hibernate解析時間只有年月日沒有時分秒

随便推点

springbbot运行无法编译成功,找不到jar包报错:Error:(3, 46) java: 程序包org.springframework.context.annotation不存在-程序员宅基地

文章浏览阅读1k次,点赞2次,收藏2次。文章目录问题描述:解决方案:问题描述:提示:idea springbbot运行无法编译成功,找不到jar包报错E:\ideaProject\demokkkk\src\main\java\com\example\demo\config\WebSocketConfig.javaError:(3, 46) java: 程序包org.springframework.context.annotation不存在Error:(4, 46) java: 程序包org.springframework.conte_error:(3, 46) java: 程序包org.springframework.context.annotation不存在

react常见面试题_recate面试-程序员宅基地

文章浏览阅读6.4k次,点赞6次,收藏36次。1、redux中间件中间件提供第三方插件的模式,自定义拦截 action -&gt; reducer 的过程。变为 action -&gt; middlewares -&gt; reducer 。这种机制可以让我们改变数据流,实现如异步 action ,action 过滤,日志输出,异常报告等功能。常见的中间件:redux-logger:提供日志输出redux-thunk:处理异步操作..._recate面试

交叉编译jpeglib遇到的问题-程序员宅基地

文章浏览阅读405次。由于要在开发板中加载libjpeg,不能使用gcc编译的库文件给以使用,需要自己配置使用另外的编译器编译该库文件。/usr/bin/ld:.libs/jaricom.o:RelocationsingenericELF(EM:40)/usr/bin/ld:.libs/jaricom.o:RelocationsingenericELF(EM:40)...._jpeg_utils.lo: relocations in generic elf (em: 8) error adding symbols: file

【办公类-22-06】周计划系列(1)“信息窗” (2024年调整版本)-程序员宅基地

文章浏览阅读578次,点赞10次,收藏17次。【办公类-22-06】周计划系列(1)“信息窗” (2024年调整版本)

SEO优化_百度seo resetful-程序员宅基地

文章浏览阅读309次。SEO全称为Search Engine Optimization,中文解释为搜索引擎优化。一般指通过对网站内部调整优化及站外优化,使网站满足搜索引擎收录排名需求,在搜索引擎中提高关键词排名,从而把精准..._百度seo resetful

回归预测 | Matlab实现HPO-ELM猎食者算法优化极限学习机的数据回归预测_猎食者优化算法-程序员宅基地

文章浏览阅读438次。回归预测 | Matlab实现HPO-ELM猎食者算法优化极限学习机的数据回归预测_猎食者优化算法

推荐文章

热门文章

相关标签