首页 > 其他 > 详细

树pao(雾)

时间:2019-02-12 00:30:25      阅读:259      评论:0      收藏:0      [点我收藏+]

今天难得睡醒了,大概睡了13个小时,打算今晚把树剖学了。

嘿嘿嘿雀魂真好玩,这篇文章咕了

有学上了,找到了水平远高于我的队友,好开心的说,,,

卡空间感觉治不了了。。。

板子题是洛谷P3384

实在不知道怎么优化了,感觉应该是哪里写错了,但是是MLE又不是WA,明天再说惹

参考博客:

https://www.cnblogs.com/pks-t/p/9194322.html#_label0

https://www.cnblogs.com/George1994/p/7821357.html

下面那篇要好懂一些,码风也不错。所以为什么还要粘上面的呢?

其实树剖真不难。。就是码量大一些,如果你会 线段树,dfs序,lca。那么应该一看就懂了。

树剖其实就是把一个树剖成很多条链,然后在链上搞来搞去。

先来了解一下概念 

  • 重结点:子树结点数目最多的结点;
  • 轻节点:父亲节点中除了重结点以外的结点;
  • 重边:父亲结点和重结点连成的边;
  • 轻边:父亲节点和轻节点连成的边;
  • 重链:由多条重边连接而成的路径;
  • 轻链:由多条轻边连接而成的路径;

技术分享图片

 

然后有几个结论,不会证真的不会证,·········  记不得了,可以看上面blog

反正就是告诉你我们这样搞不会t。。。

那么我们怎么维护呢

比方说我们要  求 x 到 y 的路径长。

当x与y不在一条重链上时,我们可以 往上跳,跳到 top[x]这个位置,然后累加这一段答案,因为x与top[x]在一条重链所以很简单。

这样跳呀跳呀就跳到一条重链上去了。 然后再算一下就阔以了。

为什么不会T呢?复杂度是什么呢?这些我都不会。

同理 更新操作也是一样的。

如果对子树呢? 显然子树内的序号都是连续的。所以直接对 tid[x],tid[x]+siz[x]-1   这段区间搞来搞去就行了。

484很简单qwq

嗯敲代码的时候可就爽死你了

哦还没有说怎么算top,hson这些东西,可以先dfs一下求siz和hson。再dfs一下求top。具体看代码。

然后粘一个适合我自己码风的板子。vector的改天再说。因为自己平时喜欢vector。。。

要板子的话可以去看上面blog,真的讲得很好,每一步的代码什么的。

 

技术分享图片
  1 #include <bits/stdc++.h>
  2 typedef long long ll;
  3 using namespace std;
  4 inline int read() {
  5     int X=0,w=1; char c=getchar();
  6     while (c<0||c>9) { if (c==-) w=-1; c=getchar(); }
  7     while (c>=0&&c<=9) X=(X<<3)+(X<<1)+c-0,c=getchar();
  8     return X*w;
  9 }
 10 const int N = 1e5+1;
 11 struct Edge{
 12     int to,next;
 13 }edge[N];int cnt=0;int head[N];
 14 void addEdge(int u,int v){
 15     edge[++cnt].to=v;
 16     edge[cnt].next=head[u];
 17     head[u]=cnt;
 18 }
 19 struct Node{
 20     int lz,val;
 21 }tree[N];
 22 int siz[N],top[N],hson[N],dep[N],fa[N],tid[N],ran[N];
 23 int n,m,r,p,val[N],tot=0;
 24 void dfs1(int u,int par,int d){
 25     dep[u]=d;
 26     fa[u]=par;
 27     siz[u]=1;
 28     for(int i=head[u];i;i=edge[i].next){
 29         int v = edge[i].to;
 30         if(v!=fa[u]){
 31             dfs1(v,u,d+1);
 32             siz[u]+=siz[v];
 33             if(hson[u]==-1||siz[v]>siz[hson[u]]){
 34                 hson[u]=v;
 35             }
 36         }
 37     }
 38 }
 39 void dfs2(int u,int t){
 40     top[u]=t;
 41     tid[u]=++tot;
 42     ran[tot]=u;
 43     if(!hson[u])return;
 44     dfs2(hson[u],t);
 45     for(int i=head[u];i;i=edge[i].next){
 46         int v = edge[i].to;
 47         if(v!=hson[u]&&v!=fa[u]){
 48             dfs2(v,v);
 49         }
 50     }
 51 }
 52 void push_up(int x){
 53     tree[x].val=(0ll+tree[x<<1].val+tree[x<<1|1].val)%p;
 54 }
 55 void push_down(int x,int l,int r){
 56     if(!tree[x].lz)return;
 57     int mid=l+r>>1;
 58     tree[x<<1].val=(0ll+tree[x<<1].val+(mid-l+1)*tree[x].lz)%p;
 59     tree[x<<1|1].val=(0ll+tree[x<<1|1].val+(r-mid)*tree[x].lz)%p;
 60     tree[x<<1].lz=(0ll+tree[x<<1].lz+tree[x].lz)%p;
 61     tree[x<<1|1].lz=(0ll+tree[x<<1|1].lz+tree[x].lz)%p;
 62     tree[x].lz=0;
 63 }
 64 void build(int k,int l,int r){
 65     if(l==r){
 66         tree[k].val=val[ran[l]]%p;
 67     } else{
 68         int mid=l+r>>1;
 69         build(k<<1,l,mid);
 70         build(k<<1|1,mid+1,r);
 71         push_up(k);
 72     }
 73 }
 74 int query(int l,int r,int ql,int qr,int rt){
 75     int res = 0;
 76     if(ql<=l&&qr>=r){
 77         res=(0ll+res+tree[rt].val)%p;
 78         return res;
 79     }
 80     int mid = (l+r)>>1;
 81     push_down(rt,l,r);
 82     if(ql<=mid) res=(0ll+res+query(l,mid,ql,qr,rt<<1))%p;
 83     if(qr>mid) res=(0ll+res+query(mid+1,r,ql,qr,rt<<1|1))%p;
 84     return res;
 85 }
 86 void upd(int l,int r,int ul,int ur,int rt,int k){
 87     if(ul<=l&&ur>=r) {
 88         tree[rt].val = (tree[rt].val + k * (r - l + 1))%p;
 89         tree[rt].lz+=k;
 90         return;
 91     }
 92     int mid=l+r>>1;
 93     push_down(rt,l,r);
 94     if(ul<=mid)upd(l,mid,ul,ur,rt<<1,k);
 95     if(ur>mid)upd(mid+1,r,ul,ur,rt<<1|1,k);
 96     push_up(rt);
 97 }
 98 int query_p(int u,int v){//x
 99     int ans=0;
100     while (top[u]!=top[v]){
101         if(dep[top[u]]<dep[top[v]])swap(u,v);
102         ans=(0ll+ans%p+query(1,n,tid[top[u]],tid[u],1))%p;
103         u=fa[top[u]];
104     }
105     if(dep[u]>dep[v])swap(u,v);
106     ans=(0ll+ans+query(1,n,tid[u],tid[v],1))%p;
107     return ans;
108 }
109 void upd_p(int u,int v,int z){
110     z%=p;
111     while (top[u]!=top[v]){
112         if(dep[top[u]]<dep[top[v]])swap(u,v);
113         upd(1,n,tid[top[u]],tid[u],1,z);
114         u=fa[top[u]];
115     }
116     if(dep[u]>dep[v])swap(u,v);
117     upd(1,n,tid[u],tid[v],1,z);
118 }
119 int main(){
120     n=read();m=read();r=read();p=read();
121     for(int i=1;i<=n;i++)val[i]=read();
122     int op,x,y,z;
123     for(int i=1;i<n;i++){
124         x=read();y=read();
125         addEdge(x,y);
126         addEdge(y,x);
127     }
128     dfs1(r,0,1);
129     dfs2(r,r);
130     build(1,1,n);
131     /**while (m--){
132         op=read();
133         if(op==1){
134             x=read();y=read();z=read();
135             upd_p(x,y,z);
136         } else if(op==2){
137             x=read();y=read();
138             printf("%d\n",query_p(x,y));
139         } else if(op==3){
140             x=read();z=read();
141             upd(1,n,tid[x],tid[x]+siz[x]-1,1,z);
142         } else{
143             x=read();
144             printf("%d\n",query(1,n,tid[x],tid[x]+siz[x]-1,1));
145         }
146     }*/
147 }
View Code

 

树pao(雾)

原文:https://www.cnblogs.com/MXang/p/10363669.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!