#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<map>
#include<vector>
#include<queue>
#include<set>
#define IL inline
#define x first
#define y second
typedef long long ll;
using namespace std;
const int N=100010,M=400010;
int h[N],e[M],ne[M],idx;
int fa[N][18],depth[N];
int n,m;
int d[N];
int q[N];
int ans;
void add(int a,int b)
{
e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}
void bfs()
{
int tt=0,hh=0;
memset(depth,0x3f,sizeof depth);
depth[0]=0,depth[1]=1;
q[tt++]=1;
while(hh<=tt)
{
int t=q[hh++];
for(int i=h[t];~i;i=ne[i])
{
int j=e[i];
if(depth[j]>depth[t]+1)
{
depth[j]=depth[t]+1;
q[tt++]=j;
fa[j][0]=t;
for(int k=1;k<=17;k++)
fa[j][k]=fa[fa[j][k-1]][k-1];
}
}
}
}
int lca(int a,int b)
{
if(depth[a]<depth[b])
{
swap(a,b);
}
for(int k=17;k>=0;k--)
{
if(depth[fa[a][k]]>=depth[b])
a=fa[a][k];
}
if(a==b) return a;
for(int k=17;k>=0;k--)
if(fa[a][k]!=fa[b][k])
a=fa[a][k],b=fa[b][k];
return fa[a][0];
}
int dfs(int u,int fa)
{
int t=d[u];
for(int i=h[u];~i;i=ne[i])
{
int j=e[i];
if(j==fa) continue;
int tt=dfs(j,u);
if(!tt) ans+=m;
else if(tt==1) ans++;
t+=tt;
}
return t;
}
int main()
{
cin>>n>>m;
memset(h,-1,sizeof h);
for(int i=0;i<n-1;i++)
{
int a,b;
cin>>a>>b;
add(a,b),add(b,a);
}
bfs();
for(int i=0;i<m;i++)
{
int a,b;
cin>>a>>b;
d[a]++,d[b]++;
d[lca(a,b)]-=2;
}
dfs(1,-1);
cout<<ans<<endl;
return 0;
}
树上差分
原创black_hole6 ©著作权
©著作权归作者所有:来自51CTO博客作者black_hole6的原创作品,请联系作者获取转载授权,否则将追究法律责任
上一篇:174. 受欢迎的牛
下一篇:差分约束(专题)

提问和评论都可以,用心的回复会被更多人看到
评论
发布评论
相关文章
-
树上差分总结
需要注意的是,修改的方法和点差分
差分 数组 i++ -
树上差分与LCA
树上差分与LCA可以看b站的这个视频讲的很不错强烈推荐LCA什么是LCALCA 全称 Least Co
深度学习 #include 差分 权值 -
[LCA 树上差分 点差分 板子] Max Flow P
[LCA 树上差分 点差分] Max Flow P题目题目链接思路比松鼠的新家更板的题目,可以
算法 #include 差分 i++ -
树上差分的两种操作
近几年来,树上差分也考地多了,这篇文章中讲一讲树上差分的两种基本操作。一、已知路径求被所有路径覆盖的边
#include 权值 i++ -
python雷达回波数据处理
本文主要针对TI 毫米波雷达的测距、测速、测角的基本方法——3DFFT处理进行简单介绍,并提供MATLAB处理程序,分析3DFFT处理结果。1、测距、测速  
python雷达回波数据处理 3DFFT 测距测速测角 数据 搜索