【BZOJ】1833 [ZJOI2010]count 数字计数
【算法】数位DP
【题解】
记忆化搜索
#include<cstdio>
#include<algorithm>
#include<cstring>
#define ll long long
using namespace std;
ll A[],B[],f[][],a[],p[];
ll dfs(ll* A,ll h,bool limit,bool pre)
{
if(h==)return ;
if(!limit&&f[h][]!=-&&!pre)
{
for(int i=;i<=;i++)A[i]+=f[h][i];
return f[h][];
}
int end=limit?a[h]:;ll ans=;
for(int i=;i<=end;i++)
{
if(limit&&i==end)p[i]=dfs(A,h-,,pre&(!i)),A[i]+=p[i],ans+=p[i];
else p[i]=dfs(A,h-,,pre&(!i)),A[i]+=p[i],ans+=p[i];
if(pre&(!i))A[i]-=p[i];
}
if(!limit&&!pre)
{
for(int i=;i<=;i++)f[h][i]=f[h-][i]*+p[i];
f[h][]=ans;
}
return ans;
}
void solve(ll* A,ll x)
{
if(x==){return;}
int n=;
while(x>)a[++n]=x%,x/=;
dfs(A,n,,);
}
int main()
{
ll l,r;
scanf("%lld%lld",&l,&r);
for(int i=;i<=;i++)f[i][]=-;
solve(B,l-);solve(A,r);
for(int i=;i<;i++)printf("%lld ",A[i]-B[i]);
printf("%lld\n",A[]-B[]);
return ;
}
递推:下面只考虑单一数字数量统计,其它一样。
第一步,预处理。(计数)
规定最低位为第1位,最高位为第len位。
f[i][j]表示前i位,最高位数字为j的答案数(不考虑前导零有前导零的数字只要待会再最高位附上数字就是有效的了)
f[i][j]=∑f[i-1][k]+10^(i-1)。
第二步,第len位为0。(前导零)
将所有len位为0的或数字长度不足len的数字先统计进来,ans+=∑f[i][j]+1,1<=i<=len-1,1<=j<=9。其中+1是数字0。
第三步,第len位不为0逐位确定。(限位)
强制确定len位不为0,然后加入每一位枚举到顶-1的答案就可以了。
ans+=∑f[i][j],1<=i<len,0<=j<=a[i],其中i=len时j从1开始。
还要计算当前数位的顶之后会出现的次数,为后面数字大小+1。
例如1211这个数字,整个数位DP的过程是:0,1~9,10~99,100~999,1000~1199,1200~1209,1210~1211。
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const int N=,M=;
int a[M];
ll fac[M];
struct cyc{
ll a[];
}f[N][N];
cyc operator + (cyc a,cyc b){
cyc c;
for(int i=;i<=;i++)c.a[i]=a.a[i]+b.a[i];
return c;
}
cyc dp(ll num){
int len=;
cyc ans;
for(int i=;i<=;i++)ans.a[i]=;
ans.a[]=;
if(!num)return ans;
ll number=num;
while(num){
a[++len]=num%;
num/=;
}
for(int i=;i<len;i++)for(int j=;j<=;j++)ans=ans+f[i][j];
for(int i=len;i>=;i--){
for(int j=(i==len);j<a[i];j++){
ans=ans+f[i][j];
}
number%=fac[i-];
ans.a[a[i]]+=number+;
}
//for(int k=0;k<=9;k++)printf("%d ",ans.a[k]);puts("");
return ans;
}
int main(){
fac[]=;for(int i=;i<=N;i++)fac[i]=fac[i-]*;
for(int j=;j<=;j++)f[][j].a[j]=;
for(int i=;i<=N;i++){
for(int j=;j<=;j++){
for(int k=;k<=;k++){
f[i][j]=f[i][j]+f[i-][k];
f[i][j].a[j]+=fac[i-];
}//printf("[%d][%d]",i,j);
//for(int k=0;k<=9;k++)printf("%d ",f[i][j].a[k]);puts("");
}
}
ll A,B;scanf("%lld%lld",&A,&B);
cyc cA=dp(A-),cB=dp(B);
for(int i=;i<;i++)printf("%lld ",cB.a[i]-cA.a[i]);
printf("%lld",cB.a[]-cA.a[]);
return ;
}
【BZOJ】1833 [ZJOI2010]count 数字计数的更多相关文章
- BZOJ 1833: [ZJOI2010]count 数字计数( dp )
dp(i, j, k)表示共i位, 最高位是j, 数字k出现次数. 预处理出来. 差分答案, 对于0~x的答案, 从低位到高位进行讨论 -------------------------------- ...
- [BZOJ 1833] [ZJOI2010] count 数字计数 【数位DP】
题目链接:BZOJ - 1833 题目分析 数位DP .. 用 f[i][j][k] 表示第 i 位是 j 的 i 位数共有多少个数码 k . 然后差分询问...Get()中注意一下,如果固定了第 i ...
- BZOJ 1833: [ZJOI2010]count 数字计数
Description 问 \([L,R]\) 中0-9的个数. Sol 数位DP. 预处理好长度为 \(i\), 最高位为 \(j\) 的数位之和. 然后从上往下计算,不要忘记往下走的同时要把高位的 ...
- bzoj 1833 [ZJOI2010]count 数字计数(数位DP)
[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=1833 [题意] 统计[a,b]区间内各数位出现的次数. [思路] 设f[i][j][k ...
- BZOJ 1833 ZJOI2010 count 数字计数 数位DP
题目大意:求[a,b]间全部的整数中0~9每一个数字出现了几次 令f[i]为i位数(算前导零)中每一个数出现的次数(一定是同样的,所以仅仅记录一个即可了) 有f[i]=f[i-1]*10+10^(i- ...
- bzoj 1833: [ZJOI2010]count 数字计数【数位dp】
非典型数位dp 先预处理出f[i][j][k]表示从后往前第i位为j时k的个数,然后把答案转换为ans(r)-ans(l-1),用预处理出的f数组dp出f即可(可能也不是dp吧--) #include ...
- 1833: [ZJOI2010]count 数字计数
1833: [ZJOI2010]count 数字计数 Time Limit: 3 Sec Memory Limit: 64 MBSubmit: 2951 Solved: 1307[Submit][ ...
- 1833: [ZJOI2010]count 数字计数 - BZOJ
Description给定两个正整数a和b,求在[a,b]中的所有整数中,每个数码(digit)各出现了多少次.Input输入文件中仅包含一行两个整数a.b,含义如上所述.Output输出文件中包含一 ...
- 【BZOJ】1833: [ZJOI2010] count 数字计数(数位dp)
题目 传送门:QWQ 分析 蒟蒻不会数位dp,又是现学的 用$ dp[i][j][k] $ 表示表示长度为i开头j的所有数字中k的个数 然后预处理出这个数组,再计算答案 代码 #include < ...
随机推荐
- Debian 7 amd64问题
Debian 7 发布了有1段时间,最近才在自己的电脑硬盘安装,用户体验还算可以.在安装Debian的过程中,有问题还是要记录一下的. 注意:遇到的问题跟硬件体系相关,可能在个别电脑没法重现. 1.默 ...
- Perfmon - Windows 自带系统监控工具
一. 简述 可以用于监视CPU使用率.内存使用率.硬盘读写速度.网络速度等. Perfmon提供了图表化的系统性能实时监视器.性能日志和警报管理,系统的性能日志可定义为二进制文件.文本文件.SQLSE ...
- opencv里vector的内存的申请和释放的问题
改成: ); CvSeq * m_contour=; cvFindContours( &IPlImage(img), m_storage, &m_contour, ,)); //释放内 ...
- thinkPHP框架单一入口文件解析
一.index.php (可参考ThinkPHP学习手册http://document.thinkphp.cn/manual_3_2.html#entrance_file) index.php单入口 ...
- binlog2sql数据恢复
牛叉的工具有好几个,包括MyFlash.binlog2Sql.mysqlbinlog_flashback,还有一些收费的等等,各有优劣,具体使用可自行百度 1.安装binlog2sql shell&g ...
- ICommand接口
WPF 中的命令是通过实现 ICommand 接口创建的.ICommand 的 WPF 实现是 RoutedCommand 类. WPF 中的主要输入源是鼠标.键盘.墨迹和路由命令.更加面向设备的输入 ...
- set(gcf,'DoubleBuffer','on')以及sort
设置的目的是为了防止在不断循环画动画的时候会产生闪烁的现象,而这样便不会了.在动画的制作比较常用. Matlab排序函数-sort sort函数的调用格式: sort(X) 功能:返回对向量X中的元素 ...
- BZOJ 1407 Savage(拓展欧几里得)
这题的时间复杂度真玄学... O(m*n^2).1e8也能过啊... 首先题目保证m<=1e6. 这启发我们枚举或者二分答案? 但是答案不满足单调性,考虑从小到大枚举m. 对于每一个m,枚举两个 ...
- 【bzoj4641】基因改造 特殊匹配条件的KMP
题目描述 如果两个长度相等的字符串,如果存在一种字符的一一映射,使得第一个字符串的所有字符经过映射后与第二个字符串相同,那么就称它们“匹配”.现在给出两个串,求第一个字符串所有长度等于第二个字符串的长 ...
- POJ2689:Prime Distance——题解
http://poj.org/problem?id=2689 题目大意,给不超过int的l,r,其中r-l+1<=1000000,筛出其中的素数,并且求出相邻素数差值最大和最小的一对. ———— ...