首页 > 其他分享 >无源汇有上下界可行流

无源汇有上下界可行流

时间:2024-02-20 13:33:34浏览次数:32  
标签:sx return int 流量 下界 无源 TP st 汇有

Describe:

\(n\) 个点,\(m\) 条边,每条边 \(e\) 有一个流量下界 \(\text{lower}(e)\) 和流量上界 \(\text{upper}(e)\),求一种可行方案满足流量守恒的同时满足每条边的限制条件。

Solution:

可以先考虑满足所有边的最低条件,获得一个初始的流量网络。

image

这时各个点是不一定满足流量守恒的。可能会有多余的流量无处可走,也可能会有流量不足以流出。因此可以建立一个源点,补充不足的流量,同时也建立一个汇点,汇入多余的流量。

image

如果能在满足所有边的条件下从源点全部流入汇点,显然就是一个可行流。

因此可以建立一个残量网络,每条边的流量改为 \(\text{upper}(e) - \text{lower}(e)\),跑一个最大流。判断一下最大流是否等于源点的流量即可。每条边的实际流量就是最低流量加残量网络中每条边流出的流量。

Code:

bool _Start;
#include<deque>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
namespace IO
{
	#define TP template<typename T>
	#define TP_ template<typename T,typename ... T_>
	#ifdef DEBUG
	#define gc() (getchar())
	#else
	char buf[1<<20],*p1,*p2;
	#define gc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<20,stdin),p1==p2)?EOF:*p1++)
	#endif
	#ifdef DEBUG
	void pc(const char &c)
	{
		putchar(c);
	}
	#else
	char pbuf[1<<20],*pp=pbuf;
	void pc(const char &c)
	{
		if(pp-pbuf==1<<20)
			fwrite(pbuf,1,1<<20,stdout),pp=pbuf;
		*pp++=c;
	}
	struct IO{~IO(){fwrite(pbuf,1,pp-pbuf,stdout);}}_;
	#endif
	TP void read(T &x)
	{
		x=0;static int f;f=0;static char ch;ch=gc();
		for(;ch<'0'||ch>'9';ch=gc())ch=='-'&&(f=1);
		for(;ch>='0'&&ch<='9';ch=gc())x=(x<<1)+(x<<3)+(ch^48);
		f&&(x=-x);
	}
	TP void write(T x)
	{
		if(x<0)
			pc('-'),x=-x;
		static T sta[35],top;top=0;
		do
			sta[++top]=x%10,x/=10;
		while(x);
		while(top)
			pc(sta[top--]^48);
	}
	TP_ void read(T &x,T_&...y){read(x);read(y...);}
	TP void writeln(const T x){write(x);pc('\n');}
	TP void writesp(const T x){write(x);pc(' ');}
	TP_ void writeln(const T x,const T_ ...y){writesp(x);writeln(y...);}
	TP void debugsp(const T x){fprintf(stderr,"%d ",x);}
	TP void debug(const T x){fprintf(stderr,"%d\n",x);}
	TP_ void debug(const T x,const T_...y){debugsp(x);debug(y...);}
	TP inline T max(const T &a,const T &b){return a>b?a:b;}
	TP_ inline T max(const T &a,const T_&...b){return max(a,max(b...));} 
	TP inline T min(const T &a,const T &b){return a<b?a:b;}
	TP_ inline T min(const T &a,const T_&...b){return min(a,min(b...));}
	TP inline void swap(T &a,T &b){static T t;t=a;a=b;b=t;}
	TP inline T abs(const T &a){return a>0?a:-a;}
	#undef TP
	#undef TP_
}
using namespace IO;
using std::cerr;
using LL=long long;
constexpr int inf=1e9,N=2e2+5,M=1e4+2*N;
struct edge
{
	int y,f,pre;
}a[M<<1];int alen=1,last[N];
void ins(int x,int y,int f)
{
	a[++alen]=edge{y,f,last[x]};
	last[x]=alen;
	a[++alen]=edge{x,0,last[y]};
	last[y]=alen;
}
int n,m,st,ed;
int w[M],in[N];
int h[N],cur[N];
bool pd()
{
	static std::deque<int>q;q.clear();
	q.push_back(st);
	memset(h,0,sizeof(h));h[st]=1;
	while(q.size())
	{
		int x=q.front();q.pop_front();
		for(int k=last[x];k;k=a[k].pre)
		{
			int y=a[k].y;
			if(a[k].f&&!h[y])
			{
				h[y]=h[x]+1;
				q.push_back(y);
			}
		}
	}
	return h[ed];
}
int findflow(int x,int f)
{
	if(x==ed)
		return f;
	int sx=0,sy;
	for(int &k=cur[x];k;k=a[k].pre)
	{
		int y=a[k].y;
		if(a[k].f&&h[y]==h[x]+1)
		{
			sy=findflow(y,min(a[k].f,f-sx));
			a[k].f-=sy;sx+=sy;
			a[k^1].f+=sy;
			if(sx==f)
				return sx;
		}
	}
	if(!sx)
		h[x]=-1;
	return sx;
}
int dinic()
{
	int s=0;
	while(pd())
	{
		memcpy(cur,last,sizeof(cur));
		s+=findflow(st,inf);
	}
	return s;
}
bool _End;
int main()
{
	// fprintf(stderr,"%.2 MBlf\n",(&_End-&_Start)/1048576.0);
	read(n,m);st=n+1;ed=st+1;
	for(int i=1;i<=m;i++)
	{
		int x,y,up,dn;
		read(x,y,dn,up);
		w[i]=dn;ins(x,y,up-dn);
		in[x]-=dn,in[y]+=dn;
	}
	int sum=0;
	for(int i=1;i<=n;i++)
		if(in[i]<0)
			ins(i,ed,-in[i]);
		else if(in[i]>0)
			ins(st,i,in[i]),sum+=in[i];
	if(dinic()==sum)
	{
		puts("YES");
		for(int i=1;i<=m;i++)
			writeln(w[i]+a[i<<1|1].f);
	}
	else
		puts("NO");
	return 0;
}

标签:sx,return,int,流量,下界,无源,TP,st,汇有
From: https://www.cnblogs.com/lofty2007/p/18022884

相关文章

  • 上下界 可行/最大/最小 网络流/费用流(有/无源汇)
    对网络的定义进行扩展,我们可以得出一堆奇奇怪怪的网络。上下界令\(Max_e\)为边\(e\)的流量上界,\(Min_e\)为边\(e\)的流量下界,一条边的流量\(f_e\)要满足\(Min_e\lef_e\leMax_e\),除此之外和普通网络流定义相同,可以发现,普通的网络就是下界为\(0\)的网络。无源汇......
  • 不搜索,无问题。冗余、上下界剪枝
    公众号:编程驿站1.搜索算法本文和大家聊聊搜索算法,计算机解决问题的抽象流程是,先搜索,或完全搜索后得到答案,或边搜索边找答案。所以,对给定的数据集进行搜索是解决问题的前置条件。不搜索,无问题。搜索算法无非就是线性、二分、深度、广度搜索算法。其它的搜索算法的底层逻辑也是建......
  • SYN5305型晶振测试仪、晶振测试设备、晶体测试仪、无源晶振测试常用仪器
    负载电容(Load Capacitance)CL:在电路中跨接晶振两端总的外界有效电容;要等于或接近晶振数据手册给出的数值才能使晶振按照预期工作;比较通常使用的24MHz晶振负载电容是20pF静态电容C0:以石英晶片为介电材料和两个电极之间的电容;很多晶振规格里是<3pF的;SYN5305型晶振测试仪集合有源和......
  • 二级交错指数时间的电路下界
    \(\newcommand{\NP}{\mathsf{NP}}\newcommand{\PP}{\mathsf{P}}\newcommand{\PPoly}{\mathsf{P/_{poly}}}\newcommand{\EXPSPACE}{\mathsf{EXPSPACE}}\newcommand{\SigmaTwo}{\mathsf{\Sigma_2}}\newcommand{\E}{\mathsf{E}}\newcommand{\F......
  • STM32实现无源蜂鸣器唱歌
    记录学习stm32中实现小demo所涉及到的知识点一、蜂鸣器发声原理蜂鸣器分为有源和无源两种。所谓的源,指的是其中内部的振荡源,有源蜂鸣器中的振荡器一般是[[多谐振荡器]],其原理就是模拟电路中RC振荡器的一般原理(放大电路、正反馈、相位差90°、稳压电路),有源蜂鸣器内部的振荡源频......
  • 上下界网络流
    学一次忘一次,搞笑。规定\(s\)和\(t\)为原图的源汇点,\(S\)和\(T\)为新建的虚拟源汇点。无源汇上下界可行流考虑先把每条边的下界流满,然后网络的边权改为\(r-l\)。但这样每个点的流量平衡不能保证,我们建源点\(S\)和汇点\(T\),如果一个点的入量大于出量,就从\(S\)到它......
  • 【玩转鲲鹏 DevKit系列】如何快速迁移无源码应用?
    本文分享自华为云社区《【玩转鲲鹏DevKit系列】如何快速迁移无源码应用?》,作者:华为云社区精选。为了帮助广大用户和开发者快速将无源码应用从x86迁移到鲲鹏,鲲鹏DevKit提供了动态二进制翻译工具ExaGear,它能在运行时将x86二进制指令翻译成鲲鹏二进制指令,使得大部分x86应用无需......
  • U型无源测温标签7535U有哪些特性
    产品名称: 测温传感器产品型号: CW7535U标签产品特性: 无源测温/陶瓷标签/U型底板/ 75mm*35mm/无源芯片1.概述CW7535U标签是一款应用在金属表面的测温标签。该测温标签采用了RFID无源无线传输技术,并在芯片内集成了超低功耗温度传感器,可用来监测接触点温度,是无源无线RFID温度检测......
  • 【雕爷学编程】Arduino动手做(49)---有源和无源蜂鸣器模块5
    37款传感器与模块的提法,在网络上广泛流传,其实Arduino能够兼容的传感器模块肯定是不止37种的。鉴于本人手头积累了一些传感器和执行器模块,依照实践出真知(一定要动手做)的理念,以学习和交流为目的,这里准备逐一动手试试多做实验,不管成功与否,都会记录下来——小小的进步或是搞不掂的问题......
  • 无源测向技术介绍——振幅法测向
    振幅法测向即为对天线接收到的信号幅度信息进行分析处理,以获得信号角度信息的方法,主要分为波束搜索法(最大振幅法)和全向振幅单脉冲测向法。**121.波束搜索法**波束搜索法的原理图如1所示。接收天线从以v的速度在一定搜索角度内顺时针旋转,当天线接收到的信号幅度高......