bzoj 1070 费用流

//可以网络流，但是要怎么分配每辆车让谁维修以及维修顺序呢。可以考虑每辆车维修时间对总结果的贡献，把每个修车人拆成n个点共n*m个点，
//n辆车连向这n*m个点，流量1，费用k*修车时间，其中k(1=<k<=n)表示这辆车是这个技术人员修的倒数第k俩车，他后面k-1辆被同一个人维修的车
//都要等待他修好，因此他的贡献值就是k*维修时间。然后源点连向n辆车，n*m个点连向汇点。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
const int INF=0x7fffffff;
const int MAXN=;
const int MAXM=;
int N,a[][],tot,head[MAXN+],pre[MAXN],dis[MAXN];
bool vis[MAXN+];
struct Edge
{
    int to,next,cap,flow,cost;
}edge[MAXM+];
void init(int x)
{
    N=x;
    tot=;
    memset(head,-,sizeof(head));
}
void add(int x,int y,int w,int c)
{
    edge[tot].to=y;
    edge[tot].cap=w;
    edge[tot].cost=c;
    edge[tot].flow=;
    edge[tot].next=head[x];
    head[x]=tot++;
    edge[tot].to=x;
    edge[tot].cap=;
    edge[tot].cost=-c;
    edge[tot].flow=;
    edge[tot].next=head[y];
    head[y]=tot++;
}
bool Spfa(int s,int t)
{
    queue<int>q;
    for(int i=;i<=N;i++){
        vis[i]=;
        dis[i]=INF;
        pre[i]=-;
    }
    vis[s]=;
    dis[s]=;
    q.push(s);
    while(!q.empty()){
        int u=q.front();q.pop();
        vis[u]=;
        for(int i=head[u];i!=-;i=edge[i].next){
            int v=edge[i].to;
            if(edge[i].cap>edge[i].flow&&dis[v]>dis[u]+edge[i].cost){
                dis[v]=dis[u]+edge[i].cost;
                pre[v]=i;
                if(!vis[v]){
                    vis[v]=;
                    q.push(v);
                }
            }
        }
    }
    if(pre[t]==-) return ;
    return ;
}
int MaxCostFlow(int s,int t)
{
    int Flow=,Cost=;
    while(Spfa(s,t)){
        int Min=INF;
        for(int i=pre[t];i!=-;i=pre[edge[i^].to])
            Min=min(Min,edge[i].cap-edge[i].flow);
        for(int i=pre[t];i!=-;i=pre[edge[i^].to]){
            edge[i].flow+=Min;
            edge[i^].flow-=Min;
            Cost+=edge[i].cost*Min;
        }
        Flow+=Min;
    }
    return Cost;
}
int main()
{
    int n,m;
    while(scanf("%d%d",&m,&n)==){
        init(n+n*m+);
        for(int i=;i<=n;i++){
            add(,i,,);
            for(int j=;j<=m;j++)
                scanf("%d",&a[i][j]);
        }
        for(int i=;i<=n;i++){
            for(int j=;j<=m;j++){
                for(int k=;k<=n;k++){
                    add(i,j*n+k,,k*a[i][j]);
                }
            }
        }
        for(int i=n+;i<=n+n*m;i++)
            add(i,n+n*m+,,);
        printf("%.2lf\n",1.0*MaxCostFlow(,n+n*m+)/(1.0*n));
    }
    return ;
}