特殊组合—错排

/*  2.计算排错方案
  有一种非常特殊的排序问题——错排
  n个元素依次标号1,2,3...n，求每个元素都不在自己原来位置的排序数。
  Ai为数i在第i个位置的全体排序（1<=i<=n)。因数字i不动，因此
  |Ai|=(n-1)!;
  同理
  |Ai交Aj|=(n-2)!(1<=i,j<=n,i!=j);
  每个元素都不在自己原来位置的排列数
  |（非A1）交（非A2）...交（非A3）|=n!-c(n,1)(n-1)!+c(n-2)(n-2)!-...(-1)^(n+1)c(n+1)1!
                                  =n!(1-1/(1!)+1/(2!)-...(-1)^(n+1)(1/(n!))
  递推方法：
  设1,2.....n错排的数目为D(n);
  情况1：i与其他n-1个数之一互换，其余n-2个错排，根据乘法原理，共产生(n-1)D(n-2)个错排
  情况2：以外的n-1个数先错排，然后i与其中每个互换。根据乘法原理，共产生(n-1)D(n-2)个错排
  使用加法原理综合上述情况，可得出递推式：
  D(1)=0,D(2)=1
  D(n)=(n-1)(D(n-1)+D(n-2),n>2;
  注意：当n较大时，要用高精度
*/
/*UVA 10497输入-1结束*/
#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
#include<algorithm>
using namespace std;
__int64 m=10;
struct Bigint
{
	__int64 s[20000];int l;
	Bigint(){l=0;memset(s,0,sizeof(s));}
	void operator*=(int x)
	{
		__int64 d=0;
		for(int i=0;i<=l;i++)
		{
			d+=s[i]*x;
			s[i]=d%m;
			d/=m;
		}
		while(d)
		{
			s[++l]=d%m;
			d/=m;
		}
	}
	void print()
	{
		printf("%I64d",s[l]);
		for(int i=l-1;i>=0;i--)
			printf("%I64d",s[i]);
	}
	void set(__int64 a)
	{
		s[l]=a%m;a/=m;
		if(a)
		{
			l++;s[l]=a%m;
		}
	}
}dp[1000];
Bigint operator+(Bigint b,Bigint &a)
{
	__int64 d=0;
	b.l=b.l>a.l?b.l:a.l;
	for(int i=0;i<=b.l;i++)
	{
		b.s[i]+=d+a.s[i];
		d=b.s[i]/m;b.s[i]%=m;
	}
	if(d) b.l++,b.s[b.l]=d;
	return b;
}
int n;
int main()
{
	dp[1].set(0);dp[2].set(1);
	for(int i=3;i<=800;i++)
	{
		dp[i]=dp[i-2]+dp[i-1];
		dp[i]*=(i-1);
	}
	while(scanf("%d",&n)==1)
	{
		if(n==-1) break;
		dp[n].print();printf("\n");
	}
	return 0;
}