首页 > 代码库 > 【bzoj3251】树上三角形 朴素LCA+暴力
【bzoj3251】树上三角形 朴素LCA+暴力
题目描述
给定一大小为n的有点权树,每次询问一对点(u,v),问是否能在u到v的简单路径上取三个点权,以这三个权值为边长构成一个三角形。同时还支持单点修改。
输入
第一行两个整数n、q表示树的点数和操作数
第二行n个整数表示n个点的点权
以下n-1行,每行2个整数a、b,表示a是b的父亲(以1为根的情况下)
以下q行,每行3个整数t、a、b
若t=0,则询问(a,b)
若t=1,则将点a的点权修改为b
输出
对每个询问输出一行表示答案,“Y”表示有解,“N”表示无解。
样例输入
5 5
1 2 3 4 5
1 2
2 3
3 4
1 5
0 1 3
0 4 5
1 1 4
0 2 5
0 2 3
样例输出
N
Y
Y
N
题解
朴素LCA+暴力
一开始想到了一个$O(n\log^3n)$的数据结构算法,然后发现自己太naive了= =
由于点权是int范围内的,所以如果想让尽量多的边不构成三角形,那么它们的边权应该为1、1、2、3、5、8、...
这显然是斐波那契数列,而斐波那契数列是指数增长的,到第50项左右就爆int了。
所以可以直接拿出两个点之间的路径,当拿出的超过50个时直接判定能构成三角形,否则排序,暴力。
时间复杂度为$O(q·50·\log 50)$
#include <cstdio>#include <algorithm>#define N 100010using namespace std;int w[N] , head[N] , to[N << 1] , next[N << 1] , cnt , fa[N] , deep[N] , a[100] , tot;void add(int x , int y){ to[++cnt] = y , next[cnt] = head[x] , head[x] = cnt;}void dfs(int x){ int i; for(i = head[x] ; i ; i = next[i]) if(to[i] != fa[x]) fa[to[i]] = x , deep[to[i]] = deep[x] + 1 , dfs(to[i]);}bool judge(int x , int y){ int i; tot = 0; if(deep[x] < deep[y]) swap(x , y); while(deep[x] > deep[y]) { a[++tot] = w[x] , x = fa[x]; if(tot > 50) return 1; } while(x != y) { a[++tot] = w[x] , a[++tot] = w[y] , x = fa[x] , y = fa[y]; if(tot > 50) return 1; } a[++tot] = w[x] , sort(a + 1 , a + tot + 1); for(i = 3 ; i <= tot ; i ++ ) if(a[i] - a[i - 1] < a[i - 2]) return 1; return 0;}int main(){ int n , m , i , opt , x , y; scanf("%d%d" , &n , &m); for(i = 1 ; i <= n ; i ++ ) scanf("%d" , &w[i]); for(i = 1 ; i < n ; i ++ ) scanf("%d%d" , &x , &y) , add(x , y) , add(y , x); dfs(1); while(m -- ) { scanf("%d%d%d" , &opt , &x , &y); if(opt) w[x] = y; else if(judge(x , y)) puts("Y"); else puts("N"); } return 0;}
【bzoj3251】树上三角形 朴素LCA+暴力
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。