P1066 [NOIP2006 提高组] 2^k进制数

题目描述

设 �r 是个 2�2k 进制数，并满足以下条件：

• �r 至少是个 22 位的 2�2k 进制数。

• 作为 2�2k 进制数，除最后一位外，�r 的每一位严格小于它右边相邻的那一位。

• 将 �r 转换为二进制数 �q 后，则 �q 的总位数不超过 �w。

在这里，正整数 �,�k,w 是事先给定的。

问：满足上述条件的不同的 �r 共有多少个？

我们再从另一角度作些解释：设 �S 是长度为 �w 的 0101 字符串（即字符串 �S 由 �w 个 00 或 11 组成），�S 对应于上述条件三中的 �q。将 �S 从右起划分为若干个长度为 �k 的段，每段对应一位 2�2k 进制的数，如果 �S 至少可分成 22 段，则 �S 所对应的二进制数又可以转换为上述的 2�2k 进制数 �r。

例：设 �=3,�=7k=3,w=7。则 �r 是个八进制数（ 23=823=8 ）。由于 �=7w=7，长度为 77 的 0101 字符串按 33 位一段分，可分为 33 段（即 1,3,31,3,3，左边第一段只有一个二进制位），则满足条件的八进制数有：

22 位数：
高位为 11：66 个（即 12,13,14,15,16,1712,13,14,15,16,17 ）,
高位为 22：55 个，
…，
高位为 66：11 个（即 6767 ）。
共 6+5+…+1=216+5+…+1=21 个。

33 位数：
高位只能是 11，
第 22 位为 22：55 个（即 123,124,125,126,127123,124,125,126,127 ），
第 22 位为 33：44 个，
…，
第 22 位为 66：11 个（即 167167 ）。
共 5+4+…+1=155+4+…+1=15 个。

所以，满足要求的 �r 共有 3636 个。

输入格式

一行两个正整数 �,�k,w 用一个空格隔开：

输出格式

一行一个个正整数，为所求的计算结果。
即满足条件的不同的 �r 的个数（用十进制数表示），要求不得有前导零，各数字之间不得插入数字以外的其他字符（例如空格、换行符、逗号等）。

（提示：作为结果的正整数可能很大，但不会超过 200200 位）

输入输出样例

输入 #1复制

3 7

输出 #1复制

36

说明/提示

【数据范围】
1≤�≤91≤k≤9
1≤�≤3×1041≤w≤3×104

NOIP 2006 提高组 第四题

#等会？？我A了？？

这题正解不是组合数学吗

我看着机房的大佬一个个都在想正解真是瑟瑟发抖

然后我用了个递推结果还跑得贼快

思路大概是这样的

用a[i][j]表示组成一个i位数，第i位选j的所有情况

就等于a[i-1][j+1~maxx]的和

用ans表示每个a[i][j]的和，就是最后的结果

然后可以把第一维压掉，用前缀和进行优化，再加上高精度

代码里面有注释

#include<iostream>
#include<cstdio>
#define rint register int
#define ini inline int 
using namespace std;
int a[30001][201],c[201],ans[201],k,w;
ini jia(int *a,int *b)
{
    int lenc=0,x=0;
    while (lenc<a[0] || lenc<b[0])
    {
        a[++lenc]=a[lenc]+b[lenc]+x;
        x=a[lenc]/10;
        a[lenc]%=10;    
    }
    if (x>0) a[++lenc]=x;
    a[0]=lenc;
}
int main()
{
    scanf("%d%d",&k,&w);
    int kk=w%k;int hh=w/k;
    int y=0;
    for (rint i=1;i<=kk;i++)
        y+=1<<(i-1);     //最高位最大是几
    int minn=(1<<k)-1;//第一位最大是多少
    if (hh==1 || (hh==2 && kk==0))
    {
        if (kk==0) y=minn;int tot=0;
        for (rint i=1;i<=y;i++)
            tot+=minn-i;
        cout<<tot;
        return 0;
    }//特判两位以内的情况
    for (rint i=1;i<=minn-1;i++)
    {
        a[i][1]=i;a[i][0]=1;
        jia(ans,a[i]);
    }//第二位是minn-1到1各有多少种可能
    for (rint i=3;i<=hh;i++)
        for (rint j=1;j<=minn-i+1;j++)
//第i位最大是minn-i+1
//第i位是1到minn-i+1各有多少种可能
        {
            jia(a[j],a[j-1]);
            jia(ans,a[j]);
        }
//以下是单独处理最高位，因为最高位有限制
    for (rint j=1;j<=minn-hh;j++)
        jia(a[j],a[j-1]);
    for (rint j=minn-hh;j>=max(minn-hh-y+1,1);j--)
        jia(ans,a[j]);
//我们的前缀和是倒着存的
    for (rint i=ans[0];i>=1;i--)
        printf("%d",ans[i]);
}