尝试着用二分图最大完美匹配来求解此题。
把每个商品店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;
}