最长公共子序列

算法原理参考 http://blog.csdn.net/v_july_v/article/details/6695482

ZOJ1733 基础题无疑

DP[i][j]存放的是对于A字符串和B字符串，当A到达 i 和B到达j 时，最大公共字符串的长度。我们定义temp，如果A的第i个字符和B的第j个字符相等，则temp为1，否则为0。则对于DP[i][j]，到达DP[i][j]的可能有三个，即从DP[i-1][j], DP[i-1][j-1], dp[i][j-1]到达。所以DP[i][j]=max(DP[i-1][j-1]+temp,DP[i][j-1],DP[i-1][j])。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
#define MAXN 1001
char a[MAXN],b[MAXN];
int dp[MAXN][MAXN];
int max3(const int &a,const int &b,const int &c)
{
	int temp=max(a,b);
	temp=max(temp,c);
	return temp;
}
int main()
{
    while(cin>>a>>b)
    {
        int t,i,j;
        memset(dp,0,sizeof(dp));
        int len1=strlen(a);
        int len2=strlen(b);
        for(i=1;i<=len1;i++)
        {
            for(j=1;j<=len2;j++)
            {
                t=(a[i-1]==b[j-1])?1:0;
                dp[i][j]=max3(dp[i-1][j-1]+t,dp[i][j-1],dp[i-1][j]);
            }
        }
        printf("%d\n",dp[len1][len2]);
    }
    return 0;
}

ZOJ1170

一开始用最长公共连续子序列写发现第三个本地案例错了= =

最长公共连续子序列参考： http://www.cnblogs.com/ErinCodeMM/archive/2012/10/30/2747042.html

我的代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int gcd(int n,int m)
{
    if(n%m==0)return m;
    else gcd(m,n%m);
}
void yuefen(int a,int b)//a<b
{
    if(a==0)printf("0\n");
    else if(a==b)printf("1\n");
    else{
        int t=gcd(a,b);
        printf("%d/%d\n",a/t,b/t);
    }
}
int main()
{
   // freopen("test.in","r",stdin);
    //freopen("test.out","w",stdout);
    char a[1001],b[1001];
    while(~scanf("%s",a))
    {
        if(a[0]=='-')break;
        scanf("%s",b);
      //  printf("%s %s\n",a,b);
        int i,j,k;
        int lena=strlen(a);
        int lenb=strlen(b);
       // printf("%d %d\n",lena,lenb);
        int ma=0;
        int s;
        for(i=0;i<lena;i++)
        {
            k=i,j=0;
            s=0;
            while(k<lena&&j<lenb)
            {
                if(a[k]==b[j])s++;
                k++;
                j++;
            }
            if(s>ma)ma=s;
        }
        for(i=0;i<lenb;i++)
        {
            k=i;j=0;
            s=0;
            while(k<lenb&&j<lena)
            {
                if(a[j]==b[k])s++;
                k++;
                j++;
            }
            if(s>ma)ma=s;
        }
        printf("appx(%s,%s) = ",a,b);
        yuefen(ma*2,lena+lenb);
    }
    return 0;
}

但是网上搜了一下题解是可以用最长公共连续子序列写的：

参考：http://blog.163.com/software_ark/blog/static/1756145942011330101437599/

关键代码：

for(i=0;i<lena;i++)  //LCS表格的改进表生成
              {
                     for(j=0;j<lenb;j++)
                     {
                            if(a[i]==b[j])
                            {
                                   if(i>0&&j>0)num[i][j]=num[i-1][j-1]+1;//可以有间隔的连续
                                   else num[i][j]=1;//i，j一旦有一个是0就要重新匹配

                            }
                            else {
                                    if(i>0&&j>0)num[i][j]=num[i-1][j-1];//可以有间隔的连续
                                    else num[i][j]=0;//i，j一旦有一个是0就要重新匹配
                            }
                            max=max>num[i][j]?max:num[i][j];
                     }

              }

学习了