POJ 3080 查找公共最长子串 kmp+枚举

找出所有串的最长公共子串,我枚举第一个串的每一个子串,然后使用kmp在剩余子串中查找,看是否存在。这里使用了一个减枝,如果枚举的串的长度j-i+1小于已经得出的串,则直接跳过。

//11347502	c00h00g	3080	Accepted	384K	16MS	G++	2503B	2013-03-14 11:54:08
#include<stdio.h>
#include<stdlib.h>
#include<string.h>

int n,m;
char str[12][65];
int next[65]; 
char output[65];

void get_next(const char* ptrn,int plen,int* next){
    int i=0;
    next[i]=-1;
    int j=-1;
    while(i<plen-1){
        if(j==-1||ptrn[i]==ptrn[j]){
            ++i;
            ++j;
            next[i]=j;
        }else{
            j=next[j];
        }
    }
}

int kmp_search(char const* src,int slen,char const* patn,int plen,int const* next,int pos){
    int i=pos;  
    int j=0;  
    while(i<slen&&j<plen){  
        if(j==-1||src[i]==patn[j]){  
            ++i;  
            ++j;  
        }else{  
            j=next[j];  
        }  
    }  
    if(j>=plen)  
        return i-plen;  
    else  
        return -1; 
}

int main(){
    scanf("%d",&n);
    while(n--){
        scanf("%d",&m);
        for(int i=0;i<m;i++)
            scanf("%s",str[i]);
        //枚举第一个中的每一个,然后在其他字符串使用kmp查找
        int len0=strlen(str[0]);
        int res=-1;
        char tmp[65],tmp1[65];
        for(int i=0;i<len0;i++){

            for(int j=i;j<len0;j++){
                if(j-i+1<res) continue;
                int num=0;
                for(int k=i;k<=j;k++) tmp[num++]=str[0][k];
                tmp[num]='\0';
                get_next(tmp,strlen(tmp),next);
                //剩下的每一个string
                bool flag=true;
                for(int ii=1;ii<m;ii++){
                    int ans=kmp_search(str[ii],strlen(str[ii]),tmp,strlen(tmp),next,0);
                    if(ans==-1){
                        flag=false;
                        break;
                    }
                }
                if(flag&&res<j-i+1){
                    res=j-i+1;
                    int pos=0;
                    for(int jj=i;jj<=j;jj++)
                        output[pos++]=str[0][jj];
                    output[pos]='\0';
                }
                if(flag&&res==j-i+1){
                    int pos=0;
                    for(int jj=i;jj<=j;jj++)
                        tmp1[pos++]=str[0][jj];
                    tmp1[pos]='\0';
                    if(strcmp(tmp1,output)<0)
                        strcpy(output,tmp1);
                }
            }
		}
		if(res<3)
			printf("no significant commonalities\n");
		else{
            for(int i=0;i<strlen(output);i++)
	            printf("%c",output[i]);
		    printf("\n");
		}
    }
    //system("pause");
    return 0;
}


  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值