题目
学习链接:
https://www.cnblogs.com/zhj5chengfeng/archive/2013/08/23/3278557.html
http://www.matrix67.com/blog/archives/682
一定是个整数。。
//难道数据范围有问题???1e3+5为什么不可以
#include<cstdio>
#include<iostream>
#include<cstring>
#define en '\n'
#define m(a,b) memset(a,b,sizeof a)
using namespace std;
typedef long long ll;
const int N=1e5+5;
int is[N],prime[N],tot;//质因子分解
int d[N],tmp[N],pri[N],cnt;//储存prufer需要的
int ans[N];//高精*低精 ans的空间需要>1e3;
void getprime()
{
for(int i=2;i<=N-5-2;++i)
{
if(!is[i]) prime[++tot]=i;
for(int j=1;j<=tot&&i*prime[j]<=N-5-2;++j)
{
is[i*prime[j]]=1;
if(i%prime[j]==0) break;
}
}
}
inline void cheng(int x)
{
for(int i=1;i<=ans[0];++i) ans[i]*=x;
for(int i=1;i<=ans[0];++i) ans[i+1]+=ans[i]/10,ans[i]%=10;
while(ans[ans[0]+1])
{
++ans[0];
ans[ans[0]+1]=ans[ans[0]]/10,ans[ans[0]]%=10;
}
}
int main()
{
getprime();
int n;scanf("%d",&n);
if(n==1)
{
scanf("%d",&d[1]);
puts(d[1]<=0?"1":"0");
return 0;
}
int sum=0;
for(int i=1;i<=n;++i)
{
scanf("%d",&d[i]);
if(!d[i]||d[i]==n) {puts("0");return 0;}//不可能存在无度的点
if(d[i]!=-1) tmp[++cnt]=d[i]-1,sum+=tmp[cnt];
}
if(n-2<sum) {puts("0");return 0;}//要求过高 不得已实现
for(int i=1;i<=cnt;++i)
{
for(int j=2;j<=tmp[i];++j)
{
int jj=j;
for(int k=1;prime[k]<=jj;++k)
while(jj%prime[k]==0) --pri[k],jj/=prime[k];
}
}
for(int i=n-2-sum+1;i<=n-2;++i)
{
int jj=i;
for(int k=1;prime[k]<=jj;++k)
while(jj%prime[k]==0) ++pri[k],jj/=prime[k];
}
ans[0]=1,ans[1]=1;
for(int i=1;i<=tot;++i)
{
while(pri[i])// 一定为正值!!!
{
--pri[i];
cheng(prime[i]);
}
}
for(int i=1;i<=n-2-sum;++i) cheng(n-cnt);
for(int i=ans[0];i>=1;--i) printf("%d",ans[i]);
puts("");
}