【NOIP 校内模拟】T3 忘了是啥名字了(dfs序+树状数组)

【NOIP 校内模拟】T3 忘了是啥名字了(dfs序+树状数组)

对于当前新加入的一条路径 他产生的贡献分为两种

1.另一条路径的LCA在当前路径上
2.当前路径的LCA在另一条上

对于情况1:

可以维护当前点到根节点有多少个LCA,查询只需查询u,v,-2*lca(u,v),修改需要对lca的子树+1

对于情况2:

显然的树上差分,查询就是lca子树的前缀和,修改u++,v++,lca-2

即开两个树状数组,一个支持单点查询+区间修改,一个支持单点修改+区间查询,不嫌麻烦的话可以尝试线段树。

需要开栈,某OJ栈空间感人。

#include<bits/stdc++.h>
#define N 1000005
#define M 1000005
#define ll long long
using namespace std;
template<class T>
inline void read(T &x)
{
    x=0; int f=1;
    static char ch=getchar();
    while((!isdigit(ch))&&ch!='-')  ch=getchar();
    if(ch=='-') f=-1,ch=getchar();
    while(isdigit(ch))  x=x*10+ch-'0',ch=getchar();
    x*=f;
}
//1e6
struct Edge
{
    int to,next;
}edge[2*N];
int n,m,tot,first[N];
inline void addedge(int x,int y)
{
    tot++;
    edge[tot].to=y; edge[tot].next=first[x]; first[x]=tot;
}
int up[N][27],depth[N],st[N],sign,ed[N];
ll con[N];
void dfs(int now,int fa)
{
    up[now][0]=fa;
    depth[now]=depth[fa]+1;
    st[now]=++sign;
    for(int i=1;i<=25;i++)  up[now][i]=up[up[now][i-1]][i-1];
    for(int u=first[now];u;u=edge[u].next)
    {
        int vis=edge[u].to;
        if(vis==fa) continue;
        dfs(vis,now);
    }
    ed[now]=sign;
}
inline int getlca(int x,int y)
{
    if(depth[x]<depth[y])   swap(x,y);
    for(int i=25;i>=0;i--) if(depth[up[x][i]]>=depth[y]) x=up[x][i];
    if(x==y)    return x;
    for(int i=25;i>=0;i--) if(up[x][i]!=up[y][i])   x=up[x][i],y=up[y][i];
    return up[x][0];
}
inline int lowbit(int x)
{
    return x&(-x);
}
struct BIT
{
    int n;
    ll tree[N];
    inline void getn(int x)
    {
        n=x;
    }
    inline void update(int x,ll del)
    {
        for(int i=x;i<=n;i+=lowbit(i))  tree[i]+=del;
    }
    inline ll query(int x)
    {
        ll ans=0;
        for(int i=x;i;i-=lowbit(i)) ans+=tree[i];
        return ans;
    }
}bit1,bit2; //区间加单点查  单点加区间查    其实就是差分,普通 bit
int main()
{
    ll size=40<<20;//40M
    __asm__ ("movq %0,%%rsp\n"::"r"((char*)malloc(size)+size));//提交用这个 
    read(n),read(m);
    for(register int i=1;i<n;i++)
    {
        int x,y;
        read(x),read(y);
        addedge(x,y); addedge(y,x);
    }
    dfs(1,0);
    bit1.getn(n); bit2.getn(n);
    ll ans=0;
    //需要分两种情况讨论:其他的lca在这条路径上  自己的lca在其他路径上 
    for(int i=1,u,v,lca;i<=m;i++)
    {
        read(u); read(v); lca=getlca(u,v);
        ans=ans+bit1.query(st[u])+bit1.query(st[v])-2*bit1.query(st[lca]);
        ans=ans+bit2.query(ed[lca])-bit2.query(st[lca]-1);
        ans=ans+con[lca];
        con[lca]++;
        bit1.update(st[lca],1); bit1.update(ed[lca]+1,-1);
        bit2.update(st[u],1); bit2.update(st[v],1); bit2.update(st[lca],-2);
    }
    cout<<ans;
    exit(0);
}

转载于:https://www.cnblogs.com/Patrickpwq/articles/9911341.html