POJ-3686 The Windy's 最小权匹配

　　题目链接：http://poj.org/problem?id=3686

　　这个题目容易想到拆点，开始我是直接拆点在求的，考虑到直接建图复杂度有点高，然后就在找增广路的同时再加点，因为并不是每个点都会达到50个，这样实际复杂度会好很多。敲完之后，样例过了，然后直接叫，WA！后来看了discuss顿然醒悟，找增光路的时候，之前的状态改变了，那么后面的状态也会改变，那么那样求下去就不一定是最优解！所有要保证所有的匹配之间是互补影响的！假设某个工作坊要处理K个命令，时间分别是t1,t2......tk，那么总共的时间就是k*t1+(k-1)*t2+......+tk。发现没有，如果k表示的是倒数第k个命令，那么总时间变为了：t1+2*t2+......+k*tk。所以我们拆点后再建图的时候，边权值w[i][j]*k,(1<=k<=50)，那么直接求最小权匹配就行了。那么我们看下时间复杂度，匹配n个点O(n)，每个点至多找n次增广路O(n)，邻接矩阵的匈牙利O(n^2)，总共就是O(50*50*2500*2500)，这是理论的时间复杂度，但实际上远远没有这么可怕，因为很难达到最坏情况，而且顶标的维护使得匈牙利算法效率大大的高，所以一般就算O(n^4)的KM算法也够用了。我的代码POJ跑了47ms，这也出乎了我的意料，总感觉POJ的题目数据很弱啊！！！

//STATUS:G++_AC_47MS_1344KB
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#include<iostream>
#include<string>
#include<algorithm>
#include<vector>
#include<queue>
#include<stack>
using namespace std;
#define LL long long
#define Max(a,b) ((a)>(b)?(a):(b))
#define Min(a,b) ((a)<(b)?(a):(b))
#define mem(a,b) memset(a,b,sizeof(a))
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
const int MAX=55,INF=0x3f3f3f3f;

int w[MAX][MAX*MAX],S[MAX],T[MAX*MAX],lx[MAX],ly[MAX*MAX],y[MAX*MAX];
int Test,n,m,slack;

int dfs(int u)
{
    int v,t;
    S[u]=1;
    for(v=0;v<m;v++){
        t=w[u][v]-lx[u]-ly[v];
        if(!t){
            if(!T[v]){
                T[v]=1;
                if(y[v]==-1 || dfs(y[v])){
                    y[v]=u;
                    return 1;
                }
            }
        }
        else if(t<slack)slack=t;
    }
    return 0;
}

double KM()
{
    int i,j;
    double s;
    mem(ly,0);
    mem(y,-1);
    for(i=0;i<n;i++){
        lx[i]=INF;
        for(j=0;j<m;j++)
            if(w[i][j]<lx[i])lx[i]=w[i][j];
    }
    for(i=0;i<n;i++){
        while(1){
            slack=INF;
            mem(S,0);mem(T,0);
            if(dfs(i))break;
            for(j=0;j<n;j++)if(S[j])lx[j]+=slack;
            for(j=0;j<m;j++)if(T[j])ly[j]-=slack;
        }
    }
    for(i=s=0;i<m;i++)
        if(y[i]!=-1)s+=w[y[i]][i];
    return s;
}

int main()
{
 //   freopen("in.txt","r",stdin);
    int i,j,k,t;
    double ans;
    scanf("%d",&Test);
    while(Test--)
    {
        scanf("%d%d",&n,&m);
        for(i=0;i<n;i++)
            for(j=0;j<m;j++)
                scanf("%d",&w[i][j]);
        for(i=0,t=m;i<t;i++)
            for(j=2;j<=50;j++,m++)
                for(k=0;k<n;k++)
                    w[k][m]=w[k][i]*j;
        ans=KM();

        printf("%.6f\n",ans/n);
    }
    return 0;
}