zoukankan      html  css  js  c++  java
  • H Minimum-cost Flow(预处理增广费用)

    题:https://ac.nowcoder.com/acm/contest/5666/H

    题意:给定网络,然后每次询问给定每条边的容量ui/vi,问从节点1到节点n在流过 1 的情况下的最小费用

    分析:把图的所有边容量设为1跑费用流,因为这样的话,在ui,vi条件下,可以等价为每条边乘上vi/ui,让每条边固定为1的容量;

       而满足条件就是要增广vi/ui(在ui的“进制”下向上取整)次,因为所有边容量相同,而所以每次增广出来的容量是一样的为ui/vi,且费用依次递增;

       所以只要预处理费用的前缀和即可(这也是预处理是设流量为1的好处);

       计算时,只要vi/ui(向下取整的部分)*ui 和余下来的部分*ui即为分子部分

    #include<bits/stdc++.h>
    using namespace std;
    #define pb push_back
    #define MP make_pair
    #define lson root<<1,l,midd
    #define rson root<<1|1,midd+1,r
    typedef long long ll;
    const int mod=1e9+7;
    const int M=105;
    const int inf=0x3f3f3f3f;
    struct node{
        int u,v,w,cost,nextt;
    }e[30000];
    int mincost;
    int head[M],cur[M],dis[M],vis[M],tmp[M*100],tot,s,t;
    void addedge(int u,int v,int w,int cost){
        e[tot].u=u;
        e[tot].v=v;
        e[tot].w=w;
        e[tot].cost=cost;
        e[tot].nextt=head[u];
        head[u]=tot++;
        e[tot].u=v;
        e[tot].v=u;
        e[tot].w=0;
        e[tot].cost=-cost;
        e[tot].nextt=head[v];
        head[v]=tot++;
    }
    bool bfs(){
        for(int i=0;i<=t;i++)
            dis[i]=inf;
        queue<int>que;
        que.push(s);
        dis[s]=0;
        while(!que.empty()){
            int u=que.front();
            que.pop();
            vis[u]=0;
            for(int i=head[u];~i;i=e[i].nextt){
                int v=e[i].v;
                if(e[i].w&&dis[u]+e[i].cost<dis[v]){
                    dis[v]=dis[u]+e[i].cost;
                    if(!vis[v]){
                        vis[v]=1;
                        que.push(v);
                    }
                }
            }
        }
        return dis[t]!=inf;
    }
    int dfs(int u,int fl){
        if(u==t)
            return fl;
        int ans=0;
        vis[u]=1;
        for(int i=cur[u];~i;i=e[i].nextt){
            int v=e[i].v;
            if(dis[v]==dis[u]+e[i].cost&&!vis[v]&&e[i].w){
                cur[u]=i;
                int x=dfs(v,min(e[i].w,fl-ans));
                e[i].w-=x;
                e[i^1].w+=x;
                ans+=x;
                mincost+=x*e[i].cost;
                if(ans==fl)
                    break;
            }
        }
        vis[u]=0;
        return ans;
    }
    int MCMF(){
        int ans=0;
        while(bfs()){
            for(int i=0;i<=t;i++)
                cur[i]=head[i];
            ans+=dfs(s,inf);
            tmp[ans]=tmp[ans-1]+mincost;
            mincost=0;
        }
        return ans;
    }
    int main(){
        int n,m;
        while(~scanf("%d%d",&n,&m)){
            memset(head,-1,sizeof(head));
            s=1,t=n;
            tot=0,mincost=0;
            for(int u,v,cost,i=1;i<=m;i++){
                scanf("%d%d%d",&u,&v,&cost);
                addedge(u,v,1,cost);
            }
            int flow=MCMF();
            int q;
            scanf("%d",&q);
            while(q--){
                ll u,v;
                scanf("%lld%lld",&u,&v);
                if(u==0){
                    puts("NaN");
                    continue;
                }
                ll need=(v+(u-1))/u;///增广次数
                if(need>flow)///v/u>flow
                    puts("NaN");
                else{
                    ll UP=1ll*tmp[need-1]*u;
                    UP+=1ll*(tmp[need]-tmp[need-1])*(v-(need-1)*u);
                    ll DOWM=v;
                    ll g=__gcd(UP,DOWM);
                    printf("%lld/%lld
    ",UP/g,DOWM/g);
                }
            }
        }
        return 0;
    }
    View Code
  • 相关阅读:
    #Leetcode# 204. Count Primes
    #Leetcode# 202. Happy Number
    #Leetcode# 8. String to Integer (atoi)
    VMWare(虚拟机) 网络配置模式 :桥接模式、nat模式、主机模式
    领域的建立
    需求获取的前期工作(不断更新)
    2.4文件查找与文件列表
    2.3录制并回放终端会话
    2命令之乐,2.1简介
    1.7比较与测试
  • 原文地址:https://www.cnblogs.com/starve/p/13330049.html
Copyright © 2011-2022 走看看