Acdream 1424 树链剖分

还是做的少，都完全不知道往树链剖分方向想，赛后听了昂神思路才恍然。

首先，石头路连接构成的是一棵树，然后，对于泥路，可以看作是在它的期末位置所在链上覆盖，被覆盖到两次的石头路都是破坏后不能达到效果的，而对于只被覆盖过1次的石头路则是可以的，这条路被断后，剩下一个要断的土路即为覆盖它的那条。所以，只要树链剖分一下，记录被覆盖次数，然后找出只被覆盖了一次的个数即为答案。

#include<iostream>#include<cstdio>#include<cstring>#include<vector>#include<queue>using namespace std;struct node{    int l,r,a;}t[88888];vector<int>g[22222];int f[22222],top[22222],w[22222],s[22222];int dep[22222],son[22222];int m,n,sz;int qx[222222],qy[222222],qnum;void dfs1(int u){    s[u]=1;    son[u]=0;    for(int i=0;i<g[u].size();i++)    {        int v=g[u][i];        if(v==f[u])continue;        dep[v]=dep[u]+1;        f[v]=u;        dfs1(v);        s[u]+=s[v];        if(!son[u]||s[v]>s[son[u]])son[u]=v;    }}void dfs2(int u,int ff){    top[u]=ff;    w[u]=++sz;    if(son[u])dfs2(son[u],ff);    for(int i=0;i<g[u].size();i++)    {        int v=g[u][i];        if(v==f[u]||v==son[u])continue;        dfs2(v,v);    }}void build(int ll,int rr,int rot){    t[rot].l=ll;    t[rot].r=rr;    t[rot].a=0;    if(ll==rr)return;    int mid=(ll+rr)/2;    build(ll,mid,rot<<1);    build(mid+1,rr,rot<<1|1);}void update(int ll,int rr,int rot){    if(ll>rr)return;    if(t[rot].l==ll&&t[rot].r==rr)    {        t[rot].a++;        return;    }    if(t[rot].a)    {        t[rot<<1].a+=t[rot].a;        t[rot<<1|1].a+=t[rot].a;        t[rot].a=0;    }    int mid=(t[rot].l+t[rot].r)/2;    if(rr<=mid)update(ll,rr,rot<<1);    else if(ll>mid)update(ll,rr,rot<<1|1);    else    {        update(ll,mid,rot<<1);        update(mid+1,rr,rot<<1|1);    }}void solve(int u,int v){    int uu=top[u];    int vv=top[v];    while(uu!=vv)    {        if(dep[u]>dep[v])        {            swap(u,v);            swap(uu,vv);        }        update(w[vv],w[v],1);        v=f[vv];        vv=top[v];    }    if(dep[u]>dep[v])swap(u,v);    update(w[son[u]],w[v],1);}int query(int rot){    if(t[rot].l==t[rot].r)    {        if(t[rot].a==1)return 1;        else return 0;    }    if(t[rot].a)    {        t[rot<<1].a+=t[rot].a;        t[rot<<1|1].a+=t[rot].a;        t[rot].a=0;    }    return query(rot<<1)+query(rot<<1|1);}int main(){    int u,v,c;    scanf("%d%d",&n,&m);    sz=0;    qnum=0;    for(int i=0;i<m;i++)    {        scanf("%d%d%d",&u,&v,&c);        if(c==0)        {            qx[qnum]=u;            qy[qnum]=v;            qnum++;        }        else        {            g[u].push_back(v);            g[v].push_back(u);        }    }    dep[1]=0;    f[1]=1;    dfs1(1);    dfs2(1,1);    build(1,sz,1);    for(int i=0;i<qnum;i++)solve(qx[i],qy[i]);    printf("%d\n",query(1));    return 0;}