hdu5909-Tree Cutting【FWT】
生活随笔
收集整理的這篇文章主要介紹了
hdu5909-Tree Cutting【FWT】
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
正題
題目鏈接:http://acm.hdu.edu.cn/showproblem.php?pid=5909
題目大意
給出nnn和mmm(m=2km=2^km=2k)。再給出一個大小為nnn的樹,每個點有點權,對于每個i∈[1,m)i\in[1,m)i∈[1,m)求有多少個聯通子圖的點權異或和為iii
1≤T≤10,1≤n≤1000,1≤m≤2101\leq T\leq 10,1\leq n\leq 1000,1\leq m\leq 2^{10}1≤T≤10,1≤n≤1000,1≤m≤210
解題思路
設fi,jf_{i,j}fi,j?表示iii的子樹中包含iii的聯通子圖里面,異或和為jjj的有多少個。那么轉移方程就是
fx,i=fx,i+∑jxork=ify,j×fy,kf_{x,i}=f_{x,i}+\sum_{j\ xor\ k=i}f_{y,j}\times f_{y,k}fx,i?=fx,i?+j?xor?k=i∑?fy,j?×fy,k?
這個是裸的FWTFWTFWT形式,所以直接做就好了
時間復雜度O(n2log?m)O(n^2\log m)O(n2logm)
比較老的題庫了,輸出格式限制是真的很嚴格
code
#include<cstdio> #include<cstring> #include<algorithm> #define ll long long using namespace std; const ll N=1030,P=1e9+7,inv2=(P+1)/2; struct node{ll to,next; }a[N<<1]; ll T,n,m,tot,ls[N],v[N]; ll f[N][N],ans[N]; void addl(ll x,ll y){a[++tot].to=y;a[tot].next=ls[x];ls[x]=tot;return; } void FWT(ll *f,ll op){for(ll p=2;p<=m;p<<=1)for(ll k=0,len=p>>1;k<m;k+=p)for(ll i=k;i<k+len;i++){ll x=f[i],y=f[i+len];f[i]=(x+y)*op%P;f[i+len]=(x-y)*op%P; }return; } void dfs(ll x,ll fa){f[x][v[x]]=1;FWT(f[x],1);for(ll i=ls[x];i;i=a[i].next){ll y=a[i].to;if(y==fa)continue;dfs(y,x);for(ll j=0;j<m;j++)f[x][j]=f[x][j]*f[y][j]%P;}FWT(f[x],inv2);for(ll j=0;j<m;j++)(ans[j]+=f[x][j])%=P;f[x][0]++;FWT(f[x],1);return; } signed main() {scanf("%lld",&T);while(T--){memset(ans,0,sizeof(ans));memset(ls,0,sizeof(ls));memset(f,0,sizeof(f));tot=0;scanf("%lld%lld",&n,&m);for(ll i=1;i<=n;i++)scanf("%lld",&v[i]);for(ll i=1;i<n;i++){ll x,y;scanf("%lld%lld",&x,&y);addl(x,y);addl(y,x);}dfs(1,1);for(ll i=0;i<m;i++){printf("%lld",(ans[i]%P+P)%P);if(i!=m-1)putchar(' '); }putchar('\n');}return 0; } 創作挑戰賽新人創作獎勵來咯,堅持創作打卡瓜分現金大獎總結
以上是生活随笔為你收集整理的hdu5909-Tree Cutting【FWT】的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 怎么在抖音上发超过15秒的长视频?
- 下一篇: 百度地图之拾取坐标系统的使用