condevs 1237 餐巾计划问题（费用流）

最新推荐文章于 2018-07-06 19:51:01 发布

wwyx2001

最新推荐文章于 2018-07-06 19:51:01 发布

阅读量349

点赞数 1

分类专栏：网络流

本文链接：https://blog.csdn.net/Blue_CuSO4/article/details/56481745

版权

网络流专栏收录该内容

38 篇文章 0 订阅

订阅专栏

题目：

我是超链接

题解：

最小费用最大流

建图

1.把每天分为二分图两个集合中的顶点Xi,Yi，建立附加源S汇T。

2.从S向每个Xi连一条容量为ri，费用为0的有向边。//每一天创造的脏餐巾只有这么些嘛

3.从每个Yi向T连一条容量为ri，费用为0的有向边。//新餐巾到货喽

4.从S向每个Yi连一条容量为无穷大，费用为p的有向边。//我可以买无数条新餐巾

5.从每个Xi向Xi+1(i+1<=N)连一条容量为无穷大，费用为0的有向边。//我的脏餐巾拖到了明天

6.从每个Xi向Yi+m(i+m<=N)连一条容量为无穷大，费用为f的有向边。//洗法1

7.从每个Xi向Yi+n(i+n<=N)连一条容量为无穷大，费用为s的有向边。//洗法2

该二分图：脏餐巾（xi）和新餐巾（yi）

代码：

#include <cstdio>
#include <iostream>
#include <cstring>
#include <queue>
#define N 30000
using namespace std;
const int INF=1e9;
int tot=-1,NN,p,m,f,n,s;
int next[N*2+5],point[N*2+5],remind[N*2+5],maxflow=0,mincost=0,v[N*2+5],last[N*2+5],dis[N*2+5],c[N*2+5],r[N*2+5];
bool vis[N*2+5];
void addline(int x,int y,int cap,int cost)
{
	++tot; next[tot]=point[x]; point[x]=tot; v[tot]=y; remind[tot]=cap; c[tot]=cost;
	++tot; next[tot]=point[y]; point[y]=tot; v[tot]=x; remind[tot]=0; c[tot]=-cost;
} 
int addflow(int s,int t)
{
	int ans=INF,now=t;
	while (now!=s)
	{
		ans=min(ans,remind[last[now]]);
		now=v[last[now]^1];
	}
	now=t;
	while (now!=s)
	{
		remind[last[now]]-=ans;
		remind[last[now]^1]+=ans;
		now=v[last[now]^1];
	}
	return ans;
}
bool bfs(int s,int t)
{
	queue <int> q;
	memset(dis,0x7f,sizeof(dis));
	memset(vis,0,sizeof(vis));	
	dis[s]=0;
	vis[s]=true;
	q.push(s);
	while (!q.empty())
	{
		int now=q.front(); q.pop();
		vis[now]=false;
		for (int i=point[now];i!=-1;i=next[i])
		  if (dis[v[i]]>dis[now]+c[i] && remind[i]){
            dis[v[i]]=dis[now]+c[i];
            last[v[i]]=i;
            if (!vis[v[i]]){
                vis[v[i]]=true;
                q.push(v[i]);
            }
          }
    }

    if (dis[t]>INF) return false;
    int flow=addflow(s,t);
    maxflow+=flow;
    mincost+=flow*dis[t];
    return true;
}
void majoy(int s,int t){while (bfs(s,t));}
int main()
{
	memset(point,-1,sizeof(point));
	memset(next,-1,sizeof(next));
	scanf("%d%d%d%d%d%d",&NN,&p,&m,&f,&n,&s);
	for (int i=1;i<=NN;i++)
	{
		scanf("%d",&r[i]);
		addline(0,i,r[i],0);
		addline(i+NN,NN*2+1,r[i],0);
		addline(0,i+NN,INF,p);
		if (i!=NN) addline(i,i+1,INF,0);
		if (i+m<=NN) addline(i,i+NN+m,INF,f);
		if (i+n<=NN) addline(i,i+NN+n,INF,s); 
	}
	majoy(0,NN*2+1);
	printf("%d",mincost); 
}