SP375 QTREE - Query on a tree

SP375 QTREE - Query on a tree

我是借这道题来说说如何从c++改到c的

1.我怕麻烦,所以把结构体拆了(忍痛割爱我封装的线段树)

2.max、swap函数进行了手写

max:

ll llmax(ll x,ll y){return x>y?x:y;}

swap:

x^=y^=x^=y

emm...这的确是一个神奇的swap,不用函数

至于为什么是对的用人类智慧法即可

c++代码:(线段树封装)

#include<bits/stdc++.h>
using namespace std;
const int N=100005;
typedef long long ll;
struct line{
    int from,to;
    ll val;
    int conec;
}s[N];
int n;
int hed[N],tal[N<<1],nxt[N<<1],cnt=0;
ll val[N<<1];
int fa[N];
int dep[N];
int size[N]={0};
int son[N]={0};
int dfn[N];
int rk[N];
int indexy=0;
int top[N];
ll v[N];
ll llmax(ll x,ll y){return x>y?x:y;}
struct Sugment_Tree{
    ll t[N<<2];
    #define mid (l+r)/2
    void push_up(int num){
        t[num]=llmax(t[num<<1],t[num<<1|1]);
    }
    void build(int l,int r,int num){
        if(l==r){
            t[num]=v[l];
            return;
        }
        build(l,mid,num<<1),build(mid+1,r,num<<1|1);
        push_up(num);
    }
    void upd(int l,int r,int num,int pos,ll SUM){
        if(l>pos||r<pos) return;
        if(l==r&&r==pos){
            t[num]=SUM;
            return;
        }
        upd(l,mid,num<<1,pos,SUM);
        upd(mid+1,r,num<<1|1,pos,SUM);
        push_up(num);
    } 
    ll ask(int l,int r,int num,int L,int R){
        if(l>R||r<L) return 0;
        if(l>=L&&r<=R) return t[num];
        return llmax(ask(l,mid,num<<1,L,R),ask(mid+1,r,num<<1|1,L,R));
    }
}T;
void addege(int x,int y,ll z){
    cnt++;
    tal[cnt]=y;
    val[cnt]=z;
    nxt[cnt]=hed[x];
    hed[x]=cnt;
}
void dfs1(int u,int pa){
    fa[u]=pa;
    size[u]=1;
    for(int i=hed[u];i;i=nxt[i]){
        int v=tal[i];
        if(v==fa[u]) continue;
        dep[v]=dep[u]+1;
        dfs1(v,u);
        size[u]+=size[v];
        if(size[v]>size[son[u]]) son[u]=v;
    }
}
void dfs2(int u,int tp){
    indexy++;
    dfn[u]=indexy;
    rk[dfn[u]]=u;
    top[u]=tp;
    if(!son[u]) return;
    dfs2(son[u],tp);
    for(int i=hed[u];i;i=nxt[i]){
        int v=tal[i];
        if(v==fa[u]||v==son[u]) continue;
        dfs2(v,v);
    }
}

ll cgn(int x,int y){
    ll ans=0;
    //LCA不能算进来!!!!!! 
    while(top[x]!=top[y]){
        if(dep[top[x]]<dep[top[y]]) swap(x,y);
        ans=llmax(ans,T.ask(1,n,1,dfn[top[x]],dfn[x]));
        x=fa[top[x]];
    }
    if(dep[x]<dep[y]) swap(x,y);
    ans=llmax(ans,T.ask(1,n,1,dfn[y]+1,dfn[x]));
    return ans;
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<n;i++){
        scanf("%d%d%lld",&s[i].from,&s[i].to,&s[i].val);    
        addege(s[i].from,s[i].to,s[i].val);
        addege(s[i].to,s[i].from,s[i].val);
    }
    dep[1]=1;
    dfs1(1,1);
    dfs2(1,1);
    for(int i=1;i<n;i++){
        if(fa[s[i].from]==s[i].to) v[dfn[s[i].from]]=s[i].val,s[i].conec=dfn[s[i].from];
        else v[dfn[s[i].to]]=s[i].val,s[i].conec=dfn[s[i].to];
    } 
    /*
    for(int i=1;i<=n;i++)
        cout<<v[i]<<" ";
    cout<<endl;
    */
    T.build(1,n,1);
    char c[10];
    while(1){
        scanf("%s",c);
        if(c[0]=='D') break;
        else if(c[0]=='Q'){
            int x,y;
            scanf("%d%d",&x,&y);
            if(x==y){
                printf("0\n");
                continue;
            }
            printf("%lld\n",cgn(x,y));
        }
        else{
            int pos;
            ll va;
            scanf("%d%lld",&pos,&va);
            T.upd(1,n,1,s[pos].conec,va);
        }
    }
    return 0;
} 

c:

#include <ctype.h>
#include <stdio.h>
#include <limits.h>
#include <stdlib.h>
#define N 100005
typedef long long ll;
int from_P[N];
int to_P[N];
ll val_P[N];
int conec_P[N];
int n;
int hed[N],tal[N<<1],nxt[N<<1],cnt=0;
ll val[N<<1];
int fa[N];
int dep[N];
int size[N]={0};
int son[N]={0};
int dfn[N];
int rk[N];
int indexy=0;
int top[N];
ll v[N];
int TE;
ll llmax(ll x,ll y){return x>y?x:y;}
    ll t[N<<2];
    #define mid (l+r)/2
    void push_up(int num)
    {
        t[num]=llmax(t[num<<1],t[num<<1|1]);
    }
    void build(int l,int r,int num)
    {
        if(l==r)
        {
            t[num]=v[l];
            return;
        }
        build(l,mid,num<<1),build(mid+1,r,num<<1|1);
        push_up(num);
    }
    void upd(int l,int r,int num,int pos,ll SUM)
    {
        if(l>pos||r<pos) return;
        if(l==r&&r==pos)
        {
            t[num]=SUM;
            return;
        }
        upd(l,mid,num<<1,pos,SUM);
        upd(mid+1,r,num<<1|1,pos,SUM);
        push_up(num);
    } 
    ll ask(int l,int r,int num,int L,int R)
    {
        if(l>R||r<L) return 0;
        if(l>=L&&r<=R) return t[num];
        return llmax(ask(l,mid,num<<1,L,R),ask(mid+1,r,num<<1|1,L,R));
    }
void addege(int x,int y,ll z)
{
    cnt++;
    tal[cnt]=y;
    val[cnt]=z;
    nxt[cnt]=hed[x];
    hed[x]=cnt;
}
void dfs1(int u,int pa)
{
    fa[u]=pa;
    size[u]=1;
    for(int i=hed[u];i;i=nxt[i])
    {
        int v=tal[i];
        if(v==fa[u]) continue;
        dep[v]=dep[u]+1;
        dfs1(v,u);
        size[u]+=size[v];
        if(size[v]>size[son[u]]) son[u]=v;
    }
}
void dfs2(int u,int tp)
{
    indexy++;
    dfn[u]=indexy;
    rk[dfn[u]]=u;
    top[u]=tp;
    if(!son[u]) return;
    dfs2(son[u],tp);
    for(int i=hed[u];i;i=nxt[i])
    {
        int v=tal[i];
        if(v==fa[u]||v==son[u]) continue;
        dfs2(v,v);
    }
}

ll cgn(int x,int y)
{
    ll ans=0;
    //LCA不能算进来!!!!!! 
    while(top[x]!=top[y])
    {
        if(dep[top[x]]<dep[top[y]]) x^=y^=x^=y;
        ans=llmax(ans,ask(1,n,1,dfn[top[x]],dfn[x]));
        x=fa[top[x]];
    }
    if(dep[x]<dep[y]) x^=y^=x^=y;
    ans=llmax(ans,ask(1,n,1,dfn[y]+1,dfn[x]));
    return ans;
}
int main()
{
    scanf("%d",&TE);
    while(TE--)
    {
        cnt=0;
        indexy=0;
        for(int i=0;i<N;i++) hed[i]=0;
        for(int i=0;i<N;i++) dep[i]=0;
        for(int i=0;i<N;i++) size[i]=0;
        for(int i=0;i<N;i++) son[i]=0;
        scanf("%d",&n);
        for(int i=1;i<n;i++)
        {
            scanf("%d%d%lld",&from_P[i],&to_P[i],&val_P[i]);    
            addege(from_P[i],to_P[i],val_P[i]);
            addege(to_P[i],from_P[i],val_P[i]);
        }
        dep[1]=1;
        dfs1(1,1);
        dfs2(1,1);
        for(int i=1;i<n;i++)
        {
            if(fa[from_P[i]]==to_P[i]) v[dfn[from_P[i]]]=val_P[i],conec_P[i]=dfn[from_P[i]];
            else v[dfn[to_P[i]]]=val_P[i],conec_P[i]=dfn[to_P[i]];
        } 
        /*
        for(int i=1;i<=n;i++)
            cout<<v[i]<<" ";
        cout<<endl;
        */
        build(1,n,1);
        char c[10];
        while(1)
        {
            scanf("%s",c);
            if(c[0]=='D') break;
            else if(c[0]=='Q')
            {
                int x,y;
                scanf("%d%d",&x,&y);
                if(x==y)
                {
                    printf("0\n");
                    continue;
                }
                printf("%lld\n",cgn(x,y));
            }
            else
            {
                int pos;
                ll va;
                scanf("%d%lld",&pos,&va);
                upd(1,n,1,conec_P[pos],va);
            }
        }
    }
    return 0;
} 

 

上一篇:BZOJ1812: [Ioi2005]riv(树形dp)


下一篇:1481F.AB Tree(树上信息统计+01背包+记录DP路径+Bitset优化时间复杂度)