bzoj4765: 普通计算姬 (分块 && BIT)
最近一直在刷分块啊
似乎感觉分块和BIT是超级棒的搭档啊
这道题首先用dfs预处理一下
得到每一个sum值
此时查询是O(1)的 (前缀和乱搞什么的
但是修改需要O(n) (需要修改该节点所有祖先的sum
复杂度就爆了呀
此时考虑分块优化
似乎弹飞绵羊也是这样思考得出分块做法的
首先分成 √n 块
sum[i]记录第i块的sum和
中间的块直接用sum数组处理 两边用树状数组暴力求
这样查询就是O(√n)的 (其实有一些常数的... 就当是 √n 好了)
修改的话在dfs时用f[i][j]表示第j个点对于第i块的贡献 (需要算几次什么的
然后直接修改块就好了 复杂度也是O (√n) 的
下面是代码
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <vector>
using namespace std;
#define isdigit(x) (x >= '0' && x <= '9')
#define lowbit(x) (x & (-x))
typedef unsigned long long ll;
const int N = 1e5 + ;
const int M = ; int n, root, cnt, sz, tot;
int d[N], b[N], f[M][N], ct[N], L[N], p[N];
ll s[N], c[N];
vector < int > E[N]; inline void read(int &ans) {
ans = ;
register int res = ;
static char buf = getchar();
for (; !isdigit(buf); buf = getchar())
if (buf == '-') res = -;
for (; isdigit(buf); buf = getchar())
ans = ans * + buf - '';
ans *= res;
} inline void addEdge(int u ,int v) {
E[u].push_back(v);
E[v].push_back(u);
} inline void add(int x, ll v) {
while (x <= n) {
c[x] += v;
x += lowbit(x);
}
} inline ll query(int x) {
ll ans = ;
while (x > ) {
ans += c[x];
x -= lowbit(x);
}
return ans;
} ll dfs(int x, int fa) {
ll sum = d[x]; p[x] = ++tot;
ct[b[x]]++; add(tot, d[x]);
for (int i = ; i <= cnt; i++) f[i][x] += ct[i];
for (int i = ; i < E[x].size(); i++) {
int u = E[x][i];
if (u == fa) continue;
sum += dfs(u, x);
}
ct[b[x]]--; L[x] = tot;
s[b[x]] += sum;
return sum;
} inline void modify(int u, int v) {
add(p[u], v - d[u]);
for (int i = ; i <= cnt; i++)
s[i] += (v - d[u]) * 1ll * f[i][u];
d[u] = v;
} inline ll query(int l ,int r) {
ll ans = ;
if (b[l] == b[r]) {
for (int i = l; i <= r; i++)
ans += query(L[i]) - query(p[i] - );
return ans;
}
for (int i = b[l] + ; i < b[r]; i++)
ans += s[i];
for (int i = l; i <= b[l] * sz; i++)
ans += query(L[i]) - query(p[i] - );
for (int i = (b[r] - ) * sz + ; i <= r; i++)
ans += query(L[i]) - query(p[i] - );
return ans;
} int main() {
int m;
read(n); read(m);
sz = sqrt(n);
for (int i = ; i <= n; i++) {
read(d[i]);
b[i] = (i - ) / sz + ;
}
cnt = b[n];
for (int i = ; i <= n; i++) {
int u, v;
read(u); read(v);
if (!u) root = v;
else addEdge(u, v);
}
dfs(root, );
while (m--) {
int op, u, v;
read(op); read(u); read(v);
if (op == )
modify(u, v);
else
printf("%llu\n", query(u, v));
}
return ;
}
bzoj4765: 普通计算姬 (分块 && BIT)的更多相关文章
- [BZOJ4765]普通计算姬(分块+树状数组)
4765: 普通计算姬 Time Limit: 30 Sec Memory Limit: 256 MBSubmit: 1725 Solved: 376[Submit][Status][Discus ...
- [bzoj4765]普通计算姬——分块
Brief Description 给定一棵n个节点的带权树,节点编号为1到n,以root为根,设sum[p]表示以点p为根的这棵子树中所有节点的权 值和.支持下列两种操作: 1 给定两个整数u,v, ...
- BZOJ4765: 普通计算姬
BZOJ4765: 普通计算姬 题目描述 传送门 题目分析 求的和非常奇怪,不具有连续性,所有上树的数据结构全死了. 考虑分块,思考对于一段连续的询问区间可以直接询问整块,零散块可以在树上dfs序暴力 ...
- 2018.06.30 BZOJ4765: 普通计算姬(dfs序+分块+树状数组)
4765: 普通计算姬 Time Limit: 30 Sec Memory Limit: 256 MB Description "奋战三星期,造台计算机".小G响应号召,花了三小时 ...
- [bzoj4765]普通计算姬(分块+树状数组+DFS序)
题意 给定一棵n个节点的带权树,节点编号为1到n,以root为根,设sum[p]表示以点p为根的这棵子树中所有节点的权值和.计算姬支持下列两种操作: 1 给定两个整数u,v,修改点u的权值为v. 2 ...
- BZOJ 4765 普通计算姬 (分块 + BIT)
4765: 普通计算姬 Time Limit: 30 Sec Memory Limit: 256 MBSubmit: 1547 Solved: 329[Submit][Status][Discus ...
- BZOJ4765 普通计算姬(分块+树状数组)
对节点按编号分块.设f[i][j]为修改j号点对第i块的影响,计算f[i][]时dfs一遍即可.记录每一整块的sum.修改时对每一块直接更新sum,同时用dfs序上的树状数组维护子树和.查询时累加整块 ...
- BZOJ 4765: 普通计算姬 [分块 树状数组 DFS序]
传送门 题意: 一棵树,支持单点修改和询问以$[l,r]$为根的子树的权值和的和 只有我这种不会分块的沙茶不会做这道题吗? 说一点总结: 子树和当然上$dfs$序了,询问原序列一段区间所有子树和,对原 ...
- BZOJ 4765: 普通计算姬 (分块+树状数组)
传送门 解题思路 树上的分块题,,对于修改操作,每次修改只会对他父亲到根这条链上的元素有影响:对于查询操作,每次查询[l,r]内所有元素的子树,所以就考虑dfn序,进标记一次,出标记一次,然后子树就是 ...
随机推荐
- LTC
LTC 即 L2C,Leads To Cash,从线索到现金的企业运营管理思想,是以企业的营销和研发两大运营核心为主线,贯穿企业运营全部流程,深度融合了移动互联.SaaS技术.大数据与企业运营智慧,旨 ...
- CF1227F2 Wrong Answer on test 233 (Hard Version)
题意 \(n\)道题,每道题有\(k\)种选项,其中第\(i\)道题正确答案是\(a_i\),但是填答案的时候填错啦,第一道题的选择填到了第二道题...第\(n\)道题的选择填到了第一道题,求在\(k ...
- tsocks代理git wget
使用clash时, 命令行的wget和git操作可能没有被代理 安装tsocks: apt-get install tsocks 修改配置文件: vi /etc/tsocks.conf 找到: ser ...
- BZOJ3473&&BZOJ3277串
BZOJ3473&&BZOJ3277串 题面 自己找去 HINT 对于所有串建立一个广义后缀自动机,对于每一个节点开一个set表示这个节点接受的子串在哪些串里出现过,然后在parent ...
- python3练习100题——049
题目:使用lambda来创建匿名函数. sum=lambda x,y:x+y from functools import reduce reduce(sum,[1,2,3,4,5])
- 回文串--Manacher算法(模板)
用途:在O(n)时间内,求出以每一个点为中心的回文串长度. 首先,有一个非常巧妙的转化.由于回文串长度有可能为奇数也有可能为偶数,说明回文中心不一定在一个字符上.所以要将字符串做如下处理:在每两个字母 ...
- .Net Core 智能提示汉化包
在.Net Core 2.x 版本,Microsoft 官方没有提供 .Net Core 正式版的多语言安装包.因此,我们在用.Net Core 2.x 版本作为框架目标编写代码时,智能提成是英文的. ...
- fastadmin选择下拉框
fastadmin中要做下拉框的效果如下: 数据库中数据: 在对应model中添加一个方法: 控制器中添加一行: 在目录lang/zh-cn中找到你控制器名称所对应的文件添加配置: 在add.html ...
- 【巨杉数据库SequoiaDB】SequoiaDB 巨杉数据库 v3.4 版本正式发布
深秋时节,SequoiaDB 巨杉数据库在深秋给大家带来了“一把火”.SequoiaDB v3.4 正式发布啦! 分布式交易场景性能大幅提升 SequoiaDB 巨杉数据库3.4版本正式发布,v3.4 ...
- jupyter CSS
.caret, .dropup > .btn > .caret { border-top-color: #000 !important; } .label { border: 1px so ...