洛谷 P3128 [USACO15DEC]最大流Max Flow (树上差分)

最新推荐文章于 2024-04-21 20:53:31 发布

c_3Po_

最新推荐文章于 2024-04-21 20:53:31 发布

阅读量326

点赞数 2

CC 4.0 BY-SA版权

分类专栏： LCA 树上差分

本文链接：https://blog.youkuaiyun.com/wzazzy/article/details/90182558

LCA 同时被 2 个专栏收录

5 篇文章

订阅专栏

树上差分

1 篇文章

订阅专栏

博客介绍了如何利用树上差分解决USACO15DEC比赛中的最大流问题。给定一棵连接N个隔间的树形网络，以及K条牛奶运输路线，每条路线会增加沿途隔间的运输压力。任务是找出压力最大的隔间的压力值。输入包括N和K，以及树的边和运输路线信息，输出是最大压力。树上差分在此问题中起到关键作用。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

FJ给他的牛棚的N(2≤N≤50,000)个隔间之间安装了N-1根管道，隔间编号从1到N。所有隔间都被管道连通了。

FJ有K(1≤K≤100,000)条运输牛奶的路线，第i条路线从隔间si运输到隔间ti。一条运输路线会给它的两个端点处的隔间以及中间途径的所有隔间带来一个单位的运输压力，你需要计算压力最大的隔间的压力是多少。

输入输出格式

输入格式：

The first line of the input contains NN and KK.

The next N-1N−1 lines each contain two integers xx and yy (x \ne yx≠y) describing a pipe

between stalls xx and yy.

The next KK lines each contain two integers ss and tt describing the endpoint

stalls of a path through which milk is being pumped.

输出格式：

An integer specifying the maximum amount of milk pumped through any stall in the

barn.

输入输出样例

输入样例#1： 复制

输出样例#1： 复制

给你一棵树给你k对点这两点之间的最短路径上的每个点权值加一问最后最大的权值是多少

考虑树上差分，这是个裸题

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int maxn=5e4+5;
std::vector<int> G[maxn];
int dep[maxn],Fa[maxn][20],power[maxn],ans,n,m;
inline void dfs(int u,int fa)
{
	dep[u]=dep[fa]+1;Fa[u][0]=fa;
	for(int i=1;(1<<i)<=dep[u];i++)Fa[u][i]=Fa[Fa[u][i-1]][i-1];
	for(auto &v:G[u]){if(v==fa)continue;dfs(v,u);}
}
inline int lca(int x,int y)
{
	if(dep[x]<dep[y])swap(x,y);
	for(int i=19;i>=0;i--)if((1<<i)<=dep[x]-dep[y])x=Fa[x][i];
	if(x==y)return x;
    for(int i=19;i>=0;i--)if(Fa[x][i]!=Fa[y][i])x=Fa[x][i],y=Fa[y][i];
    return Fa[x][0];
}
inline void solve(int u,int fa)
{
	for(auto &v:G[u])
	{
		if(v==fa)continue;
		solve(v,u);
		power[u]+=power[v];
	}
	ans=max(ans,power[u]);
}
inline int dis(int x,int y){return dep[x]+dep[y]-dep[lca(x,y)];}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<n;i++)
	{
		int u,v;scanf("%d%d",&u,&v);
		G[u].emplace_back(v);
		G[v].emplace_back(u);
	}
	dfs(1,0);
	while(m--)
	{
		int x,y;scanf("%d%d",&x,&y);
		power[x]++;power[y]++;power[lca(x,y)]--;power[Fa[lca(x,y)][0]]--;
	}
	solve(1,0);
	printf("%d\n", ans);
	return 0;
}