【算法】数位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 数字计数的更多相关文章

  1. BZOJ 1833: [ZJOI2010]count 数字计数( dp )

    dp(i, j, k)表示共i位, 最高位是j, 数字k出现次数. 预处理出来. 差分答案, 对于0~x的答案, 从低位到高位进行讨论 -------------------------------- ...

  2. [BZOJ 1833] [ZJOI2010] count 数字计数 【数位DP】

    题目链接:BZOJ - 1833 题目分析 数位DP .. 用 f[i][j][k] 表示第 i 位是 j 的 i 位数共有多少个数码 k . 然后差分询问...Get()中注意一下,如果固定了第 i ...

  3. BZOJ 1833: [ZJOI2010]count 数字计数

    Description 问 \([L,R]\) 中0-9的个数. Sol 数位DP. 预处理好长度为 \(i\), 最高位为 \(j\) 的数位之和. 然后从上往下计算,不要忘记往下走的同时要把高位的 ...

  4. bzoj 1833 [ZJOI2010]count 数字计数(数位DP)

    [题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=1833 [题意] 统计[a,b]区间内各数位出现的次数. [思路] 设f[i][j][k ...

  5. BZOJ 1833 ZJOI2010 count 数字计数 数位DP

    题目大意:求[a,b]间全部的整数中0~9每一个数字出现了几次 令f[i]为i位数(算前导零)中每一个数出现的次数(一定是同样的,所以仅仅记录一个即可了) 有f[i]=f[i-1]*10+10^(i- ...

  6. bzoj 1833: [ZJOI2010]count 数字计数【数位dp】

    非典型数位dp 先预处理出f[i][j][k]表示从后往前第i位为j时k的个数,然后把答案转换为ans(r)-ans(l-1),用预处理出的f数组dp出f即可(可能也不是dp吧--) #include ...

  7. 1833: [ZJOI2010]count 数字计数

    1833: [ZJOI2010]count 数字计数 Time Limit: 3 Sec  Memory Limit: 64 MBSubmit: 2951  Solved: 1307[Submit][ ...

  8. 1833: [ZJOI2010]count 数字计数 - BZOJ

    Description给定两个正整数a和b,求在[a,b]中的所有整数中,每个数码(digit)各出现了多少次.Input输入文件中仅包含一行两个整数a.b,含义如上所述.Output输出文件中包含一 ...

  9. 【BZOJ】1833: [ZJOI2010] count 数字计数(数位dp)

    题目 传送门:QWQ 分析 蒟蒻不会数位dp,又是现学的 用$ dp[i][j][k] $ 表示表示长度为i开头j的所有数字中k的个数 然后预处理出这个数组,再计算答案 代码 #include < ...

随机推荐

  1. java — 垃圾回收

    1. 垃圾回收的意义 在java中,当没有对象指向原先分配给某个对象的内存的时候,这片内存就变成了垃圾,JVM的一个系统级线程就会自动释放这个内存块,垃圾回收意味着程序不再需要的对象是“无用的信息”, ...

  2. 3dContactPointAnnotationTool开发日志(七)

      调了半天发现是逻辑错误,改了一下终于没那么奇怪了:   但是有的接触点很明显跑偏了.再回顾一下自己是怎么求的,我是直接用的下面的代码求解一个点是否在另一个物体内部: var bounds = us ...

  3. 3ds Max学习日记(二)

      搞了一天终于把第二章弄完了,虽然有点慢,但也是没有办法,毕竟中途学累了总得放松一下吧.   又get到了一个新知识点,是关于3ds max的捕捉工具的,详情见3DS MAX 学习笔记 一(捕捉工具 ...

  4. seaj和requirejs模块化的简单案例

    如今,webpack.gulp等构件工具流行,有人说seajs.requirejs等纯前端的模块化工具已经被淘汰了,我不这么认为,毕竟纯前端领域想要实现模块化就官方来讲,还是有一段路要走的.也因此纯前 ...

  5. C#的垃圾回收

    C#中垃圾回收 GC.Collect();强制进行内存回收.

  6. 统一日志系统 Log4Net/ExceptionLess

    一.   写在前面 本文Log4Net介绍了基础的方式,大数据量生产环境不能使用,中等日志量请日志单库. 希望爱技术的你不要错过exceptionless和ELK 第四节开始简单配置大牛们推荐的了Ex ...

  7. IPython 4.0发布:Jupyter和IPython分离后的首个版本

    IPython 4.0日前正式发布,这是IPython分离成IPython和Jupyter后的第一个重要版本. 更新Jupyter的快捷方式是: pip install --upgrade jupyt ...

  8. Java中的多线程=你只要看这一篇就够了

    如果对什么是线程.什么是进程仍存有疑惑,请先Google之,因为这两个概念不在本文的范围之内. 用多线程只有一个目的,那就是更好的利用cpu的资源,因为所有的多线程代码都可以用单线程来实现.说这个话其 ...

  9. 为windows phone listbox 添加触摸倾斜效果

    在开发windows phone程序时,经常会用到listbox或者是longlistselector等列表控件.当点击时没有触摸效果体验会稍差一些,像windows phone中的设置页面一样,点击 ...

  10. hdu 1142 A Walk Through the Forest (最短路径)

    A Walk Through the Forest Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Jav ...