[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=,memset(p,,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,,-zz); }
     int    start(const int x) { return p[x]; }
     Edge_base& operator[](const int x) { return e[x]; }
 };

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

 bool    Spfa(const int S)
 {
     int    i,t,temp;
     queue<int>    Q;
     memset(Dis,0x3f,sizeof(Dis));
     Dis[S]=;
     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^].w+=flow;
             if((rest-=flow)<=)break;
         }
     }
     if(bk==rest)Dis[S]=0x3f3f3f3f;
     return bk-rest;
 }

 int    Zkw()
 {
     while(Spfa(SSS))
     {
         do
         {
             memset(visited,,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=;i<=m;++i)
         for(j=;j<=n;++j)
             scanf("%d",&Time[i][j]);
     SSS=n*m+m+,TTT=SSS+;
     for(i=;i<=n*m;++i)e.link(SSS,i,,);
     for(i=n*m+;i<SSS;++i)e.link(i,TTT,,);
     for(k=;k<=m;++k)
         for(i=;i<=n;++i)
             for(j=;j<=m;++j)
                 e.link((i-)*m+j,n*m+k,,Time[k][i]*j);
     printf("%.2lf\n",(double)Zkw()/m);
     return ;
 }