LG2053/BZOJ1070 「SCOI2007」修车 费用流

问题描述

LG2053

BZOJ1070

---

题解

将(m)个修理工拆为(n imes m)个，将修理工和车辆看做二分图，连出一个完全二分图。

边流量为(1)，费用为时间，费用流即可。

---

(mathrm{Code})

#include<bits/stdc++.h>
using namespace std;

template <typename Tp>
void read(Tp &x){
	x=0;char ch=1;int fh;
	while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
	if(ch=='-'){
		fh=-1;ch=getchar();
	}
	else fh=1;
	while(ch>='0'&&ch<='9'){
		x=(x<<1)+(x<<3)+ch-'0';
		ch=getchar();
	}
	x*=fh;
}

const int maxn=63;
const int maxm=13;
const int INF=0x3f3f3f3f; 

int n,m;
int a[maxn][maxm];

int Head[100000],Next[1000000],to[1000000],tot=1,w[1000000],co[1000000];
int S,T;

void add(int a,int b,int c,int d){
	to[++tot]=b,Next[tot]=Head[a],Head[a]=tot,w[tot]=c,co[tot]=d;
}

bool vis[100000];
int now[1000000],pre[1000000];
int dis[100000];
int ans,flow;

bool spfa() {
    memset(vis,0,sizeof(vis));memset(dis,0x3f,sizeof(dis));
    queue<int>q;q.push(S);vis[S]=1;dis[S]=0;now[S]=INF;
    while(!q.empty()){
        int x=q.front();q.pop();vis[x]=0;
        for (int i=Head[x];i;i=Next[i]){
            int y=to[i],z=w[i],ww=co[i];
            if (!z||dis[y]<=dis[x]+ww) continue;
            dis[y]=dis[x]+ww;
            now[y]=min(now[x],z);
            pre[y]=i;
            if(!vis[y]){
                q.push(y);
                vis[y] = 1;
            }
        }
    }
    return dis[T]!=INF;
}

void upd(){
	flow+=now[T],ans+=dis[T]*now[T];
	int p=T;
	while(p!=S){
		int k=pre[p];
		w[k]-=now[T],w[k xor 1]+=now[T];
		p=to[k xor 1];
	}
}

int main(){
	read(m);read(n);tot=1;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
			read(a[i][j]);
		}
	}
	S=n*(m+1)+1,T=S+1;
	for(int i=1;i<=n*m;i++){
		add(S,i,1,0);add(i,S,0,0);
	}
	for(int i=1;i<=m;i++){
		for(int j=1;j<=n;j++){
			for(int k=1;k<=n;k++){
				int p=(i-1)*n+k;
				add(p,j+n*m,1,a[j][i]*k);add(j+n*m,p,0,-a[j][i]*k);
			}
		}
	}
	for(int i=1;i<=n;i++){
		add(n*m+i,T,1,0);add(T,n*m+i,0,0);
	}
	while(spfa())
		upd();
	double ret=(double)ans/(double)n*1.0;
	cout<<fixed<<setprecision(2)<<ret<<endl;
	return 0;
}