P3899 [湖南集训]谈笑风生 主席树解决二维数点

传送门

文章目录

  • 题意:
  • 思路:

题意:

在这里插入图片描述

思路:

由于a,ba,ba,b都比ccc厉害,那么a,ba,ba,b一定是某个是某个的祖先。那么就分为两种情况了:
(1)(1)(1) bbbaaa上面,约定depth[1]=1depth[1]=1depth[1]=1,此时答案显然为min(depth[a]−1,k)∗(se[a]−1)min(depth[a]-1,k)*(se[a]-1)min(depth[a]1,k)(se[a]1)
(2)(2)(2)bbbaaa的下下面,这个时候就不是那么容易搞了,问题转化成我们要求aaa这个子树中depthdepthdepth范围在[depth[a]+1,depth[a]+k][depth[a]+1,depth[a]+k][depth[a]+1,depth[a]+k]内的所有点的se[i]−1se[i]-1se[i]1,先考虑暴力怎么写呢?显然我们可以暴力对这颗子树的深度建线段树,这样就变成了区间查询的问题了。而这样复杂度是肯定不行的,所以我们考虑建可持久化线段树,根据树的dfsdfsdfs序建主席树,以节点深度为下标,那么查询就变成了query(root[dfn[p]],root[dfn[p]+se[p]−1],1,n,depth[p]+1,min(depth[p]+k,n))query(root[dfn[p]],root[dfn[p]+se[p]-1],1,n,depth[p]+1,min(depth[p]+k,n))query(root[dfn[p]],root[dfn[p]+se[p]1],1,n,depth[p]+1,min(depth[p]+k,n)),再加上(1)(1)(1)的答案即可。

我们还可以将其转换成二维数点的问题。
通过以上分析不难发现我们要找的点就是
深度范围是[depth[p]+1,depth[p]+k][depth[p]+1,depth[p]+k][depth[p]+1,depth[p]+k]dfsdfsdfs序范围是[dfn[p],dfn[p]+se[p]−1][dfn[p],dfn[p]+se[p]-1][dfn[p],dfn[p]+se[p]1],那么我们以深度建立xxx轴,以dfsdfsdfs序建立yyy轴,让后统计就好啦。

主席树O(nlogn)O(nlogn)O(nlogn)

//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid (tr[u].l+tr[u].r>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;const int N=1000010,mod=1e9+7,INF=0x3f3f3f3f;
const double eps=1e-6;int n,m;
int root[N],tot,idx;
int dfn[N],se[N],depth[N];
vector<int>v[N];
struct Node
{int l,r;LL sum;
}tr[N*40];void insert(int p,int &q,int l,int r,int pos,int x)
{q=++tot; tr[q]=tr[p];tr[q].sum+=x;if(l==r) return;int mid=l+r>>1;if(pos<=mid) insert(tr[p].l,tr[q].l,l,mid,pos,x);else insert(tr[p].r,tr[q].r,mid+1,r,pos,x);
}LL query(int p,int q,int l,int r,int ql,int qr)
{if(l>=ql&&r<=qr) return tr[q].sum-tr[p].sum;LL ans=0;int mid=l+r>>1;if(ql<=mid) ans+=query(tr[p].l,tr[q].l,l,mid,ql,qr);if(qr>mid) ans+=query(tr[p].r,tr[q].r,mid+1,r,ql,qr);return ans;
}void dfs1(int u,int fa)
{se[u]=1; dfn[u]=++idx;depth[u]=depth[fa]+1;for(auto x:v[u]) if(x!=fa) dfs1(x,u),se[u]+=se[x];
}void dfs2(int u,int fa)
{insert(root[dfn[u]-1],root[dfn[u]],1,n,depth[u],se[u]-1);for(auto x:v[u]) if(x!=fa) dfs2(x,u);
}int main()
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);scanf("%d%d",&n,&m);for(int i=1;i<=n-1;i++){int a,b; scanf("%d%d",&a,&b);v[a].pb(b); v[b].pb(a);}dfs1(1,0); dfs2(1,0);while(m--){int p,k; scanf("%d%d",&p,&k);LL ans=1ll*min(depth[p]-1,k)*(se[p]-1);ans+=query(root[dfn[p]],root[dfn[p]+se[p]-1],1,n,depth[p]+1,min(depth[p]+k,n));printf("%lld\n",ans);}return 0;
}
/**/

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

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

相关文章

多层图,dj+dp Gym 102501A Environment-Friendly

一般求最短路&#xff0c;限制某个条件 Gym 102501A Environment-Friendly Gym 102501A Environment-Friendly 题意&#xff1a;求最小的co2消耗量&#xff08;最短路可&#xff09; &#xff0c;有一个限制条件&#xff0c;路途的距离 不能超过B 思路&#xff1a;djdp 代码&…

打造 .NET Core 链接转发服务

我最近使用 .NET Core 2.2 造了个名为"Link Forwarder" &#xff08;链接转发器&#xff09;的 URL 转发服务&#xff0c;并已开源。目前预览版已部署到我的子域"go.edi.wang"。本文将分享我如何构建这个项目&#xff0c;以及我学到的东西。为了帮助大家了…

P3293 [SCOI2016]美味 主席树 + 伪01trie

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 看到异或的话&#xff0c;很容易想到用01trie来贪心的搞&#xff0c;但是这个题涉及区间问题&#xff0c;直接搞的话需要将[l,r][l,r][l,r]的数都插入trie里面&#xff0c;这样的复杂度显然…

Abp 0.18.0 正式发布! -ABP CLI,新模板和其他功能

ABP CLI, v0.18版本新模板和其他功能ABP v0.18已发布, 包含解决的80个issue, 550次提交.网站更改abp.io网站完全更新以突出ABP框架的目标和重要功能.文档和博客网址也会更改&#xff1a;abp.io/documents移至docs.abp.io.abp.io/blog转移到blog.abp.io.ABP CLIABP CLI(命令行界…

牛客练习赛74 E CCA的期望(算概率的技巧+floyd处理)

牛客地址 题目描述 是否经常有艺术创作的冲动&#xff0c;但却限于水平无法描绘&#xff1f;那就交给随机吧&#xff01; 给定一张 n 个点 m 条边的无向带边权连通图&#xff0c;点有颜色&#xff0c;为黑或白&#xff0c;保证无自环和重边。 定义一次操作为&#xff1a;随机选…

Codeforces Round #712 (Div. 2) D. 3-Coloring 交互 构造

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给一个n∗nn*nn∗n的格子染色&#xff0c;一共可以染三种颜色&#xff0c;每次都会给一种颜色&#xff0c;代表当前这种颜色不可以使用&#xff0c;染色要求相邻的颜色不能相同&#xff0c;让你给出一种染色…

短信验证码“最佳实践”

1、背景年初&#xff0c;从外地转移阵地到西安&#xff0c;转眼已两个多月。很久不写业务代码了&#xff0c;到了新公司&#xff0c;条件恶劣到前所未有&#xff0c;从需求&#xff0c;设计&#xff0c;架构&#xff0c;实现&#xff0c;实施&#xff0c;测试&#xff0c;bug修…

Codeforces Round #712 (Div. 2) E. Travelling Salesman Problem 思维转换

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你nnn个点&#xff0c;从iii到jjj的花费是max(ci,aj−ai)max(c_i,a_j-a_i)max(ci​,aj​−ai​)&#xff0c;求从111开始经过每个点再回到111的最小花费。 思路&#xff1a; 首先可以发现我们的起点在哪…

HDU - 5877 Weak Pair (dfs序+树状数组+离散化)

VJ地址 题意&#xff1a;给一个有根树给你&#xff0c;计算一下满足下列条件的序列对的数目 &#xff08;1&#xff09;u是v的祖先&#xff08;不能是它自己&#xff09; &#xff08;2&#xff09;a[v]*a[u]<k 思路&#xff1a;用DFS序分裂每一条链&#xff0c;使链上的点…

Codeforces Round #614 (Div. 2) D. Aroma‘s Search 暴力 + 思维

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你x0,y0,ax,ay,bx,byx_0,y_0,a_x,a_y,b_x,b_yx0​,y0​,ax​,ay​,bx​,by​&#xff0c;让后根据[ax∗xi−1bx,ay∗yi−1by][a_x*x_{i-1}b_x,a_y*y_{i-1}b_y][ax​∗xi−1​bx​,ay​∗yi−1​by​]构造出…

Insider Dev Tour 2019巡演中国站

Insider Dev Tour 2019全球巡演中国站落下帷幕&#xff0c;在线直播加六大城市会场&#xff0c;作为四十余位讲师的一份子&#xff0c;我在线分享了Microsoft Graph及Microsoft Teams开发平台的内容&#xff0c;很高兴与几千位开发者共同度过了一个特别的星期天&#xff0c;很不…

HDU - 6267 (概论/找规律/递推)

VJ地址 题目大意&#xff1a; 有n个节点 从0-&#xff08;n-1&#xff09;&#xff0c;连边的规律为 即i点的父亲只能是比i小的数&#xff0c;而且是随机的&#xff0c;现在随机选择应该一个节点作为根&#xff0c;求这子树的和的期望是多少。 思路&#xff1a;可以知道总共有…

Educational Codeforces Round 81 (Rated for Div. 2) D. Same GCDs 欧拉函数\莫比乌斯

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给定a,ma,ma,m&#xff0c;求满足gcd(a,m)gcd(ax,m)gcd(a,m)gcd(ax,m)gcd(a,m)gcd(ax,m)的xxx的个数&#xff0c;且0<x<m0<x<m0<x<m。 思路&#xff1a; 由辗转相除法得&#xff1a;gcd…

Microsoft Graph Toolkit 初探

在今年的Build大会上面, Microsoft Graph 产品组公开宣布了一套新的Web组件&#xff0c;Microsoft Graph Toolkit&#xff08;简称mgt&#xff09;&#xff0c;这套组件可以与任何前端开发平台无缝整合&#xff0c;通过几行代码就能实现基于Microsoft Graph的应用&#xff0c;而…

CodeCraft-21 and Codeforces Round #711 (Div. 2) 题解

先上链接CodeCraft-21 and Codeforces Round #711 (Div. 2) A&#xff1a; 从n开始往后找&#xff0c;不出几十个 一定能找到的&#xff0c;所以暴力就好了 void sovle(){cin>>n;while(1){ll kn;ll sum0;while(k) sumk%10,k/10;if(gcd(n,sum)>1) {cout<<n<&…

Educational Codeforces Round 81 (Rated for Div. 2) B. Infinite Prefixes 数学

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你个串sss&#xff0c;让后把它重复无限次得到ttt&#xff0c;定义前缀的价值为cnt0−cnt1cnt_0-cnt_1cnt0​−cnt1​&#xff0c;求ttt的前缀价值为xxx的前缀个数&#xff0c;若有无限多输出−1-1−1。 …

终于等到你!微软正式上线 Windows Terminal 预览版

前一段时间&#xff0c;一直在知乎、技术社区收到技术小伙伴们的终极拷问&#xff1a;微软Build 大会上提到的「6月中旬」要上Windows store 的 Windows Terminal 到底啥时候可以用到呀&#xff1f;有一次&#xff0c;我跑去我给 Windows Terminal 提的某个 issue 里问 Windows…

牛客练习赛76 E 牛牛数数(线性基加二分)

牛客地址 思路&#xff1a;全部组合异或&#xff0c;很容易想到使用线性基&#xff0c;正好线性基中有一个求第k小的用法&#xff0c;那我们可以二分来找 K是第几小的数&#xff0c;然后用总数减去。 #include <iostream> #include <cstdio> #include <fstream…

Educational Codeforces Round 81 (Rated for Div. 2) E. Permutation Separation 线段树 + dp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个打乱的排列&#xff0c;每个位置都各有一个价值&#xff0c;让你选择一个分界点&#xff0c;分成p1,p2,...,prp_1,p_2,...,p_rp1​,p2​,...,pr​和pr1,...,pn−1,pnp_{r1},...,p_{n-1},p_{n}pr1​,…

几种设计良好结构以提高.NET性能的方法

写在前面设计良好的系统&#xff0c;除了架构层面的优良设计外&#xff0c;剩下的大部分就在于如何设计良好的代码&#xff0c;.NET提供了很多的类型&#xff0c;这些类型非常灵活&#xff0c;也非常好用&#xff0c;比如List&#xff0c;Dictionary、HashSet、StringBuilder、…