[HAOI2010]计数 数位DP+组合数
题面:
你有一组非零数字(不一定唯一),你可以在其中插入任意个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+组合数的更多相关文章
- BZOJ2425:[HAOI2010]计数(数位DP)
Description 你有一组非零数字(不一定唯一),你可以在其中插入任意个0,这样就可以产生无限个数.比如说给定{1,2},那么可以生成数字12,21,102,120,201,210,1002,1 ...
- BZOJ 2425 [HAOI2010]计数:数位dp + 组合数
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=2425 题意: 给你一个数字n,长度不超过50. 你可以将这个数字: (1)去掉若干个0 ( ...
- 【BZOJ-1833】count数字计数 数位DP
1833: [ZJOI2010]count 数字计数 Time Limit: 3 Sec Memory Limit: 64 MBSubmit: 2494 Solved: 1101[Submit][ ...
- bzoj1833: [ZJOI2010]count 数字计数(数位DP+记忆化搜索)
1833: [ZJOI2010]count 数字计数 题目:传送门 题解: 今天是躲不开各种恶心DP了??? %爆靖大佬啊!!! 据说是数位DP裸题...emmm学吧学吧 感觉记忆化搜索特别强: 定义 ...
- 【题解】P2602 数字计数 - 数位dp
P2602 [ZJOI2010]数字计数 题目描述 给定两个正整数 \(a\) 和 \(b\) ,求在 \([a,b]\) 中的所有整数中,每个数码(digit)各出现了多少次. 输入格式 输入文件中 ...
- BZOJ 1833 数字计数 数位DP
题目链接 做的第一道数位DP题,听说是最基础的模板题,但还是花了好长时间才写出来..... 想深入了解下数位DP的请点这里 先设dp数组dp[i][j][k]表示数位是i,以j开头的数k出现的次数 有 ...
- bzoj1833: [ZJOI2010]count 数字计数 数位dp
bzoj1833 Description 给定两个正整数a和b,求在[a,b]中的所有整数中,每个数码(digit)各出现了多少次. Input 输入文件中仅包含一行两个整数a.b,含义如上所述. O ...
- NEUQ OJ 2004:追梦之人 (计数数位dp)
2004: 追梦之人 描述 题目描述: 为了纪念追梦人,粉丝们创造了一种新的数——“追梦数”.追梦数要满足以下两个条件:1.数字中不能出现“7”2.不能被7整除.比如:777和4396就不是追梦数,而 ...
- LightOJ 1140 计数/数位DP 入门
题意: 给出a,b求区间a,b内写下过多少个零 题解:计数问题一般都会牵扯到数位DP,DP我写的少,这道当作入门了,DFS写法有固定的模板可套用 dp[p][count] 代表在p位 且前面出现过co ...
随机推荐
- Python-内置函数4
import time # 返回时间戳 t=time.time() print(t) name="one" ''' bin() oct() hex() bytes() ascii( ...
- Docker - 容器中的tomcat如何使用startup.sh启动
网上大多介绍的catalina.sh启动,因为docker容器中,无法直接启动startup.sh. 解决方法: 编辑catalina.sh,找到 >> "$CATALINA_O ...
- 初试Docker on Debian on VirtualBox
一直以来都对Docker如雷贯耳,很想尝试一下但都被各种忙给耽误了,最近由于项目调试,需要安装 Oracle 和 SQL Server 数据库,但又不想安装到本机系统里,于是下决心啃一下docker这 ...
- 初识java atomic
2018-8-19 昨天看到java.util.concurrent.atomic相关的文章,之前有过留意但并未去了解,正好有空学习一下.本人理解atomic包是concurrent子包,当是为并发所 ...
- Linux内核设计笔记13——虚拟文件系统
虚拟文件系统 内核在它的底层文件系统系统接口上建立一个抽象层,该抽象层使Linux可以支持各种文件系统,即便他们在功能和行为上存在很大差异. VFS抽象层定义了各个文件系统都支持的基本的.概念上的接口 ...
- C# 生成行和列
private DataTable GetListBind() { DataTable dt = new DataTable(); try { dt.Columns.Add("1" ...
- wwnjld第二轮迭代测试报告
1.引言 1.1测试报告目的 被测试报告为wwnjld小组我们的时间管理软件的第二轮迭代所写的软件测试报告.在经过本小组大家不懈的努力之下,我们小组第二轮迭代的产品终于新鲜出炉了.这次测试小组的主要成 ...
- win7 个人电脑 IIS7服务器(web服务器) 同一局域网下均可访问本机网页
建立web服务器: 1.控制面板-->程序-->打开或关闭windows功能-->internet信息服务全部打钩,确定即可. 访问网页: 1.C:\inetpub\wwwroot\ ...
- 计算器软件实现系列(五)策略模式+asp.net
一 策略模式代码的编写 using System; using System.Collections.Generic; using System.Linq; using System.Web; /// ...
- TCP 的有限状态机
TCP 有限状态机的图中每一个方框都是 TCP 可能具有的状态. 每个方框中的大写英文字符串是 TCP 标准所使用的 TCP 连接状态名. 状态之间的箭头表示可能发生的状态变迁. 箭头旁边的字,表明引 ...