算法竞赛---day9(最长公共子序列「LCS」)

世间上本没有路,走的人多了,也便就有了路。。。
不要想一次就把一样东西搞得通透和全面,在处理某些事情的时候,我们只能尽自己的所能,如果错过了最佳补救的时间,那就下次,或者接着补全,不要重新翻到过去重新再来一遍,那样会适得其反,总之,一次不行,那就多来几轮,总会越来越好的,路就是这样慢慢走出来的。。。
最经典的dp问题之一

  • 最长公共子序列(LCS)
//#include<bits/stdc++.h>
#include <iostream>
#include<cstring>
#include<cstdio>
#include<stack>
using namespace std;
void LCS(string s1,string s2){
	int m=s1.length()+1;
	int n=s2.length()+1;
	int **c;
	int **b;
	c=new int* [m];
	b=new int* [m];
	for(int i=0;i<m;i++){
		c[i]=new int [n];
		b[i]=new int [n];
		for(int j=0;j<n;j++){
			b[i][j]=0;
		}
	}
	for(int i=0;i<m;i++){
		c[i][0]=0;
	}
	for(int i=0;i<n;i++){
		c[0][i]=0;
	}
	for(int i=0;i<m-1;i++){
		for(int j=0;j<n-1;j++){
			if(s1[i]==s2[j]){
				c[i+1][j+1]=c[i][j]+1;
				b[i+1][j+1]=1;//1表示箭头为 左上 
			}else if(c[i][j+1]>=c[i+1][j]){
				c[i+1][j+1]=c[i][j+1];
				b[i+1][j+1]=2;//2表示箭头为 向上 
			}else{
				c[i+1][j+1]=c[i+1][j];
				b[i+1][j+1]=3;//3表示箭头为 向左 
			} 
		}
	}
	for(int i=0;i<m;i++){//输出c数组 
		for(int j=0;j<n;j++){
			printf("%d ",c[i][j]); 
		}
		printf("\n"); 
	} 
	stack<char>same;//存LCS字符 
	stack<int>same1,same2;//存LCS字符在字符串1和字符串2中对应的下标,方便显示出来 
	for(int i=m-1,j=n-1;i>=0&&j>=0;){
		if(b[i][j]==1){
			i--;
			j--;
			same.push(s1[i]);
			same1.push(i);
			same2.push(j);
		}else if(b[i][j]==2){
			i--;
		}else{
			j--;
		}
	} 
	printf("%s\n",s1.c_str());//输出字符串 1
	for(int i=0;i<m&&!same1.empty();i++){//输出字符串1的标记 
		if(i==same1.top()){
			printf("1");
			same1.pop();
		}else{
			printf(" ");
		}	
	} 
	printf("\n%s\n",s2.c_str());//输出字符串 2
	for(int i=0;i<n&&!same2.empty();i++){//输出字符串2的标记 
		if(i==same2.top()){
			printf("1");
			same2.pop();
		}else{
			printf(" ");
		}
	}
	printf("\n最长公共子序列为:");
	while(!same.empty()){
		printf("%c",same.top());
		same.pop();
	} 
	printf("\n长度为:%d\n",c[m-1][n-1]);
	for(int i=0;i<m;i++){
		delete[]c[i];
		delete[]b[i];
	}
	delete[]c;
	delete[]b;
}
int main(){
	string s1="ABCPDSFJGODIHJOFDIUSHGD";
    string s2="OSDIHGKODGHBLKSJBHKAGHI";
    LCS(s1,s2);
	return 0;
}

Tips:printf如何输出string:
在string类型后加上c_str()就可以了。例如:printf("%s\n",s1.c_str());

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值