学习笔记：prufer编码

概念

$prufer$编码是一种将无根树转换成序列的一种方式。$n$个结点的树（结点编号为$1...n$）可以唯一对应一个$n-2$位的$prufer$序列，$n-2$位的$prufer$序列也可以唯一对应一棵$n$个结点的树（结点编号为$1...n$）。

方法

每一次找到一个编号最小的度数为$1$的结点，把它的父节点记下来。其实最后会记下$n-1$个数，但是第$n-1$位是确定的$n$，所以就只保留$n-2$位。

例题

AcWing 2419

$TreetoPrufer$：先记录所有点的度数，用一个$j$从小到大枚举，如果遇到度数为$0$的点（因为这次输入的是父亲序列）就记录并将该点的父亲的度数减$1$。如果父亲的度数变成了$0$且比$j$小就记录并把父亲的父亲的度数减$1$，一直这样递归直到不满足条件。
$PrufertoTree$：一个点在$Prufer$序列中出现了几次就会有几个子节点，注意把固定的$n-1$位的$n$加上。然后用一个$j$从小到大枚举，如果遇到了没有子节点的结点就说明它是目前最小的叶节点，父亲就是目前的$Prufer$编码。同$TreetoPrufer$，它父亲的子节点会变少，如果它变成了没有子节点的结点就会变成一个叶节点，若此时编号比$j$小就是目前最小的，要记录并迭代。

#include<bits/stdc++.h>
using namespace std;
const int NN=1e5+4;
int d[NN],p[NN],f[NN];
int main()
{
    int n,m;
    scanf("%d%d",&n,&m);
    if(m==1)
    {
        for(int i=1;i<n;i++)
        {
            scanf("%d",&f[i]);
            d[f[i]]++;
        }
        for(int i=0,j=1;i<=n-2;j++)
        {
            while(d[j])
                j++;
            p[++i]=f[j];
            while(i<=n-2&&!(--d[p[i]])&&p[i]<j)
                p[++i]=f[p[i]];
        }
        for(int i=1;i<=n-2;i++)
            printf("%d ",p[i]);
    }
    else
    {
        for(int i=1;i<=n-2;i++)
        {
            scanf("%d",&p[i]);
            d[p[i]]++;
        }
        p[n-1]=n;
        d[n]++;
        for(int i=1,j=1;i<n;i++,j++)
        {
            while(d[j])
                j++;
            f[j]=p[i];
            while(i<=n-2&&!(--d[p[i]])&&p[i]<j)
                f[p[i]]=p[i+1],i++;
        }
        for(int i=1;i<n;i++)
            printf("%d ",f[i]);
    }
    return 0;
}

AcWing 2418

题目的两个性质说明了图是联通的且不存在有公共点的环。令$f_{i,j}$表示$i$个点且有$j$个联盟（一个点也看成是一个联盟），边界显然$f_{0,0}=1$。然后可以枚举环的大小$k$，现在要加入点$i$。则与$i$构成环的方案数有$C_{i-1}^{k-1}$，剩下的点都有$f_{i-k,j-1}$种方案。需要注意的是，环内排列的方案也有$\frac{(k-1)!}{2}$种（圆排列忽略旋转和翻转），最后还要选一个点与其他的联盟相连，乘上$k$。所以状态转移方程就是：$$f_{i,j}=\sum _{k=1}^{i-j+1}{f}_{i-k,j-1}\times C_{i-1}^{k-1}\times \frac{k!}{2}$$

那我们的主题Prufer编码去哪里了呢？不要慌，因为我们还漏了一个东西：树的形态。把原图缩点以后，就会是一棵树，但是不同的形态会带来不同的方案。由$Prufer$序列唯一对应树的性质可以得知，树的形态数量有$n^{n-2}$种。最后，因为我们在状态转移方程中算了一个与其他点相连的点，但是如果整个世界都是一个联盟就没有这个插头，所以当$j$为$1$时要特殊计算，就是$n$的圆排列数。

#include<bits/stdc++.h>
using namespace std;
const int NN=204;
int C[NN][NN],g[NN],f[NN][NN];
int main()
{
    int n,m;
    scanf("%d%d",&n,&m);
    for(int i=0;i<=n;i++)
        for(int j=0;j<=i;j++)
            if(!j)
                C[i][j]=1;
            else
                C[i][j]=(C[i-1][j-1]+C[i-1][j])%m;
    g[1]=1;
    g[3]=3;
    for(int i=4;i<=n;i++)
        g[i]=g[i-1]*i%m;
    f[0][0]=1;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=i;j++)
            for(int k=1;k<=i-j+1;k++)
                f[i][j]=(f[i][j]+1ll*f[i-k][j-1]*C[i-1][k-1]%m*g[k]%m)%m;
    int res=g[n-1],p=1;
    for(int i=2;i<=n;i++)
    {
        res=(res+1ll*f[n][i]*p%m)%m;
        p=p*n%m;
    }
    printf("%d",res);
    return 0;
}