bzoj 1005: [HNOI2008]明明的烦恼

#include<cstdio>
#include<iostream>
#define M 10005
int a[M],n,d[M],f[M],b[M],tot,m;
void jia(int a1)
{
for(int i=1;i<=d[0];i++)
for(;a1%d[i]==0;)
{
b[i]++;
a1/=d[i];
}
return;
}
void jian(int a1)
{
for(int i=1;i<=d[0];i++)
for(;a1%d[i]==0;)
{
b[i]--;
a1/=d[i];
}
return;
}
void cheng(int a1)
{
for(int i=1;i<=a[0];i++)
a[i]*=d[a1];
for(int i=1;i<=a[0];i++)
{
a[i+1]+=a[i]/10;
a[i]%=10;
}
for(;a[a[0]+1];)
{
a[0]++;
a[a[0]+1]=a[a[0]]/10;
a[a[0]]%=10;
}
}
int main()
{
scanf("%d",&n);
for(int i=2;i<=n;i++)
if(!f[i])
for(int j=2*i;j<=n;j+=i)
f[j]=1;
for(int i=2;i<=n;i++)
if(!f[i])
{
d[0]++;
d[d[0]]=i;
}
for(int i=2;i<=n-2;i++)
jia(i);
for(int i=1;i<=n;i++)
{
int a1;
scanf("%d",&a1);
if(a1==-1)
m++;
else
{
for(int j=1;j<a1;j++)
jian(j);
tot+=a1-1;
}
}
for(int i=1;i<=n-2-tot;i++)
{
jia(m);
jian(i);
}
a[0]=1;
a[1]=1;
for(int i=1;i<=d[0];i++)
for(;b[i]>0;b[i]--)
cheng(i);
for(int i=a[0];i;i--)
printf("%d",a[i]);
return 0;
}