hdu5791 Two(dp求公共子序列个数)

求两个序列的公共子序列个数,dp[i][j]表示a序列前i个和b序列前j个公共子序列个数。由容斥原理dp[i][j]=dp[i-1][j]+dp[i][j-1]-dp[i-1][j-1]。特别的，如果a[i]==b[j]，会创造更多的公共子序列，即dp[i][j]+=dp[i-1][j-1]+1.(不带a[i],b[j]的公共子序列都在加上a[i],b[j]构成新的公共子序列，再加上1个a[i],b[j]自己构成的公共子序列)
tips：这题动不动就超int。。。还是用ll的省心。

#include <cstdio>
#include <cstring>
#define N 1005
int const mod=1e9+7;
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,m,a[N],b[N],dp[N][N];
int main(){
//  freopen("a.in","r",stdin);
    while(~scanf("%d%d",&n,&m)){
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=n;++i) a[i]=read();
        for(int i=1;i<=m;++i) b[i]=read();
        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]+=mod;//不小心就会超。。不如ll 
                dp[i][j]%=mod;
            }
        printf("%d\n",dp[n][m]);
    }
    return 0;
}