poj 3080 Blue Jeans ( 二分 + kmp)

本文详细介绍了如何使用二分查找和KMP算法解决POJ 3080问题。通过代码展示这两种算法在实际问题中的结合运用,帮助读者理解其工作原理。
摘要由CSDN通过智能技术生成
题目:http://poj.org/problem?id=3080
题目大意:给你m个由A、T、C、G组成的字符串,然你找出他们字典序最小的最长的公共子串,如果长度<3,就输出那句话。
思路:先找一个长度最小的,二分长度,然后枚举子串,然后再枚举m个字符串,kmp匹配,就这样。感觉kmp的初学者来做这道题很不错。。

代码如下:

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

const int MAXN = 66;

int fail[MAXN];

void get_fail(char *s)
{
    int len =strlen(s);
    fail[0] = -1;
    int j = -1;
    for(int i = 1;i<len;i++)
    {
        while(j >= 0 && s[i]!=s[j+1]) j = fail[j];
        if(s[i] == s[j+1]) j++;
        fail[i] = j;
    }
}

int kmp(char *s,char *t)
{
    get_fail(s);
    int len1 = strlen(t);
    int len2 = strlen(s);
    for(int i = 0,j = -1;i<len1;i++)
    {
        while(j >=0 && s[j+1] != t[i]) j = fail[j];
        if(s[j+1] == t[i]) j++;
        if(j == len2 - 1) return 1;
    }
    return 0;
}

char str[11][MAXN];

char p[MAXN];

char ans[MAXN];

int main()
{
    int _;
    scanf("%d",&_);
    while(_--)
    {
        int m;
        scanf("%d",&m);
        int min_len = 111;
        int x;
        for(int i = 0;i<m;i++)
        {
            scanf("%s",str[i]);
            int len = strlen(str[i]);
            if(len < min_len)
            {
                min_len = len;
                x = i;
            }
        }
        ans[0] = '\0';
        int l = 1,r = min_len;
        while(l <= r)
        {
            int mid = (l+r) >>1;
            int ok = 0;
            int first = 1;
            for(int i = 0;i<min_len;i++)
            {
                int j = i + mid - 1;
                if(j >= min_len) break;
                int cnt = 0;
                for(int k = i;k <= j; k++)
                    p[cnt++] = str[x][k];
                p[cnt] = '\0';
                int flag = 1;
                for(int k = 0;k<m;k++)
                    if(k != x)
                    {
                        if(!kmp(p,str[k]))
                        {
                            flag = 0;
                            break;

                        }
                    }
                if(flag)
                {
                    ok = 1;
                    if(first || strcmp(ans,p) > 0)
                    {
                        strcpy(ans,p);
                        first = 0;
                    }
                }
            }
            if(ok)
            {
                l = mid + 1;
            }
            else r = mid - 1;
        }
        int ans_len = strlen(ans);
        if(ans_len >= 3) puts(ans);
        else puts("no significant commonalities");
    }
    return 0;
}


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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值