国产探花免费观看_亚洲丰满少妇自慰呻吟_97日韩有码在线_资源在线日韩欧美_一区二区精品毛片,辰东完美世界有声小说,欢乐颂第一季,yy玄幻小说排行榜完本

首頁(yè) > 學(xué)院 > 開(kāi)發(fā)設(shè)計(jì) > 正文

POJ 1509 環(huán)形后綴數(shù)組

2019-11-09 19:43:07
字體:
來(lái)源:轉(zhuǎn)載
供稿:網(wǎng)友
////  main.cpp//  后綴數(shù)組////  Created by liuzhe on 17/2/5.//  Copyright © 2017年 my_code. All rights reserved.//#include <iostream>#include <cstdio>#include <cstring>#include <string>#include <algorithm>using namespace std;
//分析是學(xué)長(zhǎng)的/*題意:給個(gè)串問(wèn)最小的后綴,但是這個(gè)串是個(gè)環(huán)思路:是個(gè)環(huán)也好處理,那么就將串復(fù)制一遍加到后面就行了,
然后求出后綴數(shù)組和高度數(shù)組,剛開(kāi)始我寫(xiě)的是就求出后綴數(shù)
組sa,然后從排名第1開(kāi)始找,比如找到了I,排名為I的從第
sa[I]開(kāi)始,那么如果這個(gè)值小于輸入串的長(zhǎng)度,直接輸出就
行,提交WA,看了看討論,有這樣的數(shù)組,aaaa那么應(yīng)該輸出
1,而那個(gè)程序輸出4,然后想到用高度數(shù)組,剛開(kāi)始找到4的
位置,然后找下一個(gè)排名的高度數(shù)組,如果lcp的值大于串的長(zhǎng)
度,那么更新結(jié)果,以此類推直到lcp的值不再大于原串的長(zhǎng)度,
而為什么大于原串長(zhǎng)度就可以繼續(xù)呢,打個(gè)比方,abab,變成
abababab,sa數(shù)組為4 8 3 7 2 6 1 5,lcp數(shù)組為
0 2 4 6 0 1 3 5 ,然后我們找到的第一個(gè)pos為3,然后
它的下一個(gè)就是pos為1的位置,如果lcp[3]的值大于原串長(zhǎng)度4,
就說(shuō)明從1開(kāi)始的串,可以代替從3開(kāi)始的串,因?yàn)樽罱K的長(zhǎng)度肯
定為4嘛,在說(shuō)的具體一點(diǎn),看的時(shí)候拿紙畫(huà)一下,從3開(kāi)始的串
代表的是3 4 1 2 ,體現(xiàn)在復(fù)制一遍的串上就是3 4 5 6 ,而
從1開(kāi)始的是1 2 3 4,所以如果lcp[3]大于4,代表3 4 5 6
與1 2 3 4是相同的串,懂了嗎~~~~*/// poj 1509const int maxn = 50010;int wa[maxn],wb[maxn],wv[maxn],ww[maxn];int sa[maxn],lcp[maxn],rank[maxn],rank1[maxn];char str[maxn];int t,nn,n,k;int tmp[maxn];bool compare_sa(int i,int j){    if(rank1[i]!=rank1[j]) return rank1[i]<rank1[j];    else{        int ri=i+k<=n? rank1[i+k]:-1;        int rj=j+k<=n? rank1[j+k]:-1;        return ri<rj;    }}void construct_sa(){    for(int i=0;i<=n;i++){        sa[i]=i;        rank1[i]=i<n? str[i]:-1;    }    for(k=1;k<=n;k*=2){        sort(sa,sa+1+n,compare_sa);        tmp[sa[0]]=0;        for(int i=1;i<=n;i++){            tmp[sa[i]]=tmp[sa[i-1]]+(compare_sa(sa[i-1],sa[i])?1:0);        }        for(int i=0;i<=n;i++) rank1[i]=tmp[i];    }}/*inline bool cmp(int *r,int a,int b,int len){    return r[a]==r[b]&&r[a+len]==r[b+len];}void construct_sa(){    int m = 20010;    int i,j,p,*x=wa,*y=wb,*t;    //n++;//why?    for(i=0;i<m;i++) ww[i] = 0;    for(i=0;i<n;i++) ww[x[i]=str[i]]++;    for(i=1;i<m;i++) ww[i] += ww[i-1];    for(i=n-1;i>=0;i--) sa[--ww[x[i]]] = i;    for(j=p=1;p<n;j<<=1,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<m;i++)            ww[i] = 0;        for(i=0;i<n;i++)            ww[wv[i]=x[y[i]]] ++;        for(i=1;i<m;i++)            ww[i]+=ww[i-1];        for(i=n-1;i>=0;i--)            sa[--ww[wv[i]]] = y[i];        for(t=x,x=y,y=t,x[sa[0]]=0,p=i=1;i<n;i++)            x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;    }}*/void construct_lcp(){    for(int i=0;i<=n;i++)        rank1[sa[i]] = i;    int h = 0;    lcp[0] = 0;    for(int i=0;i<n;i++)    {        int j = sa[rank1[i]-1];        if(h>0)            h--;        for(;j+h<n && i+h<n;h++)            if(str[i+h]!=str[j+h])                break;        lcp[rank1[i]-1] = h;    }}int main(){    int T;    scanf("%d",&T);    while(T--){        scanf("%s",str);        int nn=strlen(str);        for(int i=0;i<nn;i++) str[nn+i]=str[i];        n=2*nn;        construct_sa();        construct_lcp();        int pos=n+1;        for(int i=1;i<=n;i++){            if(sa[i]<nn){                pos=sa[i++];                while(i<=n&&lcp[i-1]>=nn){                    pos=min(pos,sa[i]);i++;                }                break;            }        }        PRintf("%d/n",pos+1);    }    return 0;}/*int main(){    int T;    cin>>T;    while(T--)    {        scanf("%s",str);        int nn = strlen(str);        for(int i=0;i<nn;i++)            str[nn+i] = str[i];        int n = nn<<1;        construct_sa(n+1,200010);        construct_lcp(n);        int pos = n+1;        for(int i=1;i<=n;i++)        {            if(sa[i]<nn)            {                pos = sa[i++];                while(i<=nn&&lcp[i-1]>=nn)                {                    pos = min(pos,sa[i]);                    i++;                }                break;            }        }        printf("%d/n",pos+1);    }    return 0;}*/
上一篇:poj2181

下一篇:poj2141

發(fā)表評(píng)論 共有條評(píng)論
用戶名: 密碼:
驗(yàn)證碼: 匿名發(fā)表
主站蜘蛛池模板: 晋中市| 那坡县| 宣武区| 永仁县| 德惠市| 翁牛特旗| 蓬莱市| 车险| 文成县| 连州市| 成都市| 镇雄县| 南宁市| 樟树市| 罗源县| 忻州市| 独山县| 铜梁县| 肃宁县| 海丰县| 沙湾县| 鲁山县| 隆尧县| 来凤县| 镇雄县| 阿拉善盟| 荆门市| 永安市| 曲沃县| 西乌珠穆沁旗| 岳阳市| 吴川市| 乐平市| 虹口区| 双牌县| 丰原市| 津市市| 清新县| 靖边县| 宁河县| 建宁县|