• hdu 4165 dp


    可以用卡特兰数做

    以下分析转自:http://www.cnblogs.com/kevinACMer/p/3724640.html?utm_source=tuicool

    这道题之前自己做的时候并没有反应过来是求卡特兰数,当时是按递推来想的。后来查了下HDU4165,结果一看大标题就说是卡特兰数,自己想了想, 还真是那么回事。主要还是对于卡特兰数用的不多,也就当时没立马反应过来了。下面介绍这道题我的思路,然后对卡特兰数再做一些补充。

    本题题意:罐子里有N片相同的药片,开始的时候药片都是完整的一整片。然后一个每天从中任意取一片,如果该药片是完整的一整片,那么他会吃点半片, 然后将剩余半片扔回罐里,如果恰好是半片,那他直接吃掉拿出来的半片。显然2N天后,他吃完全部N片药片,总共有多少种不同的吃法过程。

    我当时的思路顺序是这样的:发现对于一个状态,它的下一个状态会有两种可能(当然有些特殊情况是只有一种可能),那么对于这两种可能显然容易想到只 要递推下去就行了。然后再考虑,这两种可能无非就是下个状态我可能吃完整的药片,也可能吃只有半片的药片。所以递推显然设计两个参数,完整药片的个数和只 有半片药片的个数。于是开了个二维数组dp[i][j],表示当罐子里有i片完整的药片和j片只有一半的药片时距离吃完总共有多少种不同的吃法。

    然后递推过程应该比较容易就出来了:dp[i][j]=dp[i-1][j+1]+dp[i][j-1],递推式说明:对于当前dp[i][j]这 个状态来说,它接下来会出现两个分叉,一种可能,接下来拿出的是一枚完整的药片,然后吃一半,扔回一半,于是变成了dp[i-1][j+1]这个状态,还 一种可能就是正好拿出来的是只有半片的药片,于是变成了dp[i][j-1]这个状态。然后由于DP定义的是方案数,显然这两个之后的状态加起来的方案数 就是之前那个状态的方案数了。

    当然,接下来就是给这个二维dp赋初始值了,显然dp[1][0]=1,然后发现其实dp[1][i]=i+1(可以这样理解:当罐子里只有i片半 状的药片和1片完整药片时,对于求方案数,其实主要就是看那个完整的药片在第几天吃,总共i+1天,那么显然就是C(i+1)(1)种方案数)。然后 dp[i][j]中的i相当于就从2开始递推了。对于j,当j为0时,它对应的下一个状态其实只有一种可能,就是dp[i][0]=dp[i-1] [1],于是j就从1开始递推了。具体代码如下:

    2015-05-11

     1 #include<cstdio>
     2 #include<iostream>
     3 #include<algorithm>
     4 #include<cstring>
     5 #include<cmath>
     6 #include<queue>
     7 #include<map>
     8 using namespace std;
     9 #define MOD 1000000007
    10 const int INF=0x3f3f3f3f;
    11 const double eps=1e-5;
    12 typedef long long ll;
    13 #define cl(a) memset(a,0,sizeof(a))
    14 #define ts printf("*****
    ");
    15 const int MAXN=35;
    16 ll dp[MAXN][MAXN];
    17 int n,m,tt;
    18 void init()
    19 {
    20     for(int i=0;i<MAXN;i++)    dp[1][i]=i+1;
    21     for(int i=2;i<MAXN;i++)
    22     {
    23         dp[i][0]=dp[i-1][1];
    24         for(int j=1;j<MAXN;j++)
    25         {
    26             dp[i][j]=dp[i-1][j+1]+dp[i][j-1];
    27         }
    28     }
    29 }
    30 int main()
    31 {
    32     int i,j,k;
    33     #ifndef ONLINE_JUDGE
    34     freopen("1.in","r",stdin);
    35     #endif
    36     int n;
    37     init();
    38     while(scanf("%d",&n)!=EOF)
    39     {
    40         if(n==0)break;
    41         printf("%I64d
    ",dp[n][0]);
    42     }
    43     return 0;
    44 }

    本题思路就是这样,比较简单的递推。但之后看到一些题解,对于本题,其实就是求卡特兰数。
    我们这样来思考:本题的题意转化一下可以这样理解,有N个红球,N个白球,问这2N个小球的排列数,且从第一个球到后面任意个球中红球的个数必须不小于白球的个数(理解为白球是红球产生的,一个红球产生一个白球)。那么这个就变成了卡特兰数的一个典型应用了。

    卡特兰数:

    令h(0)=1,h(1)=1,catalan数满足递归式:

    h(n)= h(0)*h(n-1) + h(1)*h(n-2) + + h(n-1)h(0) (其中n>=2)   

    该递推关系的解为:   h(n)=C(2n,n)/(n +1) (n=1,2,3,)      

    另类递归式:  h(n)=((4*n-2)/(n+1))*h(n-1);

    组合数学中有非常多.的组合结构可以用卡塔兰数来计数。以下用Cn=3和Cn=4举若干例:

    • Cn表示长度2n的dyck word的个数。Dyck word是一个有n个X和n个Y组成的字串,且所有的部分字串皆满足X的个数大于等于Y的个数。以下为长度为6的dyck words:
    XXXYYY XYXXYY XYXYXY XXYYXY XXYXYY
    • 将上例的X换成左括号,Y换成右括号,Cn表示所有包含n组括号的合法运算式的个数:
    ((())) ()(()) ()()() (())() (()())
    • Cn表示有n+1个叶子的二叉树的个数。

                                                                          

    • Cn表示所有不同构的含n个分枝结点的满二叉树的个数。(一个有根二叉树是满的当且仅当每个结点都有两个子树或没有子树。)
    • 证明:

      令1表示进栈,0表示出栈,则可转化为求一个2n位、含n个1、n个0的二进制数,满足从左往右扫描到任意一位时,经过的0数不多于1数。显然含n个1、n个0的2n位二进制数共有{2n choose n}个,下面考虑不满足要求的数目.

      考虑一个含n个1、n个0的2n位二进制数,扫描到第2m+1位上时有m+1个0和m个1(容易证明一定存在这样的情况),则后面的0-1排列中必有n-m个1和n-m-1个0。将2m+2及其以后的部分0变成1、1变成0,则对应一个n+1个0和n-1个1的二进制数。反之亦然(相似的思路证明两者一一对应)。

      从而C_n = {2n choose n} - {2n choose n + 1} = frac{1}{n+1}{2n choose n}。证毕。

      • Cn表示所有在n×n格点中不越过对角线的单调路径的个数。一个单调路径从格点左下角出发,在格点右上角结束,每一步均为向上或向右。计算这种路径的个数等价于计算Dyck word的个数: X代表“向右”,Y代表“向上”。下图为n= 4的情况:
      •                                                                        
      • Cn表示通过连结顶点而将n + 2边的凸多边形分成三角形的方法个数。下图中为n= 4的情况:

                                                                                      

      • Cn表示对{1, ...,n}依序进出置换个数。一个置换w是依序进出栈的当S(w) = (1, ...,n), 其中S(w)递归定义如下:令w=unv,其中nw的最大元素,uv为更短的数列;再令S(w) =S(u)S(v)n,其中S为所有含一个元素的数列的单位元。
      • Cn表示用n个长方形填充一个高度为n的阶梯状图形的方法个数。下图为n = 4的情况:

                                                                                               

  • 相关阅读:
    ArcEngine连接表join
    坐标系统与投影变换及在ARCGIS中的应用
    ArcEngine 最短路径分析(源码)【转载】
    TOC图例符号转换到ListView
    Dijkstra 最短路径算法的一种高效率实现
    添加、删除连接(Join)和关联(Relate)VBA实例
    C# 实现在线软件自动升级程序
    MapX八个技巧
    希腊字母的读法
    经纬度点与投影坐标系下坐标点的互相转换
  • 原文地址:https://www.cnblogs.com/cnblogs321114287/p/4310279.html
Copyright © 2020-2023  润新知