P4930「FJ2014集训」采药人的路径
题目:P4930「FJ2014集训」采药人的路径
思路:
这篇不算题解,是让自己复习的,什么都没说清楚。
很久没有写点分治了,以前为了赶课件学的太急,板子都没打对就照着题解写题,导致学得很不扎实。
这道题差不多是在郭老师的指导下一点点凑出来的,还是没能自己完整写出一道题,惭愧。
这道题大意是:给出一棵边权为0/1的树,求满足以下条件的路径总数:0的个数等于1的个数,且路径上存在一点到路径两端也满足该条件。
这种求路径总数的题,可以想到用点分治。
把0看作-1,就可以转化为路径边权和为0。
如果没有休息点这个条件是很容易统计的,加上这个条件我们就要多记一些信息了。
点分治的容斥写法应用范围有限,用另一种写法会更方便。
类似树形背包,我们在统计经过u点的路径条数时,维护了一个之前遍历过的子树上的路径集合,每次进入一棵新子树,记录一下现在走的链,和之前集合里的链拼在一起形成答案,最后把现在的链加入集合中。
所以对这道题,我们维护现在走的路径和之前的路径集合时都分作两类:有休息点和没有休息点。
现在有休息点的路径某段后缀和为0,有一段和为0等价于有两点前缀和相等,开桶记录前缀和,类似P3941 入阵曲那道题。
没有休息点的路径用路径集合中有休息点的路径匹配,有休息点的路径之前有没有休息点都能匹配。
现在走的有休息点的路径可能本身已经满足条件了,要记入答案。
注意:休息点也可以不在路径上某点,而在根节点,即两边各有一条“1、-1”的路径,要加上这种情况。
Code:
#include <bits/stdc++.h>
using namespace std;
const int N=3e5+5,inf=0x3f3f3f3f;
int n,rt,rtsiz,tot[2],d[2][N],siz[N];
long long ans;
bool del[N];
struct bucket{//桶+时间戳
int tim,vis[N],c[N];
inline void add(int pos,int val=1){
if(vis[n+pos]==tim) c[n+pos]+=val;
else vis[n+pos]=tim,c[n+pos]=val;
}
inline int ask(int pos){
if(vis[n+pos]==tim) return c[n+pos];
else return 0;
}
}cnt[2],vis;
int Top,ver[N],val[N],nxt[N],head[N];
inline void add(int x,int y,int z){
ver[++Top]=y;val[Top]=z;nxt[Top]=head[x];head[x]=Top;
}
void findrt(int u,int fa,int all){//找重心
siz[u]=1;
int mxsiz=0;
for(int i=head[u];i;i=nxt[i]){
int v=ver[i];
if(v==fa||del[v]) continue;
findrt(v,u,all);
siz[u]+=siz[v];
mxsiz=max(mxsiz,siz[v]);
}
mxsiz=max(mxsiz,all-siz[u]);
if(mxsiz<rtsiz) rtsiz=mxsiz,rt=u;
}
void getdis(int u,int fa,int dis){
int op=1;
if(vis.ask(dis)==0) vis.add(dis),op=0;//之前没有出现相同路径长度 不存在休息点
d[op][++tot[op]]=dis;//记录路径长
for(int i=head[u];i;i=nxt[i]){
int v=ver[i];
if(v==fa||del[v]) continue;
getdis(v,u,dis+val[i]);
}
if(!op) vis.add(dis,-1);//如果是第一次出现 就退栈 倒出桶
}
void calc(int u){
++cnt[0].tim;++cnt[1].tim;++vis.tim;//更新时间戳 相当于memset
for(int i=head[u];i;i=nxt[i]){
int v=ver[i];
if(del[v]) continue;
tot[0]=tot[1]=0;
getdis(v,u,val[i]);
for(int j=1;j<=tot[0];++j) ans+=cnt[1].ask(-d[0][j]);//没有休息点的路径 用以前有休息点的路径匹配
for(int j=1;j<=tot[1];++j) ans+=cnt[0].ask(-d[1][j])+cnt[1].ask(-d[1][j]);//有休息点的路径 之前有没有休息点都能匹配
for(int j=1;j<=tot[0];++j) if(d[0][j]==0) ans+=cnt[0].ask(0);//特判
for(int j=1;j<=tot[0];++j) cnt[0].add(d[0][j]);//把当前路径加入路径集合
for(int j=1;j<=tot[1];++j) {
if(d[1][j]==0) ++ans;//注意当前路径也可能贡献答案
cnt[1].add(d[1][j]);
}
}
}
void dfs(int u){
del[u]=true;//删掉u点
calc(u);//计算过u点路径条数
for(int i=head[u];i;i=nxt[i]){
int v=ver[i];
if(del[v]) continue;
rtsiz=inf;
findrt(v,u,siz[v]);
dfs(rt);
}
}
int main(){
scanf("%d",&n);
for(int i=1,u,v,w;i<=n-1;++i){
scanf("%d%d%d",&u,&v,&w);
if(!w) w=-1;
add(u,v,w);
add(v,u,w);
}
rtsiz=inf;
findrt(1,0,n);
dfs(rt);
printf("%lld\n",ans);
return 0;
}
P4930「FJ2014集训」采药人的路径的更多相关文章
- 「FJ2014集训」采药人的路径
啦啦啦 来写一篇题解 洛谷链接: P4930 「FJ2014集训」采药人的路径 统计路径?嗯往点分治上想. 把0和1转化为-1和1,求和完dis为0的路径就是阴阳平衡的路径了. 如果题目没有限制要有中 ...
- 「BZOJ3694」「FJ2014集训」最短路
「BZOJ3694」「FJ2014集训」最短路 首先树剖没得说了,这里说一下并查集的做法, 对于一条非树边,它会影响的点就只有u(i),v(i)到lca,对于lca-v的路径上所有点x,都可通过1-t ...
- 「FJ2014集训 采药人的路径」
题目 考虑一下把\(0\)看成\(-1\),那么就是找到一条边权和为\(0\)的路径,且这条路径可以被分成两段,边权和都是\(0\) 没有第二个限制就是点分裸题了 其实有了第二个限制还是点分裸题 考虑 ...
- bzoj3697_FJ2014集训_采药人的路径_solution
小道士的矫情之路: 点分治, 对于每个子树,处理其内经过根(重心)的路径,然后递归下一层子树: 如何处理经过根的合法路径 合法有两个要求: 把输入的0改成-1后 1.len=0; 2.存在一个点i使被 ...
- [LOJ#2330]「清华集训 2017」榕树之心
[LOJ#2330]「清华集训 2017」榕树之心 试题描述 深秋.冷风吹散了最后一丝夏日的暑气,也吹落了榕树脚下灌木丛的叶子.相识数年的Evan和Lyra再次回到了小时候见面的茂盛榕树之下.小溪依旧 ...
- Loj #2331. 「清华集训 2017」某位歌姬的故事
Loj #2331. 「清华集训 2017」某位歌姬的故事 IA 是一名会唱歌的女孩子. IOI2018 就要来了,IA 决定给参赛选手们写一首歌,以表达美好的祝愿.这首歌一共有 \(n\) 个音符, ...
- Loj #2324. 「清华集训 2017」小 Y 和二叉树
Loj #2324. 「清华集训 2017」小 Y 和二叉树 小Y是一个心灵手巧的OIer,她有许多二叉树模型. 小Y的二叉树模型中,每个结点都具有一个编号,小Y把她最喜欢的一个二叉树模型挂在了墙上, ...
- Loj #2321. 「清华集训 2017」无限之环
Loj #2321. 「清华集训 2017」无限之环 曾经有一款流行的游戏,叫做 *Infinity Loop***,先来简单的介绍一下这个游戏: 游戏在一个 \(n \times m\) 的网格状棋 ...
- Loj 2320.「清华集训 2017」生成树计数
Loj 2320.「清华集训 2017」生成树计数 题目描述 在一个 \(s\) 个点的图中,存在 \(s-n\) 条边,使图中形成了 \(n\) 个连通块,第 \(i\) 个连通块中有 \(a_i\ ...
随机推荐
- 《2018年云上挖矿态势分析报告》发布,非Web类应用安全风险需重点关注
近日,阿里云安全团队发布了<2018年云上挖矿分析报告>.该报告以阿里云2018年的攻防数据为基础,对恶意挖矿态势进行了分析,并为个人和企业提出了合理的安全防护建议. 报告指出,尽管加密货 ...
- 洛谷P3306 随机数生成器
题意:给你一个数列,a1 = x,ai = (A * ai-1 + B) % P,求第一个是t的是哪一项,或者永远不会有t. 解:循环节不会超过P.我们使用BSGS的思想,预处理从t开始跳√P步的,插 ...
- (2)卷积 & 卷积和
参考资料:<信号与系统(第二版)> 杨晓非 何丰 连续信号的是卷积积分,离散信号的是卷积和. 脉冲分量 任意非周期信号,将横坐标分为若干个微小等分,得到Δτ为宽,f(kΔτ)为高的一系列微 ...
- tyvj 1423 GF和猫咪的玩具
传送门 解题思路 题目比较水,floyd求出最短路取个最小值即可.结果joyoi时限写错了..好像只有0ms才能过??突然发现加了快读就T不加就A,数据在10000以下的还是scanf快啊. 代码 # ...
- scrollbar 滚动条
滚动条样式:自定义元素的滚动条 <!DOCTYPE html> <html> <head lang="en"> <meta charset ...
- 20190807-RP-Explosion
如题,RP爆发后爆炸了. 话说在七夕这天考试? 那么? 黑暗又来临了? 没有人见过那一幕. 在如漆般胶着的黑暗中, Struggle?Dying? 用鲜血浇灌花朵,用死亡迎接明天. 考试过程: 看看三 ...
- 详解PPP模式下的产业投资基金运作【基金管理】
详解PPP模式下的产业投资基金运作[基金管理] 点击标题下「搏实资本」可快速关注 搏实资本 研究型的投资机构,实操型的专家团队 ﹌﹌﹌﹌﹌﹌﹌﹌﹌﹌﹌﹌﹌﹌﹌﹌﹌ 一.产业投资基金概述 ...
- 字符串无法分割 split无效: java split()使用“.” “\” "|" "*" "+"要转义
.是特殊字符 特殊字符需要转义. 改成split(“\\.”)
- 文本流向 layout-flow
作用与语法 文本流向layout-flow用来定义网页中的文本流向方式. 即排列方式,主要有两个属性,分别是horizonta (水平的) 和vertical-ideographic (垂直的). 文 ...
- Linux预习
目录 linux系统和unix系统的简介 linux系统和unix系统的简介 unix是什么:和widows一样 特点:多用户,多任务 同一时刻,多用户同时执行多项程序,互不干扰 GNU项目 就是一个 ...