Lugou P3381 【模板】最小费用最大流【费用流】【EK+spfa】

题目

思路

这道题是费用流模板题。
那其实就是要让你在最大流的基础上还要使得花费最小，怎么搞？
首先最大流肯定是要找增广路，但是dfs不能找到花费最小的增广路，
所以就可以想到要用最短路找，我用了SPFA。

then?

what the fuck is that? happened?
然后调了3hours+，发现tot没有赋初值，原地爆炸。

tot一定要赋初值为1啊啊啊！！！

tot一定要赋初值为1啊啊啊！！！

tot一定要赋初值为1啊啊啊！！！

关于skw，我不会

代码

#include<iostream>
#include<cstdio>
#include<cmath>
#include<queue>
using namespace std;
const int INF=2147483646;
int n,m,s,t,u,vv,w,c,dis[5010],cur[5010],v[5010];
int tot=1,hd[5010],cost;
queue<int>q;
struct node
{
	int y,next,w,val;
}a[101010];
struct node1
{
	int x,edge;
}pre[10010];
void add(int x,int y,int val,int w)
{
	a[++tot]=(node){y,hd[x],val,w};
	hd[x]=tot;
}
void ycl()
{
	for(register int i=1; i<=n; i++)
	 {
	   dis[i]=INF,v[i]=0;
	   pre[i].edge=0,pre[i].x=0;
	 }
}
bool spfa()
{
	ycl();
	v[s]=0,dis[s]=0;
	q.push(s);
	while(!q.empty())
	 {
	 	int x=q.front();
	 	v[x]=0;
	 	q.pop();
	 	for(register int i=hd[x]; i; i=a[i].next)
	 	 {
	 	 	int yy=a[i].y;
	 	 	if(a[i].w&&dis[yy]>dis[x]+a[i].val)
	 	 	 {
	 	 	    dis[yy]=dis[x]+a[i].val;
	 	 	    pre[yy].edge=i;
	 	 	    pre[yy].x=x;
	 	 	    if(v[yy]==0)
	 	 	     {
	 	 	     	q.push(yy);
	 	 	     	v[yy]=1;
	 	 	     }
	 	 	 }
	 	 }
	 }
	return dis[t]!=INF;
}
int EK()
{
	int ans=0,mi=INF;
	while(spfa())
	 {
	 	mi=INF;
	 	for(register int i=t; i!=s; i=pre[i].x)
	 	   mi=min(mi,a[pre[i].edge].w);
	 	for(register int i=t; i!=s; i=pre[i].x) 
		   a[pre[i].edge].w-=mi,a[pre[i].edge^1].w+=mi;
	    ans+=mi;
	    cost+=mi*dis[t];
	 }
	return ans;
}
int main()
{
	scanf("%d%d%d%d",&n,&m,&s,&t);
	for(register int i=1; i<=m; i++)
	 {
	 	scanf("%d%d%d%d",&u,&vv,&w,&c);
	 	add(u,vv,w,c);
	 	add(vv,u,0,-c);
	 }
	printf("%d ",EK());
	printf("%d",cost);
	return 0;
}