bzoj1070: [SCOI2007]修车(最小费用最大流)

来源:互联网 发布:网络零售的理解 编辑:程序博客网 时间:2024/05/29 07:28

题目传送门
继续膜一发费用流。
考构图的神题。
然而我审错题了读入的时候m,n读反了,WA了半天美滋滋。。

解法:
一个技术员可能修多辆车。
那么他每修一辆车只能对后面需要修的车产生影响。
然而我们并不知道某技术员修多少辆车。
所以我们采用倒数。
那么这个技术员倒数修的车是谁就会对后面的车都有影响(包括自己)
那么就差点咯。
一个技术员最多修n辆车。所以每个技术员拆成n个点。
第j个技术员的第i个点表示倒数第i个修的车。
那么话费就是time*i了。
美滋滋这题真是考构图然而我并不会。

代码实现:

#include<cstdio>#include<cstdlib>#include<cstring>#include<iostream>#include<algorithm>using namespace std;struct node {    int x,y,c,d,next,other;}a[210000];int len,last[11000];void ins(int x,int y,int c,int d) {    int k1,k2;    len++;k1=len;    a[len].x=x;a[len].y=y;a[len].c=c;a[len].d=d;    a[len].next=last[x];last[x]=len;    len++;k2=len;    a[len].x=y;a[len].y=x;a[len].c=0;a[len].d=-d;    a[len].next=last[y];last[y]=len;    a[k1].other=k2;    a[k2].other=k1;}int st,ed;int head,tail,list[11000];int d[11000],fbian[11000];bool v[11000];bool spfa() {    head=1;tail=2;list[1]=st;    memset(v,false,sizeof(v));v[st]=true;    for(int i=1;i<=ed;i++)        d[i]=2147483647;    d[st]=0;    memset(fbian,0,sizeof(fbian));    while(head!=tail) {        int x=list[head];        for(int k=last[x];k;k=a[k].next) {            int y=a[k].y;            if(a[k].c>0&&d[y]>d[x]+a[k].d) {                d[y]=d[x]+a[k].d;                fbian[y]=k;                if(v[y]==false) {                    v[y]=true;                    list[tail++]=y;                    if(tail==ed+1)                        tail=1;                }            }        }        head++;        if(head==ed+1)            head=1;        v[x]=false;    }    if(d[ed]==2147483647)        return false;    return true;}int map[1100][1100];int main() {    int n,m;scanf("%d%d",&m,&n); //然而这里我输反了。    memset(map,0,sizeof(map));    st=n*m+n+1;ed=st+1;    for(int i=1;i<=n;i++)        for(int j=1;j<=m;j++)            scanf("%d",&map[i][j]);    len=0;memset(last,0,sizeof(last));     for(int i=1;i<=n;i++)        ins(st,i,1,0); //每辆车只修一次。    for(int i=n+1;i<=n+n*m;i++)        ins(i,ed,1,0); //每次只能修一辆车    for(int i=1;i<=n;i++)         for(int j=1;j<=m;j++)            for(int k=1;k<=n;k++)                ins(i,j*n+k,1,k*map[i][j]); //构图    int ans=0;    while(spfa()==true) {        int flow=2147483647,x=ed;        while(x!=st) {            int k=fbian[x];            flow=min(flow,a[k].c);            x=a[k].x;        }        x=ed;        while(x!=st) {            int k=fbian[x];            a[k].c-=flow;a[a[k].other].c+=flow;            x=a[k].x;        }        ans+=d[ed]*flow;    }    printf("%.2lf\n",double(ans)/n);    return 0;}

神题。