2729: HNOI2012 排队

题目
美其名曰月考题。
对于多重相邻限制，不能直接使用多重插空法，得考虑形成相邻再被分开的情况。
显然先放男生，再放老师，最后放女生。
如果放老师后老师不相邻，答案为：$n!A_{n+1}^2{A}_{n+3}^m$
如果放老师后相邻，得需要一个女生来分开他们，答案为：$n!2A_{n+1}^{1}m{A}_{n+2}^{m-1}$
相加化简后：$\frac{(n+1)!(n+2)!}{n-m+3!}[n(n+3)+2m]$
可以只需要高精度乘$int$.

AC Code:

#include<bits/stdc++.h>
#define maxn 10005
#define mod 100000000
#define LL long long
using namespace std;

int n,m;
int ct[maxn],len1;
LL a1[maxn];

void mul(LL *ans,int a,int &len)
{
	for(int i=0;i<len;i++) ans[i] = ans[i] * a;
	for(int i=0;i<len;i++) ans[i+1] += ans[i] / mod , ans[i] %= mod;
	for(;ans[len]>0;)
	{
		ans[len+1] += ans[len] / mod , ans[len] %= mod;
		len++;
	}
}

int main()
{
	scanf("%d%d",&n,&m);
	if((n==0 && m==0)) 
	{
		puts("0");
		return 0;
	}
	a1[0] = 1 , len1 = 1;
	for(int i=1;i<=n+1;i++) ct[i]+=2;
	ct[n+2] ++;
	for(int i=1;i<=n+3-m;i++) ct[i]--;
	for(int i=1;i<=n+3;i++)
		for(;ct[i];ct[i]--)
			mul(a1,i,len1);
	mul(a1,n*(n+3)+2*m,len1);
	for(int i=len1-1;i>=0;i--)
	{
		if(i == len1-1)
		{
			printf("%lld",a1[i]);
		}
		else
		{
			printf("%08lld",a1[i]);
		}
	}
}