hdu6188&&百度之星初赛（B） T5

度度熊的交易计划

Problem Description

度度熊参与了喵哈哈村的商业大会，但是这次商业大会遇到了一个难题：

喵哈哈村以及周围的村庄可以看做是一共由n个片区，m条公路组成的地区。

由于生产能力的区别，第i个片区能够花费a[i]元生产1个商品，但是最多生产b[i]个。

同样的，由于每个片区的购买能力的区别，第i个片区也能够以c[i]的价格出售最多d[i]个物品。

由于这些因素，度度熊觉得只有合理的调动物品，才能获得最大的利益。

据测算，每一个商品运输1公里，将会花费1元。

那么喵哈哈村最多能够实现多少盈利呢？

Input

本题包含若干组测试数据。 每组测试数据包含： 第一行两个整数n,m表示喵哈哈村由n个片区、m条街道。 接下来n行，每行四个整数a[i],b[i],c[i],d[i]表示的第i个地区，能够以a[i]的价格生产，最多生产b[i]个，以c[i]的价格出售，最多出售d[i]个。 接下来m行，每行三个整数,u[i],v[i],k[i]，表示该条公路连接u[i],v[i]两个片区，距离为k[i]

可能存在重边，也可能存在自环。

满足: 1<=n<=500, 1<=m<=1000, 1<=a[i],b[i],c[i],d[i],k[i]<=1000, 1<=u[i],v[i]<=n

Output

输出最多能赚多少钱。

Sample Input

2 1
5 5 6 1
3 5 7 7
1 2 1

Sample Output

23
———————————————————————————————————————————————
这道题是道很明显的网络流 建一波图就好了 注意是最小费用流不是最小费用最大流 当d【T】>=0时就没有必要继续跑下去了

#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long
using namespace std;
const int M=,inf=0x3f3f3f3f;
int read(){
    int ans=,f=,c=getchar();
    while(c<''||c>'') {if(c=='-') f=-; c=getchar();}
    while(c>=''&&c<='') {ans=ans*+(c-''); c=getchar();}
    return ans*f;
}
int n,m,S,T,ans;
int first[M],cur[M],cnt,d[M],q[],vis[M];
struct node{int to,next,flow,cost;}e[];
void insert(int a,int b,int flow,int cost){
    e[++cnt]=(node){b,first[a],flow,cost}; first[a]=cnt;
    e[++cnt]=(node){a,first[b],,-cost}; first[b]=cnt;
}
int spfa(){
    memset(d,0x3f,sizeof(d));
    int head=,tail=;
    q[]=S;  d[S]=;
    while(head!=tail){
        int x=q[head++]; vis[x]=; if(head>) head=;
        for(int i=first[x];i;i=e[i].next){
            int now=e[i].to;
            if(e[i].flow&&d[now]>d[x]+e[i].cost){
                d[now]=d[x]+e[i].cost;
                if(!vis[now]){q[tail++]=now; vis[now]=; if(tail>) tail=;}
            }
        }
    }
    return d[T]<;
}
int dfs(int x,int a){
    if(x==T||a==) return a;
    vis[x]=;
    int f,flow=;
    for(int &i=cur[x];i;i=e[i].next){
        int now=e[i].to;
        if(!vis[now]&&d[now]==d[x]+e[i].cost&&(f=dfs(now,min(a,e[i].flow)))>){
            e[i].flow-=f; e[i^].flow+=f;
            flow+=f; ans+=e[i].cost*f;
            a-=f; if(a==) break;
        }
    }
    vis[x]=;
    return flow;
}
int main()
{
    int x,y,w;
    while(scanf("%d %d",&n,&m)==){
        memset(first,,sizeof(first));
        ans=; cnt=;
        S=; T=n+;
        for(int i=;i<=n;i++){
            x=read(); w=read(); insert(S,i,w,x);
            x=read(); w=read(); insert(i,T,w,-x);
        }
        for(int i=;i<=m;i++){
            x=read(); y=read(); w=read();
            if(x==y) continue;
            insert(x,y,inf,w);
            insert(y,x,inf,w);
        }
        memset(vis,,sizeof(vis));
        while(spfa()) {for(int i=S;i<=T;i++) cur[i]=first[i]; dfs(S,inf);}
        printf("%d\n",-ans);
    }
    return ;
}