亚洲香蕉成人av网站在线观看_欧美精品成人91久久久久久久_久久久久久久久久久亚洲_热久久视久久精品18亚洲精品_国产精自产拍久久久久久_亚洲色图国产精品_91精品国产网站_中文字幕欧美日韩精品_国产精品久久久久久亚洲调教_国产精品久久一区_性夜试看影院91社区_97在线观看视频国产_68精品久久久久久欧美_欧美精品在线观看_国产精品一区二区久久精品_欧美老女人bb

首頁 > 學院 > 開發設計 > 正文

bzoj4539 [Hnoi2016]樹

2019-11-14 08:55:14
字體:
來源:轉載
供稿:網友

Description

  小A想做一棵很大的樹,但是他手上的材料有限,只好用點小技巧了。開始,小A只有一棵結點數為N的樹,結點的編號為1,2,…,N,其中結點1為根;我們稱這顆樹為模板樹。小A決定通過這棵模板樹來構建一顆大樹。構建過程如下:(1)將模板樹復制為初始的大樹。(2)以下(2.1)(2.2)(2.3)步循環執行M次(2.1)選擇兩個數字a,b,其中1<=a<=N,1<=b<=當前大樹的結點數。(2.2)將模板樹中以結點a為根的子樹復制一遍,掛到大樹中結點b的下方(也就是說,模板樹中的結點a為根的子樹復制到大樹中后,將成為大樹中結點b的子樹)。(2.3)將新加入大樹的結點按照在模板樹中編號的順序重新編號。例如,假設在進行2.2步之前大樹有L個結點,模板樹中以a為根的子樹共有C個結點,那么新加入模板樹的C個結點在大樹中的編號將是L+1,L+2,…,L+C;大樹中這C個結點編號的大小順序和模板樹中對應的C個結點的大小順序是一致的。下面給出一個實例。假設模板樹如下圖:

根據第(1)步,初始的大樹與模板樹是相同的。在(2.1)步,假設選擇了a=4,b=3。運行(2.2)和(2.3)后,得到新的大樹如下圖所示現在他想問你,樹中一些結點對的距離是多少。

Input

  第一行三個整數:N,M,Q,以空格隔開,N表示模板樹結點數,M表示第(2)中的循環操作的次數,Q 表示詢問數量。接下來N-1行,每行兩個整數 fr,to,表示模板樹中的一條樹邊。再接下來M行,每行兩個整數x,to,表示將模板樹中 x 為根的子樹復制到大樹中成為結點to的子樹的一次操作。再接下來Q行,每行兩個整數fr,to,表示詢問大樹中結點 fr和 to之間的距離是多少。N,M,Q<=100000

Output

  輸出Q行,每行一個整數,第 i行是第 i個詢問的答案。

Sample Input

5 2 3 1 4 1 3 4 2 4 5 4 3 3 2 6 9 1 8 5 3

Sample Output

633

HINT

經過兩次操作后,大樹變成了下圖所示的形狀:結點6到9之間經過了6條邊,所以距離為6;類似地,結點1到8之間經過了3條邊;結點5到3之間也經過了3條邊。

正解:dfs序主席樹+lca

惡心數據結構題+碼農題。。前后總共調了6個多小時。。WA顯示成RE也無語至極。。

具體做法:

  首先對于模板樹進行預處理,dfs一遍得到dfs序,為了維護子樹第k小編號的查詢操作,構主席樹。

 ?。╬s:在主席樹上查找子樹第k小編號是一個經典問題,按照dfs序依次把每個編號相應的位置+1,然后就是常規的區間第k小問題)

  對于復制操作,復制操作的話對于每次復制我只需要記錄這次復制的這棵子樹的根,根接到了哪個點的下方,當然上述記錄的都是在原樹中的相應編號。同時為了方便之后計算這棵子樹內到根的距離,不妨記錄一下根在原樹到1的距離,這樣以來當我想查詢這棵子樹內的某個點到根的距離時,直接用在原樹中與1的距離之差即可。

  對于最后的查詢操作,需要仔細考慮了,有很多細節。當x和y處在同一塊中(也就是新樹的同一結點上時),直接查詢。否則,先求出在x、y各自所處的塊(也就是新樹上的兩個結點,不妨設為Rx、Ry,且deep[Ry]<deep[Rx])在新樹上的lca,如果Ry=lca,則不用處理;否則就把x、y先跳到所在的塊的頂端,記錄貢獻,再一直跳到lca的兒子結點(實際上就是兒子結點所表示的塊中的根結點),接著往上走一步即可進入lca的塊中。對于同一塊中的直接查詢即可。細節的話太多不贅述了,最重要的一點就是因為最大情況下,點數肯定是超int的,所以最好是都開longlong。

//It is made by wfj_2048~#include <algorithm>#include <iostream>#include <cstring>#include <cstdlib>#include <cstdio>#include <vector>#include <cmath>#include <queue>#include <stack>#include <map>#include <set>#define N (100010)#define il inline#define RG register#define ll long long#define File(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout)using namespace std;il ll gll(){    ll x=0,q=1; char ch=getchar(); while ((ch<'0' || ch>'9') && ch!='-') ch=getchar();    if (ch=='-') q=-1,ch=getchar(); while (ch>='0' && ch<='9') x=x*10+ch-48,ch=getchar(); return q*x;}struct edge{ ll nt,to,dis; };struct tree{    ll head[N],dfn[N],tid[N],dis[N],id[N],anc[N],up[N],size[N],top[N],fa[N],son[N],dep[N],n,num,cnt,tot;    edge g[2*N];        il void insert(ll from,ll to,ll dis){ g[++num]=(edge){head[from],to,dis},head[from]=num; return; }        il void dfs1(ll x,ll p){	dep[x]=dep[p]+1,fa[x]=p,size[x]=1; ll mx=0,v;	for (RG ll i=head[x];i;i=g[i].nt){	    v=g[i].to; if (v==p) continue;	    dis[v]=dis[x]+g[i].dis; dfs1(v,x);	    size[x]+=size[v]; if (size[mx]<=size[v]) mx=v;	}	son[x]=mx; return;    }        il void dfs2(ll x,ll p,ll a){	dfn[++cnt]=x,tid[x]=cnt,top[x]=a;	if (son[x]) dfs2(son[x],x,a); ll v;	for (RG ll i=head[x];i;i=g[i].nt){	    v=g[i].to; if (v==p || v==son[x]) continue;	    dfs2(v,x,v);	}	return;    }        il ll lca(ll u,ll v){	while (top[u]!=top[v]){	    if (dep[top[u]]<dep[top[v]]) swap(u,v);	    u=fa[top[u]];	}	return dep[u]>dep[v] ? v : u;    }        il ll jump(ll u,ll v){	ll x=u; while (top[u]!=top[v]) x=top[u],u=fa[top[u]];	return u==v ? x : son[v];    }    }tr1,tr2;struct chairtree{    ll root[N],sum[20*N],ls[20*N],rs[20*N],sz;        il void build(ll x,ll &y,ll l,ll r,ll v){	sum[y=++sz]=sum[x]+1,ls[y]=ls[x],rs[y]=rs[x];	if (l==r) return; ll mid=(l+r)>>1;	if (v<=mid) build(ls[x],ls[y],l,mid,v);	else build(rs[x],rs[y],mid+1,r,v);	return;    }        il ll query(ll x,ll y,ll l,ll r,ll k){	if (l==r) return l; ll mid=(l+r)>>1;	if (k<=sum[ls[y]]-sum[ls[x]]) return query(ls[x],ls[y],l,mid,k);	else return query(rs[x],rs[y],mid+1,r,k-sum[ls[y]]+sum[ls[x]]);    }    }ch;il ll find(ll x){    if (x<=tr1.n) return 1;    ll l=1,r=tr2.n,mid,ans=1;    while (l<=r){	mid=(l+r)>>1;	if (x<tr2.id[mid]) r=mid-1;	else ans=mid,l=mid+1;    }    return ans;}il ll query(ll blv,ll v){    if (v<=tr1.n) return v; v-=tr2.id[blv]-1;    ll x=tr2.anc[blv],l=tr1.tid[x],r=l+tr1.size[x]-1;    return ch.query(ch.root[l-1],ch.root[r],1,tr1.n,v);}il void work(){    tr1.n=gll(); ll m=gll(),q=gll();    for (RG ll i=1;i<tr1.n;++i){ ll u=gll(),v=gll(); tr1.insert(u,v,1),tr1.insert(v,u,1); }    tr1.dfs1(1,0),tr1.dfs2(1,0,1); tr2.id[1]=1,tr2.anc[1]=1,tr2.n=1,tr2.tot=tr1.n;    for (RG ll i=1;i<=tr1.n;++i) ch.build(ch.root[i-1],ch.root[i],1,tr1.n,tr1.dfn[i]);    for (RG ll i=1;i<=m;++i){	ll u=gll(),v=gll(),blv=find(v),idv=query(blv,v);	tr2.n++,tr2.id[tr2.n]=tr2.tot+1,tr2.tot+=tr1.size[u],tr2.anc[tr2.n]=u,tr2.up[tr2.n]=idv;	tr2.insert(blv,tr2.n,tr1.dis[idv]-tr1.dis[tr2.anc[blv]]+1);    }    tr2.dfs1(1,0),tr2.dfs2(1,0,1);    for (RG ll i=1;i<=q;++i){	ll u=gll(),v=gll(),blu=find(u),blv=find(v);	ll idu=query(blu,u),idv=query(blv,v),lca,ans=0;	if (blu==blv){	    lca=tr1.lca(idu,idv);	    ans=tr1.dis[idu]+tr1.dis[idv]-2*tr1.dis[lca];	    PRintf("%lld/n",ans); continue;	}	if (tr2.dep[blu]>tr2.dep[blv]) swap(u,v),swap(blu,blv),swap(idu,idv);	lca=tr2.lca(blu,blv);	if (blu!=lca){	    ans=tr1.dis[idu]-tr1.dis[tr2.anc[blu]];	    u=tr2.jump(blu,lca); ans+=tr2.dis[blu]-tr2.dis[u]+1;	    u=tr2.up[u];	}else u=query(blu,u);	ans+=tr1.dis[idv]-tr1.dis[tr2.anc[blv]];	v=tr2.jump(blv,lca); ans+=tr2.dis[blv]-tr2.dis[v]+1;	v=tr2.up[v],lca=tr1.lca(u,v);	ans+=tr1.dis[u]+tr1.dis[v]-2*tr1.dis[lca];	printf("%lld/n",ans);    }    return;}int main(){    File("tree");    work();    return 0;}//鬼畜dfs序+主席樹+lca沃日


上一篇:c#復制文件

下一篇:1048. 數字加密(20)

發表評論 共有條評論
用戶名: 密碼:
驗證碼: 匿名發表
亚洲香蕉成人av网站在线观看_欧美精品成人91久久久久久久_久久久久久久久久久亚洲_热久久视久久精品18亚洲精品_国产精自产拍久久久久久_亚洲色图国产精品_91精品国产网站_中文字幕欧美日韩精品_国产精品久久久久久亚洲调教_国产精品久久一区_性夜试看影院91社区_97在线观看视频国产_68精品久久久久久欧美_欧美精品在线观看_国产精品一区二区久久精品_欧美老女人bb
2019亚洲日韩新视频| 亚洲a在线观看| 国产精品99久久久久久www| 国产精品美女www| 国产精品视频午夜| 国产精品av免费在线观看| 国产精品一区专区欧美日韩| 91免费国产网站| 精品久久久久久| 国产一区玩具在线观看| 91免费精品国偷自产在线| 亚洲国产欧美在线成人app| 992tv成人免费影院| 日韩av有码在线| 国产一区二区三区三区在线观看| 亚洲电影在线看| 久久久国产精品x99av| 国产精品一区二区久久| 久久天天躁狠狠躁夜夜爽蜜月| 精品欧美一区二区三区| 96精品久久久久中文字幕| 国产91久久婷婷一区二区| 日韩精品久久久久久福利| 欧美激情区在线播放| 久久精品国产亚洲精品2020| 国产精品麻豆va在线播放| 九九热精品视频国产| 欧美激情中文字幕在线| 国产欧美日韩精品丝袜高跟鞋| 亚洲欧美制服另类日韩| 欧美国产极速在线| 精品国产乱码久久久久久虫虫漫画| 日韩女在线观看| 91chinesevideo永久地址| 51ⅴ精品国产91久久久久久| 欧美人与性动交a欧美精品| 国产精品wwww| 欧美在线亚洲一区| 欧美日韩一区二区免费在线观看| 亚洲欧美日韩天堂一区二区| 中文字幕日韩精品在线观看| 久久久久国色av免费观看性色| 国产精品日韩在线一区| 九九热这里只有在线精品视| 91香蕉嫩草神马影院在线观看| 欧美一乱一性一交一视频| 91日本视频在线| 国产欧美日韩精品在线观看| 91极品女神在线| 欧美裸体视频网站| 久久影视电视剧免费网站| 日韩成人中文字幕在线观看| 欧美疯狂做受xxxx高潮| 欧美有码在线观看视频| 91久久久久久久久| 色狠狠久久aa北条麻妃| 91国产中文字幕| 九九热最新视频//这里只有精品| 国产精品看片资源| 久久久久久久久久久久久久久久久久av| 国产精品香蕉国产| 欧美精品日韩www.p站| 日本一区二区三区在线播放| 国产精品27p| 国产欧美精品一区二区三区介绍| 久久久久日韩精品久久久男男| 成人精品一区二区三区| 狠狠色噜噜狠狠狠狠97| 亚洲欧美日韩中文在线| 久久久亚洲影院你懂的| 国产91精品久久久久久| 日韩免费在线视频| 国产精品成人va在线观看| 国产精品777| 91av在线看| 91av免费观看91av精品在线| 亚洲第一国产精品| 欧美精品一本久久男人的天堂| 午夜精品久久久久久久99黑人| 欧美日韩国产精品一区| 国产香蕉97碰碰久久人人| 精品亚洲国产视频| 久久久999成人| 91国偷自产一区二区三区的观看方式| 精品国产一区二区三区久久| 亚洲剧情一区二区| 在线看片第一页欧美| 羞羞色国产精品| 久久综合久中文字幕青草| 欧美午夜激情视频| 欧美成人免费视频| 91视频免费网站| 久久久女人电视剧免费播放下载| 亚洲电影免费观看高清完整版| 久久久久久久影院| 欧美最猛性xxxx| 91久久精品国产91性色| 国产精品久久久久久亚洲影视| 8090理伦午夜在线电影| 欧美日韩国产黄| 国产69精品久久久| 成人久久18免费网站图片| 亚洲伊人久久综合| 亚洲天堂第一页| 国产精品一区二区久久国产| 日本亚洲欧美成人| 欧洲s码亚洲m码精品一区| 色天天综合狠狠色| 久久精品亚洲94久久精品| 欧美精品久久久久久久免费观看| 久久精品视频在线观看| 日韩综合视频在线观看| 日韩三级影视基地| 欧美视频在线免费看| 26uuu国产精品视频| 人人做人人澡人人爽欧美| 国产精品视频最多的网站| 欧美日本中文字幕| 欧美日韩国产在线看| 欧洲精品在线视频| 国产精品久久久久久久一区探花| 韩日欧美一区二区| 亚洲xxxxx电影| www.亚洲一区| 成人免费视频xnxx.com| 亚洲黄页视频免费观看| 亚洲欧洲在线播放| 91免费看视频.| 欧美精品精品精品精品免费| 亚洲精品国产美女| 亚洲一区免费网站| 国产日韩欧美在线视频观看| 欧美一级大片视频| 久久在线视频在线| 91牛牛免费视频| 韩国视频理论视频久久| 国产在线日韩在线| 中文字幕视频在线免费欧美日韩综合在线看| 91九色视频在线| 国产日产亚洲精品| 韩国三级日本三级少妇99| 亚洲a成v人在线观看| 庆余年2免费日韩剧观看大牛| 欧美激情免费在线| 亚洲第一精品电影| 韩国国内大量揄拍精品视频| 国产精品视频不卡| 亚洲美女免费精品视频在线观看| 国产97人人超碰caoprom| 中文字幕日韩视频| 国产精品永久免费视频| 97碰碰碰免费色视频| 国产精品草莓在线免费观看| 亚洲精品第一页| 国产一区二区三区在线| 91精品国产色综合久久不卡98口| 色狠狠av一区二区三区香蕉蜜桃| 国产午夜精品久久久| 国产精品亚洲аv天堂网| 精品国产91久久久久久| 色妞色视频一区二区三区四区| 91高潮精品免费porn| 91超碰中文字幕久久精品| 欧美成人精品三级在线观看|