HDU4507吉哥系列故事――恨7不成妻

 

　　单身! 
　　依然单身！ 
　　吉哥依然单身！ 
　　DS级码农吉哥依然单身！ 
　　所以，他生平最恨情人节，不管是214还是77，他都讨厌！ 
　　 
　　吉哥观察了214和77这两个数，发现： 
　　2+1+4=7 
　　7+7=7*2 
　　77=7*11 
　　最终，他发现原来这一切归根到底都是因为和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

 

 

 

 

代码：

 

/*总状态转移思路
(1) ans.cnt += tmp.cnt
(2) ans.s += tmp.s + [ i*10^p ]*tmp.cnt
(3) ans.ss += tmp.ss + 2*(i*10^p)*tmp.s + [(i*10^p)^2]*tmp.cnt（平方和公式）
*/
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstring>
using namespace std;
typedef long long ll;

const int N = 20;
const int mod = 1e9 + 7;

struct Node
{
	ll cnt,s,ss;                        //与 7 无关的数的个数,与 7 无关的数的和,与 7 无关的数的平方和
	Node() { cnt = -1; s = 0, ss = 0; } //默认的构造函数，初始化cnt顺便标记是否访问过
}dp[N + 2][10][10];


ll p[N + 2],dig[N + 2];                 //保存10的次方 PS：记得取模,保存数字的每一个数位上的数

void init()
{
	p[1] = 1;                           //因为我的数字的个位用dig[1]表示，所以我用p[1]记录10的0次方
	for (int i = 2; i < N; ++i) p[i] = (p[i - 1] * 10) % mod;
}

inline ll EX(ll x)                      //求平方
{
	x %= mod;
	return (x*x) % mod;
}

Node dfs(int len, int r1, int r2, bool up)
{
	if (len == 0)
	{
		Node tmp;
		tmp.cnt = (r1&&r2);             //数位和不是7的倍数且数字不是7的倍数,因为这里已经完成了一组排列，只要判断这里的就行了所以&&
		tmp.s = tmp.ss = 0;
		return tmp;
	}
	if (!up&&dp[len][r1][r2].cnt != -1)
        return dp[len][r1][r2];
	int n =up?dig[len]:9;
	Node ans; ans.cnt = 0;
	for (int i = 0; i <= n; ++i)
	{
		if (i == 7) continue;
		Node tmp = dfs(len - 1, (i + r1) % 7, (i + r2 * 10) % 7, up&&i == n);               //得到子状态
		ans.cnt = (ans.cnt + tmp.cnt) % mod;                                                //更新个数
		ans.s = (ans.s + (tmp.s + (i*p[len]) % mod*tmp.cnt%mod) % mod) % mod;               //更新数字
		ans.ss = (ans.ss + (tmp.ss + ((2LL * i*p[len] % mod)*tmp.s) % mod) % mod) % mod;    //式子太长分开写
		ans.ss = (ans.ss + (EX(i*p[len])*tmp.cnt) % mod) % mod;
	}
	if (!up) dp[len][r1][r2] = ans;
	return ans;
}

ll cal(ll x)
{
	int len = 0;
	while (x)
	{
		dig[++len] = x % 10;
		x /= 10;
	}
	return dfs(len, 0, 0, 1).ss;
}

int main()
{
    init();
	int T;
	scanf("%d", &T);
	while (T--)
	{
		ll l, r;
		scanf("%lld %lld", &l, &r);
		printf("%lld\n", ((cal(r) - cal(l - 1)) % mod + mod) % mod);
	}
	return 0;
}