5.15 牛客挑战赛40 E 小V和gcd树 树链剖分 主席树 树状数组 根号分治
LINK:小V和gcd树
时限是8s 所以当时好多nq的暴力都能跑过.
考虑每次询问暴力 跳父亲 这样是nq的 4e8左右 随便过。
不过每次跳到某个点的时候需要得到边权 如果直接暴力gcd的话 nqlogn就过不了了。
这里有两种解决办法:
一种是比赛的时候队友想的 一种是比较容易想到的方法。
前者套用树链剖分 只对重儿子的边进行修改 这样每次修改的复杂度为qlog.
考虑查询的时候沿着重链向上跳 这样重边可以O(1)得到答案 轻边暴力。
那么就得到了一个 nq+qlog^2的做法了。
当然可以进行根号分治 对于度数>sqrt的点不管 对于度数<sqrt的暴力修改。
前者在暴力跳的时候暴力查询 这样前者最多遇到sqrt(n)个这样的点 后者也是要修改sqrt次。
所以总复杂度为 nq+qsqrt(n)log.
两种暴力都是基于nq的暴力的。
题解给了一个树上带修莫队的做法 不过不是很优秀复杂度接近nq的暴力 而且写起来估计也比较难写..
对于这类问题还是考虑主席树 修改带上树状数组。
不过这个修改是多个点的点权修改 考虑树链剖分。
只修改重边上的信息 轻边暴力 那么利用主席树套树状数组 就可以了。
一个比较naive的想法是 只对重链开树状数组 这样是不对的。
考虑每次查询 是笔直的链查询 所以我们直接利用树状数组维护dfs序即可解决这个问题。
总复杂度 qlog^3+qlog^2.这个做法还算是比较优美的。
const int MAXN=20010;
int n,Q,maxx,len,cnt,id,top1,top2;
int f[MAXN],d[MAXN],a[MAXN],w[MAXN],son[MAXN],dfn[MAXN],sz[MAXN],pos[MAXN];
int lin[MAXN],ver[MAXN<<1],nex[MAXN<<1],root[MAXN],ql[MAXN],qr[MAXN],top[MAXN];
inline int gcd(int a,int b){return b?gcd(b,a%b):a;}
inline void add(int x,int y)
{
ver[++len]=y;nex[len]=lin[x];lin[x]=len;
ver[++len]=x;nex[len]=lin[y];lin[y]=len;
}
struct wy
{
int l,r;
int sum;
}t[MAXN*900];
inline void dfs(int x,int fa)
{
d[x]=d[fa]+1;f[x]=fa;sz[x]=1;
go(x)if(tn!=fa)
{
dfs(tn,x);
sz[x]+=sz[tn];
if(sz[tn]>sz[son[x]])son[x]=tn;
}
}
inline void dp(int x,int fa)
{
top[x]=fa;dfn[x]=++cnt;pos[cnt]=x;
if(!son[x])return;
dp(son[x],fa);
go(x)if(tn!=f[x]&&tn!=son[x])dp(tn,tn);
}
inline void insert(int &p,int l,int r,int x,int y)
{
if(!p)p=++id;sum(p)+=y;
if(l==r)return;
int mid=(l+r)>>1;
if(x<=mid)insert(l(p),l,mid,x,y);
else insert(r(p),mid+1,r,x,y);
}
inline void insert(int x,int v,int y)
{
while(x<=n)
{
insert(root[x],1,maxx,v,y);
x+=x&(-x);
}
}
inline void modify(int x,int v)
{
if(son[x])
{
insert(dfn[son[x]],w[son[x]],-1);
w[son[x]]=gcd(v,a[son[x]]);
insert(dfn[son[x]],w[son[x]],1);
}
if(f[x])
{
insert(dfn[x],w[x],-1);
w[x]=gcd(v,a[f[x]]);
insert(dfn[x],w[x],1);
}
a[x]=v;
}
inline int query(int l,int r,int x)//查<=x的个数.
{
if(l==r)
{
int cnt=0;
rep(1,top1,i)cnt-=sum(ql[i]);
rep(1,top2,i)cnt+=sum(qr[i]);
return cnt;
}
int mid=(l+r)>>1;
if(x>mid)
{
int cnt=0;
rep(1,top1,i)cnt-=sum(l(ql[i]));
rep(1,top2,i)cnt+=sum(l(qr[i]));
rep(1,top1,i)ql[i]=r(ql[i]);
rep(1,top2,i)qr[i]=r(qr[i]);
return cnt+query(mid+1,r,x);
}
rep(1,top1,i)ql[i]=l(ql[i]);
rep(1,top2,i)qr[i]=l(qr[i]);
return query(l,mid,x);
}
inline int ask(int s1,int s2,int k)
{
if(s1==s2)return 0;
top1=top2=0;
while(s1)ql[++top1]=root[s1],s1-=s1&(-s1);
while(s2)qr[++top2]=root[s2],s2-=s2&(-s2);
return query(1,maxx,k);
}
inline int Task(int x,int y,int k)
{
int cnt=0;
int fx=top[x],fy=top[y];
while(fx!=fy)
{
if(d[fx]<d[fy])swap(x,y),swap(fx,fy);
cnt+=ask(dfn[fx],dfn[x],k);
int ww=gcd(a[fx],a[f[fx]]);
cnt+=(ww<=k);
x=f[fx];fx=top[x];
}
if(d[x]<d[y])swap(x,y);
cnt+=ask(dfn[y],dfn[x],k);
return cnt;
}
int main()
{
//freopen("1.in","r",stdin);
get(n);get(Q);maxx=1000000;
rep(1,n,i)get(a[i]);
rep(2,n,i)add(read(),read());
dfs(1,0);dp(1,1);
rep(2,n,i)insert(dfn[i],w[i]=gcd(a[i],a[f[i]]),1);
rep(1,Q,i)
{
int op,x;
get(op);get(x);
if(op==1)modify(x,read());
else
{
int get(v);
put(Task(x,v,read()));
}
}
return 0;
}
5.15 牛客挑战赛40 E 小V和gcd树 树链剖分 主席树 树状数组 根号分治的更多相关文章
- 5.15 牛客挑战赛40 C 小V和字符串 数位dp 计数问题
LINK:小V和字符串 容易想到只有1个数相同的 才能有贡献. 知道两个01串 那么容易得到最小步数 大体上就是 第一个串的最前的1和第二个串最前的1进行匹配. 容易想到设f[i][j]表示 前i位1 ...
- 5.15 牛客挑战赛40 B 小V的序列 关于随机均摊分析 二进制
LINK:小V的序列 考试的时候 没想到正解 于是自闭. 题意很简单 就是 给出一个序列a 每次询问一个x 问序列中是否存在y 使得x^y的二进制位位1的个数<=3. 容易想到 暴力枚举. 第一 ...
- 牛客练习赛40 A 小D的剧场 (思维dp)
链接:https://ac.nowcoder.com/acm/contest/369/A 题目描述 若你摘得小的星星 你将得到小的幸福 若你摘得大的星星 你将得到大的财富 若两者都能摘得 你将得到 ...
- 牛客练习赛40 C 小A与欧拉路(树的直径)
链接:https://ac.nowcoder.com/acm/contest/369/C 题目描述 小A给你了一棵树,对于这棵树上的每一条边,你都可以将它复制任意(可以为0)次(即在这条边连接的两个点 ...
- 牛客挑战赛 30 A 小G数数
题目链接:https://ac.nowcoder.com/acm/contest/375/A 分析:我写的时候竟然把它当成了DP....... 还建了个结构体DP数组,保存一二位,不知道当时脑子在抽啥 ...
- 牛客挑战赛40 VMware和基站 set 二分 启发式合并 区间覆盖
LINK:VMware和基站 一道 做法并不常见的题目 看起来很难写 其实set维护线段就可以解决了. 容易想到 第二个操作借用启发式合并可以得到一个很不错的复杂度 不过利用线段树维护这个东西 在区间 ...
- 【牛客挑战赛30D】小A的昆特牌(组合问题抽象到二维平面)
点此看题面 大致题意: 有\(S\)张无编号的牌,可以将任意张牌锻造成\(n\)种步兵或\(m\)种弩兵中的一种,求最后步兵数量大于等于\(l\)小于等于\(r\)的方案数. 暴力式子 首先我们来考虑 ...
- 牛客小白月赛13-J小A的数学题 (莫比乌斯反演)
链接:https://ac.nowcoder.com/acm/contest/549/J来源:牛客网 题目描述 小A最近开始研究数论题了,这一次他随手写出来一个式子,∑ni=1∑mj=1gcd(i,j ...
- 牛客小白月赛13 小A的回文串(Manacher)
链接:https://ac.nowcoder.com/acm/contest/549/B来源:牛客网 题目描述 小A非常喜欢回文串,当然我们都知道回文串这种情况是非常特殊的.所以小A只想知道给定的一个 ...
随机推荐
- 组合 a 标签与 canvas 实现图片资源的安全下载的方法与技巧
普通用户下载图片时只需一个「右键另存为」操作即可完成,但当我们做在线编辑器.整个 UI 都被自定义实现时,如何解决不同域问题并实现页面中图片资源的安全下载呢?本文就解决该问题过程中所涉及的正则表达式. ...
- sun.reflect.generics.reflectiveObjects.TypeVariableImpl cannot be cast to java.lang.Class
1.首先,问题出现的点是在泛型 我出现问题的原因是,和泛型有关系,要调整泛型 2.我把问题出现的过程描述一哈子 1.基础类 @tk.mybatis.mapper.annotation.Register ...
- 最大的位或 HDU - 5969 简单思维题
题目描述 B君和G君聊天的时候想到了如下的问题. 给定自然数l和r ,选取2个整数x,y满足l <= x <= y <= r ,使得x|y最大. 其中|表示按位或,即C. C++. ...
- 前端07 /jQuery初识
前端07 /jQuery初识 目录 前端07 /jQuery初识 1.jquery介绍 1.1 jquery的优势 1.2 jquery是什么? 1.3 jquery的导入 2.jQuery的使用 2 ...
- Elasticsearch备份数据
Elasticsearch备份数据 1.建立备份目录 POST _snapshot/my_backup/ { "type": "fs", "setti ...
- 安装调试mosquitto日志
1.下载mosquitto安装包https://mosquitto.org/download/2.安装其他依赖包 yum install gcc gcc-c++ libstdc++-devel yum ...
- ASP.NET Core静态文件处理源码探究
前言 静态文件(如 HTML.CSS.图像和 JavaScript)等是Web程序的重要组成部分.传统的ASP.NET项目一般都是部署在IIS上,IIS是一个功能非常强大的服务器平台,可以直接 ...
- JVM调优工具Arthas的使用
Arthas 是Alibaba开源的Java诊断工具,深受开发者喜爱.在线排查问题,无需重启:动态跟踪Java代码:实时监控JVM状态. Arthas 支持JDK6+,支持Linux/Mac/Wind ...
- 【软件安装】CentOS7安装Tengine_2_3_2(Nginx 1_17_0)
简单比较一下Tengine 和Nginx 背景 使用最新的软件,可以处理一些bug,文章对CentOS6不做介绍(不会用) 推荐使用Tengine,理由是淘宝再用,兼容 Nginx 可以随时切换 Te ...
- P.SDA1.DEV - 一个没有服务器的图床
图床特色 P.SDA1.DEV的愿景是为大家提供一个免费.长期稳定外链分享图片的选择. P.SDA1.DEV的主要特点有: 完全建构在Serverless云服务上,致力于提供(墙外)可用性99.9%的 ...