网络流之Dinic(板子总结)

思路来源

https://blog.csdn.net/qq_37352710/article/details/88037456

Dinic复杂度

Dinic算法的时间复杂度的理论上界是O(N^2*M)(N是结点数,M是边数)。

如果所有边容量均为1,那么时间复杂度是O(min(N^{0.67},M^{0.5}) * M)

对于二分图最大匹配这样的特殊图,时间复杂度是O(N^{0.5} * M)

板子来源

航神+涛神

心得

下午涛神给航神例会讲的dinic板子加了一行代码

本来1996ms的代码瞬间跑到了170ms

卧槽这dinic板子碉堡了我要总结一下

基本可过大部分卡Dinic的板子,如hdu4280

hdu4280实测大概跑7700ms(时限10000ms),

而那道题用ISAP跑的话大概5000多ms

回头再总结个ISAP的板子叭

网络流理解了原理之后当小黑盒使用就好

代码

以hdu6214 Smallest Minimum Cut 为例 

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<map> 
using namespace std;
typedef long long ll;
const int INF=0x3f3f3f3f;
const int maxn=210;
const int maxm=2e3+10;
int level[maxn];
int head[maxn],cnt;
int t,n,m;
int ss,ee;
struct edge{int v,nex;ll w;}e[maxm];
void init()
{
	cnt=0;
	memset(head,-1,sizeof head);
}
void add(int u,int v,ll w)
{
	e[cnt].v=v;
	e[cnt].w=w;
	e[cnt].nex=head[u];
	head[u]=cnt++;
}
void add2(int u,int v,ll w,bool op)//是否为有向图 
{
	add(u,v,w);
	add(v,u,op?0:w);
}
bool bfs(int s,int t)
{
	queue<int>q;
	memset(level,0,sizeof level);
	level[s]=1;
	q.push(s);
	while(!q.empty())
	{
		int x=q.front();
		q.pop();
		if(x==t)return 1;
		for(int u=head[x];~u;u=e[u].nex)
		{
			int v=e[u].v;ll w=e[u].w;
			if(!level[v]&&w)
			{
				level[v]=level[x]+1;
				q.push(v);
			}
		}
	}
	return 0;
}
ll dfs(int u,ll maxf,int t)
{
	if(u==t)return maxf;
	ll ret=0;
	for(int i=head[u];~i;i=e[i].nex)
	{
		int v=e[i].v;ll w=e[i].w;
		if(level[u]+1==level[v]&&w)
		{
			ll MIN=min(maxf-ret,w);
			w=dfs(v,MIN,t);
			e[i].w-=w;
			e[i^1].w+=w;
			ret+=w;
			if(ret==maxf)break;
		}
	}
	if(!ret)level[u]=-1;//优化,防止重搜,说明u这一路不可能有流量了 
	return ret;
}
ll Dinic(int s,int t)
{
	ll ans=0;
	while(bfs(s,t))
	ans+=dfs(s,INF,t);
	return ans;
}
int main()
{ 
    scanf("%d",&t);
    //n个点 m条边 
    while(t--)
    {
    	init();
    	scanf("%d%d",&n,&m);
    	scanf("%d%d",&ss,&ee);
    	ss--,ee--;
    	for(int j=0;j<m;++j) 
    	{
    		int u,v;ll w;
    		scanf("%d%d%lld",&u,&v,&w);
    		u--,v--;
    		w=w*10000+1;
    		add2(u,v,w,1);//op==1代表为有向图 
    	}
    	ll ans=Dinic(ss,ee);
    	printf("%lld\n",ans%10000);
    }
	return 0;
}

 

  • 1
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 3
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Code92007

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值