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

算法學社
記錄難忘的征途
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>
            国产日韩欧美制服另类| 亚洲一区久久久| 久久网站免费| 亚洲国产成人精品女人久久久 | 欧美暴力喷水在线| 午夜精品久久久久久久99黑人| 亚洲黄一区二区| 亚洲美女在线视频| 99国产精品99久久久久久| 亚洲啪啪91| 一区二区免费在线播放| 一区二区三区精品视频在线观看| 亚洲精品网址在线观看| 一本久久综合| 亚洲综合欧美| 免费久久99精品国产自| 亚洲人成人77777线观看| 亚洲人成在线免费观看| 亚洲一区二区黄| 久久久国产精品一区| 欧美日韩一区二区三区在线 | 亚洲一区二区三区在线视频| 在线亚洲电影| 久久五月天婷婷| 99热在这里有精品免费| 久久久久久久久岛国免费| 欧美激情视频一区二区三区不卡| 欧美视频免费在线| 亚洲国产99| 美女图片一区二区| 欧美一区二区三区日韩视频| 欧美日韩在线电影| 亚洲人成欧美中文字幕| 美日韩在线观看| 欧美自拍丝袜亚洲| 久久精品91久久久久久再现| 久久中文久久字幕| 亚洲调教视频在线观看| 欧美国产一区二区三区激情无套| 国产一区美女| 久久久精品一区| 亚洲制服av| 国产精品久久久久久久久久免费 | 亚洲二区在线观看| 欧美中在线观看| 国产一区二区日韩| 国产视频欧美| 国产精品私房写真福利视频| 99视频热这里只有精品免费| 亚洲黄色免费网站| 欧美天天在线| 久久久久久久国产| 你懂的成人av| 欧美一区二区精品| 中文久久乱码一区二区| 国产精品一区二区三区成人| 久久岛国电影| 欧美日韩国产综合新一区| 国产日韩综合| 亚洲精品中文在线| 国产在线日韩| 亚洲一区二区毛片| 香蕉免费一区二区三区在线观看| 久久久一本精品99久久精品66| 亚洲国产精品女人久久久| 99精品久久久| 亚洲欧洲一区二区三区在线观看 | 狠狠色综合日日| 欧美高清影院| 久久在线观看视频| 欧美国产亚洲视频| 欧美视频官网| 亚洲精品一二| 亚洲国产乱码最新视频| 亚洲综合大片69999| 精品电影在线观看| 欧美韩国在线| 久久成人av少妇免费| 欧美国产国产综合| 亚洲国产精品va在线看黑人| 99国产精品国产精品久久| 精品成人a区在线观看| 亚洲欧美日韩精品| 欧美亚洲视频一区二区| 国产午夜精品美女视频明星a级| 99国产一区| 亚洲欧美日韩爽爽影院| 国产精品日本欧美一区二区三区| 亚洲美女诱惑| 欧美一区二区三区日韩| 极品av少妇一区二区| 欧美成人精品影院| av成人手机在线| 亚洲自拍高清| 91久久亚洲| 一本到12不卡视频在线dvd| 欧美精品在线视频| 亚洲欧美日韩网| 国产精品日韩| 久久精品国产99国产精品| 久久夜色精品国产欧美乱| 韩国精品一区二区三区| 欧美专区中文字幕| 久久综合中文| 日韩网站在线观看| 国产精品久久久久影院亚瑟| 欧美小视频在线| 亚洲精品欧美日韩专区| 一区二区三区产品免费精品久久75 | 欧美成年人网| 正在播放欧美一区| 久久久九九九九| 国产乱码精品一区二区三区五月婷| 午夜精品久久久久久久99水蜜桃 | 亚洲欧美变态国产另类| 久久久夜色精品亚洲| 中文av一区特黄| 亚洲大胆在线| 久久综合中文色婷婷| 香蕉久久夜色精品国产使用方法| 国内精品免费在线观看| 国产精品亚洲成人| 欧美理论电影在线观看| 久久久久一区二区| 久久国产精品久久精品国产| 欧美xxx成人| 亚洲电影免费在线观看| 亚洲精品一区二区三区福利| 狠狠色狠狠色综合日日tαg| 国产精品入口日韩视频大尺度| 欧美成人亚洲| 亚洲天堂成人在线视频| 老牛国产精品一区的观看方式| 亚洲一区激情| 亚洲在线视频一区| 久久精品国产99精品国产亚洲性色 | 亚洲乱码视频| 1000部精品久久久久久久久| 激情丁香综合| 亚洲国产精品成人va在线观看| 欧美精品尤物在线| 国产精品久久久久久久久久久久| 国产精品成人一区二区| 巨胸喷奶水www久久久免费动漫| 亚洲影视在线播放| 久久久美女艺术照精彩视频福利播放 | 日韩视频―中文字幕| 亚洲亚洲精品三区日韩精品在线视频| 亚洲欧洲日产国产网站| 亚洲综合色丁香婷婷六月图片| 欧美一区二区在线免费观看| 99视频国产精品免费观看| 欧美一区二区三区的| 欧美成人免费网站| 国产精品一区二区久久久| 亚洲人成亚洲人成在线观看| 亚洲永久免费视频| 亚洲福利视频一区| 蜜桃av一区二区三区| 久久男人资源视频| 这里是久久伊人| 欧美日韩色一区| 99re热这里只有精品免费视频| 香蕉精品999视频一区二区| 欧美xx视频| 久久三级福利| 午夜一区不卡| 国产精品一区一区三区| 欧美激情视频一区二区三区在线播放 | 欧美日韩美女在线| 日韩亚洲精品视频| 91久久久国产精品| 欧美激情精品久久久| 99精品免费| 亚洲欧美激情视频在线观看一区二区三区 | 亚洲经典在线看| 欧美成年人视频网站| 久久久久久久久久久久久女国产乱 | 国产欧美精品| 亚洲一区二区三区午夜| 亚洲一区二区三区免费在线观看| 小黄鸭精品aⅴ导航网站入口| 亚洲午夜精品福利| 欧美香蕉大胸在线视频观看| 久久精品国语| 午夜精品久久久久久久99黑人| 国产自产女人91一区在线观看| 欧美高清在线视频观看不卡| 欧美专区在线观看| 亚洲一区中文字幕在线观看| 午夜老司机精品| 免费亚洲一区| 久久动漫亚洲| 欧美三区在线观看| 欧美福利电影网| 国产亚洲综合在线| 一区二区三区四区五区视频| 136国产福利精品导航| 欧美激情aaaa| 一区二区三区www|