Train Problem II(卡特兰数以及高精度运算)

版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/sinat_37668729/article/details/72808832

 
Problem Description
As we all know the Train Problem I, the boss of the Ignatius Train Station want to know if all the trains come in strict-increasing order, how many orders that all the trains can get out of the railway.
 
Input
The input contains several test cases. Each test cases consists of a number N(1<=N<=100). The input is terminated by the end of file.
 
Output

            For each test case, you should output how many ways that all the trains can get out of the railway.
 
Sample Input
1
2
3
10
 
Sample Output
1
2
5
16796

<div style="" font-size:="" 14px;="" border-top:="" #b7cbff="" 1px="" dashed;="" font-family:="" times"="">
Hint
The result will be very large, so you may not process it by 32-bit integers.
 
Author
Ignatius.L

前面补充中以说过卡特兰数,不再赘述!再次,有两个重要公式:

Catalan数的组合公式为 Cn=C(2n,n) / (n+1);

递推公式为 h(n ) = h(n-1)*(4*n-2) / (n+1)

代码:

#include<iostream>  
using namespace std;  
int a[101][101]={0};  
int main()  
{  
    int n,i,j,len,r,temp,t;  
    int b[101];  
    a[1][0] = 1;  
    len = 1;  
    b[1] = 1;  
    for(i=2;i<=100;i++)  
    {  
        t = i-1;  
        for(j=0;j<len;j++) //乘法  
            a[i][j] = a[i-1][j]*(4*t+2);  
        for(r=j=0;j<len;j++)  //处理相乘结果  
        {  
            temp = a[i][j] + r;  
            a[i][j] = temp % 10;  
            r = temp / 10;  
        }  
        while(r) //进位处理  
        {  
            a[i][len++] = r % 10;  
            r /= 10;  
        }  
  
        for(j=len-1,r=0;j>=0;j--) //除法  
        {  
            temp = r*10 + a[i][j];  
            a[i][j] = temp/(t+2);  
            r = temp%(t+2);  
        }  
        while(!a[i][len-1]) //高位零处理  
            len --;  
        b[i] = len;  
    }  
    while(cin>>n)  
    {     
        for(j=b[n]-1;j>=0;j--)  
            printf("%d",a[n][j]);  
        printf("/n");  
    }  
    return 0;  
}  



另外运用递归:

//输出卡特兰数   
//首先需要肯定,程序是正确的   
//这算是大数乘除法!记住他们是如何处理的!由于数据很大,用基本数据类型根本无法满足要求,只能用数组来表示!  
#include <iostream>  
#include<cstdio>  
#include<memory.h>  
using namespace std;  
#define MAX 101  
#define BASE 10000//base只是一个基度,对最终取值并没有影响,相反,base取值愈大,计算量愈小  
//base发生改变的时候,下面的输出也要相应地做出调整,否则也会输出错误答案!除非当base取10!  
void multiply(int a[],int len,int b)//乘法   
{  
    for(int i=len-1,carry=0;i>=0;--i)//从最后一位开始相乘,依次向前与每一位相乘   
    {//问题在于,为什么BASE=10000?   
        carry+=b*a[i];  
        a[i]=carry%BASE;  
        carry/=BASE;  
    }   
}  
void divide(int a[],int len,int b)//除法,很妙的!这种除法可能想不到,仔细体会!   
{//应当如何除呢?  
    for(int i=0,div=0;i<len;++i)//从高位除起  
    {  
        div=div*BASE+a[i];  
        a[i]=div/b;//b为除数  
        div%=b;  
    }  
}  
int main()  
{  
    int i,j,h[101][MAX];  
    memset(h[1],0,MAX*sizeof(int));//赋值,每一个都置为0   
    for(i=2,h[1][MAX-1]=1;i<=100;++i)//运用递归,并且h[1]=1;   
    {  
        memcpy(h[i],h[i-1],MAX*sizeof(int));//h[i]=h[i-1];按字节拷贝,保证了h[i]和h[i-1]指向数组的一致性   
        multiply(h[i],MAX,4*i-2);//h[i]*=(4*i-2);  
        divide(h[i],MAX,i+1);//h[i]/=(i+1);          
    }//递归得到前100项的卡特兰数!  
    while(cin>>i && i>=1 && i<=100)//输入i的值   
    {  

        for(j=0;j<MAX && h[i][j]==0;++j);//从0位开始搜索,找到不为0的第一个数   
        //在c语言中,EOF=-1;  
        printf("%d",h[i][j++]);//像是这个输出,就很妙了,第一位可能不足四位,就地输出!  
              for(;j<MAX;++j)  
        {  
        printf("%04d",h[i][j]);//处在中间的值也可能没有四位,这时候要注意了,往左边加0,凑足4位,不然答案会出错!  
          
       }  
      
        printf("\n");  
    }  
    //system("pause");  
      
    return 0;  
}  
模板:::::
//h( n ) = ( ( 4*n-2 )/( n+1 )*h( n-1 ) );


#include<stdio.h>

//*******************************
//打表卡特兰数
//第 n个 卡特兰数存在a[n]中,a[n][0]表示长度;
//注意数是倒着存的,个位是 a[n][1] 输出时注意倒过来。 
//*********************************
int a[105][100];
void ktl()
{
    int i,j,yu,len;
    a[2][0]=1;
    a[2][1]=2;
    a[1][0]=1;
    a[1][1]=1;
    len=1;
    for(i=3;i<101;i++)
    {
        yu=0;
        for(j=1;j<=len;j++)
        {
            int t=(a[i-1][j])*(4*i-2)+yu;
            yu=t/10;
            a[i][j]=t%10;
        }    
        while(yu)
        {
            a[i][++len]=yu%10;
            yu/=10;
        }
        for(j=len;j>=1;j--)
        {
            int t=a[i][j]+yu*10;
            a[i][j]=t/(i+1);
            yu = t%(i+1);
        }        
        while(!a[i][len])
        {
            len--;
        }    
        a[i][0]=len;
    }    
    
}    
int main()
{
    ktl();
    int n;
    while(scanf("%d",&n)!=EOF)
    {
        for(int i=a[n][0];i>0;i--)
        {
            printf("%d",a[n][i]);
        }    
        puts("");
    }    
    return 0;
}



没有更多推荐了,返回首页