CodeForces 903G Yet Another Maxflow Problem 题解

最新推荐文章于 2019-09-30 20:09:17 发布

原创最新推荐文章于 2019-09-30 20:09:17 发布 · 212 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

题解_杂专栏收录该内容

304 篇文章

订阅专栏

本文探讨了一种解决二分图中最大流问题的算法，通过转化为最小割问题进行求解。介绍了如何利用线段树维护边的流量，并在修改操作下快速更新最大流，适用于动态变化的图结构。

题目传送门

题目大意： 有一张二分图，左边有 $n$ 个点组成集合 $A$ ，右边有 $n$ 个点组成集合 $B$ ， $A, B$ 满足 $A_i$ 向 $A_{i+1}$ 连边， $B_i$ 向 $B_{i+1}$ 连边，流量都给出，然后还有 $m$ 条 $A_i$ 连向 $B_j$ 的边，流量也给出，最后还有 $q$ 组修改，每次修改 $A_x$ 到 $A_{x+1}$ 的流量，要求每次修改后输出最大流。

题解

因为最大流等于最小割，所以这题可以转化为求最小割。

为了方便，设连接集合 $A$ 里面的点的边的集合为 $E A$ ，连接集合 $B$ 里面的点的边的集合为 $E B$ ，连接集合 $A$ 和集合 $B$ 里面的点的集合为 $E C$ 。

假如使 $EA_x$ 和 $EB_y$ 成为割边，那么要使得 $A_1$ 和 $B_n$ 不连通，那么需要把所有 $EC_{(i,j)}~~(1\leq i <x,y<j\leq n)$ 都给割掉。

设 $a n s (x, y)$ 表示使 $EA_x$ 和 $EB_y$ 成为割边的贡献，那么显然有：
$ans(x,y)=EA_{x}+EB_{y}+\sum_{1\leq i <x,y<j\leq n} EC_{(i,j)}$

对于每一个 $x$ ，我们要找到一个最优的 $y$ ，使得 $a n s (x, y)$ 最小（因为是最小割嘛），然而修改操作只会修改 $EA_x$ ，所以对于每个 $x$ ，其对应的 $y$ 是不会变的。

于是我们可以枚举 $i = 1$ ~ $n$ ，对于每个 $i$ ，将他连向的 $B$ 集合里面的点造成的贡献加入到线段树里面，显然对于一个点 $j$ ，他会对 $1$ ~ $j - 1$ 的所有点造成贡献。然后在所有点里面取一个 $\min$ ，加上 $A_i$ ，就是这个点的贡献了。

然后再对所有点的贡献取一个 $\min$ 就是最终的答案，可以用线段树维护，修改也可以用线段树轻易解决。

看到这里可能会有人产生疑问：对于一对 $(x, y)$ ，可能并不存在让 $EA_x$ 和 $EB_y$ 做割边的方案啊，就算有，这个方案里面也不一定要把所有 $EC_{(i,j)}~~(1\leq i <x,y<j\leq n)$ 都给割掉啊？

对于这样的一对不合要求的 $(x, y)$ ，其实是不会成为最小解的，因为显然它算多了很多，肯定存在另一个 $(x, z)$ 或 $(z, y)$ 比他们小，所以并不会对答案产生影响。

然后还有一个问题，可能最小割根本不需要割 $E A$ 或 $E B$ 里面的边，所以我们考虑往 $E A$ 的末尾放一条长度为 $0$ 的边，往 $E B$ 的开头放一条长度为 $0$ 的边，那么就可以将不割 $E A$ 里面的边和不割 $E B$ 里面的边这两种情况考虑进去。

代码如下：

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define maxn 200010
#define ll long long

int n,m,k;
int a[maxn],b[maxn];
struct nod{int y,z,next;};
nod e[maxn];
int first[maxn];
void buildroad(int x,int y,int z)
{
	static int len=0;
	e[++len]=(nod){y,z,first[x]};
	first[x]=len;
}
ll ans[maxn];
struct node{
	int l,r;
	ll z,lazy;
	node *zuo,*you;
	template<class T>
	node(int x,int y,T *s)
	{
		l=x,r=y;lazy=0;
		if(l<r)
		{
			int mid=l+r>>1;
			zuo=new node(l,mid,s);
			you=new node(mid+1,y,s);
			z=min(zuo->z,you->z);
		}
		else z=s[l],zuo=NULL,you=NULL;
	}
	void pushdown()
	{
		if(lazy!=0)
		{
			z+=lazy;
			if(zuo!=NULL)zuo->lazy+=lazy,you->lazy+=lazy;
			lazy=0;
		}
	}
	void change(int x,int y,int c)
	{
		if(l==x&&r==y)
		{
			lazy+=(ll)c;
			pushdown();
			return;
		}
		pushdown();
		int mid=l+r>>1;
		if(y<=mid)zuo->change(x,y,c),you->pushdown();
		else if(x>=mid)you->change(x,y,c),zuo->pushdown();
		else zuo->change(x,mid,c),you->change(mid+1,y,c);
		z=min(zuo->z,you->z);
	}
};
node *root=NULL;
void work()
{
	root=new node(1,n,b);
	for(int i=1;i<=n;i++)
	{
		for(int j=first[i];j;j=e[j].next)
		root->change(1,e[j].y,e[j].z);
		ans[i]=root->z+a[i];
	}
}

int main()
{
	scanf("%d %d %d",&n,&m,&k);
	for(int i=1;i<n;i++)
	scanf("%d %d",&a[i],&b[i+1]);
	for(int i=1,x,y,z;i<=m;i++)
	scanf("%d %d %d",&x,&y,&z),buildroad(x,y,z);
	work();
	root=new node(1,n,ans);
	printf("%lld\n",root->z);
	for(int i=1,x,y;i<=k;i++)
	{
		scanf("%d %d",&x,&y);
		root->change(x,x,y-a[x]);
		a[x]=y;
		printf("%lld\n",root->z);
	}
}