POJ 2084 Game of Connections (思路+卡特兰数+大数处理)

传送门：POJ 2084

Sample Input

2
3
-1

Sample Output

2
5

题目大意：

由数字1~2*n形成的一个圆，然后用直线连接一对数字，每个数字都必须连接到另外一个数字上（注意这里一个数字只能和另一个数字连接，比如1连2，2连3是错误的）。如果没有两个线段是相交的，问有多少种连接方式。

思路：先写出前几个的结果，显然，是卡特兰数啦~开个玩笑~

当时做的时候完全没想到卡特兰数。这个题如果想到了是卡特兰数当然会节省很多时间，但是如果没想到也完全可以自己把结果推出来，下面就跟我一起推一下吧。

我们来推一下N=3，也就是有6个节点的时候的情况。单看1的情况来说可以与1连接的点有3个，分别是2，4，6。对这3种情况分类讨论并把结果相加则就是最后的答案。当1与2相连的时候剩下的4个点又可以随意连接，这4个点连接能产生的情况也就是N=2的时候的连接可能，即2种。说到这不知道你有没有发现什么，如果没有，没关系，我们继续往下看。

上面说到与1连接的只能有2，4，6，为什么呢？如上图所示，当1与3连接的时候，2就被孤立了，它与任何节点连接都会与1、3连接的边有交叉。也就是说，连线将剩余点分成两部分，两部分点的个数必须都是偶数才可以。请注意，2与1连接或者2与3连接是不合法的。

再来看1与4连接的时候，连接线把节点分成了两部分，两边的节点个数都为2，那么这时候的连接方式种类就是左边2个节点的连接方式 乘以 右边2个节点的连接方式，即1*1=1种。而两边的连接方式都是N=1的时候的种类。

我们可以发现，当对节点1进行讨论时，只有1和相邻节点（2，6节点）连接的时候会把剩下的节点分成1部分，其他情况都会分成2部分。我们只需要把每种情况下产生的连接方式数加起来就好了。例如N=3时，结果有：2+1*1+2=5种。

明白了这个思路我们就可以把N拓展到更大，每次计算的时候先算2种只把剩余点分成1组的情况，然后循环计算分成2组的情况。我们注意到共有2*N个节点，减去1和与1相连的节点还有2*N-2个节点。分成两组的情况时，一组最少可以有2个节点，每次可以增加2个，最多可以有2*N-4个节点。这样将结果相加就好了。

思路有了，但是我第一次提交却WA了，为什么呢？我发现卡特兰数的数值太大了，当N=35的时候就已经超了long long 型的范围，题目给出的N最大为100，当然会溢出了。所以我又自己实现了大数的乘法。当然如果嫌麻烦这里用java应该是比较容易写的。

我们发现代码实现的时候只用到了大数的加法和乘法，只要实现一下就可以了。为了简化运算，我规定1个变量可以存储10^9的数，比如大数12345678987654321，就可以用两个变量存储，一个存高位12345678，一个存低位987654321。最多可以存9位。这样一来需要开辟的数组就不必那么大了。计算原理和一个变量只存1位的时候是一样的。

具体实现：

我们用数组的低位来存大数的低位，数组的高位来存大数的高位。大数加法的时候只要对应位相加再加上低位的进位就好了。当然本位加完记得 mod 1e9 。至于大数的乘法，对于一个大数的低 i 位和另一个大数的低 j 位（i和j都是从0开始），则它们的结果会对低 i+j 位产生影响。这样先打表然后输出结果就好了。

注意：

输出结果的时候如果本位不足9位数则需要在前面补0。具体实现看代码吧，我给出来了较为详细的注释。

代码：

//卡特兰数，需要做大数处理 
#include<stdio.h>
#include<string.h>
typedef long long LL;

LL mod=1e9;    //mod用于计算本位存储9位的结果和计算进位 
LL a[101][10]; //a[i][j]用于保存N=i的时候第j位的结果 

void init()
{ //初始化函数打表 
	int i,j,x,y,c;
	//注意sum数组的大小应该是a[i]的2倍 
	LL sum[20]; //sum为临时数组，先将结果保存到sum再转存回a[i] 
	memset(a,0,sizeof(a));
	a[1][0]=1; //初始化N=1,N=2,N=3时候的结果 
	a[2][0]=2;
	a[3][0]=5;
	for(i=4;i<101;i++)
	{ //计算N从1~100的结果 
		memset(sum,0,sizeof(sum));
		for(j=0;j<10;j++)
		{ //这里是计算当连线将剩余点分成1组时候的结果 
			sum[j]+=(2*a[i-1][j]);  //本位值 
			sum[j+1]+=sum[j]/mod;   //本位对更高位的进位 
			sum[j]%=mod;			//本位%1e9,保证个变量保存9位数值 
		}
		for(j=1;j<i-1;j++) //j是分成2组的时候，两边剩余点对应的N的值 
			for(x=0;x<10;x++) //x是a[j]指向大数的第x位 
				for(y=0;y<10;y++) y是a[i-1-j]指向大数的第y位
				{
					sum[x+y]+=a[j][x]*a[i-1-j][y]; //本位是两个分组的种类相乘 
					sum[x+y+1]+=sum[x+y]/mod;      //对更高位的影响 
					sum[x+y]%=mod;                 //本位%1e9,保证个变量保存9位数值 
				}
		for(j=0;j<10;j++) a[i][j]=sum[j];     //转存回a[i] 
	} 
}

int main()
{
	int i,n;	
	init(); //先打表 
	while(~scanf("%d",&n))
	{
		if(n==-1) break;
		i=9;
		while(a[n][i]==0) i--; //从最高位开始寻找第一个不为0的数，则i就是结果的最高位 
		printf("%lld",a[n][i]); //输出最高位 
		i--;
		while(i>=0)
		{
			printf("%0.9lld",a[n][i]); //注意是%0.9lld，如果本位不足9位则前面补0 
			i--;
		}
		printf("\n");
	}
	return 0;
}