uva 10069 Distinct Subsequences

动态规划问题加上大数处理，f[i][j]表示在X的前i个字符中寻找子串Z[1~j]的个数，当X[i]==Z[j]的时候，f[i][j] = f[i-1][j-1] + f[i-1][j]，当X[1]!=Z[j]的时候，f[i][j] = f[i-1][j]，用一维滚动数组来模拟二维dp迭代过程就行了。注意迭代一开始的时候，在空串里面是可以找到一个空的子串的，否则后续的递推也就不成立了。

#include <stdio.h>
#include <string.h>

char buf1[100001];
char buf2[110];

struct dp_node{
	long long value[12];
};

void clear_dp_node(struct dp_node &node_ins){
	int i;
	for(i=0; i<=11; i++)
		node_ins.value[i] = 0;
}

void add_dp_node(struct dp_node &node_ins1, struct dp_node &node_ins2, struct dp_node &node_sum){
	int i;
	
	node_sum.value[0] = 0;
	for(i=0; i<11; i++){
		node_sum.value[i] += (node_ins1.value[i] + node_ins2.value[i]);
		if(node_sum.value[i] >= 10000000000){
			node_sum.value[i+1] = node_sum.value[i]/10000000000;
			node_sum.value[i] %= 10000000000;
		}
		else{
			node_sum.value[i+1] = 0;
		}
	}
}

void print_dp_node(const struct dp_node &node_ins){
	int i, j;
	for(i=11; i>=0; i--){
		if(node_ins.value[i] != 0)
			break;
	}

	if(i == -1){
		printf("0\n");
		return;
	}

	for(j=i; j>=0; j--){
		if(j == i)
			printf("%lld", node_ins.value[j]);
		else
			printf("%010lld", node_ins.value[j]);
	}
	printf("\n");
}

struct dp_node dp[120];

void func(){
	int i, j, len1, len2;
	struct dp_node new_node;

	len1 = strlen(buf1+1);
	len2 = strlen(buf2+1);
	for(j=0; j<=len2; j++) clear_dp_node(dp[j]);
	dp[0].value[0] = 1; //空串是可以找到一个空的子串的

	for(i=1; i<=len1; i++){
		for(j=len2; j>=1; j--){
			if(buf1[i] == buf2[j]){
				add_dp_node(dp[j], dp[j-1], new_node);
				dp[j] = new_node;
			}
		}
	}
	//dp[len2].value[2] = 99;
	print_dp_node(dp[len2]);
}

int main(void){
	int n;

	//freopen("input.dat", "r", stdin);
	scanf("%d", &n);
	while(n--){
		scanf("%s %s", buf1+1, buf2+1);
		func();
	}
	return 0;
}