bzoj 2286: [Sdoi 2011 소모 전] 허수 구조

우선 우 리 는 매번 트 리 DP 를 한 번 씩 하 는 것 을 고려 합 니 다.
분명히 이것 은 TLE 의 의심 할 바 없 는 것 이다.
변경 할 노드 의 개수 가 매우 적 기 때문에 우 리 는 가상 나무 로 그림 을 다시 만 들 고 범 위 를 좁 히 는 것 을 고려 할 수 있다. 이런 나무 에는 최대 m * 2 개의 노드 가 있다.
그리고 매번 허수 수 에 대해 트 리 dp 를 진행 할 때마다 전체적인 복잡 도 는 n * logn 이 고 logn 은 lca 의 복잡 도이 다.
그리고 여러 가지 디 테 일 한 곳 이에 요...

#include 
#include 
#include 
#include 
#include 
#include 
using namespace std;
typedef long long sint;
#define pii pair
#define mp make_pair
#define deg 20
#define maxn 501000
#define inf 0x3f3f3f3f
#define INF (1ll<<40)
int getint()
{
    int res;char c;
    while(c=getchar(),c'9');
    res=c-'0';
    while(c=getchar(),c>='0'&&c<='9')
        res=res*10+c-'0';
    return res;
}
struct node
{
    int u,v,next;
}tree[maxn];
int first[maxn],next[maxn],to[maxn],fa[maxn][deg];
sint dp[maxn];
int indexs,en,en2,pre[maxn],n,dfn[maxn],dep[maxn],num,T;
pii g[maxn];
int vis[maxn],meet[maxn],val[maxn],minn[maxn][deg];
void build(int a,int b,int c)
{
    en++;
    to[en]=b;
    val[en]=c;
    next[en]=first[a];
    first[a]=en;
}
void add(int u,int v)
{
    en2++;
    tree[en2].v=v;
    if(vis[u]!=T)
    {
        vis[u]=T;
        tree[en2].next=0;
    }
    else
    {
        tree[en2].next=pre[u];
    }
    pre[u]=en2;
}
void dfs(int now)
{
    int v;
    indexs++;
    dfn[now]=indexs;
    for(int i=first[now];i;i=next[i])
    {
        v=to[i];
        if(v==fa[now][0]) continue;
        fa[v][0]=now;
        minn[v][0]=val[i];
        for(int j=1;j=0;i--)
    {
        if(dep[fa[u][i]]>=dep[v])
            u=fa[u][i];
    }
    if(u==v)return v;
    for(int i=19;i>=0;i--)
    if(fa[u][i]!=fa[v][i])
    {
        u=fa[u][i];
        v=fa[v][i];
    }
    return fa[u][0];
}
int sta[maxn],top;
inline long long getmin(int x,int y)
{
    int ans=inf;
    if(dep[x]=0;i--)
        if(dep[fa[x][i]]>=dep[y])
        {
            ans=min(ans,minn[x][i]);
            x=fa[x][i];
        }
    if(x==y)return ans;
    for(int i=19;i>=0;i--)
        if(fa[x][i]!=fa[y][i])
        {
            ans=min(ans,min(minn[x][i],minn[y][i]));
            x=fa[x][i],y=fa[y][i];
        }
    if(x!=y)ans=min(ans,min(minn[x][0],minn[y][0]));
    return ans;
}
void dfs2(int x)
{
    int v;
    dp[x]=0;
    if(vis[x]==T)
    for(int i=pre[x];i;i=tree[i].next)
    {
        dfs2(v=tree[i].v);
        dp[x]+=min(getmin(v,x),meet[v]==T?inf:dp[v]);
    }
}
void virtree()
{
    sta[top=1]=1;
    for(int i=1;i<=num;i++)
    {
        int j=g[i].second;
        int anc=lca(sta[top],j);
        for(;dep[anc]