hdu4507(数位dp，较复杂）

这题麻烦就是要求数的平方和。

我们先考虑求和的问题，一个区间，数位dp能在一些约束下计数，现在要这些数的和。其实组合数学搞搞就可以了：如 现在枚举的某一位pos,我统计了这一位枚举i的满足条件的个数cnt，其实只要算i对总和的贡献就可以了，对于一个数而已第pos位是i，那么对求和贡献就是i*10^pos,就是十进制的权值，然后有cnt个数都满足第pos位是i，最后sum=cnt*i*10^pos.原理就是这样平方和可以看做(a*10^pos+b)^2,a是你当前pos位要枚举的，b其实是个自问题，就是pos之后的位的贡献值，把这个平方展开就可以了！

具体见代码计算部分：（有的地方要乘上cnt的，有的地方不用）

#include <bits/stdc++.h>
#define ll long long 
using namespace std;
const int maxn = 20;
const ll MOD = 1e9 + 7;
struct CAL{
	ll sum;
	ll sqr;
	ll cnt;
	CAL() {sum =sqr=cnt=0;}
}dp[20][10][10];
ll a[maxn];
ll pow10[maxn];
void init()
{
	pow10[0] = 1;
	for(int i = 1;i <= 18;i++){
		pow10[i] = pow10[i - 1] * 10;
	}
	memset(dp,-1,sizeof dp);
}
ll pow2(ll x){
	return (x % MOD)* (x % MOD) % MOD;
}
CAL dfs(int pos, int sum, int val,bool lead,bool limit){
	if(pos == -1) {
		CAL p;
		if(sum != 0 && val != 0) p.cnt = 1;
		else p.cnt = 0;
		return p;
	}
	if(!limit && dp[pos][sum][val].cnt != -1) return dp[pos][sum][val];
	int up = limit ? pos[a] : 9;
	CAL ans ,temp;
	for(int i = 0;i <= up;i++){
		if(i == 7) continue;
		temp = dfs(pos - 1,(sum + i) % 7,(val * 10 + i) % 7,lead && (i == 0),limit && (i == up));
		ll tmp = ((temp.cnt * pow2(i * pow10[pos])) % MOD  + (2 * (i * pow10[pos]) % MOD * temp.sum) % MOD + temp.sqr) % MOD;
		ans.cnt = (ans.cnt + temp.cnt) % MOD;
		ans.sum = (ans.sum + temp.sum + ((i * pow10[pos]) % MOD * temp.cnt) % MOD)%MOD;
		ans.sqr = (ans.sqr + tmp) % MOD;
	}
	if(!limit) dp[pos][sum][val] = ans;
	return ans;
}
ll solve(ll x){
	int pos = 0;
	while(x > 0){
		a[pos++] = x % 10;
		x /= 10;
	}
	CAL ans = dfs(pos - 1,0,0,true,true);
	return ans.sqr;
}
int main()
{
	int kase;scanf("%d",&kase);
	init();
	while(kase--){
		ll aa,bb;
		scanf("%I64d%I64d",&aa,&bb);
		printf("%I64d\n",(solve(bb) - solve(aa - 1) + MOD) % MOD);
	}
	return 0;
}