解题报告 NOIP2015 子串

这是一道dp题，然后来想想怎么表示状态,对答案有影响的就是A串的第i个字符，B串的第j个字符，和k个子串,简单来说就是和选取的字符和子串的数量有关.
那么设 f[i][j][kk]表示在A串的前i个字符中选kk个子串匹配B串的前j个字符的方案数.求方案数可以采用加法原理，考虑A串的第i个字符,那么这个字符的决策只有取或不取，很明显，加法原理，把不取的方案数和取的方案数加起来就可以，但是状态的定义并不能看出这个字符到底取不取，或者说并不能推出结果来，怎么办呢？

那么就用一个数组s[i][j][kk]来表示在A串的前i个字符中选kk个子串匹配B串的前j个字符的方案数，A串的第i个字符会被取到.那么这个s数组该怎么推出来呢？可以发现，如果取第i个字符也有2种可能，因为kk是一定的，第i个字符可能和第i-1个字符合并成一个子串，那么从s[i-1][j][kk]转移过来，也可能不和第i-1个字符合并成一个子串，那么就要新开一个子串，故kk一定从kk-1转移过来，根据加法原理，
那么s[i][j][kk] = s[i-1][j-1][kk] + f[i-1][j-1][kk-1].

还有一个问题，f数组该怎么推呢？根据之前分析的，把不取和取的加上来即可,即f[i][j][kk] = s[i][j][kk] + f[i-1][j][kk],其实可以发现，f数组就相当于前缀和数组.

还有一个问题：这是一个三维的状态转移方程！空间不一定开的下，再看数据范围，这绝对MLE,怎么办?注意到i只能从i或i-1转移过来，可以想到滚动数组，滚动数组有一个比较好些的写法就是把第一维全部加上&1,或者干脆弄两个变量处理完一组数据就交换即可.

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>

using namespace std;

int n, m, k;

long long f[2][205][205], s[2][205][205];
char a[1005], b[205];
const int mod = 1000000007;

int main()
{
    scanf("%d%d%d", &n, &m, &k);
    scanf("%s", a + 1);
    scanf("%s", b + 1);
    int now = 1, last = 0;
    f[0][0][0] = 1;
    for (int i = 1; i <= n; i++)
    {
        f[now][0][0] = 1;
        for (int j = 1; j <= m; j++)
            for (int kk = 1; kk <= k; kk++)
            {
                if (a[i] == b[j])
                    s[now][j][kk] = (s[last][j - 1][kk] + f[last][j - 1][kk - 1]) % mod;
                else
                    s[now][j][kk] = 0;
                f[now][j][kk] = (f[last][j][kk] + s[now][j][kk]) % mod;
            }
        swap(now, last); //一定要变，否则就是一直自己跟自己转移.
    }
    printf("%lld\n", f[last][m][k]);

    return 0;
}