51Nod1601 完全图的最小生成树计数

传送门

我居然忘写题解啦！(记忆废)

不管怎么说，这题还算是一道好题啊……你觉得敦爷出的题会有水题么

……

这题比较容易把人误导到Boruvka算法之类的东西上去(我们机房去刚D题的人一开始大多也被误导了)，但仔细思考之后是可以发现问题的特殊性质的。

听说很多人是从Kruskal算法想到这道题的做法的？好吧我并不是，那我就写写我的思考过程好了……

记得算导上有一道思考题，判断一个最小生成树算法的正确性。那个算法是这样的：把当前图的点集随意划分成两半，递归两半后选出连接两个点集的边中权值最小的一条，得到最后的最小生成树。

这个算法显然是错的，因为最终的最小生成树中可能有两条连接当前层两个点集的边。但本题有特殊性，边权都是端点点权的异或值，也许可以把这个算法改造一下用到这道题中。

考虑对一个点集求最小生成树，由于边权是端点点权的异或，因此可以把所有点按照最高位划分成两半，一半最高位为0，另一半最高位为1。这样递归两半之后只选一条连接两半的最小权边就可以得到一个生成树。可以证明它是最小生成树，并且证明起来并不难：如果某一组解中有两条边都连接了两半，由于这两条边边权的最高位一定是1，而位于两半内的边边权最高位一定是0，因此把这两条边中的一条替换成两半内的边得到的解一定比当前优。

有了正确性，算法也就成型了：每次把当前点集按最高位划分为两半后递归两半，然后若两边均非空则把任意一条连接两半的最小权边加入最小生成树即可。可能有多条边都是最小权边，显然方案数应该是每层的方案数之积。

实现的时候对所有权值建一棵01-Trie，那么选边的过程就相当于对每个点计算它的左子树和右子树之间的贡献，再dfs一遍左右子树即可。每个数最多被dfs到$\frac{32^2}2=512$次，因此复杂度为$O(512n)$。

注意一个细节：在划分过程中如果递归到了叶子节点且此处点数$>1$，则说明有多个点权值相同，显然这些点随便连就行了，那么答案就应该乘上对应点数的无向完全图生成树的数量。这个在OEIS上可以找到，通项是$n^{n-2}$。

 #include<cstdio>
 #include<cstring>
 #include<cassert>
 #include<algorithm>
 using namespace std;
 const int maxn=,maxm=maxn<<,p=1e9+;
 void insert(int,int&);
 void solve(int,int);
 void dfs(int,int,int,int);
 int qpow(int,int);
 long long sum=;
 int sm[maxm]={},ch[maxm][]={{}},root=,cnt=;
 int n,a[maxn],x,ans=,mn,tmp;
 signed main(){
     scanf("%d",&n);
     for(int i=;i<=n;i++){
         scanf("%d",&x);
         insert(,root);
     }
     solve(,root);
     //assert(ans==1ll);
     printf("%lld\n%d",sum,ans);
     return ;
 }
 void insert(int k,int &rt){
     if(!rt)rt=++cnt;
     sm[rt]++;
     if(k==-)return;
     insert(k-,ch[rt][(x>>k)&]);
 }
 void solve(int k,int x){
     if(sm[x]<=)return;//printf("solve(%d,%d)\n",k,x);
     if(k==-){
         ans=(long long)ans*qpow(sm[x],sm[x]-)%p;
         return;
     }
     solve(k-,ch[x][]);
     solve(k-,ch[x][]);
     if(sm[ch[x][]]&&sm[ch[x][]]){
         mn=;
         tmp=;
         dfs(k-,ch[x][],ch[x][],<<k);//printf("solve(%d,%d)\n",k,x);printf("mn=%d tmp=%d\n",mn,tmp);
         sum+=mn;
         ans=(long long)ans*tmp%p;
     }
 }
 void dfs(int k,int x,int y,int now){
     if(!sm[x]||!sm[y])return;
     if(k==-){
         if(now<mn){
             mn=now;
             tmp=(long long)sm[x]*sm[y]%p;
         }
         else if(now==mn)tmp=(tmp+(long long)sm[x]*sm[y]%p)%p;
         return;
     }
     if(sm[ch[x][]]){
         if(sm[ch[y][]])dfs(k-,ch[x][],ch[y][],now);
         else dfs(k-,ch[x][],ch[y][],now|(<<k));
     }
     if(sm[ch[x][]]){
         if(sm[ch[y][]])dfs(k-,ch[x][],ch[y][],now);
         else dfs(k-,ch[x][],ch[y][],now|(<<k));
     }
 }
 int qpow(int a,int b){
     int ans=;
     for(;b;b>>=,a=(long long)a*a%p)if(b&)ans=(long long)ans*a%p;
     return ans;
 }