hdu3652：B-number_hdu ["b:-1", "a:1", "b:3", "a:5"]-CSDN博客

本文链接：https://blog.csdn.net/xiaobei1a2b3c/article/details/41912457

本文介绍了一种使用数位动态规划方法解决特定B-Number问题的算法。通过对比两个数列的前缀，根据不同的条件进行分类讨论，利用记忆化搜索减少重复计算，最终高效地求解出所有可能的B-Number数量。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=3652

参考思路：如果对数位dp有一些了解，那么这道题一看就知道用数位dp做，现在就看怎么做。

对于两个前缀，s1s2s3...si（表示以s1s2s3...si开头的所有数）与d1d2d3...dj（与前一个同理），假设si与dj后面分别都还有leng位未确定，并且这些位可以随便取数字，那么有三种情况，以s序列开头的B-number的数量等于以d序列开头的B-number的数量：

1. s1s2...si和d1d2...dj都含有13子序列，并且s1s2...si%13 = d1d2....dj%13。

2. s1s2...si和d1d2...dj都不含有13子序列，并且si = dj = 1，s1s2...si%13 = d1d2....dj%13。

3. s1s2...si和d1d2...dj都不含有13子序列，并且si = dj ，si != 1，s1s2...si%13 = d1d2....dj%13。

那么这样就可以用记忆性优化来做了，对于第一种情况，用dp[i][j]表示，j表示还有j位不确定，就是上面的length，i表示前缀对13的余数，对于第二种情况，用dp[i][j][k]表示，k表示还有k位不确定，j表示前缀是否已1结尾，0为否，1为是，i表示前缀对13的余数。

C++源代码：

#include <iostream>
#include <cstdio>

using namespace std;

int ds[11];//保存数n
int dp1[13][11];//对应第一种情况，做记忆性优化
int dp2[13][2][11];//对应第二、三种情况，做记忆性优化

//index表示还有还有多少位未确定，即解题思路里面的length
//m表示前缀对13的余数
//lastD表示前缀最后一位数字
//f表示该前缀是否是n的前缀
//have13表示前缀是否有13子序列

int work(int index, int m, int lastD, bool f, bool have13) {
	if(index<1) {
		if(m==0 && have13) {
			return 1;
		} else {
			return 0;
		}
	}
	if(have13 && !f && dp1[m][index]!=-1) {
		return dp1[m][index];
	}
	int temp = (lastD==1 ? 1 : 0);
	if(!have13 && !f && dp2[m][temp][index]!=-1) {
		return dp2[m][temp][index];
	}
	int i, j, result = 0;
	if(f) {
		for(i=0; i<=ds[index]; i++) {
			result += work(index-1, (m*10+i)%13, i, i==ds[index], (have13 || (lastD==1 && i==3)));
		}
	} else {
		for(i=0; i<=9; i++) {
			result += work(index-1, (m*10+i)%13, i, false, (have13 || (lastD==1 && i==3)));
		}
	}
	if(have13 && !f) {
		dp1[m][index] = result;
	}
	if(!have13 && !f) {
		dp2[m][temp][index] = result;
	}
	return result;
}

int main()
{
	int n;
	while(scanf("%d", &n)==1) {
		int len = 0;
		while(n != 0) {
			ds[++len] = n%10;
			n /= 10;
		}
		int i, j, k;
		for(i=0; i<13; i++) {
			for(j=0; j<11; j++) {
				dp1[i][j] = -1;
			}
		}
		for(i=0; i<13; i++) {
			for(j=0; j<2; j++) {
				for(k=0; k<11; k++) 
					dp2[i][j][k] = -1;
			}
		}
		printf("%d\n", work(len, 0, 0, true, false));
	}

return 0;
}