后缀数组 + Hash + 二分 or Hash + 二分 + 双指针 求 LCP ---- 2017icpc 青岛 J Suffix (假题!!)

news/2024/7/2 23:03:01

题目链接


题目大意:

就是给你n个串每个串取一个后缀,要求把串拼起来要求字典序最小!!
sum_length_of_n≤5e5sum\_length\_of\_n\leq 5e5sum_length_of_n5e5


MY Slove :

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

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

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

  4. 我们考虑优化:

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

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

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

3
ababac
b
drl
ans:ababacbdrl

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

  1. 但是如果前面都一样怎么办?我们就要看红色重叠的部分了,但是如果暴力比较红色部分复杂度还是将不下来!!因为最坏就是红色的长度了那么有nnn个后缀就是O(n2)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]); // 逆着Hashreverse(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
*/

http://lihuaxi.xjx100.cn/news/285925.html

相关文章

第2关:实现一个链接存储的栈

#if !defined(LINKED_STACK_H_985552) #define LINKED_STACK_H_985552 typedef int T; //数据元素类型 struct LNode {T data;LNode* next; };struct LinkStack {LNode* top; // 栈顶指针int len; // 栈的长度 };LinkStack* LS_Create(); void LS_Free(LinkStack* ls); void LS…

最常用 150 个Linux命令汇总(建议收藏)

点击上方[视学算法]→右上角[...]→[设为星标⭐]来源&#xff1a;banana 童www.cnblogs.com/bananaaa/p/7774467.htmllinux 命令是对 Linux 系统进行管理的命令。对于 Linux 系统来说&#xff0c;无论是中央处理器、内存、磁盘驱动器、键盘、鼠标&#xff0c;还是用户等都是文件…

牛客小白月赛25 补题+题解[A-J]

加油加油加油&#xff01; 文章目录A.AOE还是单体&#xff1f;B.k-size字符串C.白魔法师D.抽卡E.点击消除F.疯狂的自我检索者G.解方程H.神奇的字母&#xff08;二&#xff09;I.十字爆破J.异或和之和A.AOE还是单体&#xff1f; 思路&#xff1a;这题数据范围2e5&#xff0c;如…

Python | 一万多条拼车数据,看春运的迁徙图

作者 | 白苏&#xff0c;医疗健康领域产品经理一枚&#xff0c;Python&R爱好者来源 | InThirty编辑 | Jane今天是腊月二十八&#xff0c;你们都到家了吗&#xff1f;这篇文章&#xff0c;作者对北京、上海、广州、深圳、杭州等地 1万多条出行数据进行分析&#xff0c;得出了…

Unity3d多线程

为什么80%的码农都做不了架构师&#xff1f;>>> &#xff08;一&#xff09;多线程的创建 Thread t new Thread(new ThreadStart(Go)); Thread t1 new Thread(Go); 两种创建方式没有区别&#xff1b; &#xff08;二&#xff09;多线程的状态控制和优先级 多线程…

第1关:利用栈实现整数的十进制转八进制

#ifndef stack__h #define stack__h#include <stdio.h> #include <stdlib.h>typedef int T; // 数据元素的数据类型struct Stack{T* data; // 数据元素存储空间的开始地址int top; // 栈顶的位置int max; // 栈的最大长度 };Stack* Stack_Create(int maxlen)…

数据可视化[python-pyecharts]制作中国各省份近三个月新型冠状病毒肺炎变化图

大体思路&#xff1a; 通过pyecharts等库一个for循环批量绘制近几个月每天的图&#xff0c;最后通过pr将图片合成 先看一下某一天的样图&#xff0c;用pr组合起来之后就是个动态的了 文章目录安装pyecharts库数据来源代码部分1.导入库2.将路径中所有文件找出保存至列表3.处理导…

树上启发式合并问题 ---- 2019icpc南昌 K. Tree (树上启发式合并 + 动态开点线段树)

题目链接 题目大意&#xff1a; 就是给你一颗树&#xff0c;每个点有个权值viv_ivi​&#xff0c;问你有多少对(x,y)(x,y)(x,y)满足&#xff1a; xxx不是yyy的祖先yyy也不是xxx的祖先xxx和yyy的距离不超过kkkxxx和yyy最近公共祖先&#xff1a;zzz,满足vxvy2vzv_xv_y2v_zvx​vy…