POJ 3686 The Windy's （二分图最小权匹配 拆点 构图）

题意：

有n个订单m个车间，每个车间均可以单独完成任何一个订单。每个车间完成不同订单的时间是不同的。不会出现两个车间完成同一个订单的情况。给出每个订单在某个车间完成所用的时间。问订单完成的平均时间是多少。

思路：

1、这个题在建图上有一些需要思考很长时间的地方。因为每个订单所消耗的时间是车间完成订单的时间加上订单等待的时间。我们设在车间A需要完成k个订单，消耗的总时间是t1+(t1+t2)+(t1+t2+t3)……转换一下就是t1*k+t2*(k-1)+t3*(k-3)……我们就找到了规律：当第i个订单在第j个车间是倒数第k个任务时，总消耗时间需要加上订单i在车间对应消耗时间的k倍。

2、建图的时候记得用负值来建图，这样顶标数组初始化成0就行了。因为不会有比0大的值出现。

3、结果出来后取相反数，然后除以n来得到最后的结果。

（图片来源）

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<algorithm>
using namespace std;
const int maxn = 2505;
const int INF = 0x3f3f3f3f;
int t[maxn][maxn];
//
int nx, ny;
int match[maxn], lx[maxn], ly[maxn], slack[maxn];
int visx[maxn], visy[maxn], w[maxn][maxn];

int Hungary(int x)
{
    visx[x] = 1;
    for(int y = 1; y <= ny; y++)
    {
        if(visy[y]) continue;
        int temp = lx[x]+ly[y]-w[x][y];
        if(temp == 0)
        {
            visy[y] = 1;
            if(match[y] == -1 || Hungary(match[y]))
            {
                match[y] = x;
                return 1;
            }
        }
        else if(slack[y] > temp)
            slack[y] = temp;
    }
    return 0;
}

int KM()
{
    memset(match, -1, sizeof(match));
    memset(ly, 0, sizeof(ly));
    for(int i = 1; i <= nx; i++)
        lx[i] = -INF;
    for(int i = 1; i <= nx; i++)
        for(int j = 1; j <= ny; j++)
            if(w[i][j] > lx[i])
                lx[i] = w[i][j];
    for(int x = 1; x <= nx; x++)
    {
        for(int i = 1; i <= ny; i++)
            slack[i] = INF;
        while(1)
        {
            memset(visx, 0, sizeof(visx));
            memset(visy, 0, sizeof(visy));
            if(Hungary(x)) break;
            int d = INF;
            for(int i = 1; i <= ny; i++)
                if(!visy[i] && d > slack[i])
                    d = slack[i];
            for(int i = 1; i <= nx; i++)
                if(visx[i])
                    lx[i] -= d;
            for(int i = 1; i <= ny; i++)
                if(visy[i])
                    ly[i] += d;
                else
                    slack[i] -= d;
        }
    }
    int res = 0;
    for(int i = 1; i <= ny; i++)
        if(match[i] != -1)
            res += w[match[i]][i];
    return res;
}

int main(void)
{
    int n, m, T;
    cin >> T;
    while(T--)
    {
        scanf("%d%d", &n, &m);
        for(int i = 1; i <= n; i++)
            for(int j = 1; j <= m; j++)
                scanf("%d", &t[i][j]);
        nx = n, ny = n*m;
        for(int i = 1; i <= n; i++)
            for(int j = 1; j <= m; j++)
                for(int k = 1; k <= n; k++)
                    w[i][(j-1)*n+k] = -t[i][j]*(k);
        printf("%.6f\n", -KM()*1.0/n);
    }

    return 0;
}

Sample Input
3


3 4
100 100 100 1
99 99 99 1
98 98 98 1


3 4
1 100 100 100
99 1 99 99
98 98 1 98


3 4
1 100 100 100
1 99 99 99
98 1 98 98
Sample Output
2.000000
1.000000
1.333333