[R语言] hdu 5791 Two dp

#include <bits/stdc++.h>
using namespace std;
const int pr=1000000007;
int dp[2010][2010];
int n,m;
int a[2010],b[2010];
int main()
{
  while(scanf("%d%d",&n,&m)!=EOF)
  {
    for(int i=1;i<=n;i++)
      for(int j=1;j<=m;j++) dp[i][j]=0;
    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++) dp[i][0]=1;
   //   for(int i=1;i<=m;i++) dp[0][i]=1;
    for(int i=1;i<=n;i++)
    {
      for(int j=1;j<=m;j++)
      {
        dp[i][j]=dp[i-1][j]+dp[i][j-1]-dp[i-1][j-1];
        if(a[i]==b[j]) dp[i][j]+=dp[i-1][j-1]+1;
        if(dp[i][j]<0) dp[i][j]+=pr;
        if(dp[i][j]>=pr) dp[i][j]%=pr;
       //   cout<<i<<" "<<j<<" "<<dp[i][j]<<endl;
      }
    }
    cout<<dp[n][m]<<endl;
  }
  return 0;
}