费用流,也叫作最小费用最大流,是指在普通的网络流图中,每条边的流量都有一个单价,求出一组可行解,使得在满足它是最大流的情况下,总的费用最小
链接:https://www.luogu.org/problem/P3381
最小费用最大流流模板:
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#define MAXN 5005
#define MAXM 50005
using namespace std;
//dis最小花费;pre每个点的前驱;last每个点的所连的前一条边;flow源点到此处的流量
bool vis[MAXN];
int dis[MAXN],pre[MAXN],last[MAXN],flow[MAXN];
int maxflow,mincost;
int to[MAXM<<1],nxt[MAXM<<1],water[MAXM<<1],val[MAXM<<1],head[MAXM<<1];
int tot;
queue <int> que;
void init()
{
memset(head,-1,sizeof(head));
tot=1;//初始化
}
void add(int u,int v,int f,int w)
{
to[++tot]=v;
nxt[tot]=head[u];
head[u]=tot;
water[tot]=f;
val[tot]=w;
}
bool spfa(int s,int t)
{
memset(dis,0x7f,sizeof(dis));
memset(flow,0x7f,sizeof(flow));
memset(vis,0,sizeof(vis));
que.push(s); vis[s]=1; dis[s]=0; pre[t]=-1;
while (!que.empty())
{
int now=que.front();
que.pop();
vis[now]=0;
for (int i=head[now]; i!=-1; i=nxt[i])
{
int v=to[i];
if (water[i]>0 && dis[v]>dis[now]+val[i])
{
dis[v]=dis[now]+val[i];
pre[v]=now;
last[v]=i;
flow[v]=min(flow[now],water[i]);
if (!vis[v])
{
vis[v]=1;
que.push(v);
}
}
}
}
return pre[t]!=-1;
}
void MCMF(int s,int t)
{
maxflow=0;
mincost=0;
while (spfa(s,t))
{
int now=t;
maxflow+=flow[t];
mincost+=flow[t]*dis[t];
while (now!=s)
{
water[last[now]]-=flow[t];
water[last[now]^1]+=flow[t];
now=pre[now];
}
}
}
int main()
{
int n,m,s,t,u,v,f,w;
init();
scanf("%d%d%d%d",&n,&m,&s,&t);
for (int i=1; i<=m; i++)
{
scanf("%d%d%d%d",&u,&v,&f,&w);
add(u,v,f,w);
add(v,u,0,-w);
}
MCMF(s,t);
printf("%d %d",maxflow,mincost);
return 0;
}
https://vjudge.net/problem/POJ-3680
题意:给定n个带权的开区间,i号区间覆盖(ai,bi),权重为wi,现在要从中选取一些区间,要求任意点都不被超过K个区间覆盖,目标是最大化总的权重。
思路:将区间端点看成图中的点,在区间之间和区间内部建边,保证每个区间只能选一次就要在同一个区间的两个端点之间建一条权值为-wi,容量为1的边,剩余的点之间建边容量为inf,权值为0,至于每个点至多选k次,用源点的出边限制一下就好了
ac:
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#define MAXN 1005
#define MAXM 2005
using namespace std;
//dis最小花费;pre每个点的前驱;last每个点的所连的前一条边;flow源点到此处的流量
bool vis[MAXN];
int dis[MAXN],pre[MAXN],last[MAXN],flow[MAXN];
int maxflow,mincost;
int to[MAXM<<1],nxt[MAXM<<1],water[MAXM<<1],val[MAXM<<1],head[MAXM<<1];
int tot;
queue <int> que;
void init()
{
memset(head,-1,sizeof(head));
tot=1;//初始化
}
void add(int u,int v,int f,int w)
{
to[++tot]=v;
nxt[tot]=head[u];
head[u]=tot;
water[tot]=f;
val[tot]=w;
}
bool spfa(int s,int t)
{
memset(dis,0x7f,sizeof(dis));
memset(flow,0x7f,sizeof(flow));
memset(vis,0,sizeof(vis));
que.push(s); vis[s]=1; dis[s]=0; pre[t]=-1;
while (!que.empty())
{
int x=que.front();
que.pop();
vis[x]=0;
for (int i=head[x]; i!=-1; i=nxt[i])
{
int v=to[i];
if (water[i]>0 && dis[v]>dis[x]+val[i])
{
dis[v]=dis[x]+val[i];
pre[v]=x;
last[v]=i;
flow[v]=min(flow[x],water[i]);
if (!vis[v])
{
vis[v]=1;
que.push(v);
}
}
}
}
return pre[t]!=-1;
}
void MCMF(int s,int t)
{
maxflow=0;
mincost=0;
while (spfa(s,t))
{
int now=t;
maxflow+=flow[t];
mincost+=flow[t]*dis[t];
while (now!=s)
{
water[last[now]]-=flow[t];
water[last[now]^1]+=flow[t];
now=pre[now];
}
}
}
int b[MAXN<<2];
int x[MAXN],y[MAXN],f[MAXN];
int main()
{
int q,k,m,s,t;
scanf("%d",&q);
while(q--)
{
init();
scanf("%d%d",&m,&k);
int num=0;
for (int i=1; i<=m; i++)
scanf("%d%d%d",&x[i],&y[i],&f[i]),b[++num]=x[i],b[++num]=y[i];
sort(b+1,b+1+num);
num=unique(b+1,b+1+num)-b-1;
for(int i=1;i<=m;i++)
{
x[i]=lower_bound(b+1,b+1+num,x[i])-b;
y[i]=lower_bound(b+1,b+1+num,y[i])-b;
add(x[i],y[i],1,-f[i]);
add(y[i],x[i],0,f[i]);
}
add(0,1,k,0);
add(1,0,0,0);
add(num,num+1,k,0);
add(num+1,num,0,0);
for(int i=2;i<=num;i++)
{
add(i-1,i,k,0);
add(i,i-1,0,0);
}
s=0,t=num+1;
MCMF(s,t);
printf("%d\n",-mincost);
}
return 0;
}