码迷,mamicode.com
首页 > 其他好文 > 详细

BZOJ 1070: [SCOI2007]修车

时间:2017-12-10 11:18:18      阅读:158      评论:0      收藏:0      [点我收藏+]

标签:+=   inline   stdin   space   pac   class   desc   cpp   build   

BZOJ 1070: [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

Solution

十分精巧的构图方式。
把一个工人拆成n个点,即有n*m个点,A[i][j]表示第i个工人,倒数第j个修的车。
因为修一个车只会影响到后面该工人要修的车的时间。所以跑一个最小费用流就行了。

Code

#include<cstdio>
#include<cstring>
#include<cmath>
#include<queue>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
#define rep(i,x) for(int i=head[x];i;i=e[i].next)
#define mem(a,x) memset(a,x,sizeof(a))
typedef long long LL;
typedef double DB;
using namespace std;
inline int read() {
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9') f=(ch=='-')?-1:f,ch=getchar();
    while(ch>='0'&&ch<='9') x=x*10+(ch-'0'),ch=getchar();return f*x;
}
const int inf=0x3f3f3f3f,T=1001,S=0;
using namespace std;
int n,m,cnt=1,ans,l[61][10];
int d[1005],q[1005],from[1005],head[1005];
bool inq[1005];
struct edge {int from,to,next,c,v;} e[100001];
void ins(int u,int v,int r,int c) {e[++cnt]=(edge){u,v,head[u],c,r},head[u]=cnt;}
void insert(int u,int v,int r,int c) {ins(u,v,r,c),ins(v,u,0,-c);}
bool spfa() {
    fo(i,0,T) d[i]=inf;
    int l=0,r=1;
    d[S]=0,inq[S]=1,q[0]=S;
    while(l!=r) {
        int now=q[l++];if(l==T) l=0;
        rep(i,now) if(e[i].v&&d[e[i].to]>d[now]+e[i].c) {
            d[e[i].to]=d[now]+e[i].c;
            from[e[i].to]=i;
            if(!inq[e[i].to]) {
                inq[e[i].to]=1;
                q[r++]=e[i].to;if(r==T) r=0;
            }
        }
        inq[now]=0;
    }
    if(d[T]==inf) return 0;
    return 1;
}
void dfs() {
    int x=inf;
    for(int i=from[T];i;i=from[e[i].from]) x=min(x,e[i].v);
    for(int i=from[T];i;i=from[e[i].from]) e[i].v-=x,e[i^1].v+=x,ans+=e[i].c*x;
}
void ek_spfa() {while(spfa()) dfs();}
void build() {
    fo(i,1,n*m) insert(S,i,1,0);
    fo(i,n*m+1,n*m+m) insert(i,T,1,0);
    fo(i,1,n) fo(j,1,m) fo(k,1,m) insert((i-1)*m+j,n*m+k,1,l[k][i]*j);
}
int main() {
    freopen("1070.in","r",stdin);
    n=read(),m=read();
    fo(i,1,m) fo(j,1,n) l[i][j]=read();
    build(),ek_spfa();
    printf("%.2lf",(double)ans/m);
    return 0;
}

BZOJ 1070: [SCOI2007]修车

标签:+=   inline   stdin   space   pac   class   desc   cpp   build   

原文地址:http://www.cnblogs.com/patricksu/p/8016211.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!