hdu 5619 Jam's store (mcmf)

http://acm.hdu.edu.cn/showproblem.php?pid=5619
官方题解：
这次肯定有人ak不用看了，第五题是道水题，当然也要看你有没有做过这类的经典题和懂不懂，一道经典的费用流的构图，我们想想，这个点排在另一个点后面的话，只是被另一个点影响了罢了。那么这个构图就巧妙了，因为N不大，所以我们可以把修电脑的人都拆成NNN个点，每个点表示是倒数第几个顾客修的。那大家试想一下，比如说我第一个修电脑的，我现在是倒数第二个来找这个人，那么我只影响我后面的一个和自己，所以的话修一条费用为Tij∗2Tij*2Tij∗2，流量为111的边流过去就好了。可能还有不明白的自己好好理一下，之后每个拆点都直接连到汇点。前面的话就是起点连到每个顾客，然后修电脑的人连到每个拆点，当流量满流的时候就表示每个顾客要做的都流完了，当然这样子还有保证准确性，因为时间都是正数，所以的话缩小流就保证了在每个人都会从流量小也就是倒数第一个留起。

说白了其实就是拆点，把每个店员拆点拆开

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
using namespace std;
/** mcmf 最小费用最大流 */
const int MAX=5000;
const int maxm=500005;
const int inf=0x3f3f3f3f;
int  head[MAX],s,t,cnt,n,m;
int  flow=0,ans=0;
int  d[MAX];
int  pre[MAX];
bool vis[MAX];
int q[maxm+100005];

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

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()
{
    memset(d,inf,sizeof(d));
    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",&m,&n);
        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++)
            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);
                }
        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("%d\n",ans);
    }
    return 0;
}