URAL 1297 后缀数组+线段树
2017-01-15 23:27
225 查看
思路:
论文题……*n
倒过来接上 分奇偶讨论 求LCP 搞棵线段树即可
//By SiriusRen #include <cstdio> #include <cstring> #include <algorithm> using namespace std; #define N 5555 int cases,pos,maxans,n,cntA ,cntB ,A ,B ,rk ,sa ,tsa ,ht ; char s ; void SA(){ memset(cntA,0,sizeof(cntA)); for(int i=1;i<=n;i++)cntA[s[i]]++; for(int i=1;i<=256;i++)cntA[i]+=cntA[i-1]; for(int i=n;i;i--)sa[cntA[s[i]]--]=i; rk[sa[1]]=1; for(int i=2;i<=n;i++)rk[sa[i]]=rk[sa[i-1]]+(s[sa[i]]!=s[sa[i-1]]); for(int l=1;rk[sa ]<n;l<<=1){ memset(cntA,0,sizeof(cntA)); memset(cntB,0,sizeof(cntB)); for(int i=1;i<=n;i++) cntA[A[i]=rk[i]]++, cntB[B[i]=i+l<=n?rk[i+l]:0]++; for(int i=1;i<=n;i++)cntA[i]+=cntA[i-1],cntB[i]+=cntB[i-1]; for(int i=n;i;i--)tsa[cntB[B[i]]--]=i; for(int i=n;i;i--)sa[cntA[A[tsa[i]]]--]=tsa[i]; rk[sa[1]]=1; for(int i=2;i<=n;i++)rk[sa[i]]=rk[sa[i-1]]+(A[sa[i]]!=A[sa[i-1]]||B[sa[i]]!=B[sa[i-1]]); } for(int i=1,j=0;i<=n;i++){ j=j?j-1:0; while(s[i+j]==s[sa[rk[i]-1]+j])j++; ht[rk[i]]=j; } } int tree[N*8]; void build(int l,int r,int pos){ if(l==r){tree[pos]=ht[l];return;} int mid=(l+r)>>1,lson=pos<<1,rson=pos<<1|1; build(l,mid,lson),build(mid+1,r,rson); tree[pos]=min(tree[lson],tree[rson]); } int query(int l,int r,int pos,int L,int R){ if(l>=L&&r<=R)return tree[pos]; int mid=(l+r)>>1,lson=pos<<1,rson=pos<<1|1; if(mid<L)return query(mid+1,r,rson,L,R); else if(mid>=R)return query(l,mid,lson,L,R); else return min(query(l,mid,lson,L,R),query(mid+1,r,rson,L,R)); } int main(){ scanf("%s",s+1),n=strlen(s+1); s[n+1]='#'; for(int i=1;i<=n;i++)s[n+i+1]=s[n-i+1]; n=n+n+1; SA(),build(1,n,1); for(int i=1;i<=n/2;i++){ int tempx=min(rk[i],rk[n-i+1]),tempy=max(rk[i],rk[n-i+1]),jy; jy=query(1,n,1,tempx+1,tempy); if(maxans<jy*2-1)maxans=jy*2-1,pos=i-jy+1; else if(maxans==jy*2-1&&pos>i-jy+1)pos=i-jy+1; tempx=min(rk[i+1],rk[n-i+1]),tempy=max(rk[i+1],rk[n-i+1]); jy=query(1,n,1,tempx+1,tempy); if(maxans<jy*2)maxans=jy*2,pos=i-jy+1; else if(maxans==jy*2&&pos>i-jy+1)pos=i-jy+1; } for(int i=pos;i<maxans+pos;i++)printf("%c",s[i]); }
相关文章推荐
- URAL 1297 后缀数组+线段树
- ural 1297(求最长回文子串,后缀数组)
- [URAL1297 Palindrome]
- URAL 1297 Palindrome(后缀数组求最长回文子串)
- 【后缀数组、最长回文字串】ural1297
- ural1297 Palindrome 后缀数组
- URAL 1297. Palindrome(输出最长回文子串--后缀数组)
- Palindrome - URAL - 1297(求回文串)
- 【URAL1297】Palindrome【Manacher】
- URAL 1297 Palindrome (后缀数组)
- 最长回文子串(ural1297)
- ural 1297Palindrome(求最长回文)
- Ural1297(Palindrome)求最长回文子串(后缀数组)
- URAL1297------后缀数组
- HDU3068 URAL1297 Manacher求最长回文子串
- URAL 1297 Palindrome(最长回文子串:后缀数组)
- Ural 1297 Palindrome 【最长回文子串】
- [URAL 1297]Palindrome(后缀数组)
- Manacher Ural 1297 Palindrome
- 【URAL 1297】Palindrome 最长回文子串