[bzoj1005][HNOI2008][明明的烦恼] (高精度+prufer定理)

3
1
-1
-1

2

#include <stdio.h>
#include <memory.h>
#define ll long long
struct num {
    int top;
    ll x[500];
    num(ll a=0) {
        memset(x,0,sizeof x );
        top=1;
        x[1]=a; }
    num operator*(const num b) {
        num tmp;
        for(int i=1; i<=top; i++)
            for(int j=1; j<=b.top; j++)
                tmp.x[i+j-1]+=x[i]*b.x[j],
                tmp.x[i+j]+=tmp.x[i+j-1]/100000000,
                tmp.x[i+j-1]%=100000000;
        tmp.top=top+b.top;
        if(!tmp.x[tmp.top])
            tmp.top--;
        return tmp; } }ans(1);
void Q_pow(int x,int p) {
    num tmp(x);
    for(;p;p>>=1,tmp=tmp*tmp)
        if(p&1)
            ans=ans*tmp; }
int n,m,lef,cnt[1010];
void factorZ(int x,int d) {
    for(int i=2;i*i<=x;i++)
        while(!(x%i))
            cnt[i]+=d,
            x/=i;
    if(x^1)cnt[x]+=d; }
int main() {
    scanf("%d",&n);
    lef=n-2;
    for(int i=2;i<=lef;i++)
        factorZ(i,1);
    for(int i=1;i<=n;i++) {
        int x;
        scanf("%d",&x);
        if(~x) {
            if(x>1) {
                lef-=x-1;
                for(int i=2;i<=x-1;i++)
                    factorZ(i,-1); } }
        else ++m; }
    for(int i=2;i<=lef;i++)
        factorZ(i,-1);
    factorZ(m,lef);
    for(int i=2;i<=n;i++)
        if(cnt[i])
            Q_pow(i,cnt[i]);
    printf("%lld",ans.x[ans.top]);
    for(int i=ans.top-1;i;i--)
        printf("%08lld",ans.x[i]);
    return 0; }