bzoj 4390（LCA+差分）

传送门
题解：
树上差分：对于修改操作的每一对(x,y)，x,y的sum都加1，lca(x,y)和fa[lca(x,y]的sum都-1。虽然确实是一道水题，但是这种树上差分以前还真的只是听过，没有亲手尝试，要是现在不写写，万一考试写挂了怎办orz。。。

其它解法：
1.DFS序+线段树/树状数组，最后暴力查询每个点（或者更优的，本蒟蒻没想到的统计答案方式）
2.树链剖分，复杂度多一个log(n)，代码又长，直接淘汰。。。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int MAXN=5e4+4;
int n,q,mx=0;
int head[MAXN],edge=0;
struct EDGE {
    int v,nxt;
}e[MAXN<<1];
int f[18][MAXN],sum[MAXN],dis[MAXN];
inline int read() {
    int x=0;char c=getchar();
    while (c<'0'||c>'9') c=getchar();
    while (c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
    return x;
}
inline void adde(int u,int v) {
    e[edge].nxt=head[u],e[edge].v=v,head[u]=edge++;
    e[edge].nxt=head[v],e[edge].v=u,head[v]=edge++;
}
void dfs(int p,int fa) {
    for (int i=head[p];~i;i=e[i].nxt) {
        int v=e[i].v;
        if (v^fa) {
            dis[v]=dis[p]+1;
            f[0][v]=p;
            dfs(v,p);
        }
    }
}
inline void da() {
    for (int j=1;(1<<j)<=n;++j)
        for (int i=1;i<=n;++i)
            f[j][i]=f[j-1][f[j-1][i]];
}
inline int lca(int x,int y) {
    if (dis[x]<dis[y]) x^=y^=x^=y;
    int t=dis[x]-dis[y];
    for (int i=0;i<=17;++i)
        if (t&(1<<i)) x=f[i][x];
    if (x==y) return x;
    for (int i=17;~i;--i)
        if (f[i][x]^f[i][y]) x=f[i][x],y=f[i][y];
    return f[0][x];
}
inline void cal(int p,int fa) {
    for (int i=head[p];~i;i=e[i].nxt) {
        int v=e[i].v;
        if (v^fa) {
            cal(v,p);
            sum[p]+=sum[v];
        }
    }
    mx=max(mx,sum[p]);
}
int main() {
//  freopen("bzoj 4390.in","r",stdin);
    memset(sum,0,sizeof(sum));
    memset(head,-1,sizeof(head));
    n=read(),q=read();
    for (int i=1;i<n;++i) {
        int u=read(),v=read();
        adde(u,v);
    }
    dis[1]=0;
    dfs(1,0);
    da();//init the LCA array
    while (q--) {
        int u=read(),v=read();
        int fa=lca(u,v),fe=f[0][fa];
        ++sum[u],++sum[v],--sum[fa],--sum[fe];
    }
    cal(1,0);
    printf("%d\n",mx);
    return 0;
}