dsu on tree(无讲解)

CF741D. Arpa’s letter-marked tree and Mehrdad’s Dokhtar-kosh paths

分析:

• 最多有一个字符出现奇数次
• 维护某个状态下深度的最大值，注意是全局深度
• 写成非递归形式方便理解

代码:

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cstdlib>
using namespace std;
#define N 500050
#define M (1<<22)
#define inf 0x3f3f3f3f
int head[N],to[N<<1],nxt[N<<1],cnt,n,val[N<<1];
int mx[M],siz[N],ans[N],son[N];
int dfn[N],idf[N],enp[N];
int nowans,dep[N],len[N];
inline void add(int u,int v,int w) {
    to[++cnt]=v; nxt[cnt]=head[u]; head[u]=cnt; val[cnt]=w;
}
void df1(int x,int y) {
    int i;
    siz[x]=1; len[x]=len[y]+1;
    dfn[x]=++dfn[0]; idf[dfn[0]]=x;
    for(i=head[x];i;i=nxt[i]) if(to[i]!=y) {
        dep[to[i]]=dep[x]^(1<<val[i]);
        df1(to[i],x);
        siz[x]+=siz[to[i]];
        if(siz[to[i]]>siz[son[x]]) son[x]=to[i];
    }
    enp[x]=dfn[0];
}
void df2(int x,int y,int opt) {
    int i;
    for(i=head[x];i;i=nxt[i]) if(to[i]!=y&&to[i]!=son[x]) {
        df2(to[i],x,0); ans[x]=max(ans[x],ans[to[i]]);
    }
    if(son[x]) {
        df2(son[x],x,1);
        ans[x]=max(ans[x],ans[son[x]]);
        nowans=max(nowans,mx[dep[x]]-len[x]);
        for(i=0;i<22;i++) nowans=max(nowans,mx[dep[x]^(1<<i)]-len[x]);
    }
    mx[dep[x]]=max(mx[dep[x]],len[x]);
    for(i=head[x];i;i=nxt[i]) if(to[i]!=y&&to[i]!=son[x]) {
        int t=to[i],j,p,k;
        for(j=dfn[t];j<=enp[t];j++) {
            p=idf[j];
            nowans=max(nowans,mx[dep[p]]+len[p]-len[x]*2);
            for(k=0;k<22;k++) nowans=max(nowans,mx[ dep[p]^(1<<k) ]+len[p]-len[x]*2);
        }
        for(j=dfn[t];j<=enp[t];j++) {
            p=idf[j];
            mx[dep[p]]=max(mx[dep[p]],len[p]);
        }
    }
    ans[x]=max(ans[x],nowans);
    if(!opt) {
        nowans=-inf;
        for(i=dfn[x];i<=enp[x];i++) mx[dep[idf[i]]]=-inf;
    }
}
int main() {
    memset(mx,0xc0,sizeof(mx));
    scanf("%d",&n);
    int i,x; 
    char op[4];
    for(i=2;i<=n;i++) {
        scanf("%d%s",&x,op);
        add(x,i,op[0]-'a'); add(i,x,op[0]-'a');
    }
    df1(1,0);
    df2(1,0,1);
    for(i=1;i<=n;i++) printf("%d ",ans[i]);
}

BZOJ5457: 城市

代码:

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cstdlib>
using namespace std;
#define N 400050
int head[N],to[N<<1],nxt[N<<1],n,m,cnt;
int a[N],b[N],dfn[N],idf[N],fa[N],enp[N],son[N],siz[N],top[N];
#define nc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++)
char buf[100000],*p1,*p2;
inline int rd() {
    int x=0;char c=nc();
    while(c<48)c=nc();
    while(c>47)x=((x+(x<<2))<<1)+(c^48),c=nc();
    return x;
}
inline void add(int u,int v) {
    to[++cnt]=v; nxt[cnt]=head[u]; head[u]=cnt;
}
void df1(int x,int y) {
    int i; fa[x]=y;
    siz[x]=1; dfn[x]=++dfn[0]; idf[dfn[0]]=x;
    for(i=head[x];i;i=nxt[i]) if(to[i]!=y) {
        df1(to[i],x);
        siz[x]+=siz[to[i]];
        if(siz[to[i]]>siz[son[x]]) son[x]=to[i];
    }
    enp[x]=dfn[0];
}
int nowans,ans[N],c[N],tot[N];
void df2(int x) {
    int i,j,t,lim;
    for(i=head[x];i;i=nxt[i]) if(to[i]!=son[x]&&to[i]!=fa[x]) {
        df2(to[i]);
    }
    if(son[x]) df2(son[x]);
    c[a[x]]+=b[x];
    if(c[nowans]<c[a[x]]||(c[nowans]==c[a[x]]&&nowans>a[x])) {
        nowans=a[x];
    }
    for(i=head[x];i;i=nxt[i]) if(to[i]!=son[x]&&to[i]!=fa[x]) {
        t=to[i],lim=enp[t];
        for(j=dfn[t];j<=lim;j++) {
            int p=idf[j];
            c[a[p]]+=b[p];
    if(c[nowans]<c[a[p]]||(c[nowans]==c[a[p]]&&nowans>a[p])) {
        nowans=a[p];
    }
        }
    }
    ans[x]=nowans; tot[x]=c[nowans];
    if(son[fa[x]]!=x) {nowans=0;for(lim=enp[x],i=dfn[x];i<=lim;i++) c[a[idf[i]]]=0;}
}
char pbuf[4000000],*pp=pbuf;
int sta[30],tp;
void write(int x) {
    if(x<0) *pp++='-',x=-x;
    do {sta[++tp]=x%10,x/=10;}while(x);
    while(tp)*pp++=sta[tp--]+'0';
}
int main() {
    n=rd(),m=rd();
    int i,x,y;
    for(i=1;i<n;i++) x=rd(),y=rd(),add(x,y),add(y,x);
    for(i=1;i<=n;i++) a[i]=rd(),b[i]=rd();
    df1(1,0); df2(son[0]=1);
    for(i=1;i<=n;i++) {
        write(ans[i]); *pp++=' '; write(tot[i]); *pp++='\n';
    }
    fwrite(pbuf,1,pp-pbuf,stdout);
}