POJ 3686

尝试着用二分图最大完美匹配来求解此题。

把每个商品店j拆成n个点，分别代表倒数第p次在商品店加工，把每个order向每个点连线，权值为p*zij。然后用km求解。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
#include <algorithm>
#include <cmath>
#include <stack>
#include <vector>
#define LL long long
#define myabs(x) ((x)>0?(x):(-(x)))
using namespace std;
const int inf=0x3f3f3f3f;
const int maxn=50+10;
const int maxm=50*50+10;
int g[maxn][maxm];
int n,m;
int S[maxn],T[maxm];
int lv[maxn],rv[maxm],leftp[maxm];
int slack;
int tot1,tot2;
void update()
{
	int i;
	for(i=1;i<=tot1;i++) if(S[i]) lv[i]-=slack;
    for(i=1;i<=tot2;i++)  if(T[i]) rv[i]+=slack; 
}
int match(int u)
{
    S[u]=1;
	int i,j;
	for(i=1;i<=tot2;i++)
    {
		if(lv[u]+rv[i]==g[u][i])
        {
            if(!T[i])
            {
                T[i]=1;
                if(leftp[i]==-1||match(leftp[i]))
                {
                    leftp[i]=u;
                    return 1;
                }
            }
        }
		else slack=min(lv[u]+rv[i]-g[u][i],slack);
    }
    return 0;
}

int  solve()
{
    int i,j;
	memset(leftp,-1,sizeof(leftp));
	for(i=1;i<=tot1;i++) lv[i]=0;
	for(i=1;i<=tot2;i++) rv[i]=0;
	for(i=1;i<=tot1;i++)
    {
		for(j=1;j<=tot2;j++)
			lv[i]=max(lv[i],g[i][j]);
    }
	for(i=1;i<=tot1;i++)
    {
        for(;;)
        {
			for(j=1;j<=tot1;j++)  S[j]=0;
			for(j=1;j<=tot2;j++) T[j]=0;
            slack=inf;
            if(match(i)) break;
			if(slack==inf) return -inf; 
			update();
        }
    }
    int sum=0;
    for(i=1;i<=tot1;i++) sum+=lv[i];
	for(i=1;i<=tot2;i++) sum+=rv[i];
    return sum;
}
int main()
{
	int T;
	cin>>T;
	while(T--)
	{
		scanf("%d%d",&n,&m);
		int i,j,k,cost;
		for(i=1;i<=n;i++)
		{
			for(j=0;j<m;j++)
			{
				scanf("%d",&cost);
				for(k=1;k<=n;k++)
					g[i][j*n+k]=-k*cost;
			}
		}
		tot1=n;
		tot2=n*m;
		double ans=-solve();
		printf("%.6lf\n",ans/n);
	}
	return 0;
}