BZOJ3899 仙人掌树的同构（圆方树+哈希）

　　考虑建出圆方树。显然只有同一个点相连的某些子树同构会产生贡献。以重心为根后（若有两个任取一个即可），就只需要处理子树内部了。

　　如果子树的根是圆点，其相连的同构子树可以任意交换，方案数乘上同构子树数量的阶乘即可。而若是方点，注意到其相邻的圆点在原树中是有序地在一个环上的，要产生同构只能旋转或翻转该环。并且因为一开始我们选择了重心为根，所以对于非重心的方点，将其所在的环旋转显然是无法产生贡献的。所以对于方点的所有孩子按环上顺序存储，其哈希值应以该顺序计算，正反取较小的，算贡献时对非重心点只考虑翻转，重心特判一下。判同构当然采取哈希。

　　调了一年发现只有点双写错了。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
#define ll long long
#define N 2010
#define P 1000000003
#define M 10010
#define ull unsigned long long
#define p1 509
#define p2 923
char getc(){char c=getchar();while ((c<'A'||c>'Z')&&(c<'a'||c>'z')&&(c<''||c>'')) c=getchar();return c;}
int gcd(int n,int m){return m==?n:gcd(m,n%m);}
int read()
{
    int x=,f=;char c=getchar();
    while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
    while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
    return x*f;
}
int n,m,p[N],t,dfn[N],low[N],stk[N],top,cnt,tot,ans=;
vector<int> BCC[N];
struct data{int to,nxt;
}edge[M];
void addedge(int x,int y){t++;edge[t].to=y,edge[t].nxt=p[x],p[x]=t;}
void tarjan(int k)
{
    dfn[k]=low[k]=++cnt;
    stk[++top]=k;
    for (int i=p[k];i;i=edge[i].nxt)
    if (dfn[edge[i].to]) low[k]=min(low[k],dfn[edge[i].to]);
    else
    {
        tarjan(edge[i].to);
        low[k]=min(low[k],low[edge[i].to]);
        if (low[edge[i].to]>=dfn[k])
        {
            tot++;
            while (stk[top]!=edge[i].to) BCC[tot].push_back(stk[top]),top--;
            BCC[tot].push_back(edge[i].to);top--;
            BCC[tot].push_back(k);
        }
    }
}
namespace blocktree
{
    int p[N],t,size[N];
    ull hash[N],a[N];
    struct data{int to,nxt;}edge[M];
    void addedge(int x,int y)
    {
        //cout<<x<<' '<<y<<endl;
        t++;edge[t].to=y,edge[t].nxt=p[x],p[x]=t;
        t++;edge[t].to=x,edge[t].nxt=p[y],p[y]=t;
    }
    void make(int k,int from)
    {
        size[k]=;
        for (int i=p[k];i;i=edge[i].nxt)
        if (edge[i].to!=from)
        {
            make(edge[i].to,k);
            size[k]+=size[edge[i].to];
        }
    }
    int findroot(int k,int from,int s)
    {
        int mx=;
        for (int i=p[k];i;i=edge[i].nxt)
        if (edge[i].to!=from&&size[edge[i].to]>size[mx]) mx=edge[i].to;
        if ((size[mx]<<)>s) return findroot(mx,k,s);
        else return k;
    }
    void dfs(int k,int from)
    {
        for (int i=p[k];i;i=edge[i].nxt)
        if (edge[i].to!=from) dfs(edge[i].to,k);
        if (k<=n)
        {
            int cnt=;
            for (int i=p[k];i;i=edge[i].nxt)
            if (edge[i].to!=from) a[++cnt]=hash[edge[i].to];
            sort(a+,a+cnt+);
            for (int i=;i<=cnt;i++) hash[k]=hash[k]*p1+a[i];
            hash[k]=hash[k]*p2+size[k];
            for (int i=;i<=cnt;i++)
            {
                int t=i;
                while (t<cnt&&a[t+]==a[i]) t++;
                int fac=;
                for (int j=;j<=t-i+;j++) fac=1ll*fac*j%P;
                ans=1ll*ans*fac%P;
                i=t;
            }
        }
        else if (k!=from)
        {
            if (BCC[k-n].size()>)
            {
                int pos;
                for (int i=;i<BCC[k-n].size();i++)
                if (BCC[k-n][i]==from) {pos=i;break;}
                ull hash1=,hash2=;int x=pos,y=pos;
                for (int i=;i<BCC[k-n].size();i++)
                {
                    x--;if (x<) x+=BCC[k-n].size();
                    hash1=hash1*p1+hash[BCC[k-n][x]];
                    y++;if (y==BCC[k-n].size()) y=;
                    hash2=hash2*p1+hash[BCC[k-n][y]];
                }
                if (hash1==hash2) ans=2ll*ans%P;
                hash[k]=min(hash1,hash2)*p2+size[k];
            }
            else for (int i=p[k];i;i=edge[i].nxt)
            if (edge[i].to!=from) hash[k]=hash[edge[i].to]*p2+size[k];
        }
        else
        {
            ull hash1=;
            for (int i=;i<BCC[k-n].size();i++)
            hash1=hash1*p1+hash[BCC[k-n][i]];
            int cnt=;
            for (int i=;i<BCC[k-n].size();i++)
            {
                int x=i,y=i;ull h1=,h2=;
                for (int j=;j<BCC[k-n].size();j++)
                {
                    h1=h1*p1+hash[BCC[k-n][x]];
                    h2=h2*p1+hash[BCC[k-n][y]];
                    x++;if (x==BCC[k-n].size()) x=;
                    y--;if (y<) y+=BCC[k-n].size();
                }
                if (h1==hash1) cnt++;
                if (h2==hash1) cnt++;
            }
            if (BCC[k-n].size()==) cnt/=;
            ans=1ll*ans*cnt%P;
        }
    }
    void solve()
    {
        make(,);
        int root=findroot(,,size[]);
        make(root,root);
        dfs(root,root);
        //cout<<root<<endl;
        //for (int i=1;i<=19;i++) cout<<size[i]<<' '<<hash[i]<<endl;
    }
}
int main()
{
#ifndef ONLINE_JUDGE
    freopen("bzoj3899.in","r",stdin);
    freopen("bzoj3899.out","w",stdout);
    const char LL[]="%I64d\n";
#else
    const char LL[]="%lld\n";
#endif
    n=read(),m=read();
    for (int i=;i<=m;i++)
    {
        int x=read(),y=read();
        addedge(x,y),addedge(y,x);
    }
    tarjan();
    for (int i=;i<=tot;i++)
        for (int j=;j<BCC[i].size();j++)
        blocktree::addedge(n+i,BCC[i][j]);
    blocktree::solve();
    cout<<ans;
    return ;
}