HDU 5619 Jam's store（最小费用最大流-mcmf）

Description
Jam不好好学习，然后就去帮别人修电脑了，在一家店里，有M个店员，现在有N个顾客，给出每个顾客对应给每个店员的修电脑的时间为T​ij，问所有顾客要等待的最少时间。当然，一个顾客在某个店员那里完成之后，那个店员才会执行下一个顾客的任务
Input
第一行T(1≤T≤100)，表示T组数据。
接下来T组数据：
每组数据第一行为M,N(1≤M,N≤20)表示店员数和顾客数
接下来N行M列，每个整数表示第i个顾客找第j个店员的时间(1≤T​ij≤1000)
Output
输出一个数，表示用时总时间
Sample Input
1
4 3
4 4 1 5
8 2 5 6
4 5 10 5
Sample Output
7
Solution
因为一个店员最多会给n名顾客修电脑，而当其给某一个顾客修电脑时，后面找他修电脑的人每一个人都在等待，即如果第i名顾客是倒数第k个去找第j名店员修电脑的话，那么后k-1名顾客吗，每个人都需要等待Tij的时间，而这名顾客也需要等待Tij的时间才能修好自己的电脑，所以把m名店员每一个都拆成n个点，将m名顾客看作m个点，如果第j名顾客向第i名店员的第k个点连边就表示这名顾客是这位店员的倒数第k名顾客，那么这条边的容量就是1，费用就是k*Tij，因为不知道每名顾客会去找哪位店员去修电脑，也不知道其是作为该名店员的倒数第几名顾客，所以所有情况都要建边，之后跑一边最小费用最大流即可，最大流保证每名顾客的电脑都被修好，最小费用保证所有顾客等待的总时间最短
Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
#define maxn 555
#define maxm 22222
#define INF 0x3f3f3f3f
int head[maxn],d[maxn],s,e,no,vis[maxn],pre[maxn];//s为源点,e为汇点 
struct point
{
    int u,v,flow,next,cost;
    point(){};
    point(int _u,int _v,int _next,int _flow,int _cost)
    {
        u=_u,v=_v,next=_next,flow=_flow,cost=_cost;
    } 
}p[maxm];
void add(int x,int y,int z,int c)//从x到y建一条容量为z,花费为c的边 
{
    p[no]=point(x,y,head[x],z,c);   
    head[x]=no++;
    p[no]=point(y,x,head[y],0,-c);  
    head[y]=no++;
}
void init()//初始化 
{
    memset(head,-1,sizeof(head));
    no=0;
}
bool spfa()
{
    int i,x,y;
    queue<int>q;
    memset(d,0x3f,sizeof(d));
    memset(vis,false,sizeof(vis));
    memset(pre,-1,sizeof(pre));
    d[s]=0;   
    vis[s]=true;  
    q.push(s);
    while(!q.empty())
    {
        x=q.front();  
        q.pop();    
        vis[x]=false;
        for(i=head[x];i!=-1;i=p[i].next)
        {
            if(p[i].flow&&d[y=p[i].v]>d[x]+p[i].cost)
            {
                d[y]=d[x]+p[i].cost;   
                pre[y]=i;
                if(vis[y])  
                    continue;
                vis[y]=true;  
                q.push(y);
            }  
        }  
    }
    return d[e]!=d[e+1];
}
int mcmf()//最小费用最大流 
{
    int mincost=0,maxflow=0,minflow,i;
    while(spfa())
    {
        minflow=INF;
        for(i =pre[e];i!=-1;i=pre[p[i].u])
            minflow=min(minflow,p[i].flow);
        for(i=pre[e];i!=-1;i=pre[p[i].u])
        {
            p[i].flow-=minflow;
            p[i^1].flow+=minflow;
        }
        mincost+=d[e]*minflow; 
        maxflow+=minflow;
    }
    return mincost;//最小费用,最大流为maxflow 
}
int T,n,m,t[22][22]; 
int main()
{
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&m,&n);
        for(int i=1;i<=n;i++)
            for(int j=1;j<=m;j++)
                scanf("%d",&t[i][j]);
        init();
        s=0,e=(m+1)*n+1;
        for(int i=1;i<=n;i++)
            for(int j=1;j<=m;j++)
                for(int k=1;k<=n;k++)
                    add(i,j*n+k,1,k*t[i][j]);
        for(int i=1;i<=n;i++)add(s,i,1,0);
        for(int j=1;j<=m;j++)
            for(int k=1;k<=n;k++)
                add(j*n+k,e,1,0);
        int ans=mcmf();
        printf("%d\n",ans);
    }
    return 0;
}