poj 3686 时间片的网络流

有N个订单和M个机器，给出第i个订单在第j个机器完成的时间Mij，每台机器同一时刻只能处理一个订单，机器必须完整地完成一个订单后才能接着完成下一个订单。问N个订单完成时间的平均值最少为多少。

把每个工厂分成n个，每个代表某时刻，某时刻的工厂只能执行一个。这样解释就很合理了。

但是费用怎么求呢。

对于某个工厂，T=（z1) + (z1+z2)+ （z1+z2+z3)

变形得到 T=1*z3  + 2*z2 + 3*z1 ！！！

也就是倒数第i个进场，他的花费就是i*z

所以对于某个工厂，某个玩具在不同时刻的花费就确定了。z，z*2，z*3，z*k类推。

有人可能无法理解，觉着这样如果两个玩具都选择了同一时刻的工厂怎么办。

这是行不通过的，每一时刻的工厂还连接着汇点，流量是1，这代表着只能走一个，所以不会出现问题

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>


using namespace std;
const int N=4200;
const int MAXE=535000;
const int inf=1<<30;
int  head[N],s,t,cnt,n,m;
int  flow=0,ans=0;
int  d[N];
int  pre[N];
bool vis[N];
int q[MAXE+79000]; 


struct Edge
{
    int u,v,next;
    int c,w;
}edge[MAXE];


void addedge(int u,int v,int c,int  w)
{
    edge[cnt].u=u;
    edge[cnt].v=v;
    edge[cnt].w=w;
    edge[cnt].c=c;
    edge[cnt].next=head[u];
    head[u]=cnt++;
    edge[cnt].v=u;
    edge[cnt].u=v;
    edge[cnt].w=-w;
    edge[cnt].c=0;
    edge[cnt].next=head[v];
    head[v]=cnt++;
}


int SPFA()
{
    for(int i=0;i<=t;i++)
        d[i]=inf;
    memset(pre,-1,sizeof(pre));
    memset(vis,0,sizeof(vis));
    d[s]=0;
    int l=0,r=0;
    q[r++]=s;
    vis[s]=1;
    while(l<=r)
    {
        int u=q[l++];
        vis[u]=0;
        for(int j=head[u];j!=-1;j=edge[j].next)
        {
            int v=edge[j].v;
            if(edge[j].c>0&&d[u]+edge[j].w<d[v])
            {
                d[v]=d[u]+edge[j].w;
                pre[v]=j;
                if(!vis[v])
                {
                    vis[v]=1;
                    q[r++]=v;
                }
            }
        }
    }
    if(d[t]==inf)
        return 0;
    return 1;
}


void MCMF()
{
    flow=0;
    while(SPFA())
    {
        ans+=d[t];
        int u=t;
        int mini=inf;
        while(u!=s)
        {
            if(edge[pre[u]].c<mini)
                mini=edge[pre[u]].c;
                u=edge[pre[u]].u;
        }
        flow+=mini;
        u=t;
        while(u!=s)
        {
            edge[pre[u]].c-=mini;
            edge[pre[u]^1].c+=mini;
            u=edge[pre[u]].u;
        }
    }
}


int mp[160][160];


int main()
{
 int cas;
 scanf("%d",&cas);
 while(cas--)
 {
    scanf("%d%d",&n,&m);
    s=0;
    ans=0;
    flow=0;
    t=n+n*m+1;
    memset(head,-1,sizeof(head));
    cnt=0;
    for(int i=1;i<=n;i++)
      for(int j=1;j<=m;j++)
        scanf("%d",&mp[i][j]);
    for(int i=1;i<=n;i++) //s到每一个玩具，都制作一个
     addedge(s,i,1,0);
    for(int i=1;i<=n;i++)//对于当前玩具
      for(int j=1;j<=m;j++)//对于当前工厂
       for(int k=1;k<=n;k++)//当前工厂的每一时刻
        {
            addedge(i,n+(j-1)*n+k,1,mp[i][j]*k); //话费为v*k
        }
    for(int i=1;i<=m;i++)//工厂
     for(int j=1;j<=n;j++)//每一时刻
     {
         addedge(n+(i-1)*n+j,t,1,0);//每一工厂的每一时刻都与汇点链接
     }
     MCMF();
     printf("%.6lf\n",ans*1.0/n);
 }
 return 0;
}