【UTR #1】ydc的大树

【UTR #1】ydc的大树

全网唯一一篇题解我看不懂

所以说一下我的O(nlogn)做法:

以1号点为根节点

一个黑点如果有多个相邻的节点出去都能找到最远的黑点,那么这个黑点就是无敌的

所以考虑每个黑点x的最远距离和最远点是否仅在一个“方向”

然后这个方向的一些连续白点割掉可以使得x不高兴

1.如果都在一个方向,假设是x的子树,那就是这个子树最远黑点们的lca到x路径上的任意白点割掉,都可以使得x不高兴

2.如果都在往父亲的方向,找到最浅的点p,使得每个最远黑点到x的路径都经过p,p到x的路径上的任意白点割掉,都可以使得x不高兴

树形DP即可。

struct,记录最远距离、最远的方向个数、决策位置(1的lca或者是2的p)

转移较麻烦

 

树上差分打标记即可。

求lca,所以O(nlogn)

 

写了四个dfs。。。

#include<bits/stdc++.h>
#define reg register int
#define il inline
#define fi first
#define se second
#define mk(a,b) make_pair(a,b)
#define numb (ch^'0')
#define pb push_back
#define solid const auto &
#define enter cout<<endl
#define pii pair<int,int>
using namespace std;
typedef long long ll;
template<class T>il void rd(T &x){
    char ch;x=0;bool fl=false;while(!isdigit(ch=getchar()))(ch=='-')&&(fl=true);
    for(x=numb;isdigit(ch=getchar());x=x*10+numb);(fl==true)&&(x=-x);}
template<class T>il void output(T x){if(x/10)output(x/10);putchar(x%10+'0');}
template<class T>il void ot(T x){if(x<0) putchar('-'),x=-x;output(x);putchar(' ');}
template<class T>il void prt(T a[],int st,int nd){for(reg i=st;i<=nd;++i) ot(a[i]);putchar('\n');}
namespace Modulo{
const int mod=998244353;
int ad(int x,int y){return (x+y)>=mod?x+y-mod:x+y;}
void inc(int &x,int y){x=ad(x,y);}
int mul(int x,int y){return (ll)x*y%mod;}
void inc2(int &x,int y){x=mul(x,y);}
int qm(int x,int y=mod-2){int ret=1;while(y){if(y&1) ret=mul(x,ret);x=mul(x,x);y>>=1;}return ret;}
}
//using namespace Modulo;
namespace Miracle{
const int N=1e5+5;
const int inf=0x3f3f3f3f;
int n,m;
int b[N];
struct node{
    int nxt,to;
    int val;
}e[2*N];
int hd[N],cnt;
void add(int x,int y,int z){
    e[++cnt].nxt=hd[x];
    e[cnt].to=y;e[cnt].val=z;
    hd[x]=cnt;
}
struct po{
    int mx,cnt,pos;
    po(){mx=-inf,cnt=1,pos=0;}//warning!! -inf
    po(int v,int c,int p){mx=v;cnt=c;pos=p;}
    po friend operator +(po a,po b){
        if(a.mx==b.mx) return po(a.mx,a.cnt+b.cnt,a.pos);
        else if(a.mx>b.mx) return a;
        else return b;
    }
}pr[N],bc[N],f[N],g[N];
int sta[N],top;
int fa[N][17];
int dep[N],vf[N];
void pre(int x){
    dep[x]=dep[fa[x][0]]+1;
    for(reg i=hd[x];i;i=e[i].nxt){
        int y=e[i].to;
        if(y==fa[x][0]) continue;
        fa[y][0]=x;
        pre(y);
    }
}
void dfs(int x){
    int st=top;
    if(b[x]) f[x]=po(0,1,x);
    for(reg i=hd[x];i;i=e[i].nxt){
        int y=e[i].to;
        if(y==fa[x][0]) continue;
        vf[y]=e[i].val;
        dfs(y);
        sta[++top]=y;
        pr[y]=f[x];
        f[x]=f[x]+po(f[y].mx+e[i].val,1,f[y].pos);
    }
    if(f[x].mx<0){
        f[x].cnt=1;f[x].pos=0;
    }else{
        if(f[x].cnt>1) f[x].cnt=1,f[x].pos=x;
    }
    po now;
    while(top!=st){
        bc[sta[top]]=now;
        now=now+po(f[sta[top]].mx+vf[sta[top]],1,f[sta[top]].pos);
        --top;
    }
}
void gf(int x){
    if(fa[x][0]){
        int pa=fa[x][0];
        g[x]=g[pa]+pr[x]+bc[x];
        g[x].mx+=vf[x];
        if(g[x].mx<0){
            g[x].cnt=1;g[x].pos=0;
        }else{
            if(g[x].cnt>1) g[x].cnt=1,g[x].pos=pa;
        }
    }
    for(reg i=hd[x];i;i=e[i].nxt){
        int y=e[i].to;
        if(y==fa[x][0]) continue;
        gf(y);
    }
}
int lca(int x,int y){
    if(dep[x]<dep[y]) swap(x,y);
    for(reg j=16;j>=0;--j){
        if(dep[fa[x][j]]>=dep[y]) x=fa[x][j];
    }
    if(x==y) return x;
    for(reg j=16;j>=0;--j){
        if(fa[x][j]!=fa[y][j]) x=fa[x][j],y=fa[y][j];
    }
    return fa[x][0];
}
int tag[N];
int ans,tot;
void fin(int x){
    for(reg i=hd[x];i;i=e[i].nxt){
        int y=e[i].to;
        if(y==fa[x][0]) continue;
        fin(y);
        tag[x]+=tag[y];
    }
    if(!b[x]){
        if(tag[x]>ans) {
            ans=tag[x];tot=1;
        }else if(tag[x]==ans){
            ++tot;
        }
    }
}
int main(){
    rd(n);rd(m);
    for(reg i=1;i<=m;++i) {
        int x;rd(x);b[x]=1;
    }
    int x,y,z;
    for(reg i=1;i<n;++i){
        rd(x);rd(y);rd(z);
        add(x,y,z);add(y,x,z);
    }
    pre(1);
    dfs(1);
    gf(1);
    for(reg j=1;j<=16;++j){
        for(reg i=1;i<=n;++i){
            fa[i][j]=fa[fa[i][j-1]][j-1];
            // cout<<" fa "<<i<<" "<<j<<" : "<<fa[i][j]<<endl;
        }
    }
    for(reg i=1;i<=n;++i){
        if(b[i]){
            int x=i;
            po now=f[x]+g[x];
            if(now.cnt==1){
                // cout<<" tag? "<<x<<" "<<now.pos<<endl;
                int anc=lca(x,now.pos);
                // cout<<" anc "<<anc<<endl;
                ++tag[x];++tag[now.pos];
                --tag[anc];--tag[fa[anc][0]];
            }
        }
    }
    ans=-inf;
    fin(1);
    ot(ans);ot(tot);
    return 0;
}

}
signed main(){
    Miracle::main();
    return 0;
}

/*
   Author: *Miracle*
*/

 

posted @ 2019-06-03 21:55  *Miracle*  阅读(315)  评论(0编辑  收藏  举报