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 0x3f3f3f3fint 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;}