bzoj3306 树

来源:互联网 发布:财务评价软件 编辑:程序博客网 时间:2024/06/07 22:59

题目链接:bzoj3306
题目大意:
给定一棵大小为 n 的有根点权树,支持以下操作:
• 换根
• 修改点权
• 查询子树最小值

题解:
线段树+树剖
倍增好快啊..打了个树剖约等于卡时过orz
其实只用分类讨论就好了..先把树弄出来
如果新根与要询问的x不在一棵子树内答案不变
如果x就是新根的话答案就是整棵树的最小值
如果x是新根的祖先就求出子树含新根的与x最近的点,求除去其子树外其他的最小值。这个用dfs序来搞。
线段树维护区间最小值啊。
诶我好辣鸡啊QwQ

#include<cstdio>#include<cstdlib>#include<cstring>#include<iostream>#include<algorithm>using namespace std;#define maxn 101000struct node{    int y,next;}a[maxn];int len,first[maxn];int id,dfn[maxn],la[maxn],fa[maxn],dep[maxn];int v[maxn],siz[maxn],son[maxn],top[maxn],mn[maxn<<1];int mymin(int x,int y){return (x<y)?x:y;}void ins(int x,int y){    ++len;a[len].y=y;    a[len].next=first[x];first[x]=len;}void dfs(int x){    siz[x]=son[x]=0;    for (int k=first[x];k!=-1;k=a[k].next)    {        int y=a[k].y;        if (y==fa[x]) continue;        fa[y]=x;dep[y]=dep[x]+1;        dfs(y);        if (siz[son[x]]>siz[x]) son[x]=y;        siz[x]+=siz[y];    }    siz[x]++;}void dfs2(int x,int tp){    dfn[x]=++id;top[x]=tp;    if (son[x]!=0) dfs2(son[x],tp);    for (int k=first[x];k!=-1;k=a[k].next)    {        int y=a[k].y;        if (y==fa[x] || y==son[x]) continue;        dfs2(y,y);    }    la[x]=id;}int lca(int x,int y){    while (top[x]!=top[y])    {        if (dep[x]>dep[y]) swap(x,y);        y=fa[top[y]];    }    return (dep[x]>dep[y])?y:x;}void change(int now,int l,int r,int x,int k){    if (l==r) {mn[now]=k;return;}    int mid=(l+r)>>1,lc=now<<1,rc=(now<<1)+1;    if (x<=mid) change(lc,l,mid,x,k);    else change(rc,mid+1,r,x,k);    mn[now]=mymin(mn[lc],mn[rc]);}int query(int now,int l,int r,int L,int R){    if (l==L && r==R) return mn[now];    int mid=(l+r)>>1,lc=now<<1,rc=(now<<1)+1;    if (R<=mid) return query(lc,l,mid,L,R);    if (L>mid) return query(rc,mid+1,r,L,R);    return mymin(query(lc,l,mid,L,mid),query(rc,mid+1,r,mid+1,R));}int ffind(int x,int y){    while (dep[top[x]]>dep[y])      {        if (fa[top[x]]!=y)  x=fa[top[x]];        else return top[x];    }    if (top[x]==top[y]) return son[y];    while (fa[x]!=y) x=fa[x];    return x;}int main(){    //freopen("a.in","r",stdin);    //freopen("a.out","w",stdout);    int rt,n,q,i,f,x,y;char s[10];    scanf("%d%d",&n,&q);id=0;    len=0;memset(first,-1,sizeof(first));    for (i=1;i<=n;i++)    {        scanf("%d%d",&f,&v[i]);        if (f==0) rt=i;else ins(f,i);    }    memset(mn,63,sizeof(mn));    fa[rt]=0;dep[rt]=1;dfs(rt);dfs2(rt,rt);    for (i=1;i<=n;i++) change(1,1,n,dfn[i],v[i]);    for (i=1;i<=q;i++)    {        scanf("%s",s);        if (s[0]=='V')        {            scanf("%d%d",&x,&y);            change(1,1,n,dfn[x],y);        }        else if (s[0]=='E')        {            scanf("%d",&x);            rt=x;        }        else if (s[0]=='Q')        {            scanf("%d",&x);            if (rt==x) printf("%d\n",mn[1]);            else if (lca(rt,x)!=x) printf("%d\n",query(1,1,n,dfn[x],la[x]));            else            {                int ans,w=ffind(rt,x);                if (dfn[w]>1) ans=query(1,1,n,1,dfn[w]-1);                if (la[w]<n) ans=mymin(ans,query(1,1,n,la[w]+1,n));                printf("%d\n",ans);            }        }    }    return 0;}
0 0
原创粉丝点击