[bzoj1103][树链剖分]大都市
来源:互联网 发布:java中char 编辑:程序博客网 时间:2024/06/11 23:26
Description
在经济全球化浪潮的影响下,习惯于漫步在清晨的乡间小路的邮递员Blue Mary也开始骑着摩托车传递邮件了。
不过,她经常回忆起以前在乡间漫步的情景。昔日,乡下有依次编号为1..n的n个小村庄,某些村庄之间有一些双向的土路。从每个村庄都恰好有一条路径到达村庄1(即比特堡)。并且,对于每个村庄,它到比特堡的路径恰好只经过编号比它的编号小的村庄。另外,对于所有道路而言,它们都不在除村庄以外的其他地点相遇。在这个未开化的地方,从来没有过高架桥和地下铁道。随着时间的推移,越来越多的土路被改造成了公路。
至今,Blue Mary还清晰地记得最后一条土路被改造为公路的情景。现在,这里已经没有土路了——所有的路都成为了公路,而昔日的村庄已经变成了一个大都市。 Blue Mary想起了在改造期间她送信的经历。她从比特堡出发,需要去某个村庄,并且在两次送信经历的间隔期间,有某些土路被改造成了公路.现在Blue Mary需要你的帮助:计算出每次送信她需要走过的土路数目。(对于公路,她可以骑摩托车;而对于土路,她就只好推车了。)
Input
第一行是一个数n(1 < = n < = 2 50000).以下n-1行,每行两个整数a,b(1 < = a以下一行包含一个整数m
(1 < = m < = 2 50000),表示Blue Mary曾经在改造期间送过m次信。以下n+m-1行,每行有两种格式的若干信息
,表示按时间先后发生过的n+m-1次事件:若这行为 A a b(a若这行为 W a, 则表示Blue Mary曾经从比特堡送信到村庄a。
Output
有m行,每行包含一个整数,表示对应的某次送信时经过的土路数目。
Sample Input
5
1 2
1 3
1 4
4 5
4
W 5
A 1 4
W 5
A 4 5
W 5
W 2
A 1 2
A 1 3
Sample Output
2
1
0
1
题解
树链剖分裸题吧。。线段树维护区间内土路个数,然后跑树剖就好。。
#include<cstdio>#include<cstring>#include<cstdlib>#include<algorithm>#include<cmath>using namespace std;typedef long long LL;const int MAXN=410000;struct node{ int x,y,next;}a[1110000];int len,last[MAXN];void ins(int x,int y){ len++; a[len].x=x;a[len].y=y; a[len].next=last[x];last[x]=len;}struct trnode{ int lc,rc,l,r; int c,maxnum;}tr[1110000];int trlen;void bt(int l,int r){ int now=++trlen; tr[now].l=l;tr[now].r=r;tr[now].c=0;tr[now].maxnum=0; tr[now].lc=tr[now].rc=-1; if(l<r) { int mid=(l+r)/2; tr[now].lc=trlen+1;bt(l,mid); tr[now].rc=trlen+1;bt(mid+1,r); }}int fa[MAXN],dep[MAXN],tot[MAXN],son[MAXN];int n;void pre_tree_node(int x){ tot[x]=1;son[x]=0; for(int k=last[x];k;k=a[k].next) { int y=a[k].y; if(y!=fa[x]) { fa[y]=x; dep[y]=dep[x]+1; pre_tree_node(y); if(tot[y]>tot[son[x]])son[x]=y; tot[x]+=tot[y]; } }}int top[MAXN],ys[MAXN],z;void pre_tree_edge(int x,int tp){ ys[x]=++z;top[x]=tp; if(son[x]!=0)pre_tree_edge(son[x],tp); for(int k=last[x];k;k=a[k].next) { int y=a[k].y; if(y!=fa[x] && y!=son[x])pre_tree_edge(y,y); }}void change(int now,int pp,int k){ if(tr[now].l==tr[now].r){tr[now].maxnum=k;return ;} int lc=tr[now].lc,rc=tr[now].rc; int mid=(tr[now].l+tr[now].r)/2; if(pp<=mid)change(lc,pp,k); else change(rc,pp,k); tr[now].maxnum=tr[lc].maxnum+tr[rc].maxnum;}int findsum(int now,int l,int r){ if(tr[now].l==l && tr[now].r==r)return tr[now].maxnum; int lc=tr[now].lc,rc=tr[now].rc; int mid=(tr[now].l+tr[now].r)/2; if(r<=mid)return findsum(lc,l,r); else if(mid+1<=l)return findsum(rc,l,r); else return (findsum(lc,l,mid)+findsum(rc,mid+1,r));}int ans;int p[MAXN];int solvesum(int x,int y){ ans=0; int tx=top[x],ty=top[y]; while(tx!=ty) { if(dep[tx]>dep[ty]) { swap(tx,ty); swap(x,y); } ans+=findsum(1,ys[ty],ys[y]); y=fa[ty];ty=top[y]; } if(x==y)return ans+p[x]; else { if(dep[x]>dep[y])swap(x,y); ans+=findsum(1,ys[x],ys[y]); return ans; }}struct bian{ int x,y,c;}e[MAXN];char ss[10];int u;int v;int main(){// freopen("Test.txt","w",stdout); scanf("%d",&n); len=0;memset(last,0,sizeof(last)); for(int i=1;i<n;i++) { scanf("%d%d",&e[i].x,&e[i].y); ins(e[i].x,e[i].y); ins(e[i].y,e[i].x); } int x; scanf("%d",&x); dep[1]=0;fa[1]=0; pre_tree_node(1); z=0;pre_tree_edge(1,1); trlen=0;bt(1,z); for(int i=1;i<n;i++) { if(dep[e[i].x]>dep[e[i].y])swap(e[i].x,e[i].y); change(1,ys[e[i].y],1); } for(int T=1;T<=n+x-1;T++) { scanf("%s",ss+1); scanf("%d",&u); if(ss[1]=='A') { scanf("%d",&v); change(1,ys[v],0); } else printf("%d\n",solvesum(1,u)); } return 0;}
- [bzoj1103][树链剖分]大都市
- 【bzoj1103】大都市meg 树链剖分+线段树
- bzoj1103: [POI2007]大都市meg(树链剖分)
- bzoj1103【POI2007】大都市meg
- 【POI2007】【BZOJ1103】大都市meg
- bzoj1103[POI2007]大都市meg
- 【bzoj1103】【POI2007】【大都市】
- 【BZOJ1103】【POI2007】大都市meg
- bzoj1103: [POI2007]大都市meg
- 【BZOJ1103】[POI2007]大都市meg
- 【POI2007】bzoj1103 大都市
- bzoj1103: [POI2007]大都市meg
- bzoj1103 [POI2007]大都市meg
- BZOJ1103: [POI2007]大都市meg
- bzoj1103 [POI2007]大都市meg
- BZOJ1103: [POI2007]大都市meg
- 8.3 bzoj1103 [POI2007]大都市meg
- BZOJ1103 [POI2007]大都市meg 【树剖】
- 工具类-软键盘
- 将表达式转换为逆波兰表达式-LintCode
- MATLAB函数之imfinfo
- 不错的vim python配置
- 搜索练习 2002t2字串变换
- [bzoj1103][树链剖分]大都市
- windows下vue.js的安装与配置
- 1)Java语法基础知识点(可以借鉴)
- Java设计一个分数类。要求:1.分类包含的分数运算有:连个分数的加、减、乘、除运算。
- Quartz定时器
- 对机器学习的理解
- java开发环境搭建
- Eclipse中的Maven项目报错处理方法
- 工具类-万能Adapter