省选专练[USACO15DEC]最大流Max Flow

最新推荐文章于 2025-09-03 20:00:00 发布

原创最新推荐文章于 2025-09-03 20:00:00 发布 · 232 阅读

0 ·

CC 4.0 BY-SA版权

USACO 专栏收录该内容

14 篇文章

订阅专栏

本文介绍了一个关于树链剖分和线段树的应用问题，通过具体实例讲解了如何利用树链剖分进行LCA跳转，并使用线段树实现区间修改和单点查询，提供了一种有效解决特定类型问题的方法。

白金组的第一年，题目挺水的。

第一这不是网络流~~废话。~~

然后这明显是棵树，接着考虑对其进行树链剖分。

这个题有暴力做法，还记得货物运输吗？可以暴力跳LCA。这个也可以。

然后发现线段树维护区间修改单点查询，树状数组就可以了嘛。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#define N 60000
using namespace std;
struct Front_star{
	int u,v,nxt;
}e[4*N];
int cnt=0;
int first[N]={0};
void add(int u,int v){
	cnt++;
	e[cnt].u=u;
	e[cnt].v=v;
	e[cnt].nxt=first[u];
	first[u]=cnt;
} 
int n,k;
int siz[N]={0};
int dep[N]={0};
int son[N]={0};
int fa[N]={0};
void dfs1(int x,int fat,int depth){
	fa[x]=fat;
	dep[x]=depth;
	siz[x]=1;
	for(int i=first[x];i;i=e[i].nxt){
		int v=e[i].v;
		if(v==fat)
			continue;
		dfs1(v,x,depth+1);
		siz[x]+=siz[v];
		if(!son[x]||siz[son[x]]<siz[v]){
			son[x]=v;
		}
	}
}
int top[N]={0};
int tree[N]={0};
int pre[N]={0};
int tot=0;
void dfs2(int x,int tp){
	top[x]=tp;
	tot++;
	tree[x]=tot;
	pre[tree[x]]=x;
	if(!son[x])
		return;
	dfs2(son[x],tp);
	for(int i=first[x];i;i=e[i].nxt){
		int v=e[i].v;
		if(fa[x]==v)
			continue;
		if(v!=son[x])
			dfs2(v,v);
	}
}
//BIT
int T[4*N]={0};
inline int lowbit(int x){
	return x&(-x);
}
inline void update(int x,int val){
	while(x<=n){
		T[x]+=val;
		x+=lowbit(x);
	}
}
inline int Quary(int x){
	int ans=0;
	while(x){
		ans+=T[x];
		x-=lowbit(x);
	}
	return ans;
}
void up_sum(int u,int v,int val){
	int fu=top[u];
	int fv=top[v];
	while(fu!=fv){
//		cout<<fu<<" "<<fv<<endl;
		if(dep[fu]<dep[fv]){
			swap(fu,fv);
			swap(u,v);
		}
		update(tree[u]-1,-1);
		update(tree[fu],1);
		u=fa[fu];
		fu=top[u];
	}
	if(dep[u]>dep[v]){
		swap(u,v);
	}
	update(tree[u],1);
	update(tree[v]+1,-1);
}
int main(){
	scanf("%d%d",&n,&k);
	for(int i=1;i<n;i++){
		int u,v;
		scanf("%d%d",&u,&v);
		add(u,v);
		add(v,u);
	}
	dfs1(1,1,1);
	dfs2(1,1);
//	for(int i=1;i<=n;i++){
//		cout<<tree[i]<<" "<<top[i]<<endl;
//	}
//	cout<<"hello world"<<endl;
	for(int i=1;i<=k;i++){
		int u,v;
		scanf("%d%d",&u,&v);
		up_sum(u,v,1);
	}
	int ans=0;
//	for(int i=0;i<=n;i++){
//		cout<<T[i]<<" ";
//	}
//	cout<<Quary(1)<<endl;
	for(int i=1;i<=n;i++){
		ans=max(Quary(i),ans);
	}
	cout<<ans;
}