[SCOI2007]修车

一、题目

点此看题

二、解法

这道题和普通的费用流略有不同，一个工人可以匹配若干辆车，观察数据范围：2<=M<=9,1<=N<=60，发现$m$出奇的小，可以考虑拆点，把一个工人拆成$n$个，表示在不同时间段修车的不同工人，现在考虑花费的计算，可以考虑费用提前计算的思想，也就是算一个人的时间是加上后面会等待的人的时间，那如果$i$车交给$j$工人的$k$状态洗，那花费就是$T[i][j]\times k$，所以就得到了我们建图的方式：

• 源点连所有车，花费为$0$，容量为$1$。
• 车连所有工人的所有状态，花费为$T\times k$，容量为$1$。
• 所有工人的所有状态连汇点，花费为$0$，容量为$1$。

算出来时间总和之后除$n$就得到了平均时间qwq。

#include <cstdio>
#include <cstring>
#include <queue>
#define inf 0x3f3f3f3f
const int MAXN = 705;
using namespace std;
int read()
{
    int num=0,flag=1;char c;
    while((c=getchar())<'0'||c>'9')if(c=='-')flag=-1;
    while(c>='0'&&c<='9')num=(num<<3)+(num<<1)+(c^48),c=getchar();
    return num*flag;
}
int n,m,tot=1,S,T,f[MAXN],a[65][15];
int dis[MAXN],pre[MAXN],lst[MAXN],flow[MAXN];
struct edge
{
	int v,f,c,next;
}e[MAXN*MAXN];
struct node
{
	int u,c;
	bool operator < (const node &R) const
	{
		return c>R.c;
	}
};
priority_queue<node> q;
void add_edge(int u,int v,int c,int fl)
{
	e[++tot]=edge{v,fl,c,f[u]},f[u]=tot;
	e[++tot]=edge{u,0,-c,f[v]},f[v]=tot;
}
int Abs(int x)
{
	return x>0?x:-x;
}
int mht(int x1,int y1,int x2,int y2)
{
	return Abs(x1-x2)+Abs(y1-y2);
}
bool bfs()
{
	memset(dis,0x3f,sizeof dis);
	flow[S]=inf;dis[S]=0;pre[S]=-1;
	q.push(node{S,0});
	while(!q.empty())
	{
		int u=q.top().u,t=q.top().c;
		q.pop();
		if(dis[u]<t) continue;
		for(int i=f[u];i;i=e[i].next)
		{
			int v=e[i].v,c=e[i].c;
			if(dis[v]>dis[u]+c && e[i].f>0)
			{
				dis[v]=dis[u]+c;
				pre[v]=u;
				lst[v]=i;
				flow[v]=min(flow[u],e[i].f);
				q.push(node{v,dis[v]});
			}
		}
	}
	return dis[T]!=inf;
}
void get()
{
	int res=0,cost=0;
	while(bfs())
	{
		res+=flow[T];
		cost+=dis[T]*flow[T];
		int cur=T;
		while(cur!=S)
		{
			e[lst[cur]].f-=flow[T];
			e[lst[cur]^1].f+=flow[T];
			cur=pre[cur];
		}
	}
	printf("%.2lf\n",1.0*cost/n); 
}
signed main()
{
	m=read();n=read();
	S=0;T=m*n+n+1;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			a[i][j]=read();
	for(int i=1;i<=n;i++)
		add_edge(S,i,0,1);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
		{
			for(int k=1;k<=n;k++)
				add_edge(k,(i-1)*m+j+n,a[k][j]*i,1);
			add_edge((i-1)*m+j+n,T,0,1);
		}
	get();
}