Uva 11316

lycccccccc

于 2020-02-13 21:43:27 发布

阅读量125

点赞数

分类专栏：数位DP

本文链接：https://blog.csdn.net/lycccccccc/article/details/104302960

版权

数位DP 专栏收录该内容

1 篇文章 0 订阅

订阅专栏

Uva 11316 Investigating Div-Sum Property
题目链接：Uva 11316
题目大意：已知区间[A,B]（A<B<2^31）和k(k<=10000)，求在区间中有多少数字满足：该数字是k的整数倍且该数字的各个位数之和也是k的整数倍。
题目分析：一道数位DP的题，首先我们应该考虑到k的值肯定是不满的，个位数字之和应该大概是一个不到100的数（就算是10位全是9才只到90,2^31应该是略大于1e9）因此我们可以确定所有k大于90的答案肯定是0，接下来考虑其他的数，参考lrj的训练指南上的转移方程（其实有一个地方不太对）
f(d,m1,m2)=sum{f(d-1,(m1-x)mod k,(m2-x*10^(d-1))mod k)|x=0,1……9}（注意m2那里lrj少乘了一个x），最总答案就是
f（B的位数，0，0）-f（A的位数，0，0）
既然有了转移方程我们就可以写出代码了，实际我写的时候是向后转移的

for (int i=1;i<=cnt;i++){
        for (int j=0;j<=n;j++){
            for (int k=0;k<=n;k++){
                for(int x=0;x<10;x++)
                    dp[i][(k+x)%n][(j*10+x)%n]+=dp[i-1][k][j];
            }
        }
        for (int j=0;j<num[i];j++)//<-第二个j的for循环
            dp[i][(nt+j)%n][(ns*10+j)%n]++;
        ns=(ns*10+num[i])%n;
        nt=(nt+num[i])%n;
    }

emm，变量用的有点混乱，这里n代表的是模数也就是题面中的k，然后需要考虑怎么才能不超过我们已知的A和B呢，在转移时我们是从每个数的高位向低位转移，因此每次只有经过了第二个j的for循环的dp才是有值的，因此不必担心会有超过原数据范围的数字被统计上去。
下面是完整代码：

#include<stdio.h>
#include<algorithm>
#include<string.h>
#define LL long long
using namespace std;
int n,t,l,r,k,num[110],cnt;
LL dp[15][150][150];
void Pre(int m){
	memset(num,0,sizeof(num));
	memset(dp,0,sizeof(dp));
	cnt=1;
	while(m){
		num[cnt++]=m%10;
		m/=10;
	}
	for(int i=1;i<=cnt/2;i++)
		swap(num[i],num[cnt-i+1]);
}
LL DP(int m){
	Pre(m);
	if(m==0)return 1;
    memset(dp,0,sizeof(dp));
    int ns=0,nt=0;
    for (int i=1;i<=cnt;i++){
        for (int j=0;j<=n;j++){
            for (int k=0;k<=n;k++){
                for(int x=0;x<10;x++)
                    dp[i][(k+x)%n][(j*10+x)%n]+=dp[i-1][k][j];
            }
        }
        for (int j=0;j<num[i];j++)
            dp[i][(nt+j)%n][(ns*10+j)%n]++;
        ns=(ns*10+num[i])%n;
        nt=(nt+num[i])%n;
    }
    if (ns==0&&nt==0)
        dp[cnt][0][0]++;
    return dp[cnt][0][0];
}
int main()
{
	scanf("%d",&t);
	while(t--){
		scanf("%d %d %d",&l,&r,&n);
		if(n>85)printf("0\n");
		else printf("%lld\n",DP(r)-DP(l-1));
	}
	return 0;
 }