日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問(wèn) 生活随笔!

生活随笔

當(dāng)前位置: 首頁(yè) > 编程资源 > 编程问答 >内容正文

编程问答

SPOJ - PHRASES Relevant Phrases of Annihilation(后缀数组+二分)

發(fā)布時(shí)間:2024/4/11 编程问答 28 豆豆
生活随笔 收集整理的這篇文章主要介紹了 SPOJ - PHRASES Relevant Phrases of Annihilation(后缀数组+二分) 小編覺(jué)得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.

題目鏈接:點(diǎn)擊查看

題目大意:給出n個(gè)字符串,求出n個(gè)字符串中的最長(zhǎng)公共子串,且該字串必須在每個(gè)字符串中出現(xiàn)至少兩次且不能重疊

題目分析:結(jié)合之前整理過(guò)的:

  • 一個(gè)字符串中兩個(gè)不重疊的重復(fù)子串的最大長(zhǎng)度
  • n個(gè)字符串中的出現(xiàn)次數(shù)大于等于k次的最長(zhǎng)公共子串
  • 這兩點(diǎn)的解決方法,一結(jié)合就是本題的答案了,雖然說(shuō)得輕松,但是比較考驗(yàn)代碼能力,靜下心來(lái)慢慢寫(xiě)就好啦

    具體解決方法直接看我的模板博客最下方吧:https://blog.csdn.net/qq_45458915/article/details/104026037

    代碼:

    #include<iostream> #include<cstdio> #include<string> #include<ctime> #include<cstring> #include<algorithm> #include<stack> #include<queue> #include<map> #include<set> #include<cmath> #include<sstream> using namespace std;typedef long long LL;const int inf=0x3f3f3f3f;const int N=1e5+100;char str[N];int sa[N]; //SA數(shù)組,表示將S的n個(gè)后綴從小到大排序后把排好序的 //的后綴的開(kāi)頭位置順次放入SA中 int t1[N],t2[N],c[N];int rk[N],height[N],belong[N],len,n,sum;int s[N];bool vis[15];int mmax[15],mmin[15],cnt[15];void build_sa(int s[],int n,int m)//n為添加0后的總長(zhǎng) {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]=s[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) {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;} }void get_height(int s[],int n)//n為添加0后的總長(zhǎng) {int i,j,k=0;for(i=0;i<=n;i++)rk[sa[i]]=i;for(i=0;i<n;i++) {if(k) k--;j=sa[rk[i]-1];while(s[i+k]==s[j+k]) k++;height[rk[i]]=k;} }void solve(int base=128) {build_sa(s,len+1,base);get_height(s,len); }void init(int pos) {memset(vis,false,sizeof(vis));sum=0;for(int i=1;i<=n;i++)mmax[i]=-inf,mmin[i]=inf,cnt[i]=0;int id=belong[sa[pos]];mmax[id]=mmin[id]=sa[pos];cnt[id]++; }bool check(int mid) {init(1);for(int i=2;i<=len;i++){if(!belong[sa[i]])break;if(height[i]<mid){init(i);}else{int id=belong[sa[i]];mmax[id]=max(mmax[id],sa[i]);mmin[id]=min(mmin[id],sa[i]);cnt[id]++;if(cnt[id]>=2&&!vis[id]&&mmax[id]-mmin[id]>=mid){sum++;vis[id]=true;}}if(sum==n)return true;}return false; }int main() { // freopen("input.txt","r",stdin); // ios::sync_with_stdio(false);int w;cin>>w;while(w--){scanf("%d",&n);int pos=30,tot=0;for(int i=1;i<=n;i++){scanf("%s",str);for(int j=0;str[j];j++){s[tot]=str[j]-'a'+1;belong[tot]=i;tot++;}s[tot++]=pos++;}s[tot]=0;len=tot;solve(pos);int l=0,r=inf,ans;while(l<=r){int mid=l+r>>1;if(check(mid)){ans=mid;l=mid+1;}elser=mid-1;}printf("%d\n",ans);}return 0; }

    ?

    總結(jié)

    以上是生活随笔為你收集整理的SPOJ - PHRASES Relevant Phrases of Annihilation(后缀数组+二分)的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。

    如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò),歡迎將生活随笔推薦給好友。