挑战程序设计竞赛笔记-动态规划

2.3.1 记忆化搜索与动态规划

在计算过程中，有可能会有一些状态重复出现，对于这些重复出现的状态，不应该每次都重新计算。而应该将其记录下来。

2.3.1.1 最长公共子序列 Poj 1458

给出两个字符串，求出这样的一个最长的公共子序列的长度：子序列中的每个字符都能在两个原串中找到，而且每个字符的先后顺序和原串中的先后顺序一致。

输入有若干行，每行是两个字符串。对每一行输入的两个字符串，输出最长公共子串的长度。

Sample Input
	abcfbc abfcab
	programming contest 
	abcd mnp

Sample Output
	4
	2
	0
也就是，只要顺序相同的字符，就算是一个公共子序列。

对动态规划，我觉得首先要找到“状态”。像这道题，我们要找的是s1/s2的公共最长子串。定义一个二维数组，形状是len(s1+1) * len(s2+1)。
对于坐标 (x,y)，我们定义为字符串s1、s2中长度为x,y的子串之间的公共最长子串。也就是，我们的目标是（len(s1),len(s2))。也就是全串长的最长公共子串。
对于 (x,y)，我们可总结规律：

s1的前x字符子串和s2的前y字符子串

				1. 若s1[x]=s2[y]：则(x,y) = (x-1,y-1)
				2. 若s1[x]!=s2[y]：则(x,y) = Max( ( x-1 , y ) , (x , y-1) )

1.比较容易想。如果给s1[:x-1],s2[:y-1]各自加上一个字符，如果这两个字符相等，则说明最长公共子串的最后一个字符就是这个。(但是位置不一定，可能s1[x]对应的是s2更前面的字符)

而对于2.，s1[x] 和 s2[y] 不相等，不能直接判断最长公共子串的长度增加与否。但 s1[:x] 和 s2[:y-1] 、s1[:x-1] 和s2[:y] 却有可能公共子串增加。我们设 (x-1,y-1) = K，则 (x-1,y) 相当于给s2添上一个字符。其最长子串肯定≥K，若大于K则必定等于K+1。同理 (x,y-1) 相对于K。

#include<stdio.h>
#include<algorithm>
#include<string.h>
char String1[1000];
char String2[1000];
int DP[1001][1001];

int max(int a,int b){
    if(a>=b)return a;
    else{
        return b;
    }
}

int main(){
    while(scanf("%s %s",String1,String2)!=EOF){
        int l1,l2;
        l1=strlen(String1);
        l2=strlen(String2);
        
        for(int i=1;i<l1+1;i++){
            for(int j=1;j<l2+1;j++){
                if(String1[i-1]==String2[j-1]){
                    DP[i][j]=DP[i-1][j-1]+1;
                }
                else{
                    DP[i][j] = max(DP[i-1][j],DP[i][j-1]);
                }
            }
        }
        printf("%d\n",DP[l1][l2]);
        
        
    }
    
    
}