(数位dp)吉哥系列故事――恨7不成妻

吉哥系列故事――恨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

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn=20;
const int md=1e9+7;
ll l,r,c[maxn];
int arr[maxn];
struct node
{
    ll sum;
    ll cnt;
    ll sqsum;
    node() {}
    node(ll a,ll b,ll c):cnt(a),sum(b),sqsum(c){}
}dp[maxn][500][10];
node dfs(int pos,int sum,int mod,bool limit)
{
    if(pos==0)
    {
        return node(sum!=0&&mod!=0,0,0);
    }
    if(!limit&&dp[pos][sum][mod].cnt!=-1) return dp[pos][sum][mod];
    int up=limit?arr[pos]:9;
    node ans=node(0,0,0);
    for(int i=0;i<=up;i++)
    {
        if(i==7) continue;
        node t=dfs(pos-1,(sum+i)%7,(mod*10+i)%7,limit&&i==up);
        ans.cnt=(ans.cnt+t.cnt)%md;
        ans.sum+=(((c[pos-1]*i)%md*t.cnt)%md+t.sum)%md;
        ans.sum%=md;
        ans.sqsum+=(((2*c[pos-1]*i)%md*t.sum)%md+t.sqsum)%md;
        ans.sqsum%=md;
        ans.sqsum+=((i*c[pos-1]*i%md)*c[pos-1]%md*t.cnt)%md;
        ans.sqsum%=md;
    }
    if(!limit) dp[pos][sum][mod]=ans;
    return ans;
}
ll solve(ll x)
{
    ll pos=0;
    while(x)
    {
        arr[++pos]=x%10;
        x/=10;
    }
    return dfs(pos,0,0,1).sqsum;
}
int main()
{
    int t,num=0;
    scanf("%d",&t);
    c[0]=1;
    for(int i=1;i<=20;i++) c[i]=(c[i-1]*10)%md;
    memset(dp,-1,sizeof(dp));
    while(t--)
    {
        scanf("%lld%lld",&l,&r);
        printf("%lld\n",(solve(r)-solve(l-1)+md)%md);
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值