发现自己推得组合数好像不太一样

先把这个复杂的柿子写一遍

\[\sum_{i=2}^{\left \lfloor\frac{n}{k}\right \rfloor}C_{2^k-1}^{i}+\sum_{i=1}^{2^{n\text{ } \text{mod} \text{ }k}-1}C_{2^k-1-i}^{\left \lfloor\frac{n}{k}\right \rfloor}
\]

感觉这个柿子非常蛇皮

但是非常好求啊

由于\(2^k-1\)非常小,最大仅仅是\(511\),所以我们没有什么必要预处理阶乘,我们可以直接用组合数递推的方式来做

于是不需要打高精除或者高精乘了,一个高精加就够了

于是做法就非常无脑了,重要的是这个柿子是怎么推出来的

首先我们先考虑一个非常弱化的版本,就是\(k|n\)

如果\(k|n\)的话,**那么这个长度为\(n\)的二进制数就能被恰好分成\(n/k\)个块,而且每一个块能选择的数都是\(0\)到\(2^k-1\)这\(2^k\)个数

**

我们发现\(0\)这个非常不好考虑,于是我们可以先忽略掉\(0\)

所以现在有\(n/k\)个块,每个块内能填\(2^k-1\)种数

那么就有\(C_{2^k-1}^{n/k}\)种可能

之后我们再来考虑\(0\)的情况,首先最高位(如果不是第二位的话)是可以填\(0\)的,而剩下的\(n/k-1\)个块我们仍旧按照之前的方式来填,于是就有\(C_{2^k-1}^{n/k-1}\),之后对于次高位还是可以填\(0\)(同时最高位也填\(0\)),那么还有\(n/k-2\)个块,于是就是\(C_{2^k-1}^{n/k-2}\)

以此类推,直到对于第三个的块,我们还是可以填将这个块以及之前所有的块都填\(0\),那么就还有\(2\)个块,于是就是\(C_{2^k-1}^{2}\)

而第二个块可是不能填\(0\)了,于是就没了

所以对于\(k|n\)的时候,答案就是

\[\sum_{i=2}^{n/k}C_{2^k-1}^{i}
\]

之后我们再来考虑一下\(k\)不整除\(n\)的情况

这个样子的话一共会分成\(\left \lfloor\frac{n}{k}\right \rfloor+1\)个块,\(\left \lfloor\frac{n}{k}\right \rfloor\)个块内可以选择的数都是\(0\)到\(2^k-1\)这\(2^k\)个数,而最后一个不完整的块只有\(n\text{ } \text{mod} \text{ }k\)位,所以能选择的数只有\(0\)到\(2^{n\text{ } \text{mod} \text{ }k}-1\)

如果这个最高位选择填\(0\)那么退化成了\(k|n\)的情况,所以最高位填0的方案数为

\[\sum_{i=2}^{\left \lfloor\frac{n}{k}\right \rfloor}C_{2^k-1}^{i}
\]

之后最高位还有\(1\)到\(2^{n\text{ } \text{mod} \text{ }k}-1\)这些数可以填,如果我们选择填\(i\)的话,那么剩下的块内就不能填比\(i\)小的数,于是剩下的每个块内能选择的就有\(2^k-1-i\)个数,所以方案数就是\(C_{2^k-1-i}^{\left \lfloor\frac{n}{k}\right \rfloor}\)

所以最后的答案还应该加上

\[\sum_{i=1}^{2^{n\text{ } \text{mod} \text{ }k}-1}C_{2^k-1-i}^{\left \lfloor\frac{n}{k}\right \rfloor}
\]

代码

#include<cstring>
#include<string>
#include<cstdio>
#include<iostream>
#define re register
#define maxn 512
using namespace std;
string c[maxn][maxn];
int n,k;
int p,t;
int res;
int aa[201],bb[201],cc[201];
inline string sum(string a,string b)
{
memset(aa,0,sizeof(aa));
memset(bb,0,sizeof(bb));
memset(cc,0,sizeof(cc));
int lena=a.size();
int lenb=b.size();
for(re int i=0;i<lena;i++)
aa[i+1]=a[lena-i-1]-48;
for(re int i=0;i<lenb;i++)
bb[i+1]=b[lenb-i-1]-48;
int p=1;
for(p=1;p<=max(lena,lenb)||cc[p];p++)
{
cc[p]+=aa[p]+bb[p];
cc[p+1]+=cc[p]/10;
cc[p]%=10;
}
string C="\0";
for(re int i=p-1;i;i--)
C+=char(cc[i]+48);
return C;
}
int main()
{
scanf("%d%d",&k,&n);
p=n/k;
res=n%k;
t=(1<<k);
c[0][0]="1";
for(re int i=1;i<=t-1;i++)
c[i][0]=c[i][i]="1";
for(re int i=1;i<t;i++)
for(re int j=1;j<i;j++)
c[i][j]=sum(c[i-1][j-1],c[i-1][j]);
string ans="0";
for(re int i=2;i<=p;i++)
{
if(i>t-1) break;
ans=sum(ans,c[t-1][i]);
}
int pp=(1<<res)-1;
for(re int i=1;i<=pp;i++)
{
if(p>t-1-i) break;
ans=sum(ans,c[t-1-i][p]);
}
cout<<ans<<endl;
return 0;
}

【2^k进制数】的更多相关文章

  1. [转]as3 算法实例【输出1 到最大的N 位数 题目:输入数字n,按顺序输出从1 最大的n 位10 进制数。比如输入3,则输出1、2、3 一直到最大的3 位数即999。】

    思路:如果我们在数字前面补0的话,就会发现n位所有10进制数其实就是n个从0到9的全排列.也就是说,我们把数字的每一位都从0到9排列一遍,就得到了所有的10进制数. /** *ch 存放数字 *n n ...

  2. 1813. M进制数问题

    1813. M进制数问题 Constraints Time Limit: 1 secs, Memory Limit: 32 MB Description 试用 C++的类来表示一般进制数. 给定 2 ...

  3. CF459C Pashmak and Buses (构造d位k进制数

    C - Pashmak and Buses Codeforces Round #261 (Div. 2) C. Pashmak and Buses time limit per test 1 seco ...

  4. [codevs1157]2^k进制数

    [codevs1157]2k进制数 试题描述 设r是个2k 进制数,并满足以下条件: (1)r至少是个2位的2k 进制数. (2)作为2k 进制数,除最后一位外,r的每一位严格小于它右边相邻的那一位. ...

  5. noip2006 2^k进制数

    设r是个2k进制数,并满足以下条件: (1)r至少是个2位的2k进制数. (2)作为2k进制数,除最后一位外,r的每一位严格小于它右边相邻的那一位. (3)将r转换为2进制数q后,则q的总位数不超过w ...

  6. c++描述将一个2进制数转化成10进制数(用到初始化栈,进栈,入栈)

    /* c++描述将2进制数转化成10进制数 问题,1.初始化栈后,用new,不知道delete是否要再写一个函数释放内存, 还是在哪里可以加上delete 2.如果栈满了,我要分配多点空间,我想的办法 ...

  7. 关于不同进制数之间转换的数学推导【Written By KillerLegend】

    关于不同进制数之间转换的数学推导 涉及范围:正整数范围内二进制(Binary),八进制(Octonary),十进制(Decimal),十六进制(hexadecimal)之间的转换 数的进制有多种,比如 ...

  8. NOIP2006 2k进制数

    2^k进制数 题目描述 设r是个2^k 进制数,并满足以下条件: (1)r至少是个2位的2^k 进制数. (2)作为2^k 进制数,除最后一位外,r的每一位严格小于它右边相邻的那一位. (3)将r转换 ...

  9. 4位bcd数转换为2进制数

    DATA   SEGMENTBUF    DW  1234HBUF1   DW  ?    ;2进制数放到buf1内存单元DATA   ENDSCODE   SEGMENTASSUME CS:CODE ...

  10. P1066 2^k进制数

    传送门 题目描述 设r是个2^k 进制数,并满足以下条件: (1)r至少是个2位的2^k 进制数. (2)作为2^k 进制数,除最后一位外,r的每一位严格小于它右边相邻的那一位. (3)将r转换为2进 ...

随机推荐

  1. 如何高效的算出2x8的值

    原文出自:https://blog.csdn.net/seesun2012 位移算法,如何高效的算出2*8的值,为什么8<<1,4<<2,2<<3,1<< ...

  2. T_SQL 语句想已有数据表添加约束

    如果向存在数据的表里添加约束,有可能会出现数据不符合检查约束而造成添加约束失败. 如: USE DEmo--指向当前操作的数据库 GO ALTER TABLE Employee ADD CONSTRA ...

  3. Spring - 几种RPC模型的使用与比较

    Spring中,用JMS搞RPC时会用到: org.springframework.jms.remoting.JmsInvokerServiceExporter org.springframework ...

  4. IE11 F12 开发人员工具 查看 Cookie

    参考网址:Using the F12 developer tools in IE11 Step1 : IE11 => F12 打开 开发人员工具 Step2:开发人员工具 => 网络F5 ...

  5. MongoDB 学习(一)安装配置和简单应用

    一.安装和部署 1.服务端安装 1.官网下载(官方网站 https://www.mongodb.org/downloads/#production),傻瓜式安装,注意修改安装路径. 安装完成后的目录结 ...

  6. 流畅的python和cookbook学习笔记(二)

    1.元组拆包和解压序列赋值 任何的序列 (或者是可迭代对象) 可以通过一个简单的赋值语句解压并赋值给多个 变量.唯一的前提就是变量的数量必须跟序列元素的数量是一样的. 1.平行赋值: >> ...

  7. 3、springboot之热部署

    我用的是idea 一.开启idea自动make功能 1.CTRL + SHIFT + A --> 查找make project automatically --> 选中 2.CTRL + ...

  8. Spring_Spring与IoC_基于注解的DI

    一.基本注解的使用 (1)导入AOP的Jar包 (2) 与set()无关 二.组件扫描器的base-package 三.@Component相关注解 四.@Scope 五.域属性的注入 (1)byTy ...

  9. HDU 4283 (第k个出场 区间DP)

    http://blog.csdn.net/acm_cxlove/article/details/7964594 http://www.tuicool.com/articles/jyaQ7n http: ...

  10. 给tomcat7w.exe改名字

    在平常使用tomcat的过程中,经常会在自己的电脑上安装两个tomcat,而且使用tomcat的tomcat7w.exe界面比较方便,即 但是两个tomcat的名字一样,也只能启动一个tomcat. ...