HYSBZ - 2157树链剖分

【题目描述】
HYSBZ - 2157树链剖分
在这里插入图片描述
【题目分析】
这道题给出的是边权而不是点权，但是我们分析这个树就会发现每个节点都只有一个父亲，也就是每条边的边权都可以存放在儿子节点上，然后在遍历路径的时候我们在从前往后遍历，但是注意最后一条链的链首不要算在内（因为我们只算边权，链首存储的值不在路径上）
接下来就是无穷无尽的代码了，我自己写的一直一直wa，我看了一天了也没有看出哪里有问题，借鉴大佬的代码：
【AC代码】

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>using namespace std;const int maxn = 120000;
#define ll long long
#define ls now<<1
#define rs now<<1|1
inline int read()
{int x = 0, f = 1;char ch = getchar();while (ch < '0' || ch > '9') { if (ch == '-') f = -1; ch = getchar(); }while (ch >= '0' && ch <= '9') { x = x * 10 + ch - 48; ch = getchar(); }return x * f;
}
struct node
{ll f, t, v;
}e[maxn << 1];
struct tree
{ll l, r, sum, maxx, minn, add, mul , num;
}tre[maxn << 2];
ll n, m, tot, cnt;
ll head[maxn], nxt[maxn << 1], used[maxn], dep[maxn], fa[maxn], sz[maxn], w1[maxn];
ll w2[maxn] , son[maxn] , top[maxn] , num[maxn];
inline void buildnode(ll a, ll b , ll c)
{tot++;e[tot].f = a;e[tot].t = b;e[tot].v = c;nxt[tot] = head[a];head[a] = tot;
}
inline void pushup(ll now)
{tre[now].sum = tre[ls].sum + tre[rs].sum;tre[now].maxx = max(tre[ls].maxx, tre[rs].maxx);tre[now].minn = min(tre[ls].minn, tre[rs].minn);
}
inline void build(ll now, ll l, ll r)
{tre[now].l = l;tre[now].r = r;tre[now].num = r - l + 1;if (l == r){tre[now].sum = w2[l];tre[now].maxx = w2[l];tre[now].minn = w2[l];return;}ll mid = (l + r) >> 1;build(ls, l, mid);build(rs, mid + 1 , r);pushup(now);
}
inline void pushdown(ll now)
{if (tre[now].mul){tre[ls].sum *= -1;tre[rs].sum *= -1;ll t = tre[ls].maxx;tre[ls].maxx = -tre[ls].minn;tre[ls].minn = -t;t = tre[rs].maxx;tre[rs].maxx = -tre[rs].minn;tre[rs].minn = -t;}tre[ls].mul = (tre[ls].mul + tre[now].mul) % 2;tre[rs].mul = (tre[rs].mul + tre[now].mul) % 2;tre[now].mul = 0;
}
inline void update1(ll now, ll x ,  ll p)
{if (tre[now].l == tre[now].r){tre[now].sum = p;tre[now].maxx = p;tre[now].minn = p;return;}pushdown(now);ll mid = (tre[now].l + tre[now].r) >> 1;if (x <= mid) update1(ls, x, p);if (x > mid) update1(rs, x, p);pushup(now);
}
inline void update2(ll now, ll l, ll r)
{if (l <= tre[now].l && tre[now].r <= r){tre[now].sum = -tre[now].sum;tre[now].mul = (tre[now].mul + 1) % 2;ll t = tre[now].maxx;tre[now].maxx = -tre[now].minn;tre[now].minn = -t;return;}pushdown(now);ll mid = (tre[now].l + tre[now].r) >> 1;if (l <= mid) update2(ls, l, r);if (r > mid) update2(rs, l, r);pushup(now);
}
inline void dfs1(ll x, ll fat)
{dep[x] = dep[fat] + 1;fa[x] = fat;sz[x] = 1;int k = -1;for (int i = head[x]; i; i = nxt[i]){int u = e[i].t;if (u == fa[x]) continue;w1[u] = e[i].v;dfs1(u, x);sz[x] += sz[u];if (sz[u] > k) k = sz[u], son[x] = u;}
}
inline void dfs2(ll x , ll topx)
{num[x] = ++cnt;w2[cnt] = w1[x];top[x] = topx;if (!son[x]) return;dfs2(son[x], topx);for (int i = head[x]; i; i = nxt[i]){int u = e[i].t;if (u == fa[x] || u == son[x]) continue;dfs2(u, u);}
}
inline int querys(ll now, ll l, ll r)
{if (l <= tre[now].l && tre[now].r <= r)return tre[now].sum;pushdown(now);ll mid = (tre[now].l + tre[now].r) >> 1, ans = 0;if (l <= mid) ans += querys(ls, l, r);if (r > mid) ans += querys(rs, l, r);return ans;
}
inline ll query_max(ll now, ll l , ll r)
{if (l <= tre[now].l &&tre[now].r <= r)return tre[now].maxx;pushdown(now);ll mid = (tre[now].l + tre[now].r) >> 1, ans = -1e9;if (l <= mid) ans = max(ans, query_max(ls, l, r));if (r > mid) ans = max(ans, query_max(rs, l, r));return ans;
}
inline ll query_min(ll now, ll l , ll r)
{if (l <= tre[now].l &&tre[now].r <= r)return tre[now].minn;pushdown(now);ll mid = (tre[now].l + tre[now].r) >> 1, ans = 1e9;if (l <= mid) ans = min(ans, query_min(ls, l, r));if (r > mid) ans = min(ans, query_min(rs, l, r));return ans;
}
inline ll qline(ll x, ll y)
{ll ans = 0;while (top[x] != top[y]){if (dep[top[x]] < dep[top[y]]) swap(x , y);ans += querys(1, num[top[x]], num[x]);x = fa[top[x]];}if (dep[x] > dep[y]) swap(x, y);ans += querys(1, num[x] + 1, num[y]);return ans;
}
inline void cline(ll x, ll y)
{while (top[x] != top[y]){if (dep[top[x]] < dep[top[y]]) swap(x , y);update2(1, num[top[x]], num[x]);x = fa[top[x]];}if (dep[x] > dep[y]) swap(x, y);update2(1, num[x] + 1, num[y]);
}
inline ll qpmax(ll x, ll y)
{ll ans = -1e9;while (top[x] != top[y]){if (dep[top[x]] < dep[top[y]]) swap(x , y);ans = max(ans, query_max(1, num[top[x]], num[x]));x = fa[top[x]];}if (dep[x] > dep[y]) swap(x, y);ans = max(ans, query_max(1, num[x] + 1, num[y]));return ans;
}
inline ll qpmin(ll x, ll y)
{ll ans = 1e9;while (top[x] != top[y]){if (dep[top[x]] < dep[top[y]]) swap(x , y);ans = min(ans, query_min(1, num[top[x]], num[x]));x = fa[top[x]];}if (dep[x] > dep[y]) swap(x, y);ans = min(ans, query_min(1, num[x] + 1, num[y]));return ans;
}
string s1 = "C", s2 = "N", s3 = "SUM", s4 = "MAX", s5 = "MIN";
int main()
{//freopen("1.in","r",stdin);n = read();for (int i = 1; i <= n - 1; i++){int a, b, c;a = read(); b = read(); c = read();buildnode(a + 1, b + 1, c);buildnode(b + 1, a + 1, c);}dfs1(1 , 0);dfs2(1 , 1);build(1 , 1 , n);m = read();while (m--){string s;int a, b, k;cin >> s >> a >> b;a += 1 , b += 1;if (s == s1){a -= 1 , b -= 1;if (dep[e[a * 2 - 1].f] > dep[e[a * 2 - 1].t])k = e[a * 2 - 1].f;else k = e[a * 2 - 1].t;update1(1, num[k], b);}if (s == s2) cline(a, b);if (s == s3) printf("%lld\n", qline(a, b));if (s == s4) printf("%lld\n", qpmax(a, b));if (s == s5) printf("%lld\n", qpmin(a, b));}return 0;
}

我自己写的也贴上，方便日后检查

#include<cstdio>
#include<cstring>
#include<cmath>
#include<climits>
#include<cstdlib>
#include<algorithm>
#include<queue>
#include<vector>
#include<set>using namespace std;typedef long long ll;const int MAXN=120005;
int fa[MAXN],A[MAXN],val[MAXN],pos[MAXN];
int siz[MAXN],son[MAXN],h[MAXN],top[MAXN];
int cnt=0,n,m;
int Sum[MAXN<<2],Max[MAXN<<2],Min[MAXN<<2],lazy[MAXN<<2];
struct node
{int u,v,w;
}e[MAXN<<1];
int head[MAXN<<1],nxt[MAXN<<1];
int tot=0;void AddEdge(int u,int v,int w)
{tot++;e[tot].u=u; e[tot].v=v; e[tot].w=w;nxt[tot]=head[u]; head[u]=tot;
}void dfs1(int u,int f)
{int i,v;siz[u]=1;son[u]=0;fa[u]=f;h[u]=h[f]+1;for(i=head[u];i;i=nxt[i]){v=e[i].v;if(v!=f){val[v]=e[i].w;dfs1(v,u);siz[u]+=siz[v];if(siz[son[u]]<siz[v]) son[u]=v;}}
}
void dfs2(int u,int f,int k)
{int i,v;top[u]=k;pos[u]=++cnt;A[cnt]=val[u];if(!son[u]) return; if(son[u]) dfs2(son[u],u,k);for(i=head[u];i;i=nxt[i]){v=e[i].v;if(v!=f&&v!=son[u]) dfs2(v,u,v);}
}void pushup(int k)
{Sum[k]=Sum[k<<1]+Sum[k<<1|1];Max[k]=max(Max[k<<1],Max[k<<1|1]);Min[k]=min(Min[k<<1],Min[k<<1|1]);
}void pushdown(int k)
{if(lazy[k]){Sum[k<<1]*=-1; Sum[k<<1|1]*=-1;swap(Max[k<<1],Min[k<<1]);Max[k<<1]*=-1; Min[k<<1]*=-1;swap(Max[k<<1|1],Min[k<<1|1]);Max[k<<1|1]*=-1; Min[k<<1|1]*=-1;lazy[k<<1]^=1; lazy[k<<1|1]^=1;lazy[k]=0;}
}void build(int k,int l,int r)
{if(l==r){Sum[k]=Max[k]=Min[k]=A[l];return;}int mid=(l+r)>>1;build(k<<1,l,mid);build(k<<1|1,mid+1,r);pushup(k);
}void PointChange(int k,int l,int r,int x,int v)
{if(l==r && l==x){Sum[k]=Max[k]=Min[k]=v;return;}pushdown(k);int mid=(l+r)>>1;if(x<=mid) PointChange(k<<1,l,mid,x,v);else PointChange(k<<1|1,mid+1,r,x,v);pushup(k);  
}void IntervalChange(int k,int l,int r,int L,int R)
{if(l>=L && r<=R){Sum[k]*=-1; swap(Max[k],Min[k]);Max[k]*=-1; Min[k]*=-1;lazy[k]^=1;return;}pushdown(k);int mid=(l+r)>>1;if(L<=mid) IntervalChange(k<<1,l,mid,L,R);if(R>mid) IntervalChange(k<<1|1,mid+1,r,L,R);pushup(k); 
}int IntervalSum(int k,int l,int r,int L,int R)
{if(L<=l && r<=R){return Sum[k];}int mid=(l+r)/2;pushdown(k);int ret=0;if(L<=mid) ret+=IntervalSum(k<<1,l,mid,L,R);if(R>mid) ret+=IntervalSum(k<<1|1,mid+1,r,L,R);return ret;
}int IntervalMax(int k,int l,int r,int L,int R)
{if(L<=l && r<=R){return Max[k];}int mid=(l+r)/2;pushdown(k);int ret=INT_MIN;if(L<=mid) ret=max(IntervalSum(k<<1,l,mid,L,R),ret);if(R>mid) ret=max(IntervalSum(k<<1|1,mid+1,r,L,R),ret);return ret;
}int IntervalMin(int k,int l,int r,int L,int R)
{if(L<=l && r<=R){return Min[k];}int mid=(l+r)/2;pushdown(k);int ret=INT_MAX;if(L<=mid) ret=min(IntervalSum(k<<1,l,mid,L,R),ret);if(R>mid) ret=min(IntervalSum(k<<1|1,mid+1,r,L,R),ret);return ret;
}int FindSum(int u,int v)
{int ans=0;while(top[u]!=top[v]){if(h[top[u]]<h[top[v]]) swap(u,v);ans+=IntervalSum(1,1,n,pos[top[u]],pos[u]);u=fa[top[u]];}if(h[u]<h[v]) swap(u,v);ans+=IntervalSum(1,1,n,pos[v]+1,pos[u]);return ans;
}int FindMax(int u,int v)
{int ans=INT_MIN;while(top[u]!=top[v]){if(h[top[u]]<h[top[v]]) swap(u,v);ans=max(IntervalMax(1,1,n,pos[top[u]],pos[u]),ans);u=fa[top[u]];}if(h[u]<h[v]) swap(u,v);ans=max(IntervalMax(1,1,n,pos[v]+1,pos[u]),ans);return ans;
}int FindMin(int u,int v)
{int ans=INT_MAX;while(top[u]!=top[v]){if(h[top[u]]<h[top[v]]) swap(u,v);ans=min(IntervalMin(1,1,n,pos[top[u]],pos[u]),ans);u=fa[top[u]];}if(h[u]<h[v]) swap(u,v);ans=min(IntervalMin(1,1,n,pos[v]+1,pos[u]),ans);return ans;
}void update(int u,int v)
{while(top[u]!=top[v]){if(h[top[u]]<h[top[v]]) swap(u,v);IntervalChange(1,1,n,pos[top[u]],pos[u]);u=fa[top[u]];}if(h[u]<h[v]) swap(u,v);IntervalChange(1,1,n,pos[v]+1,pos[u]);
}int main()
{int u,v,w,idx;char cmd[10];scanf("%d",&n);//for(int i=1;i<=n;i++) scanf("%d",&val[i]);for(int i=1;i<n;i++){scanf("%d%d%d",&u,&v,&w);AddEdge(u+1,v+1,w); AddEdge(v+1,u+1,w);}dfs1(1 , 0);dfs2(1 ,0, 1);build(1,1,n);scanf("%d",&m);while(m--){scanf("%s",cmd);if(cmd[0]=='C'){scanf("%d%d",&idx,&w);u=e[idx*2-1].u; v=e[idx*2-1].v;if(h[u]<h[v]){PointChange(1,1,n,pos[v],w);}else{PointChange(1,1,n,pos[u],w);}}else if(cmd[0]=='N') {scanf("%d%d",&u,&v);update(u+1,v+1);}else if(cmd[0]=='S'){scanf("%d%d",&u,&v);printf("%d\n",FindSum(u+1,v+1));}else if(cmd[1]=='A'){scanf("%d%d",&u,&v);printf("%d\n",FindMax(u+1,v+1));}else if(cmd[1]=='I'){scanf("%d%d",&u,&v);printf("%d\n",FindMin(u+1,v+1));}}return 0;
}

本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：http://www.mzph.cn/news/384112.shtml

如若内容造成侵权/违法违规/事实不符，请联系多彩编程网进行投诉反馈email:809451989@qq.com，一经查实，立即删除！

相关文章

C++11中的右值引用

C++11中的右值引用

http://www.cnblogs.com/yanqi0124/p/4723698.html 在C98中有左值和右值的概念&#xff0c;不过这两个概念对于很多程序员并不关心&#xff0c;因为不知道这两个概念照样可以写出好程序。在C11中对右值的概念进行了增强&#xff0c;我个人理解这部分内容是C11引入的特性中最难以…
阅读更多...
BZOJ2115XOR——线性基

BZOJ2115XOR——线性基

【题目描述】 BZOJ2115XOR——线性基 【题目分析】 这道题看完以后很懵逼&#xff0c;人家要是走的很复杂呢&#xff1f;各种绕来绕去怎么办&#xff1f; 首先我们应该注意到一个很明显的道理&#xff1a;重复的路径会和自身抵消&#xff0c;所以我们大可以随便跑&#xff0c;…
阅读更多...
单链表的相关操作

单链表的相关操作

1.冒泡排序对单链表进行排序 void LinkListBubbleSort(LinkNode* head) {if(head NULL){ return;//空链表} if(head -> next NULL){ return;//只有一个结点} LinkNode* cur head;//趟数LinkNode* tail NULL;//尾指针LinkNode* tmp head;//次数for(; cur -…
阅读更多...
socket网络编程--epoll小结

socket网络编程--epoll小结

http://www.cnblogs.com/wunaozai/p/3895860.html 以前使用的用于I/O多路复用为了方便就使用select函数&#xff0c;但select这个函数是有缺陷的。因为它所支持的并发连接数是有限的(一般小于1024)&#xff0c;因为用户处理的数组是使用硬编码的。这个最大值为FD_SETSIZE&#…
阅读更多...
进程间通信(匿名管道)

进程间通信(匿名管道)

1.进程通信的目的 (1) 数据传输: 一个进程需要将它的数据传输给另一个进程     (2) 资源共享: 多个进程之间共享同样的资源     (3) 通知事件: 一个进程需要向另一个或一组进程发送消息, 通知它们发生了什么事情 2.管道 管道是一种进程之间通信的一种方式, 我们把从…
阅读更多...
线性基入门

线性基入门

今天学习了神奇的线性基&#xff0c;主要是在解决异或问题时比较有用。 详细的解释和证明有大佬珠玉在前&#xff0c;如果感兴趣可以移步 补充一下自己的理解&#xff1a; 可以联系线性代数极大无关组进行理解&#xff0c;线性基就相当于异或的向量空间中的极大无关组&#xff…
阅读更多...
单例模式及C++实现代码

单例模式及C++实现代码

http://www.cnblogs.com/cxjchen/p/3148582.html 单例模式 单例模式&#xff0c;可以说设计模式中最常应用的一种模式了&#xff0c;据说也是面试官最喜欢的题目。但是如果没有学过设计模式的人&#xff0c;可能不会想到要去应用单例模式&#xff0c;面对单例模式适用的情况&am…
阅读更多...
UVALive - 8512——线段树维护线性基

UVALive - 8512——线段树维护线性基

【题目描述】 UVALive - 8512XOR 【题目分析】 这种区间线性基的问题我们可以考虑用线段树维护&#xff0c;线性基的合并的话就直接暴力合并 找到所在区间的线性基后再查找最大的数&#xff0c;我看网上的博客要说消除k的影响什么的&#xff0c;我觉得没有什么必要&#xff0c;…
阅读更多...
命名管道

命名管道

1.命名管道的创建 (1) 通过命令创建 mkfifo filename (2)在程序中创建 int mkfifo(const char* filename, mode_t mode); 2. 命名管道和匿名管道的区别 (1)匿名管道由pipe函数创建并且打开     (2)命名管道有mkfifo函数创建由open函数打开     (3) fifo 之间的两…
阅读更多...
HYSBZ - 1101——莫比乌斯反演

HYSBZ - 1101——莫比乌斯反演

【题目描述】 HYSBZ - 1101 【题目分析】 昨天测试出了一道差不多的题目&#xff0c;我只能想到暴力&#xff0c;各种优化&#xff0c;最后都是运行了好久TLE&#xff0c;最后才知道要用到莫比乌斯反演&#xff0c;就想着今天研究一下&#xff0c;得出的结论就是&#xff0c;我…
阅读更多...
Linux下I/O多路转接之select --fd_set

Linux下I/O多路转接之select --fd_set

http://blog.csdn.net/li_ning_/article/details/52165993 fd_set 你终于还是来了&#xff0c;能看到这个标题进来的&#xff0c;我想&#xff0c;你一定是和我遇到了一样的问题&#xff0c;一样的疑惑&#xff0c;接下来几个小时&#xff0c;我一定竭尽全力&#xff0c;写出我…
阅读更多...
BZOJ 2844 | HYSBZ - 2844albus就是要第一个出场——线性基

BZOJ 2844 | HYSBZ - 2844albus就是要第一个出场——线性基

【题目描述】 BZOJ 2844 | HYSBZ - 2844albus 【题目分析】 题目的意思大概是给一个数列&#xff0c;他有2n个子集&#xff0c;每个子集的元素的异或和构成新的一个数列&#xff0c;排序后问数字Q在这个序列里面的下标。 假如题目是求所有元素的异或和构成一个集合就好弄了&…
阅读更多...
CodeForces - 641ELittle Artem and Time Machine——map+树状数组

CodeForces - 641ELittle Artem and Time Machine——map+树状数组

【题目描述】 CodeForces - 641ELittle Artem and Time Machine 【题目分析】 题目的意思大概是有三种操作 1.在时间t加入一个数字x 2.在时间t删除一个数字x 3.询问在时间t集合里面x的个数 虽然题目描述很简单&#xff0c;但是t和x的范围都是109&#xff0c;我一开始想到的是主…
阅读更多...
I/O多路转接之poll 函数

I/O多路转接之poll 函数

http://blog.csdn.net/li_ning_/article/details/52167224 poll 一、poll()函数&#xff1a; 这个函数是某些Unix系统提供的用于执行与select()函数同等功能的函数&#xff0c;自认为poll和select大同小异&#xff0c;下面是这个函数的声明&#xff1a; [cpp] view plaincopy …
阅读更多...
链表相关笔试面试题

链表相关笔试面试题

1.判断两个链表是否相交 两个链表是否相交可分为以下几种情况     &#xff08;1&#xff09;两个链表都不带环&#xff0c;此时两个链表所对应的最后一个节点是相等的     &#xff08;2&#xff09;两个链表一个带环&#xff0c;一个不带环&#xff0c;两个链表一定…
阅读更多...
Linux经典问题—五哲学家就餐问题

Linux经典问题—五哲学家就餐问题

http://m.blog.csdn.net/aspenstars/article/details/70149038 一、问题介绍 由Dijkstra提出并解决的哲学家进餐问题(The Dinning Philosophers Problem)是典型的同步问题。该问题是描述有五个哲学家共用一张圆桌&#xff0c;分别坐在周围的五张椅子上&#xff0c;在圆桌上有五…
阅读更多...
修改之前的myshell使之支持输入输出重定向

修改之前的myshell使之支持输入输出重定向

1.open函数     这个函数是打开一个文件&#xff08;文件名叫pathname),以 flag 权限打开&#xff0c;flag 包括了以下几种 O_RDONLY&#xff08;只读&#xff09;, O_WRONLY&#xff08;只写&#xff09;, O_RDWR&#xff08;读写&#xff09;&#xff0c;当文件打开成…
阅读更多...
HDU - 6621 K-th Closest Distance——主席树+二分

HDU - 6621 K-th Closest Distance——主席树+二分

【题目描述】 HDU - 6621 K-th Closest Distance 【题目分析】 因为看到第kkk大的要求&#xff0c;刚开始的时候一直都在想怎么运用第kkk大来解决问题&#xff0c;但是后来看其他人的博客才发现并不需要用第k大&#xff0c;但是主席树维护权值线段树还是需要的&#xff0c;这…
阅读更多...
链表相关的算法题大汇总 — 数据结构之链表奇思妙想

链表相关的算法题大汇总 — 数据结构之链表奇思妙想

http://blog.csdn.net/lanxuezaipiao/article/details/22100021基本函数&#xff08;具体代码实现见后面&#xff09; 1&#xff0c;构造节点 //定义节点类型 struct Node { int value; Node*next; }; 2&#xff0c;分配节点 //之所以要分配节点原因是需要在分配函数中…
阅读更多...
CodeForces - 372CWatching Fireworks is Fun+DP+单调队列优化

CodeForces - 372CWatching Fireworks is Fun+DP+单调队列优化

【题目描述】 CodeForces - 372CWatching Fireworks is Fun 题目的大概意思就是在一个编号为1…n的街道上现在按照时间顺序放烟花&#xff0c;每个烟花获得的幸福感为b−abs(a−x)b-abs(a-x)b−abs(a−x)&#xff0c;x为观看烟花的位置&#xff0c;为了提升我们的幸福感&#x…
阅读更多...
最新文章

Copyright @ 2022~2023