Hdu 1029 Ignatius and the Princess IV

卡特兰数的模拟。

CODE：

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
using namespace std;

#define MAX 101
#define BASE 10000
int a[101][MAX];

void multiple(int a[], int size, int b)
{
    int i, j;
    int carry = 0;
    for(i = size-1; i >= 0; i--)
    {
        carry += a[i]*b;
        a[i] = carry%BASE;
        carry = carry/BASE;
    }
    return ;
}

void divide(int a[], int size, int b)
{
    int i, j;
    int div = 0;
    for(i = 0; i < size; i++)
    {
        div = div*BASE + a[i];
        a[i] = div/b;
        div = div%b;
    }
    return; 
}



void init()
{
    int i;
    memset(a, 0, sizeof(a));         //数组高位存放大数低位 
    a[1][MAX-1] = 1;
    for(i = 2; i < 101 ; i++)
    {
        memcpy(a[i], a[i-1], sizeof(a[i-1]));    //h(n) = h(n-1)
        multiple(a[i], MAX, 4*i-2);              // h(n) = h(n-1)*(4*n-2)
        divide(a[i], MAX, i+1);                     //h(n) = h(n-1)*(4*n-2)/(n+1)
    }
    return ;
}



void output(int n)
{
    int i;
    for(i = 0; a[n][i] == 0; i++);     //跳过数组前面的0 
    printf("%d", a[n][i++]);           //输出第一位非0数 
    for(; i < MAX; i++)
    {
        printf("%04d", a[n][i]);       //由于进位是10000，所以每位保持4位长度 
    }
    printf("\n");
}

int main()
{
    int n;
    init();
    while(~scanf("%d", &n))
    {
        output(n);
    }
    return 0;

}