Substrings poj1226 (后缀数组)
来源:互联网 发布:淘宝售前客服工资算法 编辑:程序博客网 时间:2024/06/05 05:35
给n个串,找到一个字符串S。
使得,每个串里要么有S,要么有S的反转。
首先构造字符串,把每个字符串和其反转连接在一起(这个中间也要用字符分割!!很重要!!),然后把所有串拼接起来,两两之间用一个字符(不会出现在串里的)分割。
然后二分答案,分段height数组,检验。
答案错误的同学们试试看
1
2
bcbcb
bcbcb
这个东西。
#include<iostream>#include<cstdio>#include<cstring>#include<vector>using namespace std;const int nMax = 100001;char str[nMax+1];int sa[nMax], rank[nMax], height[nMax];int wa[nMax], wb[nMax], wv[nMax], wd[nMax];int cmp(int *r, int a, int b, int l){ return r[a] == r[b] && r[a+l] == r[b+l];}void da(char *r, int n, int m){ // 倍增算法 r为待匹配数组 n为总长度 m为字符范围 int i, j, p, *x = wa, *y = wb, *t; for(i = 0; i < m; i ++) wd[i] = 0; for(i = 0; i < n; i ++) wd[x[i]=r[i]] ++; for(i = 1; i < m; i ++) wd[i] += wd[i-1]; for(i = n-1; i >= 0; i --) sa[-- wd[x[i]]] = i; for(j = 1, p = 1; p < n; j *= 2, 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 < n; i ++) wv[i] = x[y[i]]; for(i = 0; i < m; i ++) wd[i] = 0; for(i = 0; i < n; i ++) wd[wv[i]] ++; for(i = 1; i < m; i ++) wd[i] += wd[i-1]; for(i = n-1; i >= 0; i --) sa[-- wd[wv[i]]] = y[i]; for(t = x, x = y, y = t, p = 1, x[sa[0]] = 0, i = 1; i < n; i ++){ x[sa[i]] = cmp(y, sa[i-1], sa[i], j) ? p - 1: p ++; } }}void calHeight(char *r, int n){ // 求height数组。 int i, j, k = 0; for(i = 1; i <= n; i ++) { rank[sa[i]] = i; //cout<<sa[i]<<endl; } 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 n;int len[120];char input[120][120];bool appear[120];int find(int sa){ if(sa<len[0]-1) return 0; for(int i=0;i<=n-2;i++) { if(sa>=len[i] && sa<len[i+1]-1) return i+1; } return 110;}bool test(int ans){ int seq=0; for(int j=0;j<n;j++) appear[j]=0; for(int i=len[n-1]-1;i>=1;i--) { if(height[i]<ans) { seq=0; bool ac=true; for(int j=0;j<n;j++) { if(appear[j]==false) ac=false; } if(ac==true) return true; for(int j=0;j<n;j++) appear[j]=0; } else { if(seq==0) { seq+=2; int ia=find(sa[i]); int ib=find(sa[i-1]); appear[ia]=appear[ib]=1; } else { seq++; int ib=find(sa[i-1]); appear[ib]=1; } } } return false;}int main(){ int t; cin>>t; while(t--) { cin>>n; int pt=0; for(int i=0;i<n;i++) { scanf("%s",input[i]); int now_len=strlen(input[i]); if(i==0) len[i]=2*now_len+2; else len[i]=len[i-1]+2*now_len+2; for(int j=0;j<now_len;j++) str[pt++]=input[i][j]; str[pt++]=(i+2)%10+1;//这边一定要加一个,否则4个bcbcb出来的结果是10个。。。 for(int j=now_len-1;j>=0;j--) str[pt++]=input[i][j]; if(i==n-1) str[pt++]=0; else str[pt++]=(i+1)%10+1; } da(str,len[n-1],140); calHeight(str,len[n-1]-1); int l=1; int r=100; int ans=0; while(l<=r) { int mid=(l+r)/2; if(test(mid)==true) { ans=max(ans,mid); l=mid+1; } else r=mid-1; } printf("%d\n",ans); } return 0;}
0 0
- Substrings poj1226 (后缀数组)
- poj1226 Substrings 后缀数组
- POJ1226:Substrings(后缀数组)
- POJ1226---Substrings(后缀数组+二分)
- [POJ1226]Substrings(后缀数组+二分)
- 后缀数组poj1226
- poj1226(kmp算法,后缀数组)
- POJ1226 Substrings
- poj1226 Substrings
- POJ1226--Substrings
- poj1226--Substrings
- POJ1226-Substrings
- POJ-1226(后缀数组)(Substrings)
- spoj 694 Distinct Substrings(后缀数组)
- SPOJ 694 Distinct Substrings(后缀数组)
- CSU 1632 Repeated Substrings(后缀数组)
- Hdu 5769 Substrings (后缀数组)
- [SPOJ705]DISUBSTR - Distinct Substrings(后缀数组)
- android中的MVP模式
- HDU 1224(dp46)
- nerdTree安装详细介绍
- Android - 监听软键盘状态以及获取软键盘的高度
- ZZULIOJ 1917: E(序列的价值,思维,vector)
- Substrings poj1226 (后缀数组)
- LoadRunner11的安装与破解
- BZOJ2827 千山鸟飞绝
- 使用ZXing扫描和生成二维码
- Ubuntu双系统安装,为后面Ubuntu编译ndk做准备
- Linux 系统下原版 texlive 2016 的安装与配置
- mysql如何查询一个表中所有字段的名字
- HTTP协议形象展现
- servletconfig