青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

算法學社
記錄難忘的征途
posts - 141,comments - 220,trackbacks - 0

題目描述

   在一個點數為N(N<10,000)的帶權樹上,支持兩個操作:1. 改變一個邊權 2. 詢問u和v之間的路徑上的最大邊權

吐槽:

    1. 這就是樹鏈剖分最基本的應用了吧,寫得我要死了....
    2. 常數好大... 2.9s

算法分析:

    就是樹的輕重邊剖分了。輕邊直接修改,重邊因為是一段一段連續的,所以用一顆線段樹維護起來....
    具體怎么實現的明天再說....
    查詢兩個點u和v就是求max(path(u,lca(u,v)),path(v,lca(u,v)))。
    lca用的是clj版的....
    查詢就是重邊用RMQ,輕邊直接查詢。
    明早起來再詳細講....實驗室要關門了...
  1 #include<iostream>
  2 #include<cstdio>
  3 #include<cstring>
  4 #include<cassert>
  5 using namespace std;
  6 template <typename T> inline void chkmax(T &a, T b){ if(a < b) a = b;}
  7 // build tree
  8 const int V = 10005;
  9 const int E = V*2;
 10 int head[V], nxt[E], pnt[E], e, n, cost[E];
 11 void add_edge(int u,int v,int c){
 12     nxt[e] = head[u];
 13     head[u] = e;
 14     pnt[e] = v;
 15     cost[e] = c;
 16     e++;
 17 }
 18 // segment tree & prepare
 19 const int inf = ~0u>>2;
 20 int seg[E<<2], parent[V], M, UP[V], heavy[V], num[E], deep[V], P[V], sz[V], segsz;
 21 void ins(int pos,int val){
 22     pos += M; 
 23     seg[pos] = val;
 24     while(pos>>=1) {
 25         seg[pos] = max(seg[pos<<1], seg[pos<<1|1]);
 26     }
 27 }
 28 int find(int u){ return u==parent[u] ? u: parent[u] = find(parent[u]);}
 29 int rmq(int l,int r){
 30     int ans = -inf;
 31     for(l += M-1, r += M+1; l^r^1; l>>=1, r>>=1){
 32         if(~l&1) chkmax(ans,seg[l^1]);
 33         if(r&1) chkmax(ans,seg[r^1]);
 34     }
 35     return ans;
 36 }
 37 void dfs(int u,int f){
 38     int mx = 0 ,s = -1;
 39     sz [u] = 1;
 40     for(int i= head[u]; i!=-1; i=nxt[i]){
 41         int v = pnt[i];
 42         if(v == f) continue;
 43         P[v] = i^1;
 44         deep[v] = deep[u]+1;
 45         dfs(v,u);
 46         sz[u] += sz[v];
 47         if(sz[v] > mx){
 48             mx = sz[v];
 49             s = i;
 50         }
 51     }
 52     heavy[u] = s;
 53     if(s != -1){
 54         parent[pnt[s]] = u;
 55     }
 56 }
 57 void prepare(){
 58     for(int i =30; i;i--) if((1<<i) > n) M = 1<<i;
 59     for(int i = M; i<2*M; i++) seg[i] = -inf;
 60     memset(num,-1,sizeof(num));
 61     segsz = deep[0] = 0;
 62     for(int i=0;i<n;i++) parent[i] = i;
 63     dfs(0,0);
 64 //    for(int i=0;i<n;i++) cout<<deep[i]<<" "; cout<<endl;
 65     for(int i=0; i<n; i++) if(heavy[i] == -1){
 66         int pos = i;
 67 //        cout<<"i:"<<i<<endl;
 68         while(pos && pnt[heavy[pnt[P[pos]]]] == pos){
 69             int t = P[pos];
 70 //            cout<<"pos : "<<pos<<" "<<t<<endl;
 71             num[t^1] = num[t] = ++segsz;
 72             ins(segsz,cost[t]);
 73             pos = pnt[P[pos]];
 74             //cout<<pos<<" "<<pnt[P[pos]]<<" "<<pnt[heavy[pnt[P[pos]]]]<<endl;
 75         }
 76     }
 77 //    for(int i=0;i<e;i++) cout<<num[i]<<" ";cout<<endl;
 78 }
 79 // operator
 80 void change(int pos,int val){
 81     if(num[pos] == -1) { cost[pos] = cost[pos^1] = val; return;}
 82     ins(num[pos],val);
 83 }
 84 int query(int u,int v){
 85     //cout<<u<<" "<<v<<endl;
 86     int ans = -inf;
 87     while(u!=v){
 88         int r = P[u];
 89         if(num[r] == -1){
 90             chkmax(ans,cost[r]);
 91             u = pnt[r];
 92         }
 93         else {
 94             int p = parent[u];
 95             if(deep[p] < deep[v]) p = v;
 96             u = p;
 97             //cout<<r<<" "<<heavy[p]<<endl;
 98             int l = num[heavy[p]]; r = num[r];
 99     //        cout<<"query : "<<l<<" "<<r<<endl;
100             assert(l>0 && l>=r);
101             int t = rmq(r,l);
102             chkmax(ans,t);
103         }
104     }
105     return ans;
106 }
107 int lca(int u,int v){
108     while(1){
109         int a = find(u) , b = find(v);
110         if(a == b) return deep[u]< deep[v] ? u : v;
111         else if(deep[a] >= deep[b]) u = pnt[P[a]]; else v = pnt[P[b]];
112     }
113 }
114 int ask(int u,int v){
115     int p = lca(u,v);
116     return max(query( u, p),query( v, p));
117 }
118 int main(){
119     int t; cin >> t;
120     while(t--){
121         int u,v,c; scanf("%d",&n);
122         e = 0; memset(head,-1,sizeof(head));
123         for(int i=1;i<n;i++){
124             scanf("%d%d%d",&u,&v,&c);
125             u--;v--;
126             add_edge(u,v,c);
127             add_edge(v,u,c);
128         }
129         prepare();
130         char ch[20];
131         while(scanf("%s",ch) && strcmp(ch,"DONE")){
132             if(ch[0]=='C'){
133                 scanf("%d%d",&u,&c);
134                 change((u-1)*2,c);
135             }
136             else {
137                 scanf("%d%d",&u,&v);
138                 printf("%d\n",ask(u-1,v-1));
139             }
140         }
141     }
142 }
143 
    
posted on 2012-05-14 22:17 西月弦 閱讀(855) 評論(2)  編輯 收藏 引用 所屬分類: 解題報告經典題目

FeedBack:
# re: spoj 375 樹鏈剖分+LCA+RMQ(zkw線段樹)
2012-06-04 15:54 | wuyiqi
大神,網上這方面的講解都比較晦澀,求細講  回復  更多評論
  
# re: spoj 375 樹鏈剖分+LCA+RMQ(zkw線段樹)
2012-06-04 19:08 | 西月弦
@wuyiqi
代碼比較渣... 有問題郵件聯系吧~  回復  更多評論
  
青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            欧美人在线视频| 在线亚洲美日韩| 欧美日韩国产黄| 欧美激情综合色综合啪啪| 欧美一区二区三区久久精品| 亚洲图片欧美一区| av成人激情| 性做久久久久久久久| 香蕉精品999视频一区二区| 亚洲天堂免费观看| 亚洲欧美在线免费观看| 久久婷婷麻豆| 欧美午夜国产| 一区二区在线免费观看| 亚洲精品国产精品国自产在线| 亚洲国产91| 午夜欧美不卡精品aaaaa| 亚洲尤物在线视频观看| 欧美成年人视频网站欧美| 99精品欧美| 欧美国产亚洲另类动漫| 国产区精品在线观看| 99精品福利视频| 亚洲欧美日韩综合| 亚洲国产精品精华液2区45| 国产精品久久久久永久免费观看 | 久久久久久久综合| 美女在线一区二区| 亚洲天堂av图片| 欧美日韩一级黄| 日韩午夜在线电影| 亚洲国产日韩美| 欧美高清在线播放| 亚洲精品孕妇| 亚洲欧洲精品一区二区三区 | 免费的成人av| 国产一区在线观看视频| 久久久国产精品一区二区中文| 一区二区三区国产| 国产精品美女在线观看| 亚洲综合第一| 欧美一二区视频| 激情综合久久| 亚洲第一黄网| 欧美激情综合亚洲一二区| 国产精品99久久99久久久二8| 亚洲第一主播视频| 国产精品v日韩精品v欧美精品网站| 欧美一区二区福利在线| 久久动漫亚洲| 亚洲欧美在线一区| 欧美一区二区在线视频| 亚洲国产综合在线看不卡| 亚洲精品视频免费观看| 国产欧美一区二区精品秋霞影院| 狂野欧美激情性xxxx| 亚洲国产免费看| 国产精一区二区三区| 欧美激情在线观看| 国产亚洲网站| 一区二区三区日韩精品| 亚洲视频在线观看三级| 久久久久久久国产| 久久激情视频免费观看| 欧美精品日韩| 欧美好吊妞视频| 一区二区三区在线免费观看| 亚洲欧美在线观看| 亚洲欧美高清| 国产伦精品一区二区三区免费| 99在线精品视频| 99国内精品久久| 欧美成人精品一区二区| 免费毛片一区二区三区久久久| 国产亚洲精品激情久久| 欧美亚洲系列| 蜜桃av综合| 亚洲精品三级| 国产精品成人一区二区三区夜夜夜| 亚洲精品日韩在线观看| 亚洲视频第一页| 国内精品久久久久伊人av| 久久婷婷久久| 亚洲一区二区视频| 久久夜色精品| 亚洲毛片视频| 国产亚洲精品aa| 免费不卡中文字幕视频| 99re6热只有精品免费观看| 午夜精品久久久久久久男人的天堂 | 先锋影音久久| 国产日韩精品视频一区| 久久久综合网| 香蕉av福利精品导航| 亚洲欧洲日韩在线| 久久综合99re88久久爱| 99re成人精品视频| 亚洲福利视频二区| 国产欧美 在线欧美| 欧美精品色综合| 蜜臀a∨国产成人精品| 欧美亚洲在线播放| 一区二区久久久久| 亚洲高清视频在线| 欧美激情1区2区3区| 久久国产乱子精品免费女| 亚洲免费在线电影| 亚洲永久免费av| 亚洲影院免费观看| 欧美在线电影| 久久国产精品一区二区| 久久黄金**| 欧美国产日韩一区二区| 欧美91视频| 亚洲国产高清一区| 亚洲国产视频a| 亚洲欧洲一区二区在线观看| 亚洲黄色精品| 亚洲午夜影视影院在线观看| 亚洲男人的天堂在线| 久久精选视频| 欧美日韩视频在线第一区| 国产精品另类一区| 国产一区二区三区网站 | 国产一区二区福利| 国产欧美欧洲在线观看| 在线国产日韩| 在线一区二区三区四区五区| 久久丁香综合五月国产三级网站| 久久久久国产精品一区二区| 国产欧亚日韩视频| 合欧美一区二区三区| 亚洲精品一区中文| 久久精品水蜜桃av综合天堂| 欧美88av| 久久er99精品| 国产精品美女xx| aa亚洲婷婷| 欧美搞黄网站| 久久电影一区| 国产精品午夜在线观看| 一个色综合导航| 欧美成人中文字幕| 久久精品国产精品| 国产精品少妇自拍| 亚洲午夜一二三区视频| 亚洲国内欧美| 欧美精品久久久久久久| 亚洲国产三级网| 欧美国产日本高清在线| 久久久夜夜夜| 亚洲久色影视| 日韩亚洲国产精品| 欧美日韩精品免费观看视频完整| 亚洲三级影院| 亚洲另类春色国产| 欧美视频专区一二在线观看| 亚洲午夜一区二区三区| 亚洲精品久久久久| 欧美色欧美亚洲高清在线视频| 亚洲视频福利| 久久本道综合色狠狠五月| 亚洲国产欧美在线人成| 亚洲激情综合| 国产一区二区三区在线免费观看 | 午夜精品成人在线| 亚洲在线免费| 最新成人在线| 亚洲综合激情| 久久久午夜精品| 日韩亚洲综合在线| 亚洲一区二区三区高清| 影音先锋日韩资源| 99伊人成综合| 99在线|亚洲一区二区| 亚洲一区二区三区中文字幕在线| 91久久久亚洲精品| 午夜视频在线观看一区二区| 一本色道久久88精品综合| 欧美在线观看你懂的| 一区二区三区日韩精品视频| 久久视频国产精品免费视频在线| 欧美一区二区观看视频| 欧美日韩1080p| 欧美激情精品久久久久久免费印度| 国产精品高清在线观看| 亚洲免费大片| 亚洲视频欧洲视频| 欧美久久婷婷综合色| 亚洲国产精品第一区二区三区| 激情久久久久久| 久热国产精品| 亚洲激情偷拍| 亚洲天堂网站在线观看视频| 欧美久久视频| 亚洲欧美久久久久一区二区三区| 亚洲视频在线看| 国产欧美一区二区三区国产幕精品| 亚洲视频 欧洲视频| 欧美一区成人|