[bzoj1070][SCOI2007]修车[ 网络流]

把每个工人拆成N个点。记为A[i,j]表示第i个工人修倒数第j辆车.每个车跟所有N*M个工人拆出的点连边。流量为1，费用为$time[i,j]*k$。源和每辆车连边，N*M个点和汇连边，流量都为1，费用同为0。

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <ctime>
#include <cmath>
#include <queue>

using namespace std;

template<const int _n,const int _m>
struct Edge
{
    struct Edge_base { int to,next,w,c; }e[_m]; int    cnt,p[_n];
    void    clear() { cnt=1,memset(p,0,sizeof(p)); }
    Edge() { clear(); }
    void    insert(const int x,const int y,const int z,const int zz)
    { e[++cnt].to=y; e[cnt].next=p[x]; e[cnt].w=z; e[cnt].c=zz; p[x]=cnt; return ; }
    void    link(const int x,const int y,const int z,const int zz)
    { insert(x,y,z,zz); insert(y,x,0,-zz); }
    int    start(const int x) { return p[x]; }
    Edge_base& operator[](const int x) { return e[x]; }
};

int    n,m,Time[61][10],SSS,TTT,Cost;
int    Dis[610],cur[610];
bool    visited[610];
Edge<610,71000>e;

bool    Spfa(const int S)
{
    int    i,t,temp;
    queue<int>    Q;
    memset(Dis,0x3f,sizeof(Dis));
    Dis[S]=0;
    visited[S]=true;
    Q.push(S);
    while(!Q.empty())
    {
        t=Q.front();Q.pop();visited[t]=false;
        for(i=e.start(t);i;i=e[i].next)
        {
            temp=e[i].to;
            if(e[i].w && Dis[t]+e[i].c<Dis[temp])
            {
                Dis[temp]=Dis[t]+e[i].c;
                if(!visited[temp])
                {
                    visited[temp]=true;
                    Q.push(temp);
                }
            }
        }
    }
    return Dis[TTT]!=0x3f3f3f3f;
}

int    Dfs(const int S,const int bk)
{
    visited[S]=true;
    if(S==TTT)return bk;
    int    rest=bk;
    for(int &i=cur[S];i;i=e[i].next)
    {
        if(!visited[e[i].to] && Dis[S]+e[i].c==Dis[e[i].to] && e[i].w)
        {
            int    flow=Dfs(e[i].to,min(rest,e[i].w));
            Cost+=flow*e[i].c;
            e[i].w-=flow;
            e[i^1].w+=flow;
            if((rest-=flow)<=0)break;
        }
    }
    if(bk==rest)Dis[S]=0x3f3f3f3f;
    return bk-rest;
}

int    Zkw()
{
    while(Spfa(SSS))
    {
        do
        {
            memset(visited,0,sizeof(visited));
            memcpy(cur,e.p,sizeof(cur));
            Dfs(SSS,0x3f3f3f3f);
        }while(visited[TTT]);
    }
    return Cost;
}

int main()
{
    int    i,j,k;

    scanf("%d%d",&n,&m);
    for(i=1;i<=m;++i)
        for(j=1;j<=n;++j)
            scanf("%d",&Time[i][j]);
    SSS=n*m+m+1,TTT=SSS+1;
    for(i=1;i<=n*m;++i)e.link(SSS,i,1,0);
    for(i=n*m+1;i<SSS;++i)e.link(i,TTT,1,0);
    for(k=1;k<=m;++k)
        for(i=1;i<=n;++i)
            for(j=1;j<=m;++j)
                e.link((i-1)*m+j,n*m+k,1,Time[k][i]*j);
    printf("%.2lf
",(double)Zkw()/m);
    return 0;
}