【bzoj1070】[SCOI2007]修车 网络流

Description

　　同一时刻有N位车主带着他们的爱车来到了汽车维修中心。维修中心共有M位技术人员，不同的技术人员对不同
的车进行维修所用的时间是不同的。现在需要安排这M位技术人员所维修的车及顺序，使得顾客平均等待的时间最
小。 说明：顾客的等待时间是指从他把车送至维修中心到维修完毕所用的时间。

Input

　　第一行有两个m,n，表示技术人员数与顾客数。 接下来n行，每行m个整数。第i+1行第j个数表示第j位技术人
员维修第i辆车需要用的时间T。

Output

　　最小平均等待时间，答案精确到小数点后2位。

Sample Input

2 2

3 2

1 4

Sample Output

1.50

HINT

数据范围: (2<=M<=9,1<=N<=60), (1<=T<=1000)

Source

其实就是让总等待时间最小。

不好贪心什么的，只能网络流了。

把每个人拆成n个点，表示他修的倒数第i个车，然后就枚举每个车被哪个人第几次修就行了。
第i个车被第j个人倒数第k次修，则贡献是k*a[i][j]，因为后面的人要排队等他修完。

其他的很好说，流量为1，费用除了人和车是k*a[i][j]以外其他的都是0

我数组又开小了 自SDOI的MLE之后不敢开太大啊

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
using namespace std;

typedef long long LL;
const int SZ = 10010;
const int INF = 1000000010;

int head[SZ],nxt[300010],tot = 1;

struct edge{
    int f,t,d,c;
}l[300010];

void build(int f,int t,int d,int c)
{
    l[++ tot] = (edge){f,t,d,c};
    nxt[tot] = head[f];
    head[f] = tot;
}

void insert(int f,int t,int d,int c)
{
    build(f,t,d,c); build(t,f,0,-c);
}

deque<int> q;
bool use[SZ];
int dist[SZ],pre[SZ];

bool spfa(int s,int e)
{
    memset(dist,63,sizeof(dist));
    dist[s] = 0;
    use[s] = 1;
    q.push_front(s);
    while(q.size())
    {
        int f = q.front(); q.pop_front();
        use[f] = 0;
        for(int i = head[f];i;i = nxt[i])
        {
            int v = l[i].t;
            if(l[i].d && dist[v] > dist[f] + l[i].c)
            {
                dist[v] = dist[f] + l[i].c;
                pre[v] = i;
                if(!use[v])
                {
                    use[v] = 1;
                    if(q.empty() || dist[v] < dist[q.front()])
                        q.push_front(v);
                    else
                        q.push_back(v);
                }
            }
        }
    }
    if(dist[e] > INF) return false;
    return true;
}

int dfs(int s,int e)
{
    int maxflow = INF;
    for(int i = pre[e];i;i = pre[l[i].f])
        maxflow = min(maxflow,l[i].d);
    int ans = 0;
    for(int i = pre[e];i;i = pre[l[i].f])
        ans += maxflow * l[i].c,l[i].d -= maxflow,l[i ^ 1].d += maxflow;
    return ans;
}   

int EK(int s,int e)
{
    int ans = 0;
    while(spfa(s,e))
        ans += dfs(s,e);
    return ans;
}

int a[233][233];

int main()
{
    int n,m;
    scanf("%d%d",&m,&n);
    for(int i = 1;i <= n;i ++)
        for(int j = 1;j <= m;j ++)
            scanf("%d",&a[i][j]);
    int s = 10000,e = 10001;
    for(int i = 1;i <= n;i ++)
        insert(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 ++)
                insert(i,j * n + k,1,k * a[i][j]);
    for(int i = 1;i <= m;i ++)
        for(int j = 1;j <= n;j ++)
            insert(i * n + j,e,1,0);
    int ans = EK(s,e);
    printf("%.2lf\n",(double)ans / n);
    return 0;
}