uva 10069 Distinct Subsequences(dp,大数加法)

uva 10069 Distinct Subsequences

dp[i][j],字符串x的前i个字符包含字符串z前j个字符的个数

转移方程为(下标从1开始)：

if(x[i]==z[j]) dp[i][j]=dp[i-1][j]+dp[i-1][j-1]+(j==1?1:0)

else dp[i][j]=dp[i-1][j]

本题数据非常之大，要用大整数加法及大进制，这里用了亿进制

#include<stdio.h>
#include<string.h>
#define MAXN 10005
#define MAXM 105
#define MOD 1000000000
struct node
{
	int num[MAXM];
	int n;
};
char x[MAXN],z[1005];
node dp[MAXN][105];
int max(int a,int b) {return a>b?a:b;}
void add(node &a,node &b,node &c)//c=a+b
{
	int i,up=0;
	int len=max(a.n,b.n);
	for(i=0;i<len;i++)
	{
		int temp=(i>=a.n?0:a.num[i])+(i>=b.n?0:b.num[i])+up;
		c.num[i]=temp%MOD;
		up=temp/MOD;
	}
	if(up)
		c.num[i]=up,c.n=i+1;
	else
		c.n=i;
}
void add1(node &a)
{
	int up=1,i=0;
	while(up)
	{
		int temp=(i>=a.n?0:a.num[i])+up;
		a.num[i]=temp%MOD;
		up=temp/MOD;
		i++;
	}
	if(i>a.n) a.n=i;
}
int main()
{
	int cas;
	scanf("%d",&cas);
	while(cas--)
	{
		scanf("%s%s",&x[1],&z[1]);
		int lenx=strlen(&x[1]),lenz=strlen(&z[1]);
		for(int i=1;i<=lenx;i++) dp[i][0].num[0]=0,dp[i][0].n=0;
		for(int i=1;i<=lenx;i++)
			for(int j=1;j<=lenz;j++)
			{
				if(j>i) break;
				dp[i][j].num[0]=0,dp[i][j].n=0;
				if(x[i]==z[j])
				{
					add(dp[i-1][j],dp[i-1][j-1],dp[i][j]);
					if(j==1)
						add1(dp[i][j]);
				}
				else
					dp[i][j]=dp[i-1][j];
			}
			int len=dp[lenx][lenz].n;
			if(len==0) printf("0\n");
			else 
				printf("%d",dp[lenx][lenz].num[len-1]);
			for(int i=len-2;i>=0;i--) printf("%09d",dp[lenx][lenz].num[i]);
			printf("\n");
	}
	return 0;
}