首页 > 搜索 > DFS搜索 > HDU 4507-吉哥系列故事――恨7不成妻-动态规划-[解题报告]HOJ
2015
07-17

HDU 4507-吉哥系列故事――恨7不成妻-动态规划-[解题报告]HOJ

吉哥系列故事――恨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无关的数字的平方和。

输入:

输入数据的第一行是case数T(1 <= T <= 50),然后接下来的T行表示T个case;每个case在一行内包含两个正整数L, R(1 <= L <= R <= 10^18)。

输出:

输入数据的第一行是case数T(1 <= T <= 50),然后接下来的T行表示T个case;每个case在一行内包含两个正整数L, R(1 <= L <= R <= 10^18)。

样例输入:

3
1 9
10 11
17 17

样例输出:

236
221
0

题目链接:Click here~~

题意:

中文题不解释。

解题思路:

之前做的都是统计满足那些性质的数的 count,这次直接蹦到统计 square sum 了。。。

先考虑如何统计 sum。

统计 sum 维护两个值 count 和 sum 就可以了。

想象状态转移时,相当于在一个具有同样性质的后缀的数的集合前面,加某个数字。

那么,sum[new_state] = Σ{ sum[old_state] + (number to add at the postion) * (its base) * count[old_state] }。(语言不好描述,看式子吧。。。)

类似的,维护三个值 count & sum & square sum 可以解决这个问题了。自己想一下吧。O(∩_∩)O。

启发:原式是 a[1]^2 + a[2]^2 + … + a[n]^2,新式是 (a[1]+b)^2 + (a[2]+b)^2 + … + (a[n]+b)^2。

#include <math.h>
#include <string>
#include <stdio.h>
#include <string.h>

using namespace std;

typedef long long LL;
typedef pair< int,pair<LL,LL> > PILL;

#define F first
#define S second
#define MP make_pair

const int mod = 1e9 + 7;

int digit[20];

LL pow10[20];

PILL dp[20][7][7][2];

bool vis[20][7][7][2];

PILL dfs(int len,int sum,int remain,bool contain,bool fp)
{
    if(!len)
        if(!contain && sum && remain)
            return MP(1,MP(0LL,0LL));
        else
            return MP(0,MP(0LL,0LL));
    if(!fp && vis[len][sum][remain][contain])
        return dp[len][sum][remain][contain];
    PILL ret = MP(0,MP(0,0));
    int fpmax = fp ? digit[len] : 9;
    for(int i=0;i<=fpmax;i++)
    {
        PILL nxt = dfs(len-1,(sum + i) % 7,(remain * 10 + i) % 7,contain | (i == 7),fp && i == fpmax);
        LL pref = i * pow10[len-1] % mod;
        (ret.F += nxt.F) %= mod;
        (ret.S.F += nxt.S.F + pref * nxt.F) %= mod;
        (ret.S.S += nxt.S.S + pref * pref % mod * nxt.F + 2 * pref * nxt.S.F) %= mod;
    }
    if(!fp)
    {
        vis[len][sum][remain][contain] = true;
        dp[len][sum][remain][contain] = ret;
    }
    return ret;
}

long long f(long long n)
{
    int len = 0;
    while(n)
    {
        digit[++len] = n % 10;
        n /= 10;
    }
    return dfs(len,0,0,0,true).S.S;
}

int main()
{
    pow10[0] = 1;
    for(int i=1;i<20;i++)
        pow10[i] = pow10[i-1] * 10 % mod;
    int T;
    scanf("%d",&T);
    while(T--)
    {
        long long a,b;
        scanf("%I64d%I64d",&a,&b);
        printf("%I64d\n",(f(b) - f(a-1) + mod) % mod);
    }
    return 0;
}

版权声明:本文为博主原创文章,未经博主允许不得转载。

参考:http://blog.csdn.net/dgq8211/article/details/17469457