【模板】SPFA增广 最小费用最大流

简单的用SPFA增广进行费用流的求解
与EK求最大流类似，只是此时要求最大流的同时费用最小
所以用SPFA增广，就可以费用尽量小

模板：

#include<cstdio>#include<cstring>#include<algorithm>using namespace std;const int maxn=705,maxe=144005,INF=0x3f3f3f3f;int N,M,S,T,n,tim[65][15],id[65][15],ans=0;int tot=1,son[maxe],nxt[maxe],lnk[maxn],flw[maxe],cap[maxe],w[maxe];void add(int x,int y,int z){    son[++tot]=y;nxt[tot]=lnk[x];lnk[x]=tot;flw[tot]=0;cap[tot]=1;w[tot]=z;    son[++tot]=x;nxt[tot]=lnk[y];lnk[y]=tot;flw[tot]=0;cap[tot]=0;w[tot]=-z;}#define nc getcharinline int red(){    int tot=0,f=1;char ch=nc();    while (ch<'0'||'9'<ch) {if (ch=='-') f=-f;ch=nc();}    while ('0'<=ch&&ch<='9') tot=tot*10+ch-48,ch=nc();    return tot*f;}int que[maxn],dst[maxn],fa[maxn],ed[maxn];bool vis[maxn];bool spfa(){    memset(dst,63,sizeof(dst));    memset(vis,0,sizeof(vis));    int hed=0,til=1;    que[1]=S;dst[S]=0;fa[S]=0;    while (hed!=til){        int x=que[hed=(hed+1)%maxn];        vis[x]=0;        for (int j=lnk[x];j;j=nxt[j])         if (cap[j]>flw[j]&&dst[son[j]]>dst[x]+w[j]){            dst[son[j]]=dst[x]+w[j];            fa[son[j]]=x;ed[son[j]]=j;            if (!vis[son[j]])             vis[son[j]]=1,             que[til=(til+1)%maxn]=son[j];         }    }    if (dst[T]==INF) return 0;    return 1;}int main(){    //do something...    while (spfa()){        int Min=INF;        for (int j=T;j!=S;j=fa[j]) Min=min(Min,cap[ed[j]]-flw[ed[j]]);        for (int j=T;j!=S;j=fa[j]) flw[ed[j]]+=Min,flw[ed[j]^1]-=Min;        ans+=dst[T]*Min;    }    printf("%.2lf",(double)ans/N);    return 0;}