タイトル説明
Nの車の所有者が自分の車であると同時に、車の修理センターに来ました。Mサービスセンターは、技術スタッフを持って、使用される様々な異なる車の修理時間の技術スタッフが異なっています。今、このMを最小限に車やオーダーの修復にスキルのビット平均顧客の待ち時間なるように配置する必要があります。
説明:お客様の待ち時間はサービスが彼の車からのサービスセンターにして完了した時点です。
入力形式
最初の行は、二つの数M、Nを有し、その数は、当該技術分野における顧客の数を表します。
m個の整数の次のn行。I + 1番目の行及びjは技術者が車両Tを使用する必要があるI-メンテナンスの回数j番目のビットの数を表します。
出力フォーマット
最小平均待ち時間、2桁まで正確な回答。
#include<cmath>
#include<queue>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int N=1e4+10,M=2e5+10,inf=0x3f3f3f3f;
int n,m,s,t;
int nxt[M],head[N],go[M],edge[M],cost[M],tot=1;
inline void add(int u,int v,int o1,int o2){
nxt[++tot]=head[u],head[u]=tot,go[tot]=v,edge[tot]=o1,cost[tot]=o2;
nxt[++tot]=head[v],head[v]=tot,go[tot]=u,edge[tot]=0,cost[tot]=-o2;
}
int dis[N],ret;
bool vis[N];
inline bool spfa(){
memset(dis,0x3f,sizeof(dis)); dis[s]=0;
queue<int>q; q.push(s);
while(q.size()){
int u=q.front(); q.pop(); vis[u]=0;
for(int i=head[u];i;i=nxt[i]){
int v=go[i];
if(edge[i]&&dis[v]>dis[u]+cost[i]){
dis[v]=dis[u]+cost[i];
if(!vis[v])q.push(v),vis[v]=1;
}
}
}
return dis[t]!=inf;
}
int dinic(int u,int flow){
if(u==t)return flow;
int rest=flow,k;
vis[u]=1;
for(int i=head[u];i&&rest;i=nxt[i]){
int v=go[i];
if(edge[i]&&!vis[v]&&dis[v]==dis[u]+cost[i]){
k=dinic(v,min(rest,edge[i]));
if(!k)dis[v]=-1;
ret+=k*cost[i];
edge[i]-=k;
edge[i^1]+=k;
rest-=k;
}
}
vis[u]=0;
return flow-rest;
}
int c[65][65];
signed main(){
cin>>m>>n;
t=n*m+n+1;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
scanf("%d",&c[i][j]);
for(int i=1;i<=n*m;i++)add(s,i,1,0);
for(int i=1;i<=m;i++)
for(int j=1;j<=n;j++)
for(int k=1;k<=n;k++){
int now=(i-1)*n+k;
add(now,j+n*m,1,c[j][i]*k);
}
for(int i=1;i<=n;i++)add(n*m+i,t,1,0);
int flow=0,maxflow=0;
while(spfa())
while(flow=dinic(s,inf))maxflow+=flow;
printf("%.2lf\n",1.0*ret/n);
}