【HDU】 5791 Two

Two

---

题目链接

• Two

---

题目大意

    问A,B两个序列的子序列有多少个相同的。

---

题解

DP

    这个题刚分析的时候感觉像个LCS+容斥的东西…结果后来发现其实可以直接DP。
    我们设$dp_{i,j}$为a序列中到i，B序列中到j的相同子序列数，据此，可以很容易的写出方程：

$$dp_{i,j}=\begin{cases} dp_{i,j-1}+dp_{i-1,j}-dp_{i-1,j-1} &&& a[i]\neq b[j] （不匹配） \\ dp_{i-1,j}+dp_{i,j-1}+1 &&& a[i]= b[j] （匹配） \end{cases}$$

---

代码

#include <iostream>
#include <cstring>
#include <cstdio>
#define maxn 1005
#define LL long long
#define mod 1000000007

using namespace std;

int n,m,a[maxn],b[maxn],d[maxn][maxn];

int main()
{
    while (scanf("%d%d",&n,&m)!=EOF)
    {
        memset(d,0,sizeof(d));
        memset(a,0,sizeof(a));
        memset(b,0,sizeof(b));

        for (int i=1;i<=n;i++) scanf("%d",&a[i]);
        for (int i=1;i<=m;i++) scanf("%d",&b[i]);
        for (int i=1;i<=n;i++)
            for (int j=1;j<=m;j++)
        {
            if (a[i]==b[j]) d[i][j]=(d[i-1][j]+d[i][j-1]+1)%mod;
            else d[i][j]=(d[i-1][j]-d[i-1][j-1]+d[i][j-1])%mod;
        }
        printf("%d\n",(d[n][m]+mod)%mod);
    }
    return 0;
}