吉哥系列故事——恨7不成妻 HDU - 4507

单身!
　　依然单身！
　　吉哥依然单身！
　　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;
}