bzoj 1070: [SCOI2007]修车 费用流

时间：2023-03-08 21:28:47

1070: [SCOI2007]修车

Time Limit: 1 Sec  Memory Limit: 162 MB
Submit: 2785  Solved: 1110
[Submit][Status][Discuss]

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)

　　╮(╯▽╰)╭学了那么久费用流，还是属于见到题根本想不到的情况啊。如果省选出这道题，那我不得哭死，一方面是想不到费用流，一方面是关于题目读入中n，m的输入顺序问题。

　　还记得不久前听说有人随便的一句话：网络流的题都能用贪心水过，当时只是笑笑罢了，咋遇到了“贪心”题就想不到用网络流水过呢？

　　只要想到是费用流，图什么的直接秒建，就不赘述了。

#include<iostream>

#include<cstdio>

#include<cstring>

#include<algorithm>

using namespace std;

#define MAXN 100000

#define MAXE MAXV*20

#define MAXV MAXN

#define INF 0x3f3f3f3f

struct Edge

{

        int np,val,cost;

        Edge *next,*neg;

}E[MAXE],*V[MAXV];

int sour=,sink=;

int tope=-;

void addedge(int x,int y,int z,int c)

{

    //    cout<<"Add:"<<x<<" "<<y<<" "<<z<<" "<<c<<endl;

        E[++tope].np=y;

        E[tope].val=z;

        E[tope].cost=c;

        E[tope].next=V[x];

        V[x]=&E[tope];

        E[++tope].np=x;

        E[tope].val=;

        E[tope].cost=-c;

        E[tope].next=V[y];

        V[y]=&E[tope];

        V[x]->neg=V[y];

        V[y]->neg=V[x];

}

int dis[MAXN];

int q[MAXN*];

int vis[MAXN],vistime=;

int prv[MAXN];

Edge *prve[MAXN];

int spfa(int now)

{

        int head=-,tail=;

        Edge *ne;

        memset(dis,INF,sizeof(dis));

        dis[now]=;

        q[]=now;

        vis[now]=++vistime;

        while (head<tail)

        {

                now=q[++head];

                vis[now]=;

                for (ne=V[now];ne;ne=ne->next)

                {

                        if (ne->val && dis[ne->np]>dis[now]+ne->cost)

                        {

                                dis[ne->np]=dis[now]+ne->cost;

                                prv[ne->np]=now;

                                prve[ne->np]=ne;

                                if (vis[ne->np]!=vistime)

                                        q[++tail]=ne->np;

                        }

                }

        }

        return dis[sink]!=INF;

}

pair<int,int> cost_flow()

{

        pair<int,int> res;

        while (spfa(sour))

        {

                int mxflow=INF;

                for (int x=sink;x!=sour;x=prv[x])

                        mxflow=min(mxflow,prve[x]->val);

                for (int x=sink;x!=sour;x=prv[x])

                        prve[x]->val-=mxflow,prve[x]->neg->val+=mxflow;

                res.first+=mxflow;

                res.second+=mxflow*dis[sink];

        }

        return res;

}

int main()

{

        freopen("input.txt","r",stdin);

        int n,m;

        scanf("%d%d",&m,&n);

        int x,y,z;

        for (int i=;i<=n;i++)

        {

                for (int j=;j<=m;j++)

                {

                        scanf("%d",&x);

                        for (int k=;k<=n;k++)

                                addedge(+i,+n+j*n+k,,x*k);

                }

        }

        for (int i=;i<=n;i++)

                addedge(sour,+i,,);

        for (int i=;i<=m;i++)

                for (int j=;j<=n;j++)

                        addedge(+n+i*n+j,sink,,);

        printf("%.2lf\n",(double)cost_flow().second/n);

}

相关文章