Cod sursa(job #1970638)

Utilizator george_stelianChichirim George george_stelian Data 19 aprilie 2017 15:05:08
Problema Heavy Path Decomposition Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 2.91 kb
#include <bits/stdc++.h>

using namespace std;

vector<int> g[100010];
int niv[100010],size[100010],val[100010];
int lant[100010],poz_lant[100010],lant_size[100010],lant_tata[100010],nr_lanturi;
vector<int> arb[100010];

void dfs1(int nod,int tata,int lev)
{
    size[nod]=1;
    niv[nod]=lev;
    for(vector<int>::iterator it=g[nod].begin();it!=g[nod].end();it++)
        if(*it!=tata)
        {
            dfs1(*it,nod,lev+1);
            size[nod]+=size[*it];
        }
}

void dfs2(int nod,int tata,int l)
{
    lant[nod]=l;
    poz_lant[nod]=++lant_size[l];
    int heavy=0;
    for(vector<int>::iterator it=g[nod].begin();it!=g[nod].end();it++)
        if(*it!=tata && size[*it]>size[heavy]) heavy=*it;
    if(!heavy) return;
    dfs2(heavy,nod,l);
    for(vector<int>::iterator it=g[nod].begin();it!=g[nod].end();it++)
        if(*it!=tata && *it!=heavy)
        {
            lant_tata[++nr_lanturi]=nod;
            dfs2(*it,nod,nr_lanturi);
        }
}

void arbint_update(int nod,int st,int dr,int poz,int val,vector<int> &arb)
{
    if(st==dr) arb[nod]=val;
    else
    {
        int mid=(st+dr)/2;
        if(poz<=mid) arbint_update(nod*2,st,mid,poz,val,arb);
        else arbint_update(nod*2+1,mid+1,dr,poz,val,arb);
        arb[nod]=max(arb[nod*2],arb[nod*2+1]);
    }
}

int arbint_query(int nod,int st,int dr,int left,int right,vector<int> &arb)
{
    if(left<=st && dr<=right) return arb[nod];
    else
    {
        int mid=(st+dr)/2,ret=0;
        if(left<=mid) ret=max(ret,arbint_query(nod*2,st,mid,left,right,arb));
        if(mid<right) ret=max(ret,arbint_query(nod*2+1,mid+1,dr,left,right,arb));
        return ret;
    }
}

int main()
{
    freopen("heavypath.in", "r" ,stdin);
    freopen("heavypath.out", "w", stdout);
    int n,m,x,y,tip;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",&val[i]);
    for(int i=1;i<n;i++)
    {
        scanf("%d%d",&x,&y);
        g[x].push_back(y);
        g[y].push_back(x);
    }
    dfs1(1,0,1);
    nr_lanturi++;
    dfs2(1,0,nr_lanturi);
    for(int i=1;i<=nr_lanturi;i++) arb[i].resize(4*lant_size[i]);
    for(int i=1;i<=n;i++)
        arbint_update(1,1,lant_size[lant[i]],poz_lant[i],val[i],arb[lant[i]]);
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d%d",&tip,&x,&y);
        if(tip==0) arbint_update(1,1,lant_size[lant[x]],poz_lant[x],y,arb[lant[x]]);
        else
        {
            int ret=0;
            while(lant[x]!=lant[y])
            {
                if(niv[lant_tata[lant[x]]]<niv[lant_tata[lant[y]]]) swap(x,y);
                ret=max(ret,arbint_query(1,1,lant_size[lant[x]],1,poz_lant[x],arb[lant[x]]));
                x=lant_tata[lant[x]];
            }
            if(poz_lant[x]>poz_lant[y]) swap(x,y);
            ret=max(ret,arbint_query(1,1,lant_size[lant[x]],poz_lant[x],poz_lant[y],arb[lant[x]]));
            printf("%d\n",ret);
        }
    }
    return 0;
}