题面:

你有一组非零数字(不一定唯一),你可以在其中插入任意个0,这样就可以产生无限个数。比如说给定{1,2},那么可以生成数字12,21,102,120,201,210,1002,1020,等等。

现在给定一个数,问在这个数之前有多少个数。(注意这个数不会有前导0).

样例输入:1020     样例输出:7

题解:

刚看到这道题的时候有点懵,,,,,

其实仔细观察一下发现这题可以用组合做。

注意到0的个数是不限的,

而且如果位数小于给定n的话,肯定可以随便搭配,

所以直接加上给定数的有效位+补全的0的全排列就好了,(全排列注意不能重复)

有效位指的是不为0的位,

这样虽然有前导0,但是可以看做是位数少的数,比如

0012,0120分别可以看做12和120,

这样不仅不用去掉前导0,还可以只用加一次,非常方便

因此我们要求就只剩下位数相同的方案了,

这就要用到数位DP了

我们从前往后枚举,因为如果确定的当前位小于给定n的当前位的话,

由于是从前往后枚举,只要当前位小了,后面的位不管怎么排列都是小于n的,

因此我们加上去掉已经被确定的数后的全排列(全排列不能重复)就可以了,

如果当前位相同,那方案数就要放在后面求(注意每确定一位都要在有的数集合中减去它)

注意在第一位的时候不要枚举0就可以了

 #include<bits/stdc++.h>
using namespace std;
#define R register int
#define AC 60 #define LL long long
LL len,ans,tot;
int num[AC],have[AC],use[AC];
/*第i位是什么,在哪里?
设给定n的位数为len,那么位数小于len的显然可以随意组合,
于是利用有效位组合计算。
位数大于len的显然不合法
所以只要计算位数等于len的即可,
假设n的第i位不为0,那么如果用比第i位小的数摆在第i位上,
这个数就肯定比n小了,所以后面就可以随意搭配了
所以利用组合来算?
之前只想到位数小于len的可以组合算,
原来确定了前面的后面也可以组合算啊(废话),
如果g[i]表示i位数,自由组合的合法方案,i>=有效位
则g[i] = i! - (i-1)!
看上去怎么这么玄啊。。。。
果然是错的,,,因为如果有效位中有重复的话,这一部分就会被重复计算,
比如26637,那么对于每一种合法排列,都会计算2次,因为6可以交换
所以要去重,如果数i的出现次数为have[i],那么它会使得整个数列被多计算have[i]!次,
同理,拓展到所有数,整个数列会被多计算have[i]! * have[i+1]! + .......
因此设f[i]表示长度为i时的排列(无重复,有前导零)
g[i]表示长度为i时的排列(无重复,无前导零)
那么g[i] = f[i] - f[i-1];//相当于确定了第1位是0时,就必然有前导零了.
设比当前位i大的数有k位,
那么ans += k * newf[n-i];
因为这时前几位已经确定,所以have和tot都变了,因此要用newf,而不是f
因此可以发现,前面求的g和f到后面根本就没有用,所以还是不要求了,,,,,
直接算就好了
ans += tot!/have[1]! * have[2]!.....
*/
void pre()
{
char c=getchar();
while(c > '' || c < '') c=getchar();
while(c >= '' && c <= '')
{
num[++len]=c - '',c=getchar();
if(num[len] != ) ++have[num[len]],++tot;//统计有效位
}//读入
have[]=len - tot - ;//获取0的个数
if(!tot)
{
printf("0\n");
exit();
}
} LL cal(int cnt)
{
LL tmp=;
/* for(R i=1;i<=9;i++) tmp+=have[i],use[i]=have[i];
use[0]=cnt - tmp;//获取0的个数,因为work中获取了,所以不用再次获取*/
for(R i=;i<=;i++) use[i] = have[i];
tmp=;
for(R i=;i<=cnt;i++)
{
tmp*=i;
for(R j=;j<=;j++)
{
if(use[j] <= ) continue;
while(!(tmp % use[j]))
{
tmp /= use[j];
--use[j];
if(use[j] == ) break;
}
}
}
return tmp;
} void work()
{
LL tmp;
have[]++;//把之前为了计算小于len位的答案而减掉的0加回来
for(R i=;i<=len;i++)
{
tmp=;
for(R j=;j<=;j++) tmp+=have[j];
have[]=(len - i + ) - tmp;//获取0的个数
int b=(i == );//因为第一位不能用0代替
for(R j=b;j<num[i];j++)//枚举放哪个,但只能尝试比当前位小的,且这样还可以保证不重复
{
if(!have[j]) continue;//没有就不能放
--have[j];//先减掉
/*if(i == 1)//error!!!但是只要第一个不是0的话,,,,后面的就不是前导零啊,,,有什么关系。。。。
{//求g值,如果第一个是0的话,其他的have不会变,have[0]又是临时求的,所以没关系
ans += cal(len - i);
--have[0];
ans -= cal(len - i -1);
++have[0];
}*/
ans+=cal(len - i);//只有第一位才不能有前导0,,,,
++have[j];
}
if(num[i]) have[num[i]]--;//假设这一位也相同,那么就要减掉,0的时候不用减,因为0是临时算的(不过应该也可以不临时算?)
}
printf("%lld\n",ans);
} int main()
{
// freopen("in.in","r",stdin);
pre();
if(len - >= tot) ans+=cal(len-);//这样会有前导零,但可以看做是位数不同
work();
// fclose(stdin);
return ;
}

[HAOI2010]计数 数位DP+组合数的更多相关文章

  1. BZOJ2425:[HAOI2010]计数(数位DP)

    Description 你有一组非零数字(不一定唯一),你可以在其中插入任意个0,这样就可以产生无限个数.比如说给定{1,2},那么可以生成数字12,21,102,120,201,210,1002,1 ...

  2. BZOJ 2425 [HAOI2010]计数:数位dp + 组合数

    题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=2425 题意: 给你一个数字n,长度不超过50. 你可以将这个数字: (1)去掉若干个0 ( ...

  3. 【BZOJ-1833】count数字计数 数位DP

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

  4. bzoj1833: [ZJOI2010]count 数字计数(数位DP+记忆化搜索)

    1833: [ZJOI2010]count 数字计数 题目:传送门 题解: 今天是躲不开各种恶心DP了??? %爆靖大佬啊!!! 据说是数位DP裸题...emmm学吧学吧 感觉记忆化搜索特别强: 定义 ...

  5. 【题解】P2602 数字计数 - 数位dp

    P2602 [ZJOI2010]数字计数 题目描述 给定两个正整数 \(a\) 和 \(b\) ,求在 \([a,b]\) 中的所有整数中,每个数码(digit)各出现了多少次. 输入格式 输入文件中 ...

  6. BZOJ 1833 数字计数 数位DP

    题目链接 做的第一道数位DP题,听说是最基础的模板题,但还是花了好长时间才写出来..... 想深入了解下数位DP的请点这里 先设dp数组dp[i][j][k]表示数位是i,以j开头的数k出现的次数 有 ...

  7. bzoj1833: [ZJOI2010]count 数字计数 数位dp

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

  8. NEUQ OJ 2004:追梦之人 (计数数位dp)

    2004: 追梦之人 描述 题目描述: 为了纪念追梦人,粉丝们创造了一种新的数——“追梦数”.追梦数要满足以下两个条件:1.数字中不能出现“7”2.不能被7整除.比如:777和4396就不是追梦数,而 ...

  9. LightOJ 1140 计数/数位DP 入门

    题意: 给出a,b求区间a,b内写下过多少个零 题解:计数问题一般都会牵扯到数位DP,DP我写的少,这道当作入门了,DFS写法有固定的模板可套用 dp[p][count] 代表在p位 且前面出现过co ...

随机推荐

  1. 后续博客转移到zhylj.cc

    此博客暂不更新了 zhylj.cc

  2. js 去掉下划线,后首个字母变大写

    1.驼峰转连字符: var s = "fooStyleCss";  s = s.replace(/([A-Z])/g,"-$1").toLowerCase(); ...

  3. PS 证件照换颜色

    1.打开要修改的图片,然后先Ctrl+J备份一份 2.点击魔法棒,点击要换颜色的地方,如衣服,之后会出现虚线,如果自动选择的不全,可以按住Shift键自行选择区域 3.然后Shift+Fn+F5(由于 ...

  4. 在Linux CentOS7系统中搭建LNMP

    LNMP就是Linux+Nginx+MySQL+PHP,既然是在Linux CentOS7那么Linux就是已经安装好了.所以接下百度一下接下来的教程,整理测试如下: 教程是centos6.2的有点老 ...

  5. Anyproxy抓包工具

    1.安装Nodejs和AnyProxy以及安装模拟器(移动端抓包工具) 第一步:nodejs下载地址:http://nodejs.cn/download/ 下载Windows版本,直接运行安装即可, ...

  6. 骰子涂色 (Cube painting,UVa 253)

    题目描述:算法竞赛入门习题4-4  题目思路:1.旋转其中一个骰子进行匹配 2.进行遍历,如果匹配,就进行相对面的匹配 3.三个对立面都匹配即是一样等价的 //没有按照原题的输入输出 #include ...

  7. 理解glance

    摘要: 本节介绍 OpenStack Image 服务 Glance 的基本概念. OpenStack 由 Glance 提供 Image 服务. 理解 Image 要理解 Image Service ...

  8. HDU 4302 Holedox Eating (线段树模拟)

    题意:一个老鼠在一条长度为L的直线上跑,吃蛋糕,老鼠只能沿直线移动.开始时没有蛋糕,老鼠的初始位置是0. 有两个操作,0 x 代表在位置x添加一个蛋糕: 1 代表老鼠想吃蛋糕.老鼠每次都会选择离自己最 ...

  9. HBase 参考文档翻译之 Getting Started

    本篇是对HBase官方参考文档的大体翻译,介于本人英文水平实在有限,难免有纰漏之处.本篇不只是对官方文档的翻译,还加入了一些本人对HBase的理解.在翻译过程中,一些没有营养的废话,我就忽略了没有翻译 ...

  10. C++计算器项目的初始部分(第三次作业)

    C++计算器项目的初始部分 项目源文件地址:calculator 项目信息: * 项目名称:Calculator * 项目实现: * 对四则运算表达式进行拆解 * 对拆解的表达式进行简单的错误判断 * ...