POJ 3080 Blue Jeans(KMP)

来源：互联网发布：java http获得响应编辑：程序博客网时间：2024/05/16 02:06

题意：有多个基因碱基序列(只由AGCT组成)，在这些序列中找出一个最长的公共子串，如果长度相等，优先输出字典序较小的串。

思路：由于数据较小，所以直接枚举第一个串中的子串，然后依次对之后的串每次作一次KMP判断是否匹配即可。

注意：当串的长度相同时，要优先输出字典序较小的串。(尼玛最开始没看到好蛋疼)

# include<iostream># include<string># include<cstdio>using namespace std;char st1[1000],st2[10000],str[10000][10000];int next[10010];int num,l1,l2,i,j,k,start,end;void kmp1()//求next值{    int p,q;    q=-1;    next[0]=-1;    for (p=1;p<=l2-1;p++)    {        while (q>-1&&st2[q+1]!=st2[p])            q=next[q];        if (st2[q+1]==st2[p])            q++;        next[p]=q;    }}void kmp2()//KMP函数{    int p,q;    q=-1;    for (p=0;p<=l1-1;p++)    {        while (q>-1&&str[k][p]!=st2[q+1])            q=next[q];        if (str[k][p]==st2[q+1])            q++;        if (q==l2-1)        {            num++;//如果当前串与模式串匹配则num++            break;        }    }}int main(void){    //freopen("in.txt","r",stdin);    //freopen("out.txt","w",stdout);    int n,t;    cin>>t;    while(t--)    {        int ok=0,max=-1;        start=end=0;        scanf("%d",&n);        for(i=1;i<=n;i++)            scanf("%s",str[i]);        //枚举第一个串中的每个子串作为模式串        for(i=0;i<=57;i++)            for(j=i+2;j<=59;j++)//模式串从i开始到j            {                l2=j-i+1;//l2为模式串长度                l1=60;//l1为主串长度                int count=i;                for(k=0;k<l2;k++)//将模式串存入str数组                    st2[k]=str[1][count++];                kmp1();//求出模式串的next值                num=0;                for(k=2;k<=n;k++)//对之后每一个串求一次KMP                    kmp2();                if(num==n-1)                {                    ok=1;                    if(end-start<j-i)                    {                        start=i;                        end=j;                    }                    else if(end-start==j-i)//当两个模式串相等时要保证字典序                        for(k=0;k<=end-start;k++)                            if(str[1][start+k]>str[1][i+k])                                start=i,end=j;                }            }        if(ok)        {            for(int temp=start;temp<=end;temp++)                printf("%c",str[1][temp]);            printf("\n");        }        else            printf("no significant commonalities\n");    }    return 0;}