一道关于KMP的题求助orz 附代码
  • 板块题目总版
  • 楼主shulker31
  • 当前回复0
  • 已保存回复0
  • 发布时间2023/10/19 17:09
  • 上次更新2023/10/22 13:19:06
查看原帖
一道关于KMP的题求助orz 附代码
855922
shulker31楼主2023/10/19 17:09

纯萌新 不知道没有考虑到什么特殊情况,一直WA

求助大佬指点qaq

题目是这样的:

已知n个人的供体基因库和1个人的受体基因,请利用朴素的字符串匹配算法或KMP算法等实现基因匹配。

用例1:

3
TGAAGGTGATTAG
TTAAGGTATATTTATATT
GATTGGTTAGATATTG
TATAGA
2 7 66.67%

注:用例的第一行代表供体基因库的基因串数量n,n为≥1的整数,之后的n行代表基因库中的每个基因串。最后一行代表匹配结果。匹配结果的第一个整数表示和基因库中第几个基因串实现了最佳匹配。第二个整数表示匹配的点位,第三个数表示匹配度。匹配点位是指受体基因和基因库中某个基因串最高相似的第一次最高相似的起始位置p;匹配度是指受体基因和基因库中某个基因串的最高相似度,公式为从p开始连续相似的基因数/受体基因串的长度 *100%。相似度要求保留小数点后2位。如果匹配失败,则输出结果为:0 0 0.00%

我的代码如下:

#include<bits/stdc++.h>
using namespace std;

//定长字符串存储的声明
#define MAXSIZE 255
typedef unsigned char SString [MAXSIZE+1];//0号位置存放字符串长度

//函数提前声明
int getStr(SString T1,SString T2);
void getNext(SString T,int next[]);
void strAssign(SString SS, char *s);
int useNext(SString T1,SString T2,int next[]);
int result(char *t,char *s);


int main(){
	int n;
	cin>>n;
	string input[n];
	
	for(int i=0;i<n;i++){
		cin.clear();
		if(i==0)cin.ignore();//如何清除缓存区    cin.ignore();作用:吃掉一个缓冲区里的字符 
		getline(cin, input[i]);
	}
	
	string T;
	cin.clear();
	getline(cin, T);
    
    int best=0;//第best个供体串是答案 
    int bestLength=0;//最长长度 
	int bestResult[n];//第n个基因与受体匹配的长度 
	int bestPos[n];//第n个基因与受体匹配的初始位置
    
    for(int i=0;i<n;i++){
    	
    	
    	for(int j=0;j<T.length();j++){
    		string T_changedLength=T.substr(0,T.length()-j);
    		
    		char* s = new char[input[i].length() + 1];
		    strcpy(s, input[i].c_str());
		    
		    char* t = new char[T_changedLength.length() + 1];
		    strcpy(t, T_changedLength.c_str());
		    
		    int r=result(t,s);//如果r为0 则结束 此时 bestResult[i]里存储最长匹配段长度 
		    bestPos[i]=r;
		    if(r){
		    	bestResult[i]=T.length()-j;
//		    	cout<<i+1<<"串的最长匹配段长度为:"<< bestResult[i]<<endl;
//		    	cout<<i+1<<"串与受体匹配的初始位置:"<< bestPos[i]<<endl;
		    	if(bestLength<bestResult[i]){
		    		bestLength=bestResult[i];
		    		best=i+1;
				}
		    	break;
			}
		    delete[] s;
		    delete[] t;
		}
	}
//	cout<<best<<"串为为佳"<<endl;
//	cout<<bestLength<<"为最长长度"<<endl;
//  cout<<bestPos[best-1]<<"为最佳初始位置"<<endl;
    cout<<best<<" "<<bestPos[best-1]<<" ";
	printf("%.2f",100.0*bestLength/T.length());
	cout<<"%";
}


int result(char *t,char *s){//返回0 代表匹配失败 其余数字代表第一次匹配成功的位置 
    
    
	int next [MAXSIZE];//next数组的0号位置不管
	int i;
	SString SS,SSS;
	strAssign(SS,t);
	strAssign(SSS,s);
	getNext(SS,next);
	
	return useNext(SSS,SS,next);

}

void strAssign(SString SS,char *s){
	int i,j;
	for(i = 0;i<=MAXSIZE;++i)//清空S,下标从0到255
		SS[i]='\0';
	SS[0] = strlen(s);
	for(j = 1;*s!='\0';++j)//填充S
		SS[j] = *s++;
	
}

int getStr(SString T1, SString T2){
	int i = 1,j =1,k=i;
	while(i<=T1[0]&&j<=T2[0]){
		if(T1[i]==T2[j]){
			++i;
			++j;
		}	
		else{
			j=1;
			i=++k;
		}
	}
	if(j>T2[0])
		return k;
	else
		return -1;
}

void getNext(SString T,int next[]){
	int i = 1,j = 0;				//声明两个指针一个前指针j一个后指针i
	next[1]=0;						//next数组1号位对应模式串第一个字符,即首位置零
	while(i<T[0]){				    //遍历模式串,直到最后一个元素下标为T[0]-1	
		if(j==0||T[i]==T[j]){		//首次或当字符相同时,两指针同时后移,并将i下标的next值设为j
			++i;
			++j;
			next[i]=j;
		}
		else
			j=next[j];				//当字符不同时,将前指针去找对应的next值;
	}
}

int useNext(SString T1,SString T2,int next[]){
	int i = 1,j=1;
	while(i<=T1[0]&&j<=T2[0]){
		if(j==0||T1[i]==T2[j]){
			++i;
			++j;
		}
		else
			j=next[j];
	}
	if(j>T2[0])
		return i-T2[0];
	else
		return 0;
}
2023/10/19 17:09
加载中...