[HDU4507]吉哥系列故事——恨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)。

输出

请计算[L,R]中和7无关的数字的平方和,并将结果对10^9 + 7 求模后输出。

输入示例


输出示例


数据规模及约定

见“输入

题解

这题在状态设计上挺普通的,就是需要记好几个东西:设状态 (i, j, k, l) 表示 i 位数,最高位为 j,数字和对 7 取模等于 k,数对 7 取模等于 l。num(i, j, k, l) 记录个数,sum(i, j, k, l) 记录总和,f(i, j, k, l) 记录平方和,然后转移时推推式子就好了。

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cctype>
#include <algorithm>
using namespace std;
#define LL long long LL read() {
LL x = 0, f = 1; char c = getchar();
while(!isdigit(c)){ if(c == '-') f = -1; c = getchar(); }
while(isdigit(c)){ x = x * 10 + c - '0'; c = getchar(); }
return x * f;
} #define maxn 20
#define MOD 1000000007
LL f[maxn][10][7][7], sum[maxn][10][7][7], N[maxn][10][7][7], ten[maxn];
bool has[maxn][10][7][7]; int num[maxn];
LL cal(LL x) {
if(!x) return 0;
int cnt = 0; LL tx = x;
while(x) num[++cnt] = x % 10, x /= 10;
LL ans = 0, sm = 0, s = 0, Sum = 0; int sn = 0, aa = 0;
for(int i = cnt - 1; i; i--)
for(int j = 1; j <= 9; j++)
for(int k1 = 1; k1 < 7; k1++)
for(int k2 = 1; k2 < 7; k2++)
(ans += f[i][j][k1][k2]) %= MOD,
(aa += N[i][j][k1][k2]) %= MOD,
(Sum += sum[i][j][k1][k2]) %= MOD;
bool ok = 1;
for(int i = cnt; i; i--) {
for(int j = i < cnt ? 0 : 1; j < num[i]; j++)
for(int k1 = 1; k1 < 7; k1++)
for(int k2 = 1; k2 < 7; k2++) {
int K1 = ((k1 - sn) % 7 + 7) % 7,
K2 = (((LL)k2 - s) % 7 + 7) % 7;
ans += f[i][j][K1][K2] + (sm * sum[i][j][K1][K2] << 1) + (sm * sm % MOD) * N[i][j][K1][K2];
ans %= MOD;
(Sum += sum[i][j][K1][K2] + sm * N[i][j][K1][K2]) %= MOD;
(aa += N[i][j][K1][K2]) %= MOD;
}
(sm += ten[i-1] * num[i]) %= MOD; s += ten[i-1] * num[i]; sn += num[i];
if(num[i] == 7){ ok = 0; break; }
}
LL ttx = tx % MOD;
if(sn % 7 && s % 7 && ok) (ans += ttx * ttx) %= MOD, (Sum += ttx) %= MOD, aa++;
// printf("%lld: %d %lld %lld\n", tx, aa, Sum, ans);
return ans;
} int main() {
ten[0] = 1;
for(int i = 1; i < maxn - 1; i++) ten[i] = ten[i-1] * 10;
N[0][0][0][0] = 1; has[0][0][0][0] = 1;
for(int i = 0; i < maxn - 1; i++)
for(int j = 0; j <= 9; j++)
for(int k = 0; k < 7; k++)
for(int l = 0; l < 7; l++) if(has[i][j][k][l]) {
for(int x = 0; x <= 9; x++) if(x != 7) {
LL t = x * ten[i], tt = t % MOD;
has[i+1][x][(k+x)%7][(t+l)%7] = 1;
(N[i+1][x][(k+x)%7][(t+l)%7] += N[i][j][k][l]) %= MOD;
(sum[i+1][x][(k+x)%7][(t+l)%7] += sum[i][j][k][l] + tt * N[i][j][k][l]) %= MOD;
f[i+1][x][(k+x)%7][(t+l)%7] += f[i][j][k][l] + (tt * tt % MOD) * N[i][j][k][l] + (tt * sum[i][j][k][l] << 1);
f[i+1][x][(k+x)%7][(t+l)%7] %= MOD;
}
// printf("%d %d %d %d: %lld\n", i, j, k, l, sum[i][j][k][l]);
}
int T = read();
while(T--) {
LL l = read(), r = read();
printf("%lld\n", (cal(r) - cal(l - 1) + MOD) % MOD);
} return 0;
}

[HDU4507]吉哥系列故事——恨7不成妻的更多相关文章

  1. hdu4507吉哥系列故事——恨7不成妻 (数位dp)

    Problem Description 单身! 依然单身! 吉哥依然单身! DS级码农吉哥依然单身! 所以,他生平最恨情人节,不管是214还是77,他都讨厌! 吉哥观察了214和77这两个数,发现: ...

  2. HDU-4507 吉哥系列故事——恨7不成妻 数位DP

    题意:给定区间[L, R]求区间内与7无关数的平方和.一个数当满足三个规则之一则认为与7有关:1.整数中某一位是7:2.整数的每一位加起来的和是7的整数倍:3.这个整数是7的整数倍: 分析:初看起来确 ...

  3. hdu-4507 吉哥系列故事——恨7不成妻 数位DP 状态转移分析/极限取模

    http://acm.hdu.edu.cn/showproblem.php?pid=4507 求[L,R]中不满足任意条件的数的平方和mod 1e9+7. 条件: 1.整数中某一位是7:2.整数的每一 ...

  4. 2018.09.27 hdu4507吉哥系列故事——恨7不成妻(数位dp)

    传送门 一道比较综合的数位dp. 维护三个值:[L,R][L,R][L,R] 区间中与7无关的数的数量,与7无关的数之和,与7无关的数的的平方和. 然后可以用第一个值推第二个,第一个和第二个值推第三个 ...

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

    这题面什么垃圾玩意儿 首先看到问题格式想到数位DP,但是求的是平方和.尝试用数位DP推出. 先尝试拼出和.设$f[len][sum][mod]$表示填到$len$位,已填位置数位和$sum$,数字取余 ...

  6. HDU4507 吉哥系列故事——恨7不成妻 题解 数位DP

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4507 题目大意: 找到区间 \([L,R]\) 范围内所有满足如下条件的数的 平方和 : 不包含'7' ...

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

    吉哥系列故事——恨7不成妻 http://acm.hdu.edu.cn/showproblem.php?pid=4507 Time Limit: 1000/500 MS (Java/Others)   ...

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

    吉哥系列故事--恨7不成妻 传送门 Problem Description 单身! 依然单身! 吉哥依然单身! DS级码农吉哥依然单身! 所以,他生平最恨情人节,不管是214还是77,他都讨厌! 吉哥 ...

  9. HDU - 4507 - 吉哥系列故事——恨7不成妻(数位DP,数学)

    链接: https://vjudge.net/problem/HDU-4507 题意: 单身! 依然单身! 吉哥依然单身! DS级码农吉哥依然单身! 所以,他生平最恨情人节,不管是214还是77,他都 ...

随机推荐

  1. WinForm------关于子窗体刷新父窗体问题

    链接: http://wenwen.sogou.com/z/q242758397.htm

  2. spark 加载文件

    spark 加载文件 textFile的参数是一个path,这个path可以是: 1. 一个文件路径,这时候只装载指定的文件 2. 一个目录路径,这时候只装载指定目录下面的所有文件(不包括子目录下面的 ...

  3. mysql_query()与mysql_real_query()

    mysql_query() cannot be used ” character, which mysql_query() interprets as the end of the statement ...

  4. 对象Clone

    //================================================= // File Name : Clone_demo //-------------------- ...

  5. PyCharm 教程(五)断点 调试

    PyCharm 教程(五)断点 调试 PyCharm 作为IDE,断点调试是必须有的功能.否则,我们还真不如用纯编辑器写的快. [运行]和[调试]前的设置,详见前面的文章,helloword. 1,设 ...

  6. C#----操作应用程序配置文件App.config

    对配置文件的一些疑问: 在应用程序的目录下,有两处值得注意的地方,一个是应用程序根目录下的App.config文件,和bin\debug\name.exe.config 或者 bin\Release\ ...

  7. Cloudservie将LocalStroage中的内容通过WAD自动上传到BLOB中

    开发云服务程序,如果使用Local stroage存储我们临时生成的日志或者文件并将它们自动上传到BLOB中,可以通过WAD来实现,具体如下: 1.配置webrole,开启Local stroage功 ...

  8. 缓存 HttpContext.Current.Cache和HttpRuntime.Cache的区别

    先看MSDN上的解释: HttpContext.Current.Cache:为当前 HTTP 请求获取Cache对象. HttpRuntime.Cache:获取当前应用程序的Cache.  我们再用. ...

  9. span设为inline-block之后,未包含文字时下面会多出一条空白问题

    1.问题的引出: 产品列表页面场景: 上面是产品图片[img], 中间是提示库存信息[span](始终存在,有库存则不显示文字,但元素占位.所以设置display:inline-block), 下面是 ...

  10. ajava包的命名

    2. Package的命名 Package名的第一部分应是小写ASCII字符,并且是顶级域名之一,通常是com.edu.gov.mil.net.org或由ISO标准3166.1981定义的国家唯一标志 ...