HihoCoder 1033交错和(数位DP第三题)
描述
给定一个数 x,设它十进制展从高位到低位上的数位依次是 a0, a1, ..., an - 1,定义交错和函数:
f(x) = a0 - a1 + a2 - ... + ( - 1)n - 1an - 1
例如:
f(3214567) = 3 - 2 + 1 - 4 + 5 - 6 + 7 = 4
给定

输入
输入数据仅一行包含三个整数,l, r, k(0 ≤ l ≤ r ≤ 1018, |k| ≤ 100)。
输出
输出一行一个整数表示结果,考虑到答案可能很大,输出结果模 109 + 7。
提示
对于样例 ,满足条件的数有 110 和 121,所以结果是 231 = 110 + 121。
更多样例:
| Input |
| 4344 3214567 3 |
| Output |
| 611668829 |
| Input |
| 404491953 1587197241 1 |
| Output |
| 323937411 |
| Input |
| 60296763086567224 193422344885593844 10 |
| Output |
| 608746132 |
| Input |
| 100 121 -1 |
| Output |
| 120 |
- 样例输入
-
100 121 0
- 样例输出
-
231
#include<iostream>
using namespace std; #define ll long long int //需要用long long 来进行存放
const int mod = ;//结果的最大值,用于取模
struct node{
ll s, n;//s: 各数位交错和 n:满足条件的个数
};
node dp[][][];//[长度][以哪个数作为开头][该数的数位交错和]
int bits[];//存放一个数的各位,低到高存放在0到20
ll base[];//每一数位的基准 比如十位数=base[2]=10 //len数位长度, dig是首个数字, begin_zero表示从最高位到当前位是否全部为0, limit表示下一位枚举是否有限制(bit[len-2]或9), sum是要求的数字和
node dfs(int len, int dig, bool begin_zero, bool limit, int sum){
node t;//存放长度为 len的结果
t.s = , t.n = ;
//超过边界值
if (len <= || len >= || dig < || dig > || sum < - || sum >= )
return t;
//返回已有的DP结果,即记忆化搜索
if (!limit && dp[len][dig + (begin_zero ? : )][sum + ].n != -)
return dp[len][dig + (begin_zero ? : )][sum + ];
//长度只有一位,就不需要枚举下一位了,直接讨论返回即可
if (len == ){
if (dig != sum)
return t;
t.n = , t.s = sum;
return t;
}
//开始枚举下一位的数字
int end = limit ? bits[len - ] : ;//下一位数字的最大值
int newsum = dig - sum;
node tmp;
for (int j = ; j < end + ; j++)
{
if (begin_zero){//前面都是0,接下来的就由当前位决定 j是否为0
tmp = dfs(len - , j, j == , limit && (j == end), sum);
}
else{//前面不是全为0
tmp = dfs(len - , j, false, limit && (j == end), newsum);
}
//将tmp的值累加到t上
t.n += tmp.n;//满足条件的个数
//计算满足条件的长度为 len 的data[len]所有数的和, tmp 是长度为 i-1的data[len-1]所有数的和
//每一个: data[len] = dig * base[len] + data[len-1] 且共有n个
//t.s = t.s + tmp.n * (dig * base[len]) + tmp.s 增加 %mod 就得到下面的结果
t.s = ((t.s + tmp.s) % mod + ((tmp.n * dig) % mod * base[len]) % mod) % mod;
}
//当长度为len,且以dig开头的,数位和为sum 的所有结果都计算完成,才将其进行存储
if (!limit)
dp[len][dig + (begin_zero ? : )][sum + ] = t;
//dig + (begin_zero ? 0 : 10) 用来区分两种状态,(1. 前导都为0;2. 前导包含其它数)
return t;
} int solve(ll n, int s){
if (n <= )
return ;
int l = ;
for (int i = ; i < ; i++)
bits[i] = ;
//将n的每一位从低到高放到 bits[0]到bits[l] 中
while (n){
bits[l++] = n % ;
n /= ;
}
//从l+1开始,比n的长度大1,并且第l+1位数置为0
return dfs(l + , , true, true, s).s;
} int main(){
ll l, r, s;
node t;
t.n = -;
t.s = ;
for (int i = ; i < ; i++)//长度
for (int j = ; j < ; j++)//第i位的取值
for (int k = ; k < ; k++)// i到0 所有位上的交错和 + 200
dp[i][j][k] = t;
base[] = ;
for (int i = ; i < ; i++)//base 是作为基准数 比如两位数需要 *10, 三位数需要*100
base[i] = base[i - ] * % mod;
cin >> l >> r >> s;
cout << (solve(r, s) - solve(l - , s) + mod) % mod << endl;
return ;
}
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<memory>
#include<cstring>
#include<cmath>
using namespace std;
#define LL long long
const LL Mod=;
LL a[],cnt,L,R,S;
LL ans,dp[][][][][]; LL dfs(LL pos,LL sum,LL limit,LL sign,LL fir,LL sub)
{
if(pos==) return sum==S?sub:;
if(!limit&&sign&&dp[pos][S-sum+][][][]) {
for(int i=;i<=pos;i++) sub=sub*%Mod;
sub=(sub+dp[pos][S-sum+][][][])%Mod;
return sub;
}
LL Up=limit?a[pos]:;
LL tmp=;
for(LL i=;i<=Up;i++){
if(fir&&i==) continue;
tmp=(tmp%Mod+dfs(pos-,sum+(sign?:-)*i,limit&&i==Up,-sign,,sub*%Mod+i))%Mod;
}
dp[pos][sum][limit][sign][fir]=tmp;
return tmp;
} LL cal(int v)
{
cnt=;ans=;
while(v){
a[++cnt]=v%;
v/=;
}
for(LL i=;i<=cnt;i++){
LL tmp=dfs(i,,i==cnt,,,);//枚举第一位 ,避开前导0
ans+=tmp;
}
return ans;
} int main()
{
scanf("%lld%lld%lld",&L,&R,&S);
S+=;
printf("%lld",(cal(R)-cal(L-)+Mod)%Mod);
return ;
}
HihoCoder 1033交错和(数位DP第三题)的更多相关文章
- [hihocoder 1033]交错和 数位dp/记忆化搜索
#1033 : 交错和 时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描写叙述 给定一个数 x,设它十进制展从高位到低位上的数位依次是 a0, a1, ..., an - 1 ...
- hihoCoder 1033 : 交错和 数位dp
思路:数位dp,dp(i, j, k)表示考虑i位数,每位数可以任意取[0~9],并且这i位数的交错和为j,k=1表示前缀全是0(如000456),k=0表示前缀不为0.注意,前缀是否为0是这道题的一 ...
- hihoCoder1033 交错和 数位DP
题目:交错和 链接:http://hihocoder.com/problemset/problem/1033# 题意:对于一个十进制整数x,令a0.a1.a2.....an是x从高位到低位的数位,定义 ...
- hihoCoder 1033: 交错和
(1)题目描述: 时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 给定一个数 x,设它十进制展从高位到低位上的数位依次是 a0, a1, ..., an - 1,定义交错 ...
- hihoCoder #1033 : 交错和 (数位Dp)
题目大意: 给定一个数 x,设它十进制展从高位到低位上的数位依次是 a0, a1, ..., an - 1,定义交错和函数: f(x) = a0 - a1 + a2 - ... + ( - 1)n - ...
- HihoCoder1050 树中的最长路 树形DP第三题(找不到对象)
题意:求出的树中距离最远的两个结点之间相隔的距离. 水题一道,以前只会用路的直径来解. 代码如下: #include<cstdio> #include<cstdlib> #in ...
- hdu 3555数位dp基础入门题
#include<stdio.h> #define N 20 long long dp[N][3]; void init(){ long long i; dp[0][0]=1; for ...
- HDU 2089 不要62(数位dp入门)
题意:统计区间 [a,b] 中不含 4 和 62 的数字有多少个. 题解:这是数位DP的入门题了,首先要理解数DP的原理,DP[i][j]:代表第i位的第j值,举个栗子:如4715 数位数是从右向 ...
- bzoj3209 花神的数论题——数位dp
题目大意: 花神的题目是这样的 设 sum(i) 表示 i 的二进制表示中 1 的个数.给出一个正整数 N ,花神要问你 派(Sum(i)),也就是 sum(1)—sum(N) 的乘积. 要对1000 ...
随机推荐
- HighCharts 在IE8下饼图不显示的问题
HighCharts饼图用来做数据统计时,在IE8下发现某些图形不能正确显示出来. 在IE8下面会报 'this.renderer.gradients' 为空或不是对象 这样的错误.. 解决方法: ...
- 51Nod 1521 一维战舰
http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1521 思路:先计算出一开始最多能放多少艘战舰,然后每次输入一个点后,找到 ...
- 对dataframe中某一列进行计数
本来是一项很简单的任务...但很容易忘记搞混..所以还是记录一下 方法一: df['col'].value_counts() 方法二: groups = df.groupby('col') group ...
- IDEA.环境
1.下载: https://www.jetbrains.com/idea/download/#section=windows 安装的版本是:Ultimate:ideaIU-2017.3.3.exe 2 ...
- for each/in/of的解释and example
for-of 循环:代码示例for (var value of myArray) {console.log(value);}循环的对象需为一个数组 无法记录索引 可以相应break.continue. ...
- 整合Struts2与Spring
一.需要的JAR文件为:Spring和Struts2框架本身需要的JAR文件以及他们所依赖的JAR文件
- mybatis定义拦截器
applicationContext.xml <bean id="sqlSessionFactory" class="org.mybatis.spring.SqlS ...
- 二十、dbms_stats(用于搜集,查看,修改数据库对象的优化统计信息)
1.概述 作用:用于搜集,查看,修改数据库对象的优化统计信息. 2.包的组成 1).get_column_stats作用:用于取得列的统计信息语法:dbms_stats.get_column_stat ...
- js上传图片&预览(filereader)
fileReader HTML5定义了FileReader作为文件API的重要成员用于读取文件,根据W3C的定义,FileReader接口提供了读取文件的方法和包含读取结果的事件模型. FileRea ...
- ant+jmeter+jenkins+git持续集成以及邮件报告展示
前序准备工作: ant--下载地址:http://ant.apache.org/bindownload.cgi jmeter--下载地址:http://jmeter.apache.org/downlo ...