zoukankan      html  css  js  c++  java
  • 蓝桥杯 问题 1110: 2^k进制数 (排列组合+高精度巧妙处理)

    题目链接

    题目描述
    设r是个2^k 进制数,并满足以下条件:
    (1)r至少是个2位的2^k 进制数。
    (2)作为2^k 进制数,除最后一位外,r的每一位严格小于它右边相邻的那一位。
    (3)将r转换为2进制数q后,则q的总位数不超过w。
    在这里,正整数k(1≤k≤9)和w(k〈w≤30000)是事先给定的。

    问:满足上述条件的不同的r共有多少个?
    我们再从另一角度作些解释:设S是长度为w 的01字符串(即字符串S由w个“0”或“1”组成),S对应于上述条件(3)中的q。将S从右起划分为若干个长度为k 的段,每段对应一位2^k进制的数,如果S至少可分成2段,则S所对应的二进制数又可以转换为上述的2^k 进制数r。
    例:设k=3,w=7。则r是个八进制数(2^3=8)。由于w=7,长度为7的01字符串按3位一段分,可分为3段(即1,3,3,左边第一段只有一个二进制位),则满足条件的八进制数有:
    2位数:高位为1:6个(即12,13,14,15,16,17),高位为2:5个,…,高位为6:1个(即67)。共6+5+…+1=21个。
    3位数:高位只能是1,第2位为2:5个(即123,124,125,126,127),第2位为3:4个,…,第2位为6:1个(即167)。共5+4+…+1=15个。
    所以,满足要求的r共有36个。

    输入
    只有1行,为两个正整数,用一个空格隔开: k w

    输出
    1行,是一个正整数,为所求的计算结果,即满足条件的不同的r的个数(用十进制数表示),要求最高位不得为0,各数字之间不得插入数字以外的其他字符(例如空格、换行符、逗号等)。
    (提示:作为结果的正整数可能很大,但不会超过200位)

    样例输入
    3 7

    样例输出
    36

    分析:
    这是一个组合数学问题,

    注意这句话:作为2^k 进制数,除最后一位外,r的每一位严格小于它右边相邻的那一位。

    其实这是在暗示组合数,

    显然r中的不会有相同的位,如果每一位都不同,显然只有严格递增的排列是合法的,这便是组合,

    将 r 转化成这种形式(设k为3) 000 000 000 000

    显然除首位外每一位的取值范围为 000 to 111(2^k-1)

    在首位为0的情况下,最多可取 w/k 位,且题目要求大于2位,

    则在首位为0的合法解有 ∑ C(2^k-1,i)(2<=i<=w/k) ,

    Ps. 如果 w 模 k 等于 0 仅考虑上述情况即可。

    考虑首位不为0的情况,显然首位不为0的话,r 就有 w/k+1 位,

    除首位外还有w/k位,可以枚举首位的取值范围为 1 to 2^(w mod k)-1

    设首位取值为 val,则剩下 w/k 位 取值范围为 val+1 to 2^k-1,也就是有 2^k-1-val 个数可取,所以首位不为0的合法解有 ∑ C(2^k-1-val,w/k)(1<=val<=2^(w mod k)-1)

    所以上述两者相加便是正解(需要高精运算)

    对于高精运算的处理,我看到了一个比较巧妙的方法避开了复杂的数组运算,就是把上面那些组合数的运算

    都转换成了 C(2^k-1,i) -----> C(2^k-i,i) 然后写C函数的时候不是计算C(2^k-1,i),而是计算C(2^k-i+i-1,i)。 巧妙之处就在这里,这样可以有效的避免有溢出吧。

    注意事项:
    1、本题的关键就是看懂题目,知道这是个排列组合问题。

    2、在排列组合的计算时,用了一个巧妙的方法,希望能看懂。

    代码:

    #include <iostream>
    #include <cmath>
    #include<stdio.h>
    using namespace std;
    long C(int n,int m)                 //公式为C(n+m-1)(m)[重点]   希望结合上面说的看懂
    {
        int i;
        long sum=1;
        for (i=1;i<=m;i++)
            sum*=(n+m-i);
        for (i=1;i<=m;i++)
            sum=sum/i;
        return sum;
    }
    int main()
    {
        int k, w;
        cin >> k >> w;
        int part_num = w / k + 1;             //一共分为几段
        int maxx_num_per_part = pow(2.0,k);             //取不到,每部分最大取到maxx_num_per_part减去1,就是每一位都是1的时候
        int gaow_num_max = pow(2.0,w%k) - 1;         //最高位的数最大值,可以是0
        //开始计算,分两种情况,第一种,首段为0,那么后面n位数对应的个数符合C[maxx_num_per_part-1][n]
        long long sum = 0;
        for (int i=2; i<=part_num-1; i++)           //去掉最高位部分,还有至少两位数
        {
            sum += C(maxx_num_per_part-i, i);
        }
        //第二种情况,首段不是0,如果首段为x,解就有C[maxx_num_per_part-1-x][w/k]
        for (int i=1; i<=gaow_num_max; i++)
        {
            sum += C(maxx_num_per_part-w/k-i,w/k);
        }
        cout << sum;
    return 0;
    }
    
  • 相关阅读:
    虚拟化与KVM部署
    Jenkins部署及使用
    Tomcat及LAMT架构搭建
    版本控制与Gitlab配置
    Cobbler介绍及搭载
    MYSQL集群配置
    HAProxy——配置与实例
    LVS——配置实例
    Linux——集群
    关于IntellIJ IDEA 2016.2.4新建项目的Java Enterprise没有显示问题
  • 原文地址:https://www.cnblogs.com/cmmdc/p/8635040.html
Copyright © 2011-2022 走看看