hdu 4507 - 吉哥系列故事——恨7不成妻(数位dp)

原创 2015年07月06日 19:18:55

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

Problem Description
  单身!
  依然单身!
  吉哥依然单身!
  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

Source
2013腾讯编程马拉松初赛第一场(3月21日)

思路:刚开始我居然就拿dp保存满足条件的数的平方和,然而并不应该这样。
dp[cur][s][mod1][mod2]表示当前位置为止,是否包含7,前面的数对7的余数,前面的数位的和对7的余数,这个状态下,后面的数的平方和,和,还有有多少个这三个值
那么相当于知道对于当前这个状态a21+a22+...+a2n还有n1ai,还有有多少个数满足题意条件(假设是cnt个),现在这个位置上的数是i,那么相当于求n1(ai+b)2,于是可以化简一下相当于求n1a2i+2bn1ai+cntb2,于是可得到答案

#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<vector>
#include<cmath>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<algorithm>
using namespace std;
typedef long long LL;
const int maxn=30;
const LL MOD=1e9+7;
LL N,M;
int dig[maxn];
pair<LL,pair<LL,LL>> dp[maxn][2][10][10];
bool vis[maxn][2][10][10];
LL pow10[20];
pair<LL,pair<LL,LL>> dfs(int cur,int e,int z,int s,int mod1,int mod2)
{
    if(cur<0)
    {
        if(s||(mod1%7==0)||(mod2%7==0))
            return make_pair(0,make_pair(0,0));
        return make_pair(1,make_pair(0,0));
    }
    if(!e&&vis[cur][s][mod1][mod2])
        return dp[cur][s][mod1][mod2];
    int end=(e?dig[cur]:9);
    pair<LL,pair<LL,LL>> nxt,ans=make_pair(0,make_pair(0,0));
    for(int i=0;i<=end;i++)
    {
        int tmp1=(mod1*10+i)%7;
        int tmp2=(mod2+i)%7;

        nxt=dfs(cur-1,e&&i==end,0,s||(i==7),tmp1,tmp2);

        LL presum=i*pow10[cur]%MOD;
        (ans.first+=nxt.first)%=MOD;
        (ans.second.first+=nxt.second.first+presum*nxt.first)%=MOD;
        (ans.second.second+=nxt.second.second+presum*presum%MOD*nxt.first+2*nxt.second.first*presum)%=MOD;
    }
    if(!e)dp[cur][s][mod1][mod2]=ans,
        vis[cur][s][mod1][mod2]=1;
    return ans;
}
LL solve(LL n)
{
    int len=0;
    memset(vis,0,sizeof(vis));
    while(n)
    {
        dig[len++]=n%10;
        n/=10;
    }
    return dfs(len-1,1,1,0,0,0).second.second;
}
int main()
{
    int T;
    scanf("%d",&T);
    pow10[0]=1;
    for(int i=1;i<20;i++)pow10[i]=(pow10[i-1]*10)%MOD;
    while(T--)
    {
        scanf("%I64d%I64d",&N,&M);
        LL res=solve(M)-solve(N-1);
        printf("%I64d\n",(res+MOD)%MOD);
    }
    return 0;
}
版权声明:本文为博主原创文章,未经博主允许不得转载。

HDU 4507 吉哥系列故事――恨7不成妻(数位dp&好魔性的一道好题)

题目链接:[kuangbin带你飞]专题十五 数位DP J - 吉哥系列故事――恨7不成妻题意Time Limit:500MS Memory Limit:32768KB 64bit I...
  • to_be_better
  • to_be_better
  • 2016年02月24日 03:32
  • 2310

HDU 4507 吉哥系列故事——恨7不成妻(数位DP)

题目链接:Click here~~ 题意: 中文题不解释。 解题思路: 之前做的都是统计满足那些性质的数的 count,这次直接蹦到统计 square sum 了。。。 先考虑如何统计 s...
  • dgq8211
  • dgq8211
  • 2013年12月21日 17:57
  • 2803

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

J - 吉哥系列故事――恨7不成妻  HDU - 4507  题解:数位dp 和简单的数位dp不同,这道题要算所有合法数的平方和 考虑到一个数可以写成X=ΣAi*Pi,(其中Ai为X每一位的...
  • qq_31759205
  • qq_31759205
  • 2017年02月08日 21:17
  • 340

HDU - 4507 吉哥系列故事――恨7不成妻 (数位DP)

Description   单身!   依然单身!   吉哥依然单身!   DS级码农吉哥依然单身!   所以,他生平最恨情人节,不管是214还是77,他都讨厌!      吉哥观察了214和77这...
  • u011345136
  • u011345136
  • 2014年07月25日 16:32
  • 1443

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

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4507 题意:统计在[L,R]区间内所有与7无关数字的平方和。 什么样的数和7有关呢?   如果一个...
  • chy20142109
  • chy20142109
  • 2016年03月17日 08:27
  • 186

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

Hdu 4507 吉哥系列故事——恨7不成妻求[L,R][L,R]中不满足以下任意一条的数的平方和1 整数中某一位是72 整数的每一位加起来的和是7的整数倍3 这个整数是7的整数倍;一眼看去又是一个数...
  • a1s4z5
  • a1s4z5
  • 2016年10月09日 10:15
  • 188

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

吉哥系列故事——恨7不成妻 Time Limit: 1000/500 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others) T...
  • Zhuhuangjian
  • Zhuhuangjian
  • 2014年07月21日 13:48
  • 1871

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

吉哥系列故事――恨7不成妻 Time Limit:500MS     Memory Limit:32768KB     64bit IO Format:%I64d & %I64u   D...
  • AcIsFun
  • AcIsFun
  • 2016年05月18日 21:38
  • 112

hdu 吉哥系列故事——恨7不成妻 数位DP

唉,以前就做过了,重做了一遍。居然花了4、5个小时,才做出来。各种惨不忍睹的错误o(╯□╰)o 其实,无论是求和,还是求平方和,都是较裸的数位DP。求和还是比较简单的,这里就谈谈怎么求平方和。 比...
  • ACMmaxx
  • ACMmaxx
  • 2014年02月18日 17:03
  • 521

hdu4507 吉哥系列故事——恨7不成妻 数位DP

吉哥系列故事——恨7不成妻 Time Limit: 1000/500 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others) T...
  • corncsd
  • corncsd
  • 2015年03月19日 20:52
  • 225
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:hdu 4507 - 吉哥系列故事——恨7不成妻(数位dp)
举报原因:
原因补充:

(最多只允许输入30个字)