hdu 1133

转帖自 http://xuyemin520.is-programmer.com/posts/26317.html

还有 http://hi.baidu.com/a363310925/blog/item/aced542d719b2b5d4fc22695.html
题意：M+N个人排队买票，票的单价是50￥，每个人只能买一张。 M个人拿50的去买，N个人拿100的去买，然后悲剧的是售票处开始的时候没有钱，所以如果拿100块买票人前面的拿50块买票的人小于或者等于用100块买票的人，这种排队方式就不合法，也就是不能顺利全部都买到票（因为没零钱找了）！

卡特兰数的应用

该题的公式：(C(m+n, n)-C(m+n, m+1))*m!*n! 化简即(m+n)!*(m-n+1)/(m+1)

m个人拿50，n个人拿100 ， 所以如果 n>m，那么排序方法数为 0 这一点很容易想清楚

现在我们假设 拿50的人用 ‘0’表示， 拿100的人用 1 表示。

如果有这么一个序列 0101101001001111..........
当第K个位置出现1的个数多余0的个数时就是一个不合法序列了
假设n=4 n=3的一个序列是：0110100 显然，它不合法， 现在我们把它稍微变化一下：
把第二个1（这个1前面的都是合法的）后面的所有位0变成1，1变成0
就得到 0111011 这个序列1的数量多于0的数量， 显然不合法， 但现在的关键不是看这个序列是不是合法的
关键是：它和我们的不合法序列 0110100 成一一对应的关系
也就是说任意一个不合法序列(m个0，n个1)， 都可以由另外一个序列(n-1个0和m+1个1)得到
另外我们知道，一个序列要么是合法的，要么是不合法的
所以，合法序列数量 = 序列总数量 - 不合法序列的总量
序列总数可以这样计算m+n 个位置中， 选择 n 个位置出来填上 1， 所以是 C(m+n, n)
不合法序列的数量就是： m+n 个位置中， 选择 m+1 个位置出来填上 1 所以是 C(m+n, m+1)
然后每个人都是不一样的，所以需要全排列 m! * n!

题目还涉及到了大数乘小数，还有大数除小数

#include<iostream>
#include<cstdio>
#include<cstring>
#define MAX 201
using namespace std;
int factor[205][MAX]={0};
int sim[201]={0};
int multiply(int s[],int Max,int b)
{
int ans=0,i;
for(i=Max;i>=1;i--)
{
ans+=s[i]*b;
s[i]=ans%10000;
ans=ans/10000;
}
return 0;
}
int div(int s[],int Max,int b)
{
int ans=0,t,i;
for(i=1;i<=Max;i++)
{
t=ans*10000+s[i];
s[i]=t/b;
ans=t%b;
}
return 0;
}
int getfactor()
{
int i;
factor[0][MAX-1]=factor[1][MAX-1]=1;
for(i=2;i<=203;i++)
{
memcpy(factor[i],factor[i-1],MAX*sizeof(int));
multiply(factor[i],MAX-1,i);
}
return 0;
}
int output(int *s,int k)
{
int i=1;
printf("Test #%d:\n",k);
while(s[i]==0&&i<MAX)
i++;
printf("%d",s[i++]);
for(;i<MAX;i++)
printf("%04d",s[i]);
printf("\n");
return 0;
}
int main()
{
int m,n,i,k=1;
getfactor();
while(scanf("%d %d",&m,&n),m+n)
{
memcpy(sim,factor[m+n],sizeof(int)*MAX);
if(n>m){
printf("Test #%d:\n",k++);
printf("0\n");                         //别忘记了 判断这种情况，
continue;
}
multiply(sim,MAX-1,m-n+1);
div(sim,MAX-1,m+1);
output(sim,k);
k++;
}
return 0;
}