题目背景
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;
}