恨7不成妻 数位dp (计数转求和)

题目
题意: 求所有满足数位中不存在7、数位和不是7的倍数、数字本身不是7的倍数的数的平方和。
思路: 难就难在平方和。
核心思路是这样的，由于递归，假设第i位之后的平方和以及统计完了。那么我们可以观察到一些规律。

时间复杂度: O(能过)
代码:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<complex>
#include<cstring>
#include<cmath>
#include<vector>
#include<map>
#include<unordered_map>
#include<list>
#include<set>
#include<queue>
#include<stack>
#define OldTomato ios::sync_with_stdio(false),cin.tie(nullptr),cout.tie(nullptr)
#define fir(i,a,b) for(int i=a;i<=b;++i) 
#define mem(a,x) memset(a,x,sizeof(a))
#define p_ priority_queue
// round() 四舍五入 ceil() 向上取整 floor() 向下取整
// lower_bound(a.begin(),a.end(),tmp,greater<ll>()) 第一个小于等于的
// #define int long long //QAQ
using namespace std;
typedef complex<double> CP;
typedef pair<int,int> PII;
typedef long long ll;
// typedef __int128 it;
const double pi = acos(-1.0);
const int INF = 0x3f3f3f3f;
const ll inf = 1e18;
const int N = 2e5+10;
const int M = 1e6+10;
const int mod = 1e9+7;
const double eps = 1e-6;
inline int lowbit(int x){ return x&(-x);}
template<typename T>void write(T x)
{
    if(x<0)
    {
        putchar('-');
        x=-x;
    }
    if(x>9)
    {
        write(x/10);
    }
    putchar(x%10+'0');
}
template<typename T> void read(T &x)
{
    x = 0;char ch = getchar();ll f = 1;
    while(!isdigit(ch)){if(ch == '-')f*=-1;ch=getchar();}
    while(isdigit(ch)){x = x*10+ch-48;ch=getchar();}x*=f;
}
int n,m,k,T;
int a[20];
#define int long long
struct node{
	int cnt; //多少个数
	int sum; //数的和
	int sqr; //平方和
	node(){cnt=-1;sum=0;sqr=0;}
	node(int a,int b,int c)
	{
		cnt=a,sum=b,sqr=c;
	}
};
int p[20];
node f[20][7][7];
int dfs(int cur,int count,int val,int &cnt,int &sum,bool limit)
{
	// cout<<f[cur][count][val].cnt<<"?\n";
	if(cur==-1) {cnt=count&&val;return 0;}
	if(!limit&&f[cur][count][val].cnt>0) 
	{
		cnt = f[cur][count][val].cnt;
		sum = f[cur][count][val].sum;
		return f[cur][count][val].sqr;
	}
	int up = limit?a[cur]:9;
	int ans = 0;
	for(int i=0;i<=up;++i)
	{
		if(i==7) continue;
		int cnt2 = 0;
		int sum2 = 0;
		int base = p[cur]; //cur位应该乘10的多少次方
		ans += dfs(cur-1,(count+i)%7,(val*10+i)%7,cnt2,sum2,limit&&i==up); ans %= mod;
		(ans += cnt2*i%mod*base%mod*i%mod*base%mod + 2*i*base%mod*sum2%mod)%=mod;
		cnt = (cnt + cnt2) % mod;
		sum = (sum + cnt2*i%mod*base%mod + sum2 % mod) % mod;
	}
	if(!limit) f[cur][count][val] = node(cnt,sum,ans);
	return ans;
}
int fun(int x)
{
	int pos;
	for(pos=0;x;x/=10) a[pos++]=x%10;
	int l = 0,r = 0;
	return dfs(pos-1,0,0,l,r,1)%mod;
}
void solve()
{
   int l,r;
   cin>>l>>r;
   ll ans = fun(r)-fun(l-1); ans %= mod; ans = (ans+mod)%mod;
   cout<<ans<<"\n";
}
signed main(void)
{   
   p[0] = 1;
   for(int i=1;i<=19;++i) p[i] = p[i-1]*10%mod;
   // T = 1;
   OldTomato; cin>>T;
   // read(T);
   while(T--)
   {
   	 solve();
   }
   return 0;
}