HDU 6118 度度熊的交易计划（最小费用可行流\费用流）

思路：

其实就是个费用流的裸题。。
s=0，t=n+1；
对于每个村庄，$addedge(s,i,b[i],a[i])$ 表示增加一条s到i的边，容量为 $b[i]$ ，花费为 $a[i]$
同时 $addedge(i,t,d[i],c[i])$

对于每一条路，容量无限大，费用为距离。

自我总结：

比赛的时候为什么要拆成一条条容量为1的边啊啊啊啊啊。。反正费用都是针对每个流来说的，最后乘一下流量不就好了。。

#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
#include<map>
#define inf 0x3f3f3f3f
typedef long long int lli;
using namespace std;
const lli maxn = 550;
struct edge{
    int to,v,next,c;
}ed[4000000];
int cnte,head[maxn],pre[maxn];//ans是最小费用
int dis[maxn],ans;bool visq[maxn];
inline void ae(int x,int y,int v,int c) {
    ed[cnte].to = y;
    ed[cnte].v = v;
    ed[cnte].c = c;
    ed[cnte].next = head[x];
    head[x] = cnte++;

    ed[cnte].to = x;
    ed[cnte].v = 0;
    ed[cnte].c = -c;
    ed[cnte].next = head[y];
    head[y] = cnte++;
}

bool spfa(int s,int t){
    memset(dis,0x3f,sizeof(dis));
    memset(visq,0,sizeof(visq));
    memset(pre,-1,sizeof(pre));
    int u;queue<int> q;
    q.push(s),dis[s]=0;visq[s] = true;
    while(!q.empty()){
        u=q.front(),q.pop();
        visq[u]=false;
        for(int i=head[u];~i;i=ed[i].next){
          int v=ed[i].to;
          if(ed[i].v && dis[v] > dis[u]+ed[i].c){
            dis[v] = dis[u]+ed[i].c;
            pre[v] = i;
            if(!visq[v]){
              q.push(v);
              visq[v] = true;
            }
          }
        }
    }
    if(pre[t] == -1) return false;
    return true;
}

int n,m;
int mcmf(int s,int t){//返回的是最大流量
    int res = 0;
    while(spfa(s,t)){
        if(dis[t] >= 0) break;
        int f = 0x3f3f3f3f;
        for(lli i = pre[t];~i;i = pre[ed[i^1].to]){
            f = min(f,ed[i].v);
        }
        for(lli i = pre[t];~i;i = pre[ed[i^1].to]){
            ed[i].v -= f;ed[i^1].v += f;
        }
        ans += f*dis[t];
        res += f;
    }
    return res;
}
void ini(){
    memset(head,-1,sizeof(head));cnte = ans = 0;
}
int main(){
    while(~scanf("%d%d",&n,&m)){
        ini();
        int a,b,c,d;
        int s = 0,t = n+1;
        for(int i = 1;i <= n;i++){
            scanf("%d%d%d%d",&a,&b,&c,&d);
            ae(s,i,b,a);
            ae(i,t,d,-c);
        }
        for(int i = 1;i <= m;i++){
            scanf("%d%d%d",&a,&b,&c);
            ae(a,b,inf,c);
            ae(b,a,inf,c);
        }
        mcmf(s,t);
        printf("%d\n",-ans);
    }
    return 0;
}