动态规划——最长公共子序列

给出两个字符串,求它们的最长公共子序列。

a={ATCGATCC}

b={ATCCATCG}

#include<stdio.h>
int max(int a,int b){
	if(a>b)
		return a;
	else
		return b;
}
int main(){
	char a[9],b[9];
	gets(a);
	gets(b);
	int opt[9][9];
	int i,j;
	for(i=0;i<9;i++){
		opt[0][i]=0;
		opt[i][0]=0;
	}
	for(i=0;i<9;i++){
		for(j=0;j<9;j++){
			if(i==0||j==0){
				opt[i][j]=0;
			}
			if(i>0&&j>0&&a[i]==b[j]){
				opt[i][j]=opt[i-1][j-1]+1;
			}
			if(i>0&&j>0&&a[i]!=b[j]){
				opt[i][j]=max(opt[i][j-1],opt[i-1][j]);
			}
		}
	}
	for(i=0;i<9;i++){
		if(i==0)
			printf("    %c ",a[i]);
		else
			printf("%c ",a[i]);
	} 
	printf("\n");
	for(i=0;i<9;i++){
		if(i==0)
			printf("  ");
		for(j=0;j<9;j++){
			printf("%d ",opt[i][j]);
		}
		printf("\n");
		printf("%c ",b[i]);
	}
	return 0;
} 

可以看到最长公共子序列为:ATGATC

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值