UVa 11081 - Strings（dp）

Strings
题意：给出三个字符串，用前两个字符串的子序列构造出后一个，有多少种构造方法
思路：dp[i][j][k]表示第一个字符串前i个，第二个字符串前j个，构造第三个字符串的前k个有多少种方法，刚开始只用了一个数组，怎么都不对，后来看了题解。。UVA 11081
我们用dp的思路将它分解，例如：给三个长度为len1，len2，len3的字符串，那么我们可不可以先求出长度为a，b，c，其中$a<len1,b<len2,c<len3$.由此原问题可以分解为子问题。那么当长度为a,b,c时如何推出len1，len2，len3这个组合呢？i,j,k为三个串的的那个前位置，首先我们给第一个串加一个字符，长度i+1，如果这个字符和第k个字符相等，则当前总的可能情况为f1[i][j][k] = f1[i-1][j][k]，相等时f1[i][j][k] += f[i-1][j][k-1]；同样第二个串也是这样做，最后f[i][j][k] = f1[i][j][k]+f2[i][j][k];这里为什么要用两个数组f1，f2呢，因为这里一个数组已经不能满足题意了，因为前两个串的字符都有可能和第三个串相等，如果我们用一个数组f[i][j][k]表示第一个到i位第二个到j位，第三个到k位，如果s1[i]=s2[j]=s3[k],这种情况无法表示出来。

#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<vector>
#include<cmath>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<algorithm>
using namespace std;
typedef long long LL;
const int maxn=100;
const int maxm=1010;
const int MOD=10007;
const int INF=0x3f3f3f3f;
char s1[maxn],s2[maxn],s3[maxn];
int dp[maxn][maxn][maxn],dp1[maxn][maxn][maxn],dp2[maxn][maxn][maxn];
int main()
{freopen("in.txt","r",stdin);
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%s%s%s",s1+1,s2+1,s3+1);

        int len1=strlen(s1+1),len2=strlen(s2+1),len3=strlen(s3+1);
        memset(dp,0,sizeof(dp));
        memset(dp1,0,sizeof(dp1));
        memset(dp2,0,sizeof(dp2));
        for(int i=0;i<=len1;i++)
            for(int j=0;j<=len2;j++)
                dp[0][i][j]=dp1[0][i][j]=dp2[0][i][j]=1;
        for(int i=1;i<=len3;i++)
        {
            for(int j=0;j<=len1;j++)
            {
                for(int k=0;k<=len2;k++)
                {
                    if(j)
                    {
                        dp1[i][j][k]=dp1[i][j-1][k];
                        if(s3[i]==s1[j])dp1[i][j][k]+=dp[i-1][j-1][k];
                        dp1[i][j][k]%=MOD;
                    }
                    if(k)
                    {
                        dp2[i][j][k]=dp2[i][j][k-1];
                        if(s2[k]==s3[i])dp2[i][j][k]+=dp[i-1][j][k-1];
                        dp2[i][j][k]%=MOD;
                    }
                    dp[i][j][k]=(dp1[i][j][k]+dp2[i][j][k])%MOD;
                }
            }

        }
        printf("%d\n",dp[len3][len1][len2]);
    }
    return 0;
}