不要666升级版（数位DP，三次方和）

不要666升级版

∑(pre+suc)2n×pre2+2×pre∑suc+∑suc2\sum(pre + suc) ^ 2\\ n \times pre ^ 2 + 2 \times pre \sum suc + \sum suc ^ 2\\ ∑(pre+suc)2n×pre2+2×pre∑suc+∑suc2

∑(pre+suc)3∑(pre3+3×pre2×suc+3×pre×suc2+suc3)n×pre3+3×pre2∑suc+3×pre∑suc2+∑suc3\sum (pre + suc) ^ 3\\ \sum \left(pre ^ 3 + 3 \times pre ^ 2 \times suc + 3 \times pre \times suc ^ 2 + suc ^ 3\right)\\ n \times pre ^ 3 + 3 \times pre ^ 2 \sum suc + 3 \times pre \sum suc ^ 2 + \sum suc ^ 3\\ ∑(pre+suc)3∑(pre3+3×pre2×suc+3×pre×suc2+suc3)n×pre3+3×pre2∑suc+3×pre∑suc2+∑suc3

由此维护n,sum,sum2,sum3n, sum, sum^ 2, sum ^ 3n,sum,sum2,sum3即可，定义f[i][j][k]f[i][j][k]f[i][j][k]，为枚举到第iii为，数位和为jjj，数字余kkk的状态，记忆化搜索即可。

#include <bits/stdc++.h>using namespace std;const int mod = 1e9 + 7;inline int add(int x, int y) {return x + y < mod ? x + y : x + y - mod;
}inline int sub(int x, int y) {return x >= y ? x - y : x - y + mod;
}int cnt, num[20], p[20];struct Res {int cnt = -1, sum1, sum2, sum3;
}f[20][10][10];Res dfs(int pos, int sum, int remain, int flag) {if (!pos) {if (sum != 0 && remain != 0) {return {1, 0, 0, 0};}else {return {0, 0, 0, 0};}}if (!flag && f[pos][sum][remain].cnt != -1) {return f[pos][sum][remain];}Res ans = {0, 0, 0, 0};int nex = flag ? num[pos] : 9;for (int i = 0; i <= nex; i++) {if (i == 6) {continue;}Res cur = dfs(pos - 1, (sum + i) % 6, (remain * 10 + i) % 6, i == nex && flag);int s = 1ll * i * p[pos] % mod;ans.cnt = add(ans.cnt, cur.cnt);ans.sum1 = add(ans.sum1, add(1ll * s * cur.cnt % mod, cur.sum1));ans.sum2 = add(add(ans.sum2, 1ll * cur.cnt * s % mod * s % mod), add(cur.sum2, 2ll * s % mod * cur.sum1 % mod));ans.sum3 = add(ans.sum3, add(add(1ll * cur.cnt * s % mod * s % mod * s % mod, 3ll * s % mod * s % mod * cur.sum1 % mod), add(3ll * s % mod * cur.sum2 % mod, cur.sum3)));}if (!flag) {f[pos][sum][remain] = ans;}return ans;
}int calc(long long x) {cnt = 0;while (x) {num[++cnt] = x % 10;x /= 10;}return dfs(cnt, 0, 0, 1).sum3;
}int main() {// freopen("in.txt", "r", stdin);// freopen("out.txt", "w", stdout);p[1] = 1;for (int i = 2; i < 20; i++) {p[i] = 1ll * p[i - 1] * 10 % mod;}long long l, r;while (scanf("%lld %lld", &l, &r) != EOF) {printf("%lld\n", sub(calc(r), calc(l - 1)));}return 0;
}

不要666升级版（数位DP，三次方和）相关推荐

输入一个三位数，求个个数位数字的三次方之和，并判断和是否与该数相同
解题思路: 先把三位数(n)的3个数位值算出来然后求数位数字的三次方之和最后判断数位数字的三次方之和(sum)是否与该数(n)相同代码: # 输入一个三位数,求个个数位数字的三次方之和,并判断和 ...
牛客网暑期ACM多校训练营(第三场): C. Chiaki Sequence Reloaded（数位DP）
题目描述 Chiaki is interested in an infinite sequence a1, a2, a3, ..., which defined as follows: Chiaki ...
bzoj 4318 osu 【三次方问题】【dp】
二次方的问题可以说是经典了,三次方yy一下也就出来了. (x+1)^3=x^3+3x^2+3x+1 所以维护x^2 x的期望值.注意x^2的期望值不等于x期望值的平方. #include<cst ...
动态规划专题（三）——数位DP
前言数位DPDPDP 真的是最恶心的DPDPDP. 简介看到那种给你两个数,让你求这两个数之间符合条件的数的个数,且这两个数非常大,这样的题目一般就是数位DPDPDP 题. 数位DPDPDP一般 ...
数位 DP 入门 (不要 62+windy 数)
\[I\] 平常的做法是设 \(f_{i,j}\) 为 \(0\)~\(j \times 10^{i-1}\) 的合法个数,这里用某种神奇而快速的做法. 简化题意: 不要 \(6\ 2\) 连在一起的 ...
SPOJ - BALNUM Balanced Numbers(数位dp+进制转换)
题目链接:点击查看题目大意:给出平衡数的定义:每一个偶数出现的次数必须是奇数次,每一个奇数出现的次数必须是偶数次,求给定区间中有多少个平衡数题目分析:数位dp,这个题目就难在怎么确定状态转移,本来 ...
AcWing1081.度的数量(数位DP)题解
题目传送门求给定区间 [X, Y] 中满足下列条件的整数个数:这个数恰好等于 K 个互不相等的 B 的整数次幂之和. 例如,设 X=15,Y=20,K=2,B=2X=15,Y=20,K=2,B=2, ...
数位DP———POJ 3208 启示录
POJ3208 启示录题目传送门题目大意:我们称只要某数字的十进制表示中有三个连续的6,我们就称它为"魔鬼数",比如:666,1666,6663,16666 等.现给出一个数X ...
数位 dp 相邻位数字差值的绝对值不能超过 2_干货来了！13种行测数字推理技巧总结...
数学除了公式要记忆,解题技巧同样需要"记忆"所谓的解题技巧来源于实战,运用于实战它是在通过大量做题.充分理解和把握命题规律的基础上形成的解题方法论今天陕西京佳教育为大家盘点13种 ...

不要666升级版（数位DP，三次方和）

不要666升级版

不要666升级版（数位DP，三次方和）相关推荐

最新文章

热门文章