单身!
依然单身!
吉哥依然单身!
DS级码农吉哥依然单身!
所以,他生平最恨情人节,不管是214还是77,他都讨厌!
吉哥观察了214和77这两个数,发现:
2+1+4=7
7+7=72
77=711
最终,他发现原来这一切归根到底都是因为和7有关!所以,他现在甚至讨厌一切和7有关的数!
什么样的数和7有关呢?
如果一个整数符合下面3个条件之一,那么我们就说这个整数和7有关——
1、整数中某一位是7;
2、整数的每一位加起来的和是7的整数倍;
3、这个整数是7的整数倍;
现在问题来了:吉哥想知道在一定区间内和7无关的数字的平方和。
Input
输入数据的第一行是case数T(1 <= T <= 50),然后接下来的T行表示T个case;每个case在一行内包含两个正整数L, R(1 <= L <= R <= 10^18)。
Output
请计算[L,R]中和7无关的数字的平方和,并将结果对10^9 + 7 求模后输出。
Sample Input
3
1 9
10 11
17 17
Sample Output
236
221
0
思路:求区间内与7无关的数其实很好求,不带7在枚举每一位数时不枚举7即可保证最终的数不带7,用state表示该整数的数位和,sum表示该整数的和对7取余则向下搜索时state=(state%7+i%7) %7,sum=(sum*10+i)%7,当state!=0&&sum!=0时表示该数与7无关。
重点是求这些数的平方和。
对于枚举的一个数,第pos位是i,那么对求和贡献就是i * 10^ pos, 则对于一个pos位的与7无关的数X有X=( A*10 ^ pos+B)^2,A是你当前pos位枚举的值,B是一个子问题,展开得到A * A * 10 ^ 2 * pos+2 * A * 10 ^ pos * B+B^2,我们用一个结构体表示,结构体里有cnt,sum,qsum。cnt表示7无关的数的个数,sum表示这些数的和,qsum表示这些数的平方和。令A=A * 10 ^pos,则now.qsum=now.qsum+cnt * A * A+2 * A * sum+next.qsum^2;
#include<iostream>
#include<cstring>
using namespace std;
typedef long long ll;
ll dit[20],pp[20];
const ll mod=1000000007;
typedef struct node{
ll cnt;
ll sum;
ll qsum;
node():cnt(0),sum(0),qsum(0){}
node(ll c,ll s,ll q):cnt(c),sum(s),qsum(q) {}
}node;
node dp[20][10][10];
node dfs(int pos,int state,int sum,int limit)//state:数位和 sum:sum和
{
if(pos<0)
{
if(state!=0&&sum!=0)
return node(1,0,0);
else
return node(0,0,0);
}
if(!limit&&dp[pos][state][sum].sum!=0) return dp[pos][state][sum];
int up=limit?dit[pos]:9;
node ans;
for(int i=0;i<=up;i++)
{
if(i==7) continue;
node tmp=dfs(pos-1,(state%7+i%7)%7,(sum*10+i)%7,limit&&i==dit[pos]);
ans.cnt=(ans.cnt+tmp.cnt)%mod;
ll res=(i*pp[pos])%mod;
ans.sum=(ans.sum+(tmp.sum+(res*tmp.cnt)%mod)%mod)%mod;
ll a=((res*res)%mod*tmp.cnt)%mod;
ll b=((2*res)%mod*tmp.sum)%mod;
ans.qsum=(ans.qsum+((a+b)%mod+tmp.qsum%mod)%mod)%mod;
}
if(!limit) dp[pos][state][sum]=ans;
return ans;
}
ll solve(ll x)
{
int len=0;
while(x)
{
dit[len++]=x%10;
x/=10;
}
node ans=dfs(len-1,0,0,1);
return ans.qsum%mod;
}
int main()
{
int t;
ll n,m;
cin>>t;
pp[0]=1;
for(int i=1;i<=18;i++)
pp[i]=pp[i-1]*10;
while(t--)
{
cin>>n>>m;
cout<<(solve(m)-solve(n-1)+mod)%mod<<endl;
}
return 0;
}