[51nod1673]树有几多愁

　　lyk有一棵树，它想给这棵树重标号。
　　重标号后，这棵树的所有叶子节点的值为它到根的路径上的编号最小的点的编号。
　　这棵树的烦恼值为所有叶子节点的值的乘积。
　　lyk想让这棵树的烦恼值最大，你只需输出最大烦恼值对1e9+7取模后的值就可以了。
　　注意一开始1号节点为根，重标号后这个节点仍然为根。

　　update：数据保证叶子节点个数<=20。

　Input
　　第一行一个数n(1<=n<=100000)。
　　接下来n-1行，每行两个数ai,bi(1<=ai,bi<=n)，表示存在一条边连接这两个点。
Output
　　一行表示答案

　　显然小的编号应该丢给深度大的点，也就是说，从小到大确定编号的话，一个点子树内的所有其他点都被确定了之后 这个点才会（并且一定要）被确定。

　　但具体叶子之间谁先谁后还是有影响的。。。

　　就直接状压一波，f[i]表示已经确定编号的叶子的状态为i时的最大烦恼值（叶子只要给了编号，对烦恼值的贡献就确定下来了）。

　　先把原树上一些没用的点删掉，只保留叶子和有多个儿子的节点（其实就是虚树...）

　　每次枚举一个状态的时候，直接在虚树上暴力求出到底哪些点的编号已被确定了。这样就知道下一个叶子的编号是什么...再枚举下一个确定的是哪个叶子并转移就好了。

　　因为答案很大，比较方案优劣的时候可以用double。。

　　时间复杂度O(2^n*虚树节点数)，虚树节点数大概就40个左右吧？

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<queue>
#include<cmath>
#include<cstdlib>
#include<bitset>
//#include<ctime>
#define ll long long
#define ull unsigned long long
#define ui unsigned int
#define d double
//#define ld long double
using namespace std;
const int maxn=102333,modd=1000000007;const d eps=1e-7;
struct zs{int too,pre;}e[maxn<<1],E[maxn];int tot,last[maxn],TOT,LAST[maxn];
int sz[maxn];bool leaf[maxn],gg[maxn];
d f[(1<<20)+23333];int g[(1<<20)+23333];
int i,j,k,n,m;

int ra;char rx;
inline int read(){
    rx=getchar(),ra=0;
    while(rx<'0')rx=getchar();
    while(rx>='0')ra=ra*10+rx-48,rx=getchar();return ra;
}


void dfs(int x,int fa){
    int son=0;
    for(int i=last[x];i;i=e[i].pre)if(e[i].too!=fa)
        dfs(e[i].too,x),son++,sz[x]+=sz[e[i].too];
    leaf[x]=!son,sz[x]++;
    gg[x]=!leaf[x]&&son==1;
}

inline void insert(int a,int b){
    e[++tot].too=b,e[tot].pre=last[a],last[a]=tot,
    e[++tot].too=a,e[tot].pre=last[b],last[b]=tot;
}
inline void ins(int a,int b){
    E[++TOT].too=b,E[TOT].pre=LAST[a],LAST[a]=TOT;
}
int a[maxn],cnt;int pos[23],LEAF;int sz1[maxn],got[maxn];int num[maxn];
void dfs2(int x,int _fa,int tmp){
    if(!gg[x]){
        a[++cnt]=x,num[cnt]=tmp;
        if(leaf[x])pos[LEAF++]=cnt;
        if(_fa)ins(_fa,cnt);
        _fa=cnt,tmp=0;
    }
    for(int i=last[x];i;i=e[i].pre)if(sz[e[i].too]<sz[x])dfs2(e[i].too,_fa,tmp+1);
}
int main(){
    n=read();
    for(i=1;i<n;i++)insert(read(),read());
    dfs(1,0),dfs2(1,0,1);
    
//    for(i=1;i<=cnt;i++)printf("  %d",num[i]);puts("");
//    for(i=0;i<LEAF;i++)printf("    %d",pos[i]);puts("");
    
    for(j=0;j<LEAF;j++)sz1[pos[j]]=1;
    for(j=cnt;j;j--)for(k=LAST[j];k;k=E[k].pre)sz1[j]+=sz1[E[k].too];
    
    f[0]=g[0]=1;int mx=1<<LEAF,st,tozt,tog;d tof;register int j,k;
    for(i=0;i<mx-1;i++){
        memset(got+1,0,cnt<<2);
        for(j=0;j<LEAF;j++)got[pos[j]]=(i&(1<<j))>0;
        
        for(j=cnt,st=1;j;st+=got[j]==sz1[j]?num[j]:0,j--)
            for(k=LAST[j];k;k=E[k].pre)got[j]+=got[E[k].too];
        tof=f[i]*st,tog=1ll*g[i]*st%modd;
//        printf("zt:%d   st:%d
",i,st);
        for(j=0;j<LEAF;j++)if(!(i&(1<<j))&& f[tozt=(i|(1<<j))]<tof )f[tozt]=tof,g[tozt]=tog;
    }printf("%d
",g[mx-1]);
}