【BZOJ1070】【codevs2436】修车，费用流

传送门1
传送门2
写在前面：五一已过，五四示威
思路：分层建图，其实可以联想星际转移问题，把m个修车人拆成n*m个点，对于同一个修车人来说，他的第i个点就是说负责了第n-i+1个修车任务（车的编号不一定是n-i+1）
（以下来自Yveh’s blog）

这样从源连到这n*m个点，容量为1费用为0

后面n个点（表示b辆车）连到汇点，容量为1费用为0

n*m个工人与n辆车互相连接，容量为1，费用为cost[第i个工人][第j辆车]*k（k表示这是修的倒数第几辆车）

考虑修任意一辆车，影响的只能是在它后面修的车(和他本身)，所以它后面的车每辆都要多花的时间，就是这辆车的维修时间，所以要*k。

注意：n，m的顺序注意啊！小心数组越界
代码：

#include<bits/stdc++.h>#define inf 0x7fffff using namespace std;int m,n,tot=1,s,t,ans;int up[1000],dis[1000],first[1000],a[65][10];bool vis[1000];struct edge{    int u,v,w,next,cost;}e[150000];queue<int>q;void add(int x,int y,int z,int c){e[++tot]=(edge){x,y,z,first[x],c};first[x]=tot;}bool spfa(){    memset(dis,63,sizeof(dis));    q.push(s);dis[s]=0;    while (!q.empty())    {        int k=q.front();q.pop();vis[k]=0;        for (int i=first[k];i;i=e[i].next)            if (e[i].w&&dis[e[i].v]>dis[k]+e[i].cost)            {                up[e[i].v]=i;                dis[e[i].v]=dis[k]+e[i].cost;                if (!vis[e[i].v]) vis[e[i].v]=1,q.push(e[i].v);            }    }    return dis[t]<inf;}void flow(){    int minn=inf;    for (int i=up[t];i;i=up[e[i].u])        minn=min(minn,e[i].w);    for(int i=up[t];i;i=up[e[i].u])        ans+=minn*e[i].cost,        e[i].w-=minn,e[i^1].w+=minn;} main(){    scanf("%d%d",&m,&n);    for (int i=1;i<=n;i++)        for (int j=1;j<=m;j++) scanf("%d",a[i]+j);    t=n*m+n+1;    for (int i=1;i<=n*m;i++) add(s,i,1,0),add(i,s,0,0);    for (int i=1;i<=n;i++) add(i+n*m,t,1,0),add(t,i+n*m,0,0);    for (int i=1;i<=m;i++)        for (int j=1;j<=n;j++)            for (int k=1;k<=n;k++)                add((j-1)*m+i,n*m+k,1,j*a[k][i]),                add(n*m+k,(j-1)*m+i,0,-j*a[k][i]);    while (spfa()) flow();    printf("%.2lf",(double)ans/n); }