1070: [SCOI2007]修车 费用流

来源:互联网 发布:黑客帝国矩阵革命在线 编辑:程序博客网 时间:2024/04/30 00:55

想了一会儿还是不知道怎么处理顺序这个情况。。。于是看了题解(太弱啦qwq。。)


膜hzwer大大

N辆车,M个工人。
把每个工人拆成N个点。记为A[i,j]表示第i个工人修倒数第j辆车。
每个车跟所有N*M个工人拆出的点连边。流量为1,费用为time[i,j]*k。
源和每辆车连边,N*M个点和汇连边,流量都为1,费用同为0。
为什么这么构图呢?
考虑第i个工人,他修第j辆车只对后面要修的车有影响,而前面修过的车已经对当前没有影响了。
而这个影响就是后面每个将要修理的车都多等待了time的时间。
其他边流量都为1是显然的,每辆车修一次,每个工人一个时段只能修理一辆车。

#include<iostream>#include<cstdio>#include<cstring>#define inf 1000000007#define M 605using namespace std;int n,m,ans,cnt=1,T;int t[65][10];int head[605],dis[605],q[605],path[605];int next[100005],list[100005],from[100005],flow[100005],cost[100005];bool vis[605];inline int read(){    int a=0,f=1; char c=getchar();    while (c<'0'||c>'9') {if (c=='-') f=-1; c=getchar();}    while (c>='0'&&c<='9') {a=a*10+c-'0'; c=getchar();}    return a*f;}inline void insert(int x,int y,int z,int w){    next[++cnt]=head[x];    head[x]=cnt;    from[cnt]=x;    list[cnt]=y;    flow[cnt]=z;    cost[cnt]=w;}inline bool spfa(){    for (int i=0;i<=T;i++) dis[i]=inf;    vis[0]=1; dis[0]=0; q[1]=0;    int t=0,w=1,x;    while (t!=w)    {        t=(t+1)%M;        x=q[t];        for (int i=head[x];i;i=next[i])            if (flow[i]&&dis[list[i]]>dis[x]+cost[i])            {                dis[list[i]]=dis[x]+cost[i];                path[list[i]]=i;                if (!vis[list[i]])                {                    vis[list[i]]=1;                    w=(w+1)%M;                    q[w]=list[i];                }            }        vis[x]=0;    }    return dis[T]!=inf;}inline void mcf(){    int x=inf;    for (int i=path[T];i;i=path[from[i]]) x=min(x,flow[i]);    for (int i=path[T];i;i=path[from[i]])        ans+=x*cost[i],flow[i]-=x,flow[i^1]+=x;}int main(){    m=read(); n=read(); T=n*m+n+1;    for (int i=1;i<=n;i++)        for (int j=1;j<=m;j++)            t[i][j]=read();    for (int i=1;i<=n*m;i++)        insert(0,i,1,0),insert(i,0,0,0);    for (int i=1;i<=n;i++)        insert(i+n*m,T,1,0),insert(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++)                insert((i-1)*n+j,n*m+k,1,t[k][i]*j),insert(n*m+k,(i-1)*n+j,0,-t[k][i]*j);    while (spfa()) mcf();    printf("%.2lf",(double)(ans)/(double)(n));    return 0;}
0 0
原创粉丝点击