bzoj1005 [HNOI2008]明明的烦恼 prufer+组合数学

直接考虑prufer序列，每个度数-1就是再n-2的序列中出现的次数

对于不确定的-1，每个没有填的位置都是有可能填的，所以直接离散用乘法原理

对于确定的就是组合数选位置，注意要开高精

对于高精除法的话就使用类似exlucas的方法提取因子，最后再乘

码：

#include<iostream>
#include<cstdio>
#include<vector>
using namespace std;
int ans[1999999],daan[1999999],i,j,n,lin,x,wz,jw,gs[1005][1005],cnt,k;
vector<int>v[1005];
void C(int a,int b)
{
	int i,j;
	if(b-a>a)
	{
	for(i=b;i>b-a;i--)
	{
		for(j=0;j<v[i].size();j++)
		ans[v[i][j]]+=gs[i][v[i][j]];	
	}
	for(i=1;i<=a;i++)
	{	for(j=0;j<v[i].size();j++)
		ans[v[i][j]]-=gs[i][v[i][j]];
	}	
	}else
{
	for(i=b;i>a;i--)
	{
		for(j=0;j<v[i].size();j++)
		ans[v[i][j]]+=gs[i][v[i][j]];	
	}
	for(i=1;i<=b-a;i++)
	{	for(j=0;j<v[i].size();j++)
		ans[v[i][j]]-=gs[i][v[i][j]];
	}	
}
}
int main()
{
	scanf("%d",&n);
for(i=1;i<=n;i++)
{int lin2=i;
	for(j=2;j*j<=n;j++)
	{
		if(lin2%j==0){lin=0;while(lin2%j==0)lin++,lin2/=j;   gs[i][j]=lin;v[i].push_back(j);}	
	}
	if(lin2!=1)gs[i][lin2]=1,v[i].push_back(lin2);
}
wz=n-2;
	for(i=1;i<=n;i++)
	{
		scanf("%d",&x);
		if(x!=1){ 
	if(x==-1)++cnt;
    else C(x-1,wz),wz-=(x-1);
                }
	}
daan[0]=1;
daan[1]=1;
ans[cnt]+=wz;
for(i=1;i<=n;i++)
{
for(j=1;j<=ans[i];j++)
{
	for(k=1;k<=daan[0];k++)
	{
	long long o=jw+daan[k]*i;
	daan[k]=o%1000;
	jw=o/1000;		
	}
	while(jw)daan[++daan[0]]=jw%1000,jw/=1000;		
}			
}	
	for(i=daan[0];i>=1;i--)
	{
		if(i!=daan[0])
		{
			if(daan[i]<100)printf("0");
				if(daan[i]<10)printf("0");	
		}
	printf("%d",daan[i]);	
	}
}