我这份代码已经奇怪到一定程度了~

洛谷上一直 $TLE$,但是本地造了几个数据都过了.

简单说一下题解:

先建出来点分树.
对于每一个询问,在点分树中尽可能向上跳祖先,看是否能够处理这个询问.
找到最高点的好处就是该点的询问可以全部由那个祖先来统计.
因为祖先到 $x$ 是合法的,而那个祖先会统计子树里所有的点,当然也包括 $x$ 的所有子树.
假设现在枚举到点 $x$ ,并处理 $x$ 上面的所有询问.
$x$ 子树中的每一个点都可以用一个三元组来表示:$(l,r,c)$ 代表该点到 $x$ 路径的值域在 $[l,r],$颜色为 $c.$
那么对于一个询问,就是查询所有在 $[l_{q},r_{q}]$ 中 $c$ 的不同种类.
我们先将 $l$ 从大到小排序,依次处理点和询问.
如果有多个点,那么显然 $r$ 小的优先级会更高,即后加入的同颜色的点如果 $r$ 更小就直接替换.
这么做就能保证所有颜色的点在当前局面只出现一次.
现在的问题就是统计 $(l,r,l,r)$ 这个矩形内点的数量.
对于这个问题,可以用 $O(logn)$ 的树状数组来进行数点,总时间复杂度为 $O(nlog^2n)$.

Code:

#include <bits/stdc++.h>
#define N 1000005
#define inf 100001
#define setIO(s) freopen(s".in","r",stdin) ,freopen(s".out","w",stdout)
using namespace std;
int edges,n,flag;
int hd[N],to[N<<1],nex[N<<1],val[N],rt[N<<2];
void add(int u,int v)
{
nex[++edges]=hd[u],hd[u]=edges,to[edges]=v;
}
namespace BIT
{
int C[N];
int lowbit(int t)
{
return t&(-t);
}
void update(int x,int v)
{
for(;x<N;x+=lowbit(x)) C[x]+=v;
}
int query(int x)
{
int re=0;
for(;x>0;x-=lowbit(x)) re+=C[x];
return re;
}
void re(int x)
{
for(;x<N;x+=lowbit(x)) C[x]=0;
}
};
struct Node
{
int l,r,x;
Node(int l=0,int r=0,int x=0):l(l),r(r),x(x){}
};
vector<Node>e[N];
struct Point
{
int l,r,id,val;
Point(int l=0,int r=0,int id=0,int val=0):l(l),r(r),id(id),val(val){}
};
vector<Point>F[N];
bool cmp(Point a,Point b)
{
return a.l==b.l?a.id>b.id:a.l>b.l;
}
int root,sn;
int mx[N],size[N],vis[N],Fa[N],answer[N],lsty[N],lstx[N];
void dfs(int u,int ff)
{
size[u]=1;
for(int i=hd[u];i;i=nex[i])
if(to[i]!=ff&&!vis[to[i]])
dfs(to[i],u),size[u]+=size[to[i]];
}
void getroot(int u,int ff)
{
size[u]=1,mx[u]=0;
for(int i=hd[u];i;i=nex[i])
if(to[i]!=ff&&!vis[to[i]])
getroot(to[i],u),size[u]+=size[to[i]],mx[u]=max(mx[u],size[to[i]]);
mx[u]=max(mx[u],sn-size[u]);
if(mx[u]<mx[root]) root=u;
}
void calc(int u,int ff,int Min,int Max,int rt)
{
Min=min(Min,u),Max=max(Max,u);
F[rt].push_back(Point(Min,Max,u,val[u])), e[u].push_back(Node(Min,Max,rt));
for(int i=hd[u];i;i=nex[i])
if(!vis[to[i]]&&to[i]!=ff)
calc(to[i],u,Min,Max,rt);
}
void prepare(int u)
{
vis[u]=1;
calc(u,0,u,u,u);
for(int i=hd[u];i;i=nex[i])
if(!vis[to[i]])
dfs(to[i],u),sn=size[to[i]],root=0,getroot(to[i],u),Fa[root]=u,prepare(root);
}
void Push(int u,int l,int r,int id)
{
for(int i=0;i<(int)e[u].size();++i)
{
if(e[u][i].l>=l&&e[u][i].r<=r)
{
F[e[u][i].x].push_back(Point(l,r,-1,id));
break;
}
}
}
void solve(int u)
{
int i,j;
sort(F[u].begin(),F[u].end(),cmp);
for(i=0;i<(int)F[u].size();++i)
{
Point p=F[u][i];
if(p.id==-1)
answer[p.val]=BIT::query(p.r);
else
{
if(!lsty[p.val]||p.r<=lsty[p.val])
{
if(lsty[p.val])
{
BIT::update(lsty[p.val],-1);
}
BIT::update(p.r,1);
lstx[p.val]=p.l, lsty[p.val]=p.r;
}
}
}
for(i=0;i<(int)F[u].size();++i)
if(F[u][i].id!=-1)
{
if(lsty[F[u][i].val]) BIT::re(lsty[F[u][i].val]);
lstx[F[u][i].val]=lsty[F[u][i].val]=0;
}
}
int main()
{
int i,j,m;
// setIO("input");
scanf("%d%d",&n,&m);
for(i=1;i<=n;++i) scanf("%d",&val[i]);
for(i=1;i<n;++i)
{
int x,y;
scanf("%d%d",&x,&y),add(x,y),add(y,x);
}
root=0,mx[0]=sn=n,getroot(1,0),prepare(root);
for(i=1;i<=m;++i)
{
int l,r,x;
scanf("%d%d%d",&l,&r,&x),Push(x,l,r,i);
}
for(i=1;i<=n;++i) if(F[i].size()) solve(i);
for(i=1;i<=m;++i) printf("%d\n",answer[i]);
return 0;
}

  

  

luogu 5311 [Ynoi2011]D1T3 动态点分治+树状数组的更多相关文章

  1. HDU 4918 Query on the subtree(动态点分治+树状数组)

    题意 给定一棵 \(n\) 个节点的树,每个节点有点权.完成 \(q\) 个操作--操作分两种:修改点 \(x\) 的点权.查询与 \(x\) 距离小于等于 \(d\) 的权值总和. \(1 \leq ...

  2. 【BZOJ-3730】震波 动态点分治 + 树状数组

    3730: 震波 Time Limit: 15 Sec  Memory Limit: 256 MBSubmit: 626  Solved: 149[Submit][Status][Discuss] D ...

  3. bzoj 4372 烁烁的游戏——动态点分治+树状数组

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4372 和 bzoj 3070 震波 是一个套路.注意区间修改的话,树状数组不能表示 dis ...

  4. bzoj 3730 震波——动态点分治+树状数组

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3730 查询一个点可以转化为查询点分树上自己到根的路径上每个点对应范围答案.可用树状数组 f ...

  5. bzoj 4372 烁烁的游戏 —— 动态点分治+树状数组

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4372 本以为和 bzoj3730 一样,可以直接双倍经验了: 但要注意一下,树状数组不能查询 ...

  6. bzoj 3730 震波 —— 动态点分治+树状数组

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3730 建点分树,每个点记两个树状数组,存它作为重心管辖的范围内,所有点到它的距离情况和到它在 ...

  7. BZOJ_3295_[Cqoi2011]动态逆序对_CDQ分治+树状数组

    BZOJ_3295_[Cqoi2011]动态逆序对_CDQ分治+树状数组 Description 对于序列A,它的逆序对数定义为满足i<j,且Ai>Aj的数对(i,j)的个数.给1到n的一 ...

  8. [BZOJ 3295] [luogu 3157] [CQOI2011]动态逆序对(树状数组套权值线段树)

    [BZOJ 3295] [luogu 3157] [CQOI2011] 动态逆序对 (树状数组套权值线段树) 题面 给出一个长度为n的排列,每次操作删除一个数,求每次操作前排列逆序对的个数 分析 每次 ...

  9. BZOJ 2683 简单题 cdq分治+树状数组

    题意:链接 **方法:**cdq分治+树状数组 解析: 首先对于这道题,看了范围之后.二维的数据结构是显然不能过的.于是我们可能会考虑把一维排序之后还有一位上数据结构什么的,然而cdq分治却可以非常好 ...

随机推荐

  1. 1. centos7 的安装

    选择上海时间 我们选择桌面版 选择手动配置分区 选择标准分区 修改主机名 开始安装 设置密码 设置用户名 接下来等待安装完 同样的操作再安装3台机器!!这里不多赘述. 安装完之后就重启 接下来给不同机 ...

  2. python简单验证码识别

    在学习python通过接口自动登录网站时,用户名密码.cookies.headers都好解决但是在碰到验证码这个时就有点棘手了:于是通过网上看贴,看官网完成了对简单验证码的识别,如果是复杂的请看大神的 ...

  3. 在Docker中部署ASP.NET Core 2.2

    ⒈新建一个ASP.NET Core2.2 Web程序 因为Windows的Docker和Linux的Docker有所不同,本次测试采用的是Linux的Docker,因此没有勾选启用Docker支持. ...

  4. JS实现把一个页面层数据传递到另一个页面的两种方式

    本博客整理了两种方式从一个页面层向另一个页面层传递参数. 一. 通过cookie方式 1. 传递cookie页面的html,此处命名为a.html 请输入用户名和密码: <input id=&q ...

  5. C# 面向对象8 值类型和引用类型

    值类型和引用类型 概念 示意图: 1.值类型,在栈中开辟一块空间,存储 2.引用类型,在堆中开辟一块空间,存储数据,然在栈中开辟一块空间存储堆中的数据的地址

  6. C# 面向对象4 构造函数

    构造函数 1.构造函数用来创建对象,并且可以在构造函数中对对象进行初始化. (给对象的每个属性依次的赋值) 2.构造函数是用来创建对象的特殊方法: 1.方法名和类名一样. 2.没有返回值,连void都 ...

  7. Seeker:一款可获取高精度地理和设备信息的工具分析

    Seeker是一款可以获取高精度地理和设备信息的工具.其利用HTML5,Javascript,JQuery和PHP来抓取设备信息,以及Geolocation接口实现对设备高精度地理位置的获取. See ...

  8. NumPy 简介及安装

    NumPy(Numerical Python) 是 Python 语言的一个扩展程序库,支持大量的维度数组与矩阵运算,此外也针对数组运算提供大量的数学函数库.NumPy 的前身 Numeric 最早是 ...

  9. linux服务器硬件信息查看

    1.linux 查看服务器序列号(S/N) [root@oss20hb106 ~]# dmidecode -t 1 # dmidecode 2.11 # SMBIOS entry point at 0 ...

  10. openstack Rocky系列之keystone:(二)keystone中API注册

    主要说一下initialize_application中的application_factory def loadapp(): app = application.application_factor ...