费用流模板

inline void add(int x,int y,int f,int c){
    node[hcnt]=Node(y,head[x],f,c);head[x]=hcnt++;
    node[hcnt]=Node(x,head[y],0,-c);head[y]=hcnt++;
}
int spfa(){
    mst(d,-1);d[S]=0;
    mst(vis,0);q.push(S);
    mst(pre,-1);
    while(!q.empty()){
        int x=q.front();q.pop();
        vis[x]=0;
        for(int i=head[x];~i;i=node[i].next){
            int e=node[i].to;
            if(node[i].f&&d[e]<d[x]+node[i].c){
                d[e]=d[x]+node[i].c;
                pre[e]=i^1;
                if(!vis[e]){
                    vis[e]=1;
                    q.push(e);
                }
            }
        }
    }
    return pre[T]!=-1;
}
void work(){
    int fl=inf;
    for(int i=pre[T];~i;i=pre[node[i].to]){
        fl=min(fl,node[i^1].f);
    }
    for(int i=pre[T];~i;i=pre[node[i].to]){
        node[i].f+=fl;node[i^1].f-=fl;
       // ans+=node[i^1].c;
    }
    ans+=d[T]; ///ans是费用
}