hdu5909 Tree Cutting

假定1为本树的根，对于任意的一个联通子图，可以认为是从树上一个点向子树的儿子节点延伸产生的树。

从而考虑dp:
$f(x,j)$ 表示从x点向子树延伸而出异或和为j的连通子图的个数
考虑从 $f(p,j)$ 转移到 $f(x,j)$
用类似背包的方法：$h(i,j)$ 表示考虑前i个儿子从点x向下延伸产生的异或和为j的连通子图个数，从而有
$$h(i,j) = \sum_{k} h(i-1,k)*f(p,j \oplus k)$$
即 $h(i) = h(i-1) \oplus f(p) $
用fwt变换$O(nlogn)$实现异或卷积。fwt原理见Picks博客。

PS：

***判定叶子节点时发生了问题（我的判定是建双向边后当一个点出度为1时认为其为叶子，可实际上也有可能为根节点）***

以后在判定叶子节点时乖乖judge一下时候是否有儿子。

  1 #include <iostream>
  2 #include <cstdio>
  3 #include <cstring>
  4
  5 #define P 1000000007
  6 #define N 2010
  7 #define M (1<<12)
  8 #define LL long long
  9
 10 using namespace std;
 11
 12 struct edge{
 13     int x,to;
 14 }E[N<<1];
 15
 16 int n,m,totE,g[N],f[N][M];
 17
 18 void ade(int x,int y){
 19     E[++totE]=(edge){y,g[x]}; g[x]=totE;
 20     E[++totE]=(edge){x,g[y]}; g[y]=totE;
 21 }
 22
 23 int add(int a,int b){
 24     if(a+b>=P) return a+b-P;
 25     return a+b;
 26 }
 27
 28 int mul(int a,int b){
 29     return (int)((LL)a*(LL)b%(LL)P);
 30 }
 31
 32 int qpow(int x,int n){
 33     int ans=1;
 34     for(;n;n>>=1,x=mul(x,x)) if(n&1) ans=mul(ans,x);
 35     return ans;
 36 }
 37
 38 int inv2 = qpow(2,P-2);
 39
 40 // tf(A) = tf(tf(A0)-tf(A1) , tf(A0)+tf(A1) )
 41 void fwt(int a[],int l,int r){
 42     if(l==r) return;
 43     int mid=(l+r)>>1,len=(r-l+1)/2;
 44     fwt(a,l,mid);
 45     fwt(a,mid+1,r);
 46     for(int i=l;i<=mid;i++){
 47         int A0 = a[i],A1 = a[i+len];
 48         a[i] = add(A0, P-A1);
 49         a[i+len] = add(A0, A1);
 50     }
 51 }
 52
 53 //utf(A) = utf(utf((A0+A1)/2) , utf((A1-A0)/2))
 54 void ufwt(int a[],int l,int r){
 55     if(l==r) return;
 56     int mid=(l+r)>>1,len=(r-l+1)/2;
 57     for(int i=l;i<=mid;i++){
 58         int A0 = a[i],A1 = a[i+len];
 59         a[i] = mul(add(A0, A1),inv2);
 60         a[i+len] = mul(add(P-A0, A1),inv2);
 61     }
 62     ufwt(a,l,mid);
 63     ufwt(a,mid+1,r);
 64 }
 65
 66 #define p E[i].x
 67
 68 int ans[M],h[M],a[N];
 69
 70 void dp(int x,int fa){
 71     if(!E[g[x]].to && fa){
 72         f[x][a[x]]=1;
 73         ans[a[x]]++;
 74         return;
 75     }
 76     for(int i=g[x];i;i=E[i].to)
 77         if(p!=fa) dp(p,x);
 78     h[a[x]]=1;
 79     fwt(h,0,m-1);
 80     for(int i=g[x];i;i=E[i].to)
 81         if(p!=fa){
 82             f[p][0]=add(f[p][0],1);
 83             fwt(f[p],0,m-1);
 84             for(int j=0;j<m;j++) h[j]=mul(h[j],f[p][j]);
 85         }
 86     ufwt(h,0,m-1);
 87     for(int i=0;i<m;i++){
 88         f[x][i]=h[i];
 89         ans[i] = add(ans[i],h[i]);
 90         h[i]=0;
 91     }
 92 }
 93
 94 int main(){
 95 //    freopen("C.txt","r",stdin);
 96     int T;
 97     scanf("%d",&T);
 98     while(T--){
 99         memset(ans,0,sizeof(ans));
100         memset(f,0,sizeof(f));
101         memset(g,0,sizeof(g));
102         totE=0;
103         scanf("%d%d",&n,&m);
104         for(int i=1;i<=n;i++) scanf("%d",&a[i]);
105         for(int i=1,x,y;i<n;i++){
106             scanf("%d%d",&x,&y);
107             ade(x,y);
108         }
109         dp(1,0);
110         for(int i=0;i<m;i++)
111             printf("%d%c",ans[i],i==m-1? ‘\n‘:‘ ‘);
112     }
113     return 0;
114 }