ICode9

精准搜索请尝试: 精确搜索
首页 > 编程语言> 文章详细

kmp算法完成DNA的病毒感染检测

2021-07-04 17:02:42  阅读:297  来源: 互联网

标签:DNA int chars next maxn kmp return SString 病毒感染


运行结果:

watermark,size_14,text_QDUxQ1RP5Y2a5a6i,color_FFFFFF,t_100,g_se,x_10,y_10,shadow_20,type_ZmFuZ3poZW5naGVpdGk=

代码:

#include<cstring>
#include<iostream>

using namespace std;

#define OK 1
#define ERROR 0;
#define OVERFLOW -2
typedef int Status;
#define MAXSTRLEN 255
typedef char SString[MAXSTRLEN + 1];
const int maxn = 100;

Status StrAssign(SString T, char *chars){
    //生成一个其值等于chars的串T
    int i;
    if (strlen(chars) > MAXSTRLEN)
    {
        return ERROR;
    }
    else
    {
        T[0] = strlen(chars);
        for (i = 1; i <= T[0]; i++)
            T[i] = *(chars + i - 1);
        return OK;
    }//else
}//StrAssign
void get_next(SString T, int next[]){
    //求模式串T的next函数值并存入数组next。
    int i = 1;
    next[1] = 0;
    int j = 0;
    while (i < T[0])
    {
        if (j == 0 || T[i] == T[j])
        {
            ++i; ++j; next[i] = j;
        }//if
        else j = next[j];
    }//while
}//get_next
int Index_KMP(SString S, SString T, int pos)
{
    //利用模式串T的next函数求T在主串S中第pos个字符之后的位置的
    //KMP算法。其中,T非空,1<<pos<<StrLength(S).
    int i = pos;
    int j = 1;
    int next[maxn];
    get_next(T, next);
    while (i <= S[0] && j <= T[0]){
        if (j == 0 || S[i] == T[j]){//继续比较后继字符
            ++i; ++j;
        }//if
        else j = next[j];//模式串向右移动
    }//while
    if (j > T[0])
        return i - T[0];//匹配成功
    else return 0;
}//Index_KMP

int main()
{
    int next[maxn];
    char str[maxn], mo[maxn];
    SString S, T;

    cout << "请输入DNA序列" << endl;
    cin >> mo;
    cout << "请输入病毒序列" << endl;
    cin >> str;
    StrAssign(S, mo);
    StrAssign(T, str);
    next[0] = -1;
    get_next(T,next);
    //printf("%d\n", Index_KMP(S,T,0));
    cout << "病毒在DNA在第" << Index_KMP(S, T, 1) << "个结点处首次匹配\n";
    return 0;

}

标签:DNA,int,chars,next,maxn,kmp,return,SString,病毒感染
来源: https://blog.51cto.com/u_15098794/2973858

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有