SPOJ 3408 DNA Sequences

有点卡，一直纠结于dp数组转移相关sum数组转移的问题，看了下别人的题解发现，其实独立开来就行了。求最长公共子序列，且每一连续子段的长度不小于k。用sum数组记录前面有多少连续的序列。sum[i][j] 只能由sum[i - 1][j - 1]转移而来。sum数组的转移与dp数组的转移独立开来。当sum[i][j] > k时，可以选择从前一位连过来，也可以从前k位连过来。[i,j]序列还是表示当前位置的状态好点，若表示[i - 1, j - 1]的状态则最后会略纠结。

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <cmath>
using namespace std;
typedef long long LL;
#define N 1005
int sum[N][N];
int dp[N][N];
int k;
char a[N];
char b[N];
int main(){
    dp[0][0] = 0;
    dp[0][1] = 0;
    dp[1][0] = 0;
    while(scanf("%d", &k) , k){
        scanf("%s", a + 1);
        scanf("%s", b + 1);
        int lena = strlen(a + 1);
        int lenb = strlen(b + 1);

        for(int i = 1; i <= lena; i++){
            for(int j = 1; j <= lenb; j++){
                sum[i][j] = 0;
                if(a[i] == b[j]){
                    sum[i][j] = sum[i - 1][j - 1] + 1;
                }
                dp[i][j] = max(dp[i - 1][j], dp[i][j - 1]);
                if(sum[i][j] == k){
                    dp[i][j] = max(dp[i][j], dp[i - k][j - k] + k);
                }
                if(sum[i][j] > k){
                    dp[i][j] = max(dp[i - k][j - k] + k, dp[i - 1][j - 1] + 1);
                }

            }
        }
        printf("%d\n", dp[lena][lenb]);
    }
    return 0;
}