博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
[2019.1.15]BZOJ2152 聪聪可可
阅读量:5766 次
发布时间:2019-06-18

本文共 1554 字,大约阅读时间需要 5 分钟。

明显的点分治。

我们记\(num_{i,j}\)为点\(i\)及其子树中的节点与\(i\)之间的路径的和对3取余为\(j\)的种类数,显然\(0\le j<3\)

我们再记\(val_{u,v}\)为两个直接相连的点\(u,v\)之间的边权,我们把边权和距离全部对3取余。

我们发现这个东西很好计算。基础的树形DP就好了。

那么如何统计答案呢?

当我们将一个节点\(v\)合并到它的父节点\(u\)上时,会产生2种新的路径:

1.以\(u\)为一端,\(v\)及其子树中的节点为另一端的;

2.以\(u\)已经合并的子树中的节点为一端,\(v​\)及其子树中的节点为另一端的。

考虑第一种,种类数就是\(2\times num_{v,3-val_{u,v}}\),乘2是因为可以调换链的两端。

对于第二种,我们记\(sum_i\)为已经合并的子树中,到\(u\)的距离为\(i\)的节点的数量。

那么这一种的种类数就是\(2\times \sum_{i=0}^2num_{v,i}\times sum_{3-i-val_{u,v}}\)

最后输出的时候记得约分。

code:

#include
#define V(x) (x>2?x-3:(x<0?x+3:x))using namespace std;struct edge{ int t,v,nxt;}e[40010];int n,u,v,w,be[20010],cnt,vis[20010];long long num[20010][3],ans;void add(const int &x,const int &y,const int &val){ e[++cnt].t=y,e[cnt].v=val,e[cnt].nxt=be[x],be[x]=cnt;}void Upd(const int &x,const int &y,const int &ev,long long *sum){ num[x][ev]+=num[y][0],num[x][V(ev+1)]+=num[y][1],num[x][V(ev+2)]+=num[y][2],ans+=num[y][V(3-ev)]*2; for(int i=0;i<3;++i)ans+=num[y][i]*sum[V(V(3-ev)-i)]*2; sum[ev]+=num[y][0],sum[V(ev+1)]+=num[y][1],sum[V(ev+2)]+=num[y][2];}void dfs(const int &x){ vis[x]=1,num[x][0]=1,ans++; long long sum[3]; sum[0]=sum[1]=sum[2]=0; for(int i=be[x];i;i=e[i].nxt)!vis[e[i].t]?dfs(e[i].t),Upd(x,e[i].t,e[i].v,sum),0:0;}long long gcd(const long long &x,const long long &y){ return y?gcd(y,x%y):x;}void Print(const long long &x,const long long &y){ long long g=gcd(x,y); printf("%lld/%lld",x/g,y/g);}int main(){ scanf("%d",&n); for(int i=1;i

转载于:https://www.cnblogs.com/xryjr233/p/BZOJ2152.html

你可能感兴趣的文章
正则高级用法例子1
查看>>
iOS 时间
查看>>
26期20180719
查看>>
zookeeper 本地测试环境搭建
查看>>
apache2.4中输入网址未能指向index
查看>>
大数据是什么?华为云学院带你探索大数据之旅
查看>>
JAVA学习路线
查看>>
国内互联网公司漏洞响应平台
查看>>
百度员工揭秘8月算法变动
查看>>
13个不容错过的Java项目
查看>>
Ansible简介、自动化介绍
查看>>
银行、券商们的下一代App该往哪里走?
查看>>
【数据库】阿里云分析型数据库基本认识
查看>>
市场是衡量一切的标准
查看>>
七款android地图 [Map] 开源项目框架分类总汇
查看>>
编程珠玑:向量旋转(旋转交换)
查看>>
服务器之07源码安装LAMPMX(同台主机)
查看>>
我的友情链接
查看>>
oracle中以dba_、user_、v$_、all_、session_、index_开头的常用表和视图
查看>>
NUnit单元测试
查看>>