Tree Cutting HDU - 5909

Tree Cutting HDU - 5909

题意:

一个无根树,n个点,n-1条边,每个节点有一个权值,一棵树的权值就是其节点(包含本身及其子节点)的权值的异或和;求价值为[0,m)的树有多少颗?(所谓的树其实就是原连通图的任意子图)
n<=1000
m<=2102^{10}210

题解:

不难想到树形dp
设dp[u][i]表示以u节点为根的价值为i的树的数量
能得到转移方程:
dp[u][j⨁k]=dp[u][j⨁k]+dp[u][j]∗d[v][k]dp[u][j\bigoplus k]=dp[u][j\bigoplus k]+dp[u][j]*d[v][k]dp[u][jk]=dp[u][jk]+dp[u][j]d[v][k]
u是v的父亲节点
不过这个式子直接算会超时,复杂度为O(n∗m∗m)O(n*m*m)O(nmm)
我们将式子变形:
dp[u][i]=∑j∗k=idp[u][j]∗d[v][k]dp[u][i]=\sum_{j*k=i}dp[u][j]*d[v][k]dp[u][i]=jk=idp[u][j]d[v][k]
这个就长得很像FWT
没错,就可以用FWT优化了
每次将dp[u][]和dp[v][]卷起来,并记录当前节点u为根的答案
优化后复杂度为O(nmlogm)O(nmlogm)O(nmlogm)

代码:

#include <bits/stdc++.h>
#include <unordered_map>
#define debug(a, b) printf("%s = %d\n", a, b);
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> PII;
clock_t startTime, endTime;
//Fe~Jozky
const ll INF_ll= 1e18;
const int INF_int= 0x3f3f3f3f;
void read(){};
template <typename _Tp, typename... _Tps> void read(_Tp& x, _Tps&... Ar)
{x= 0;char c= getchar();bool flag= 0;while (c < '0' || c > '9')flag|= (c == '-'), c= getchar();while (c >= '0' && c <= '9')x= (x << 3) + (x << 1) + (c ^ 48), c= getchar();if (flag)x= -x;read(Ar...);
}
template <typename T> inline void write(T x)
{if (x < 0) {x= ~(x - 1);putchar('-');}if (x > 9)write(x / 10);putchar(x % 10 + '0');
}
void rd_test()
{
#ifdef ONLINE_JUDGE
#elsestartTime = clock ();freopen("data.in", "r", stdin);
#endif
}
void Time_test()
{
#ifdef ONLINE_JUDGE
#elseendTime= clock();printf("\nRun Time:%lfs\n", (double)(endTime - startTime) / CLOCKS_PER_SEC);
#endif
}
const int P=1e9+7;
const int mod=1e9+7;
const int maxn1=(1<<13);
#define int long long
int val[1020];
vector<int>vec[2000];
int ans[maxn1];
int tmp[maxn1];
int dp[1020][maxn1];
int n,m;
/*
设dp[i][j]:表示以i为根的子树中异或和为j的数量 
*/
int qpow(int a,int b){int ans=1;while(b){if(b&1)ans=ans*a%mod;a=a*a%mod;b>>=1;}return ans%mod;
} 
void FWT(int x[],int t1,int t2,int len)
{const ll inv2= qpow(2,mod-2);for(int i=1;i<len;i<<=1)for(int j=0;j<len;j+=(i<<1))for(int k=0;k<i;k++){ll p=x[j+k],q=x[i+j+k];if(t1==0) x[i+j+k]=(q+P+t2*p)%P;  //orelse if(t1==1) x[j+k]=(p+P+t2*q)%P;   //andelse if(t1==2)   //xor{x[j+k]=(p+q)%P*(t2<0?inv2:1)%P;x[i+j+k]=(p+P-q)%P*(t2<0?inv2:1)%P;}     }
}
void say(int a[],int len){for(int i=0;i<len;i++){printf("a[%d]=%d\n",i,a[i]);}}
void solve(int a[],int b[],int len){
//	say(a,len);
//	say(b,len);FWT(a,2,1,len);FWT(b,2,1,len);for(int i=0;i<len;i++)a[i]=a[i]*b[i]%mod;FWT(a,2,-1,len);
}
void dfs(int u,int fa){dp[u][val[u]]=1;for(auto v:vec[u]){if(v==fa)continue;dfs(v,u);for(int j=0;j<m;j++){tmp[j]=dp[u][j];}solve(dp[u],dp[v],m);for(int j=0;j<m;j++){dp[u][j]=(tmp[j]+dp[u][j])%mod;}}for(int i=0;i<m;i++){ans[i]=(ans[i]+dp[u][i])%mod;}
}
signed main()
{
//    rd_test();int t;read(t);while(t--){read(n,m);for(int i=1;i<=n;i++)read(val[i]);for(int i=1;i<n;i++){int u,v;read(u,v);vec[u].push_back(v);vec[v].push_back(u);}memset(dp,0,sizeof dp);memset(ans,0,sizeof ans);dfs(1,1);for(int i=0;i<m;i++){if(i==0)cout<<ans[i];else cout<<" "<<ans[i];}cout<<endl;for(int i=1;i<=n;i++)vec[i].clear();}//Time_test();
}
/**/

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/315859.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Redis作者antirez:开源维护者的挣扎

这两天&#xff0c;一篇名为《开源维护者的挣扎》的文章被迅速顶至 Hacker News 首页&#xff0c;这是 Redis 作者 antirez 发布的最新博客。几个月前&#xff0c;一名开源项目的维护者向 antirez 发邮件&#xff0c;倾诉自己苦心维护项目多年&#xff0c;这或多或少带来了一些…

CF662C Binary Table

CF662C Binary Table 题意&#xff1a; 有一个 n 行 m 列的表格&#xff0c;每个元素都是 0/1 &#xff0c;每次操作可以选择一行或一列&#xff0c;把 0/1 翻转&#xff0c;即把 0 换为 1 &#xff0c;把 1 换为 0 。请问经过若干次操作后&#xff0c;表格中最少有多少个 1 …

P4111 [HEOI2015]小Z的房间

P4111 [HEOI2015]小Z的房间 题目描述 详见&#xff1a;P4111 [HEOI2015]小Z的房间 solution Matrix-Tree裸题。 题意为求一个有障碍的网格图的生成树个数。 因此对于每一个点向左边和上面的节点连边&#xff0c;Matrix-Tree定理求解生成树个数即可。 #include<bits/s…

黑科技抢先尝(续) - Windows terminal中WSL Linux 终端的极简美化指南

之前&#xff0c;本人写了两篇文章 黑科技抢先尝 | Windows全新终端初体验(附代码Build全过程) 和 程会玩 | 无需自行编译也能玩转 Windows Terminal&#xff0c;介绍了玩转Windows terminal的两种方式。今天这篇文章&#xff0c;主要介绍如何美化 Windows terminal 中 WSL 的 …

P3203 [HNOI2010]弹飞绵羊

P3203 [HNOI2010]弹飞绵羊 题目描述 详见&#xff1a;P3203 [HNOI2010]弹飞绵羊 solution 这是一道LCT的裸题。 但是我并不想用LCT解决此题&#xff08;In fact 是不会LCT ~QAQ&#xff09; 于是我们开始大力分块。 考虑把弹跳装置分块&#xff0c;我们每次需要知道在一…

L - Two Ants Gym - 102823L

L - Two Ants Gym - 102823L 题意&#xff1a; 有两个线段A&#xff0c;B&#xff0c;两个线段不会超过一个公共点&#xff0c; 你站在线段B上&#xff0c;整个平面你看不到的区域的面积(如图中S所在区域) 题解&#xff1a; 计算几何&#xff0c;恶心题。调了一个小时还是…

C# 8.0 中开启默认接口实现

当你升级到 C# 8.0 和 .NET Core 3.0 之后&#xff0c;你就可以开始使用默认接口实现的功能了。从现在开始&#xff0c;你可以在接口里面添加一些默认实现的成员&#xff0c;避免在接口中添加成员导致大量对此接口的实现崩溃。要写出并且正常使用接口的默认实现&#xff0c;你需…

P4396 [AHOI2013]作业

P4396 [AHOI2013]作业 题目描述 详见&#xff1a;P4396 [AHOI2013]作业 solution 莫队树状数组的裸题&#xff08;莫队分块&#xff0c;CDQ分治都可以&#xff0c;莫队线段树大概需要卡常&#xff09;。 时间复杂度 Code #include<bits/stdc.h> using namespac…

Codeforces Round #732 (Div. 1Div. 2)

Codeforces Round #732 (Div. 2) 题号题目知识点AAquaMoon and Two Arrays贪心BAquaMoon and Stolen String异或CAquaMoon and Strange SortDAquaMoon and ChessEAquaMoon and PermutationsFAquaMoon and Wrong CoordinateG1AquaMoon and Time Stop (easy version)G2AquaMoon …

在编写异步方法时,使用 ConfigureAwait(false) 避免使用者死锁

我在 使用 Task.Wait()&#xff1f;立刻死锁&#xff08;deadlock&#xff09; 一文中站在类库使用者的角度看 async/await 代码的死锁问题&#xff1b;而本文将站在类库设计者的角度来看死锁问题。阅读本文&#xff0c;我们将知道如何编写类库代码&#xff0c;来尽可能避免类库…

P1903 [国家集训队]数颜色 / 维护队列

P1903 [国家集训队]数颜色 / 维护队列 题目描述 详见&#xff1a;P1903 [国家集训队]数颜色 / 维护队列 solution 裸的带修莫队。 多记录一维询问时已经进行多少次修改。 用莫队思想暴力即可。 &#xff08;不开O2时需要卡常&#xff0c;洛谷评测两发一样的代码&#xf…

妄想集合(牛客练习赛90)

妄想集合(牛客练习赛90) 题意&#xff1a; 开始有 n 个可重集合&#xff0c;开始时每一个集合中都有一个数&#xff0c;有 m 个操作。 Quant l r x\text{Quant l r x}Quant l r x&#xff1a;往编号在 l∼rl\sim rl∼r 的每个集合中加入一个数 x。 Ask l r\text{Ask l r}Ask …

韩国政府计划从Windows 7迁移到Linux

据《韩国先驱报》的报道&#xff0c;韩国内政部于上周四表示&#xff0c;韩国政府计划将其用于办公的计算机操作系统从 Windows 迁移至 Linux。当然不是马上全面迁移 —— 内政部将会“第一个吃螃蟹”。他们准备在其 PC 上试运行 Linux&#xff0c;如果没有出现安全问题&#x…

P3332 [ZJOI2013]K大数查询

P3332 [ZJOI2013]K大数查询 题目描述 详见&#xff1a;P3332 [ZJOI2013]K大数查询 solution 整体二分树状数组 code #include<bits/stdc.h> using namespace std; const int MAXN1e550; typedef long long ll; ll n,m,Ans[MAXN],id[MAXN],quel[MAXN],quer[MAXN],f[…

黑科技抢先尝(续2) - Windows terminal中Powershell Tab的美化全攻略

温馨提示: 原文中含有一些外部链接&#xff0c;点击全文左下角的"阅读原文"体验会更佳喔~ 接着之前的文章 黑科技抢先尝(续) - Windows terminal中WSL Linux 终端的极简美化指南&#xff0c;依然假定你安装好了windows terminal预编译版本。这次我的目标是将PowerShe…

CF455B A Lot of Games

CF455B A Lot of Games 题意&#xff1a; Andrew,Fedor和Alex是三个善于创造的人。现在&#xff0c;他们发明了一种字符串双人游戏。 给出n个非空字符串。在游戏中&#xff0c;两位玩家一起建造一个单词。开始时&#xff0c;单词是一个空字符串。两位玩家轮流操作。一位玩家…

P1712 [NOI2016]区间

P1712 [NOI2016]区间 题目描述 P1712 [NOI2016]区间 Solution 尺取法线段树 一个显然的想法是按区间长度排序。 每一次多选取一个区间相当于区间覆盖次数加1&#xff0c;每一次少选取一个区间就有区间覆盖次数减1。 可以用线段树维护区间覆盖次数的最大值。 于是转化成…

.NET开发人员如何开始使用ML.NET

随着谷歌&#xff0c;Facebook发布他们的工具机器学习工具Tensorflow 2和PyTorch &#xff0c;微软的CNTK 2.7之后不再继续更新&#xff08;https://docs.microsoft.com/zh-cn/cognitive-toolkit/releasenotes/cntk_2_7_release_notes&#xff09;&#xff0c;Build 2019 微软也…

CF1271D Portals

CF1271D Portals 题意&#xff1a; 题意选自洛谷 题解&#xff1a; 首先要先发现性质&#xff1a;对于任何一个城堡u&#xff0c;如果要往u驻军&#xff0c;则在攻占编号最大的能向u行军的城堡后再驻军&#xff0c;答案一定不会变劣。你想想&#xff0c;如果后面有编号更大…

P2605 [ZJOI2010]基站选址

P2605 [ZJOI2010]基站选址 题目描述 详见&#xff1a;P2605 [ZJOI2010]基站选址 Solution 首先不难想到一个 的DP。 表示前个村庄选择了个基站的总费用。 考虑如何优化这个转移。 对于村庄&#xff0c;我们记录它覆盖范围内最靠前的村庄 和最靠后的村庄 。 倘若在…