|hdu 3966|树链剖分|线段树|Aragorn's Story
来源:互联网 发布:什么是网络接口 编辑:程序博客网 时间:2024/05/29 03:04
hdu 3966
裸树剖+线段树维护,while写成if, 数组开小搞得我调试了好久。。静态查错真的不能快了
#pragma comment(linker, "/STACK:1024000000,1024000000") #include<cstdio>#include<cstring>#include<algorithm>#define fo(i, j, k) for (i=(j);i<=(k);i++)#define fd(i, k, j) for (i=(k);i>=(j);i--)#define fe(i, u) for (i=head[u];i!=-1;i=e[i].next)#define rd(a) scanf("%d", &a)#define rd2(a, b) scanf("%d%d", &a, &b)#define rd3(a, b, c) scanf("%d%d%d", &a, &b, &c)#define ms(i, j) memset(i, j, sizeof i)#define FN2 "hdu3966" using namespace std;const int MAXN = 50000 + 5;struct date{int to, next;}e[MAXN*2];int head[MAXN], cnt;int n, m, Q, wi[MAXN];int siz[MAXN], top[MAXN], p[MAXN], son[MAXN], fa[MAXN], dep[MAXN], pre;void dfs1(int u, int p) { int i; siz[u] = 1, fa[u] = p, dep[u] = dep[p] + 1; fe (i, u) { int v = e[i].to; if (v!=p) { dfs1(v, u); siz[u] += siz[v]; if (son[u]==-1||siz[son[u]]<siz[v]) son[u] = v; } }}void dfs2(int u, int chain) { int i; p[u] = ++pre, top[u] = chain; if (son[u]!=-1) { dfs2(son[u], chain); fe (i, u) { int v = e[i].to; if (v!=fa[u]&&v!=son[u]) dfs2(v, v); } }}int sumv[MAXN*4], addv[MAXN*4];void pushup(int o) { int lc = o*2, rc = o*2+1; sumv[o] = sumv[lc] + sumv[rc];}void pushdown(int o, int len) { int lc = o*2, rc = o*2+1; if (addv[o]) { addv[lc] += addv[o], addv[rc] += addv[o]; sumv[lc] += (len-len/2)*addv[o], sumv[rc] += (len/2)*addv[o]; addv[o] = 0; }}void update(int o, int l, int r, int x, int y, int v) { int lc = o*2, rc = o*2+1, M = (l+r)/2; if (x<=l&&r<=y) { addv[o] += v; sumv[o] += (r-l+1)*v; return ; } pushdown(o, r-l+1); if (x<=M) update(lc, l, M, x, y, v); if (M<y) update(rc, M+1, r, x, y, v); pushup(o);}int query(int o, int l, int r, int p) { int lc = o*2, rc = o*2+1, M = (l+r)/2; if (l==r) return sumv[o]; pushdown(o, r-l+1); if (p<=M) return query(lc, l, M, p); else if (M<p) return query(rc, M+1, r, p);}void add(int u, int v, int c) { int f1 = top[u], f2 = top[v]; while (f1!=f2) { if (dep[f1]<dep[f2]) swap(f1, f2), swap(u, v); update(1,1,n,p[f1],p[u],c); u = fa[f1]; f1 = top[u]; } if (dep[u]<dep[v]) swap(u, v); update(1,1,n,p[v],p[u],c);}void ins(int u, int v) { cnt++; e[cnt].to = v, e[cnt].next = head[u], head[u] = cnt; cnt++; e[cnt].to = u, e[cnt].next = head[v], head[v] = cnt;}void init() { int i; cnt = pre = 0; fo (i, 1, n) head[i] = son[i] = -1, siz[i] = top[i] = p[i] = fa[i] = dep[i] = 0; fo (i, 1, 2*n) e[i].to = 0, e[i].next = -1; fo (i, 1, n) rd(wi[i]); fo (i, 1, n*4) sumv[i] = addv[i] = 0; fo (i, 1, m) { int a, b; rd2(a, b), ins(a, b); } dfs1(1, 0), dfs2(1, 1);}void solve() { int i; fo (i, 1, n) add(i,i,wi[i]); fo (i, 1, Q) { int c1, c2, k; char ch[5]; scanf("%s", ch); if (ch[0]=='I') { rd3(c1, c2, k); add(c1, c2, k); } else if (ch[0]=='D') { rd3(c1, c2, k); add(c1, c2, -k); } else if (ch[0]=='Q') { rd(k); printf("%d\n", query(1,1,n,p[k])); } }}int main() { #ifndef ONLINE_JUDGE freopen(FN2".in","r",stdin);freopen(FN2".out","w",stdout); #endif while (rd3(n, m, Q)==3) init(), solve(); return 0;}
阅读全文
0 0
- 【HDU】3966 Aragorn's Story(树链剖分+线段树)
- HDU 3966 Aragorn's Story(树链剖分+线段树)
- 【HDU 3966】Aragorn's Story(树链剖分+线段树)
- HDU 3966 Aragorn's Story (树链剖分 + 线段树)
- |hdu 3966|树链剖分|线段树|Aragorn's Story
- HDU-3966 Aragorn's Story(树剖+线段树)
- 树链剖分+线段树 hdu3966 Aragorn's Story
- hdu3966 Aragorn's Story【树链剖分+线段树】
- 【HDU3966】Aragorn's Story(树链剖分+线段树)
- hdoj 3966 Aragorn's Story 【树链剖分+线段树||树状数组】
- HDU 3966 Aragorn's Story(树链剖分+线段树区间更新+手动扩大内存)
- HDU 3966 Aragorn's Story 树链剖分+线段树(区间更新 单点查值)
- HDU 3966Aragorn's Story(树链剖分+线段树+入点)
- HDU 3966 Aragorn's Story 树链剖分模板
- hdu 3966 Aragorn's Story(树链剖分)
- hdu 3966 Aragorn's Story (树链剖分)
- hdu 3966 Aragorn's Story 树链剖分
- Hdu 3966 Aragorn's Story (树链剖分)
- Tensorflow 学习笔记 vs code 的用户配置
- JS定时跳转页面
- Android Studio 查看签名文件中的签名
- 插入HTML标签或表格
- 【Java学习笔记】19.登录和注册案例(IO流实现)
- |hdu 3966|树链剖分|线段树|Aragorn's Story
- Codeforces 794C
- 中文技术文档写作规范
- java之多线程
- Lintcode 转换字符串到整数
- BZOJ4849 [Neerc2016]Mole Tunnels
- 蓝桥杯_地址转换
- 8090影子银行
- RegExp(正则表达式)常用知识点小结