bzoj1070: [SCOI2007]修车（最小费用最大流）

题目传送门
继续膜一发费用流。
考构图的神题。
然而我审错题了读入的时候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;
}

神题。