洛谷3171 网络吞吐量（网络流）

t开成n结果cur赋值的时候也只赋值到t令人智熄

【题目分析】

好吧我承认这个错误真的呵呵。。。。。。。。

题目有那~~~~~么长，然后画画图这道题就基本看出正解了，再一看数据范围，n<=500简直良心，好了，网络流没得跑了。

因为按最短路进行传递，所以网络流的建图肯定是在最短路的基础上，所以先进行一次SPFA。

考虑一条路如果加入网络流的图，那么这条路一定是在最短路上，dfs一次即可。

然后考虑拆点限制流量（一开始sb的写成了边权结果还跑过了样例然后一交立马咕咕），最后跑最大流即可。

PS：此题还有一个坑点，就是INF要设的很大很大，否则咕咕。

【代码~】

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const LL MAXN=1e3+10;
const LL MAXM=4e5+10;
const LL INF=0x3f3f3f3f3f;

LL n,m,cnt,s,t;
LL head[MAXN],cur[MAXN],depth[MAXN],val[MAXN];
LL nxt[MAXM],to[MAXM],w[MAXM];
LL cnt1,head1[MAXN];
LL nxt1[MAXM],to1[MAXM],w1[MAXM];
LL dis[MAXN],vis[MAXN];

LL Read()
{
	LL i=0,f=1;
	char c;
	for(c=getchar();(c>'9'||c<'0')&&c!='-';c=getchar());
	if(c=='-')
	  f=-1,c=getchar();
	for(;c>='0'&&c<='9';c=getchar())
	  i=(i<<3)+(i<<1)+c-'0';
	return i*f;
}

void sc(LL x)
{
	if(x>=10)
	  sc(x/10);
	putchar(x%10+48);
}

void Add1(LL x,LL y,LL z)
{
	cnt1++;
	nxt1[cnt1]=head1[x];
	head1[x]=cnt1;
	to1[cnt1]=y;
	w1[cnt1]=z;
}

void add1(LL x,LL y,LL z)
{
	Add1(x,y,z);
	Add1(y,x,z);
}

void SPFA()
{
	queue<LL> q;
	memset(dis,0x3f3f3f3f,sizeof(dis));
	dis[s]=0;
	q.push(s);
	while(!q.empty())
	{
		LL u=q.front();
		q.pop();
		vis[u]=0;
		for(LL i=head1[u];i!=-1;i=nxt1[i])
		{
			LL v=to1[i];
			if(dis[v]>dis[u]+w1[i])
			{
				dis[v]=dis[u]+w1[i];
				if(!vis[v])
				{
					vis[v]=1;
					q.push(v);
				}
			}
		}
	}
}

void Add(LL x,LL y,LL z)
{
	nxt[cnt]=head[x];
	head[x]=cnt;
	to[cnt]=y;
	w[cnt]=z;
	cnt++;
}

void add(LL x,LL y,LL z)
{
	Add(x,y,z);
	Add(y,x,0);
}

bool bfs()
{
	queue<LL> q;
	memset(depth,0,sizeof(depth));
	depth[s]=1;
	q.push(s);
	while(!q.empty())
	{
		LL u=q.front();
		q.pop();
		for(LL i=head[u];i!=-1;i=nxt[i])
		{
			LL v=to[i];
			if(!depth[v]&&w[i])
			{
				depth[v]=depth[u]+1;
				q.push(v);
			}
		}
	}
	if(depth[t]==0)
	  return false;
	return true;
}

LL dfs(LL u,LL dist)
{
	if(u==t)
	  return dist;
	for(LL &i=cur[u];i!=-1;i=nxt[i])
	{
		LL v=to[i];
		if(depth[v]==depth[u]+1&&w[i])
		{
			LL di=dfs(v,min(dist,w[i]));
			if(di>0)
			{
				w[i]-=di;
				w[i^1]+=di;
				return di;
			}
		}
	}
	return 0;
}

LL dinic()
{
	LL ans=0;
	while(bfs())
	{
		for(LL i=s;i<=t+n;++i)
		  cur[i]=head[i];
		while(LL d=dfs(s,INF))
		  ans+=d;
	}
	return ans;
}

void buildgraph(LL u,LL fa)
{
	for(LL i=head1[u];i!=-1;i=nxt1[i])
	{
		LL v=to1[i];
		if(v==fa)
		  continue;
		if(dis[v]==dis[u]+w1[i])
		{
			add(u+n,v,INF);
			buildgraph(v,u);
		}
	}
}

int main()
{
	memset(head1,-1,sizeof(head1));
	memset(head,-1,sizeof(head));
	n=Read(),m=Read();
	s=1,t=n;
	for(LL i=1;i<=m;++i)
	{
		LL x=Read(),y=Read(),z=Read();
		add1(x,y,z);
	}
	for(LL i=1;i<=n;++i)
	  val[i]=Read();
	SPFA();
	buildgraph(1,-1);
	add(1,n+1,INF);
	for(LL i=2;i<n;++i)
	  add(i,i+n,val[i]);
	sc(dinic());
	return 0;
}

 

转载于:https://www.cnblogs.com/Ishtar/p/10010755.html