洛谷P4999 烦人的数字游戏（数位dp）

题目背景

NOIP2018初赛结束了，但H还是要上学的啊啊。。上学肯定要完成老师布置的作业，H十分头疼。在如山的作业中，Mr.G布置的数学作业最烦人，H总是完不成~~

题目描述

Mr.G最近在看一些关于数字题的书，他每天愁同学们太聪明了，所有的作业同学们都能做到全对（拿到答案）。Mr.G蒙在鼓里（心知肚明）。为了使同学们进步，Mr.G总是创造一些简单（毒瘤）题来作为作业。以下是数学作业的最后一题题干——

给出一个区间 L~R，求L到R区间内每个数的数字和，如123这个数的数字和为1+2+3=6。

（1≤L≤R≤10的平方）

同学们纷纷做出来了，Mr.G一看这最后一题跟摆设没区别了呀，于是他迅速修改了题目，把范围定得非常非常大，且有T组数据，将最终的答案mod 10的九次方+7。

（1≤L≤R≤10的18次方） （1≤T≤20）

同学们纷纷被难住了。但H为了备战NOIP2018，没有时间完成Mr.G的数学作业~~（其实是不想做QwQ）~~，所以Ta找到了你，希望你帮助Ta和同学完成这烦人的数学作业！

输入格式

输入共T+1行，

第11行读入T。代表有T组数据；

第2~T+1行。读Li​和Ri​

输出格式

输出共T行，

每行输出Li​和Ri​的区间数字和mod 10的9次方+7。

输入输出样例

输入 

2
24 69
70 120

输出

411
498

思路：

1.首先先预处理数组， f[i][j]表示有i位，最高位数字为j

2.然后就是让ans加上比最高位数字小于a[cnt]的所有数，f[cnt][0]表示0到cnt-1位9

3.然后就开始处理最高位数字为a[cnt]的数，依次处理每一位，sum是记录前面已经走过位上数字之和，因为在后续计算中，这些该加上的数字并不会加上，所以需要额外处理

4.非常重要的一点，这个函数处理的是0到x-1区间里的数

5.在计算f[i][0]的时候一定要让i到19，否则就会被卡

ACcode

#include <iostream>
#include <cstring>
using namespace std;
#define ll long long
ll l,r;
int t;
ll f[20][20];
ll power[20];
int a[20];
const int mod=1000000007;

void chushihua(){
	power[0]=1;
	for(int i=1;i<=18;i++) power[i]=power[i-1]*10%mod;
	for(int j=0;j<=9;j++) f[1][j]=j;
	for(int k=2;k<=19;k++){ 这里必须算到第十九位，数据第十八个点就是卡10的十八次方
		for(int j=0;j<=9;j++){
			f[k][j]=j*power[k-1]%mod;
			for(int i=0;i<=9;i++){
				f[k][j]+=f[k-1][i],f[k][j]%=mod;
			}
		}
	}
}

ll dp(ll x){
	memset(a,0,sizeof a);
	ll ans=0;
	int cnt=0;
	while(x) a[++cnt]=x%10,x/=10;
    for(int i=0;i<a[cnt];i++) ans+=f[cnt][i],ans%=mod;
    int sum=a[cnt];
    for(int j=cnt-1;j>=1;j--){
    	ans+=sum*power[j-1]*a[j]%mod,ans%=mod;
    	for(int k=0;k<a[j];k++){
    		ans+=f[j][k],ans%=mod;
		}
		sum+=a[j];
	}
	return ans;
}

int main()
{
	chushihua();
	cin>>t;
	while(t--){
		cin>>l>>r;
		cout<<((dp(r+1)-dp(l))%mod+mod)%mod<<endl;
	}
	return 0; 
} 

下面这种这种不用算到第十九位

#include <iostream>
#include <cstring>
using namespace std;
#define ll long long
ll l,r;
int t;
ll f[20][20];
ll power[20];
int a[20];
const int mod=1000000007;

void chushihua(){
	power[0]=1;
	for(int i=1;i<=18;i++) power[i]=power[i-1]*10%mod;
	for(int j=0;j<=9;j++) f[1][j]=j;
	for(int k=2;k<=18;k++){
		for(int j=0;j<=9;j++){
			f[k][j]=j*power[k-1]%mod;
			for(int i=0;i<=9;i++){
				f[k][j]+=f[k-1][i],f[k][j]%=mod;
			}
		}
	}
}

ll dp(ll x){
	memset(a,0,sizeof a);
	ll ans=0;
	int cnt=0;
	while(x) a[++cnt]=x%10,x/=10;
	for(int i=1;i<cnt;i++)
		for(int j=1;j<=9;j++) ans+=f[i][j],ans%=mod;  区别在这里
    for(int i=1;i<a[cnt];i++) ans+=f[cnt][i]%mod,ans%=mod;
    int sum=a[cnt];
    for(int j=cnt-1;j>=1;j--){
    	ans+=sum*power[j-1]*a[j]%mod,ans%=mod;
    	for(int k=0;k<a[j];k++){
    		ans+=f[j][k],ans%=mod;
		}
		sum+=a[j];
	}
	return ans;
}

int main()
{
	chushihua();
	cin>>t;
	while(t--){
		cin>>l>>r;
		cout<<((dp(r+1)-dp(l))%mod+mod)%mod<<endl;
	}
	return 0; 
}