数位DP。
一般是利用DFS来求数位DP了,结合了记忆化搜索。设dp[i][j][k]为前i位,并且前i位的数位和mod7为j,前i位的数字的表示数字值mod7。为什么可以这样呢?因为继续DFS下去,必定是得到一个不是7倍数的数的,因而,k这个位只是在确定叶子结点时有用的。
然后,可以这样选,求一个一些数的平方和(A+B)^2=A^2+2AB+B^2,由于B是可以是多个,因为会乘上一个cnt个数。
至于dp,记录的是当前状态下,其后续(因为是DFS)能得到的符合要求的个数,它们的和,以及它们的平方和。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define LL __int64
using namespace std;
const LL MOD=1000000007;
struct Node{
LL cnt,msum,sqsum;
void init(){
cnt=-1,msum=0,sqsum=0;
};
}dp[35][7][7];
LL l,r;
int num[35],len;
LL mod[35];
Node dfs(int len,int prem,int presqm,bool flag){
if(len==0){
Node st;
st.init();
if(prem&&presqm){
st.cnt=1;
}
else st.cnt=0;
return st;
}
if(!flag&&dp[len][prem][presqm].cnt!=-1) return dp[len][prem][presqm];
int up=flag?num[len]:9;
Node ans;
ans.init(); ans.cnt=0;
for(int i=0;i<=up;i++){
if(i==7) continue;
Node st=dfs(len-1,(prem+i)%7,(presqm*10+i)%7,(flag&&i==up)?true:false);
ans.cnt+=st.cnt;
ans.cnt%=MOD;
ans.msum=(ans.msum+st.msum+(st.cnt%MOD*(mod[len]*i)%MOD)%MOD)%MOD;
ans.sqsum=(ans.sqsum+st.sqsum+((2*mod[len])%MOD*(i*st.msum)%MOD)%MOD)%MOD;
ans.sqsum=(ans.sqsum+((mod[len]*i)%MOD*(i*mod[len])%MOD)%MOD*st.cnt)%MOD;
}
if(!flag) dp[len][prem][presqm]=ans;
return ans;
}
LL slove(LL n){
len=0;
while(n){
num[++len]=(int)(n%10);
n/=10;
}
Node tmp=dfs(len,0,0,true);
return tmp.sqsum;
}
int main(){
int T;
mod[1]=1;
for(int i=2;i<35;i++){
mod[i]=(mod[i-1]*10)%MOD;
}
scanf("%d",&T);
while(T--){
for(int i=0;i<35;i++){
for(int k=0;k<7;k++){
for(int l=0;l<7;l++)
dp[i][k][l].init();
}
}
scanf("%I64d%I64d",&l,&r);
printf("%I64d\n",((slove(r)-slove(l-1))%MOD+MOD)%MOD);
}
return 0;
}