poj 3686 The Windy's

题意：

n个工件，m个工场，第i个工件如果在第j个工场加工的时间为z[i][j]，每个工场可以任意加工一个工件，但是同一时刻只能加工一个且一旦开始加工，且这个工件必须要加工完才能停止，求n个工件最小的平均结束时间。 (1 ≤ N,M ≤ 50)

题解：

假设一台机器J要完成K个工件，这K个工件单独在J机器完成时间为t1,t2,...,tk（顺序为1~K），则这K个工件总用时=t1*k+t2*(k-1)+...+tk。（关键点）

建图方法就是采取经典的“拆点法”。将工场拆点，拆成N个。第j个工场的第p个点表示第该工场执行倒数第p个任务。连接该点的边的权值即为t[i][j]*p，这就转换成了二分图求最小权匹配。可以采用KM算法。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
const int maxn=2510;
const int inf=1<<28;
int w[maxn][maxn];
int lx[maxn],ly[maxn];
int linky[maxn];
int visx[maxn],visy[maxn];
int slack[maxn];
int nx,ny;
bool find(int x)
{
    visx[x]=true;
    for(int y=1;y<=ny;y++)
    {
        if(visy[y])continue;
        int t=lx[x]+ly[y]-w[x][y];
        if(t==0)
        {
            visy[y]=true;
            if(linky[y]==-1||find(linky[y]))
            {
                linky[y]=x;
                return true;        //找到增广轨
            }
        }
        else if(slack[y]>t)
            slack[y] = t;
    }
    return false;                   //没有找到增广轨（说明顶点x没有对应的匹配，与完备匹配(相等子图的完备匹配)不符）
}

int KM()                //返回最优匹配的值
{
    int i,j;
    memset(linky,-1,sizeof(linky));
    memset(ly,0,sizeof(ly));

    for(i=1;i<=nx;i++)
        for(j=1,lx[i]=-inf;j<=ny;j++)
            if(w[i][j]>lx[i])
                lx[i]=w[i][j];
    for(int x=1;x<=nx;x++)
    {
        for(i=1;i<=ny;i++)
            slack[i]=inf;
        while(true)
        {
            memset(visx,0,sizeof(visx));
            memset(visy,0,sizeof(visy));
            if(find(x)) break;                  //找到增广轨，退出
            int d=inf;
            for(i=1;i<=ny;i++)          //没找到，对l做调整(这会增加相等子图的边)，重新找
            {
                if(!visy[i]&& d>slack[i])
                    d=slack[i];
            }
            for(i=1;i<=nx;i++)
            {
                if(visx[i])
                    lx[i]-=d;
            }
            for(i=1;i<=ny;i++)
            {
                if(visy[i])
                     ly[i]+= d;
                else
                     slack[i]-= d;
            }
        }
    }
    int result=0;
    for(i=1;i<=ny;i++)
    if(linky[i]>-1)
        result+=w[linky[i]][i];
    return result;
}
int z[maxn][maxn];
int main()
{
    int sec;
    scanf("%d",&sec);
    while(sec--)
    {
        int n,m;
        scanf("%d%d",&n,&m);
        for(int i=1;i<=n;i++)
         for(int j=1;j<=m;j++)
         scanf("%d",&z[i][j]);

        for(int i=1;i<=n;i++)
         for(int j=1;j<=m;j++)
          for(int k=1;k<=n;k++)
          w[i][k+(j-1)*n]=-k*z[i][j];

        nx=n;ny=n*m;
        int ans=KM();
        printf("%.6f\n",-(double)ans/n);
    }
    return 0;
}