根号均摊 ---- E. Xenia and Tree(树形dp + 暴力根号均摊)

news/2024/7/5 2:59:11

题目大意


题目大意:

你有一颗树:树开始时候1号点是红色的
现在就是你有两次操作:
1 u 把u点涂成红色
2 u 询问离u点最近红点在哪里?


  1. 数据范围是n,m∈[1,1e5]n,m\in[1,1e5]n,m[1,1e5]时间限制是5s5s5s
  2. 如果直接查询我们是可以直接树形dpO(n)O(n)O(n)去查询,跟新是O(1)O(1)O(1)
  3. 如果直接去维护每个点的答案是跟新就是O(n)O(n)O(n)查询是O(1)O(1)O(1)的?
  4. 那么我们可以考虑去均摊两个复杂度O(nm)O(n\sqrt m)O(nm)
  5. 就是我们可以这样考虑,我们每次先把要修改的点先存起来不修改,等修改个数超过m\sqrt mm的时候我们再去修改
  6. 首先我们对于询问我们要去查询这m\sqrt mm个点查询两点距离要LCA→log\rightarrow loglog那么复杂度是(m−m)∗log(n)(m-\sqrt m)*log(n)(mm)log(n)
  7. 对于跟新我们直接跑一遍树形dpO(n∗m)O(n*\sqrt m)O(nm)最多更新m\sqrt mm
  8. 总的复杂度是O(max((m−m)∗log(n),nm))O(max((m-\sqrt m)*log(n),n\sqrt m))O(max((mm)log(n),nm))
  9. cf机子很快cf机子很快cf

AC code

#include <bits/stdc++.h>
#define mid ((l + r) >> 1)
#define Lson rt << 1, l , mid
#define Rson rt << 1|1, mid + 1, r
#define ms(a,al) memset(a,al,sizeof(a))
#define log2(a) log(a)/log(2)
#define lowbit(x) ((-x) & x)
#define IOS std::ios::sync_with_stdio(0); cin.tie(0); cout.tie(0)
#define INF 0x3f3f3f3f
#define LLF 0x3f3f3f3f3f3f3f3f
#define f first
#define s second
#define endl '\n'
using namespace std;
const int N = 2e6 + 10, mod = 1e9 + 9;
const int maxn = 500010;
const long double eps = 1e-5;
const int EPS = 500 * 500;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> PII;
typedef pair<ll,ll> PLL;
typedef pair<double,double> PDD;
template<typename T> void read(T &x) {x = 0;char ch = getchar();ll f = 1;while(!isdigit(ch)){if(ch == '-')f*=-1;ch=getchar();}while(isdigit(ch)){x = x*10+ch-48;ch=getchar();}x*=f;
}
template<typename T, typename... Args> void read(T &first, Args& ... args) {read(first);read(args...);
}
int n, m;
vector<int> G[maxn];
int dp[maxn];
bool is[maxn];
int depth[maxn], fa[maxn][20];
inline void init(int u, int f) {depth[u] = depth[f] + 1;fa[u][0] = f;for(int i = 1; i < 19; ++ i) fa[u][i] = fa[fa[u][i-1]][i-1];for(auto it : G[u]) {if(it == f) continue;init(it,u);}
}inline int LCA(int u, int v) {if(depth[u] > depth[v]) swap(u,v);int delta = depth[v] - depth[u];for(int i = 19; i >= 0; -- i) if(delta >> i & 1) v = fa[v][i];if(u == v) return v;for(int i = 19; i >= 0; -- i) if(fa[u][i] != fa[v][i])u = fa[u][i], v = fa[v][i];return fa[u][0];
}inline void dfs(int u, int fa) {if(is[u]) dp[u] = 0;for(auto it : G[u]) {if(it == fa) continue;dfs(it,u);dp[u] = min(dp[it]+1,dp[u]);}
}inline void dfs1(int u, int fa) {dp[u] = min(dp[fa]+1,dp[u]);for(auto it : G[u]) {if(it == fa) continue;dfs1(it,u);}
}
vector<int> stk;
int main() {IOS;ms(dp,INF);is[1] = 1;stk.push_back(1);cin >> n >> m;for(int i = 2; i <= n; ++ i) {int u, v;cin >> u >> v;G[u].push_back(v);G[v].push_back(u);}init(1,0);int lim = sqrt(m);while(m --) {int op, u;cin >> op >> u;if(op == 1) stk.push_back(u);else {int ans = dp[u];for(auto it : stk) ans = min(ans,depth[u]+depth[it]-2*depth[LCA(u,it)]);cout << ans << endl;}if(stk.size() > lim) { // 单修改超过lim时候直接跟新for(auto it : stk) is[it] = 1;dfs(1,0);dfs1(1,0);stk.clear();}}return 0;
}

http://lihuaxi.xjx100.cn/news/286001.html

相关文章

转置卷积实现

import tensorflow as tf import numpy as np import tensorflow as keras from tensorflow.keras import layers,optimizers,losses #创建X矩阵&#xff0c;高宽各位5*5 xtf.range(16)1 #调整为合法的维度张量 xtf.reshape(x,[1,4,4,1]) xtf.cast(x,dtypetf.float32) #创建固定…

全家都是博士是一种什么样的体验?

那些在北大未名湖畔或者清华园里&#xff0c;寒窗苦读获得博士学位的学子&#xff0c;都是家长们眼中“别人的孩子”。身边的朋友以及街坊邻居听闻博士的名号&#xff0c;都会投来敬佩而崇拜的目光。虽然如今随着学历水涨船高&#xff0c;博士学位存在一定程度的贬值&#xff0…

来自程序员的福利!用Python做一款翻译软件

来源 | Ahab杂货铺&#xff08;ID&#xff1a;PythonLearningCamp&#xff09;前两天吃了平哥的一波狗粮&#xff0c;他给女朋友写了一个翻译软件&#xff0c;自己真真切切的感受到了程序员的浪漫。在学习requests请求的时候做过类似的Demo&#xff0c;给百度翻译发送一个post请…

记忆优化搜索(简单题)(洛谷P3183 [HAOI2016]食物链 )( P5635 【CSGRound1】天下第一 )

昨天做了蓝桥杯的时候&#xff0c;发现自己对于记忆优化搜索甚是不熟悉&#xff0c;所以今天随便找了几个基础题做做&#xff0c;顺便写下两片题解&#xff0c;顺便用了一下devc敲的代码&#xff0c;发现没有代码补全真的可以说是灰常难受了。。。 洛谷P3183 [HAOI2016]食物链 …

Springboot毕业设计毕设作品,纯净水商城配送系统 开题报告

本科生毕业论文 基于Java纯净水商城配送系统&#xff08;springboot框架&#xff09; 开题报告 学 院&#xff1a; 专 业&#xff1a; 计算机科学与技术 年 级&#xff1a; 学生姓名&#xff1a; …

iptables踩坑记

1&#xff1a;第一坑&#xff1a;众所周知nf_conntrack&#xff0c;下面会有介绍补坑方法。2&#xff1a;连环坑&#xff1a;要解决第一个坑&#xff0c;需要修改内核参数&#xff0c;如&#xff1a;net.netfilter.nf_conntrack_tcp_timeout_established 600 net.netfilter…

编程25年后,现实将我打回菜鸟程序员的起点

在从事了 25 年的编程工作后&#xff0c;我却发现自己没有什么具体的东西可以展示给大家。我需要克服完美主义思想和对被批评的恐惧&#xff0c;向大家展示真实的自己。作者 | DeChamp译者 | 苏本如责编 | 仲培艺出品 | CSDN&#xff08;ID&#xff1a;CSDNNews&#xff09;下面…

tensorflow中转置卷积运算例子

import tensorflow as tf import numpy as np import tensorflow as keras from tensorflow.keras import layers,optimizers,losses #创建X矩阵&#xff0c;高宽各位5*5 xtf.range(25)1 #调整为合法的维度张量 xtf.reshape(x,[1,5,5,1]) xtf.cast(x,dtypetf.float32) #创建固定…