AC日记——明明的烦恼 bzoj 1005

1005

思路：

　　prufer编码+组合数；

代码：

#include <bits/stdc++.h>
using namespace std;
#define maxn 1005
#define maxlen 4000500
struct int666 {
    int len;
    char ai[maxlen];
    void mul(int pos)
    {
        long long cnt=;
        for(int i=;i<len;i++)
        {
            cnt+=pos*(ai[i]);
            ai[i]=cnt%,cnt/=;
        }
        while(cnt)
        {
            ai[len]=cnt%;
            cnt/=,len++;
        }
        while(len>&&!ai[len-]) len--;
    }
    void print()
    {
        for(int i=len-;i>=;i--) putchar(ai[i]+'');
        putchar('\n');
    }
};
struct int666 ans;
int n,tot,ai[maxn],m,pi[maxn*],num,bi[maxn*];
bool if_p[maxn*];
inline void in(int &now)
{
    int if_z=;now=;
    char Cget=getchar();
    while(Cget>''||Cget<'')
    {
        if(Cget=='-')if_z=-;
        Cget=getchar();
    }
    while(Cget>=''&&Cget<='')
    {
        now=now*+Cget-'';
        Cget=getchar();
    }
    now*=if_z;
}
void euler(int lit)
{
    for(int i=;i<=lit;i++)
    {
        if(!if_p[i]) pi[++num]=i;
        for(int v=;pi[v]*i<=lit&&v<=num;v++)
        {
            if_p[i*pi[v]]=true;
            if(i%pi[v]==) break;
        }
    }
}
void apart(int now,int di)
{
    int temp=now;
    for(int v=;v<=num;v++)
    {
        now=temp;
        while(now)
        {
            bi[v]+=di*now/pi[v];
            now/=pi[v];
        }
    }
}
int main()
{
//    freopen("bzoj_1005.in","r",stdin);
//    freopen("bzoj_1005.out","w",stdout);
    ans.len=,ans.ai[]=;
    in(n),euler(maxn-);
    if(n==)
    {
        int pos;
        in(pos);
        if(pos<=) printf("");
        else printf("");
        return ;
    }
    for(int i=;i<=n;i++)
    {
        in(ai[i]);
        if(!ai[i])
        {
            printf("");
            return ;
        }
        if(ai[i]==-) m++;
        else tot+=--ai[i];
    }
    if(tot>n-)
    {
        puts("");
        return ;
    }
    apart(n-,),apart(n--tot,-);
    for(int i=;i<=n;i++) if(ai[i]) apart(ai[i],-);
    for(int i=;i<=num;i++)
    {
        for(int v=;v<=bi[i];v++) ans.mul(pi[i]);
    }
    for(int i=;i<=n--tot;i++) ans.mul(m);
    ans.print();
    return ;
}