bzoj 4108: [Wf2015]Catering （有源汇有上下界的费用流）

4108: [Wf2015]Catering

Time Limit: 4 Sec   Memory Limit: 128 MB
Submit: 165   Solved: 58
[ Submit][ Status][ Discuss]

Description

有一家装备出租公司收到了按照时间顺序排列的n个请求.

这家公司有k个搬运工.每个搬运工可以搬着一套装备按时间顺序去满足一些请求.一个搬运工从第i个请求的位置把东西搬到第j个请求的位置需要一些费用.公司的编号是1,请求的编号是2到n+1.所有搬运工必需从公司出发.

求满足所有请求所需的最小搬运费用.

Input

有可能有多组数据(我也不知道).

第一行两个正整数n,k.

接下来n行,第i行有n-i+1个数.第j个数表示把装备从i搬到i+j的花费.

Output

输出一行一个整数表示最小花费.

Sample Input

1 1
10

Sample Output

10

HINT

n, k <= 100;

花费 <= 1,000,000

Source

鸣谢laekov提供译文

[ Submit][ Status][ Discuss]

题解：有源汇有上下界的费用流

建出原图后重建，直接跑最小费用最大流即可。注意要加权限制总的路线数。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<queue>
#include<cstring>
#define N 50000
#define inf 1000000000
using namespace std;
int ans,n,m,tot,val[N];
int point[N],v[N],next[N],remain[N],c[N],dis[N],can[N],last[N];
void add(int x,int y,int z,int k)
{
	tot++; next[tot]=point[x]; point[x]=tot; v[tot]=y; remain[tot]=z; c[tot]=k;
	tot++; next[tot]=point[y]; point[y]=tot; v[tot]=x; remain[tot]=0; c[tot]=-k;
	//if (z)  cout<<x<<" "<<y<<" "<<z<<" "<<k<<endl;
}
int addflow(int s,int t)
{
	int now=t; int ans=inf;
	while (now!=s) {
		//cout<<now<<" ";
		ans=min(ans,remain[last[now]]);
		now=v[last[now]^1];
	}
	//cout<<now<<endl;
	now=t;
	while (now!=s) {
		remain[last[now]]-=ans;
		remain[last[now]^1]+=ans;
		now=v[last[now]^1];
	}
	return ans;
}
bool spfa(int s,int t)
{
	for (int i=0;i<=t;i++) dis[i]=inf,can[i]=0;
	dis[s]=0; can[s]=1;
	queue<int> p; p.push(s);
	while (!p.empty()){
		int now=p.front(); p.pop();
		for (int i=point[now];i!=-1;i=next[i])
		 if (remain[i]&&dis[v[i]]>dis[now]+c[i]){
		 	dis[v[i]]=dis[now]+c[i];
		 	last[v[i]]=i;
		 	if (!can[v[i]]) {
		 		can[v[i]]=1;
		 		p.push(v[i]);
			 }
		 }
		can[now]=0;
	}
	if (dis[t]==inf) return false;
	int mx=addflow(s,t);
	//cout<<dis[t]<<" "<<mx<<endl;
	ans+=mx*dis[t];
	return true;
}
void solve(int s,int t)
{
	while (spfa(s,t));
}
int main() 
{
	freopen("a.in","r",stdin);
	freopen("my.out","w",stdout);
	while (scanf("%d%d",&n,&m)!=EOF) {
		tot=-1;
		memset(point,-1,sizeof(point));
		int t=2*(n+1)+1;
		for (int i=1;i<=n;i++) {
		 int x; scanf("%d",&x);
		 add(1,i+1,1,x); 
	    }
	  	for (int i=1;i<=n;i++)
		 for (int j=i+1;j<=n;j++){
		  int x; scanf("%d",&x);
		  add(i+n+1,j+1,1,x);
	    }
	    for (int i=1;i<=n;i++) 
		 add(i+1,i+n+1,0,0),add(i+n+1,t,1,0);
		int S=t+1; int T=t+2;
	    add(t,0,inf,0); add(0,1,m,0); //add(S,1,1,0); add(0,T,1,0);
	    for (int i=2;i<=n+1;i++)
	     add(S,i+n,1,0),add(i,T,1,0);
	    solve(S,T);
	    printf("%d\n",ans);
	}
}