P2624 [HNOI2008]明明的烦恼

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;

typedef long long LL;

inline int read()
{
    char c=getchar();int num=;bool f=;
    for(;!isdigit(c);c=getchar())
        f=c=='-'?:f;
    for(;isdigit(c);c=getchar())
        num=num*+c-'';
    return f?num:-num;
}

const int N=1e4+;

int n,m;
int d[N],tot;

int prime[N],cnt;
bool flag[N];
void init()
{
    for(int i=,t;i<N;++i)
    {
        if(!flag[i])
            prime[++cnt]=i;
        for(int j=;j<=cnt&&(t=i*prime[j])<N;++j)
        {
            flag[t]=;
            if(i%prime[j]==)
                break;
        }
    }
}

inline LL ksm(LL x,LL k)
{
    LL res=;
    while(k)
    {
        if(k&)
            res=res*x;
        x*=x;
        k>>=;
    }
    return res;
}

int num[N];

inline void chai(int x,int f)
{
    for(int i=;i<=cnt&&x>;++i)
    {
        while(x%prime[i]==)
        {
            num[i]+=f;
            x/=prime[i];
        }
    }
    return;
}

int ans[N],l,x;
inline void mul(int *a)
{
    x=;
    ans[]=l=;
    for(int i=;i<=cnt;++i)
    {
        for(int j=;j<=a[i];++j)
        {
            for(int k=;k<=l;++k)
            {
                ans[k]=ans[k]*prime[i]+x;
                x=ans[k]/;
                ans[k]%=;
            }
            while(x)
            {
                ans[++l]=x%,
                x/=;
            }
        }
    }
}

int main()
{
//    freopen("3.in","r",stdin);
//    freopen("466.out","w",stdout);
    init();
    n=read();
    if(n==)
    {
        n=read();
        if(n)
            putchar('');
        else
            putchar('');
        return ;
    }
    for(int i=;i<=n;++i)
    {
        d[i]=read();
        if(d[i]==)
        {
            puts("");
            return ;
        }
        if(d[i]==-)
            ++m;
        else
        {
            tot+=(d[i]-=);
            for(int j=;j<=d[i];++j)
                chai(j,-);
        }
    }
    if(tot>n-)
    {
        puts("");
        return ;
    }
    for(int i=;i<=n--tot;++i)
        chai(i,-);
    for(int i=;i<=n-;++i)
        chai(i,);
    if(m>)
    {
        for(int i=;i<=n--tot;++i)
            chai(m,);
    }
    mul(num);
    for(;l;--l)
        cout<<ans[l];
    return ;
}