bzoj1877[SDOI2009]晨跑【费用流】

最新推荐文章于 2020-02-14 23:03:18 发布

BPM136

最新推荐文章于 2020-02-14 23:03:18 发布

阅读量453

点赞数

分类专栏：网络流文章标签：网络流费用流

本文链接：https://blog.csdn.net/BPM136/article/details/50389503

版权

网络流专栏收录该内容

12 篇文章 0 订阅

订阅专栏

我们发现每个点只能经过一次，于是点权转边权，拆点费用流。

一直不能A，调大了下大小突然A了= =这数据范围不科学

友情提示：读入优化反而会变慢= =

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<bitset>
#define LL long long
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define efo(i,x) for(int i=last[x];i!=0;i=e[i].next)
using namespace std;
inline LL read()
{
	LL d=0,f=1;char s=getchar();
	while(s<'0'||s>'9'){if(s=='-')f=-1;s=getchar();}
	while(s>='0'&&s<='9'){d=d*10+s-'0';s=getchar();}
	return d*f;
}
#define N 255
#define M 50005
#define inf 0x7fffffff
struct edge
{
	int y,w,f,next;
}e[M*2];
int last[N*2],ne=1;
int n,m;
int ans=0,anss=0;

void add(int x,int y,int f,int w)
{
	e[++ne].y=y;e[ne].f=f;e[ne].w=w;e[ne].next=last[x];last[x]=ne;
}

void add2(int u,int v,int f,int c)
{
	add(u,v,f,c);add(v,u,0,-c);
}

int q[N*20],dis[N*2],pre[N*2],lo[N*2];
bitset<N*2>inq;
bool spfa(int s,int tt)
{
	inq.reset();
	fo(i,0,N*2-1)dis[i]=inf;
	memset(lo,0,sizeof(lo));
	memset(pre,0,sizeof(pre));
	int h=0,t=1;
	q[1]=s,dis[s]=0,inq[s]=1;
	while(h<t)
	{
		int now=q[++h];
		efo(i,now)
		if(e[i].f&&dis[now]+e[i].w<dis[e[i].y])
		{
			dis[e[i].y]=dis[now]+e[i].w;
			pre[e[i].y]=now;
			lo[e[i].y]=i;
			if(inq[e[i].y]==0)
			{
				inq[e[i].y]=1;
				q[++t]=e[i].y;
			}
		}
		inq[now]=0;
	}
	if(dis[tt]==inf)return 0;
	else return 1;
}

void mcf(int s,int t)
{
	while(spfa(s,t))
	{
		int x=t,mi=inf;
		while(pre[x]!=0)
		{
			mi=min(e[lo[x]].f,mi);
			x=pre[x];
		}
		x=t;
		while(pre[x]!=0)
		{
			e[lo[x]].f-=mi;
			e[lo[x]^1].f+=mi;
			x=pre[x];
		}
		ans+=dis[t]*mi;
		anss+=mi;
	}
}
	
void init()
{
	n=read(),m=read();
	fo(i,1,m)
	{
		int x=read(),y=read(),c=read();
		add2(x+n,y,1,c);
	}
	fo(i,2,n-1)add2(i,i+n,1,0);
	add2(1,n+1,inf,0);add2(n,n+n,inf,0);
}

void ini()
{
	scanf("%d%d",&n,&m);
	fo(i,1,m)
	{
		int x,y,c;
		scanf("%d%d%d",&x,&y,&c);
		add2(x+n,y,1,c);
	}
	fo(i,2,n-1)add2(i,i+n,1,0);
	add2(1,n+1,inf,0);add2(n,n+n,inf,0);
}

int main()
{
	ini();
	mcf(1,2*n);
	printf("%d %d",anss,ans);
	return 0;
}