TOJ 3601.Longest Contiguous Subsequence 最长连续公共子序列

最长连续公共子序列问题

最长连续公共子序列，和公共子序列不一样。这里要求子序列连续，相对而言要简单一些。

---

算法设计思路是动态规划（虽然我感觉和动态规划的感觉不怎么相对）。假设两个序列s1[m]和s2[n]，循环遍历s1和s2，当出现 s1[i] == s2[j] 的时候，这一对（i，j）就是公共序列的一部分，如果（i == 0 || j == 0）则（i，j）就是这个公共序列的开始，否则，则可能是一个公共序列的中间部分或者也可能是开始。用一个二维矩阵s[m][n]来标记，s[m][n]表示s1[i]和s2[j]之间的对应关系。
如果s1[i] != s2[j]，则 s[i][j] == s[i - 1][j - 1] + 1
若s[i][j] == 0，则表明s1[i] != s2[j]，若 s[i][j] == 1，则表明（i，j）是一个公共序列的开始部分，若s[i][j] > 1，表明（i，j）是公共序列的中间部分。

程序代码

#include <stdio.h>

int main(){
    int m,n;
    scanf("%d %d",&m,&n);
    const int p = m;
    const int q = n;
    int s1[p];
    int s2[q];
    for(int i = 0;i < m;i ++)
        scanf("%d",s1 + i);
    for(int i = 0;i < n;i ++)
        scanf("%d",s2 + i);
    int s[p][q];
    int max = 0;
    for(int i = 0;i < m;i ++)
        for(int j = 0;j < n;j ++)
            s[i][j] = 0;
    for(int i = 0;i < m;i ++){
        for(int j = 0;j < n;j ++){
            if(s1[i] == s2[j]){
                if(i == 0||j == 0)
                    s[i][j] = 1;
                else{
                    s[i][j] = s[i - 1][j - 1] + 1;
                    if(s[i][j] > max)
                        max = s[i][j];
                }
            }
        }
    }
    printf("%d\n",max);
    return 0;
}