題目鏈接:點擊查看
題目大意:給出一張有向圖,每條邊都有單獨的花費,現在給出 q 次詢問,每次詢問給出一個 u 和 v 需要回答:將所有邊的流量都設置為 u / v 后,從點 1 到點 n 流量為 1 時的最小花費為多少
題目分析:因為詢問非常多,肯定不能對于每次詢問重新建邊,但因為總的邊數非常小,所以我們利用等比例的特性,初始時先將每條邊的流量都設置為 1 ,在求最小費用最大流的時候,記錄一下每個流量下的最小花費,這樣對于每次詢問,就可以貪心求解了
因為初始時每條邊的流量為 1 ,每次詢問時每條邊的流量為 u / v ,每次詢問時的總流量為 1 ,設初始時對應的流量為 x ,根據比例公式不難求解出需要在初始的圖中使用 v / u 的流量進行對應,貪心的話選前 v / u 的流量所對應的權值和就好了
代碼:
?
#include<iostream>
#include<cstdio>
#include<string>
#include<ctime>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<stack>
#include<climits>
#include<queue>
#include<map>
#include<set>
#include<sstream>
#include<cassert>
using namespace std;typedef long long LL;typedef unsigned long long ull;const int inf=0x3f3f3f3f;const int N=210;vector<LL>node,sum;struct Edge
{int to,w,cost,next;
}edge[N*N];int head[N],cnt;void addedge(int u,int v,int w,int cost)
{edge[cnt].to=v;edge[cnt].w=w;edge[cnt].cost=cost;edge[cnt].next=head[u];head[u]=cnt++;edge[cnt].to=u;edge[cnt].w=0;edge[cnt].cost=-cost;edge[cnt].next=head[v];head[v]=cnt++;
}int d[N],incf[N],pre[N];bool vis[N];bool spfa(int s,int t)
{memset(d,inf,sizeof(d));memset(vis,false,sizeof(vis));memset(pre,-1,sizeof(pre));queue<int>q;q.push(s);vis[s]=true;incf[s]=inf;d[s]=0;while(!q.empty()){int u=q.front();q.pop();vis[u]=false;for(int i=head[u];i!=-1;i=edge[i].next){int v=edge[i].to;int w=edge[i].w;int cost=edge[i].cost;if(!w)continue;if(d[v]>d[u]+cost){d[v]=d[u]+cost;pre[v]=i;incf[v]=min(incf[u],w);if(!vis[v]){vis[v]=true;q.push(v);}}}}return pre[t]!=-1;
}int update(int s,int t)
{int x=t;while(x!=s){int i=pre[x];edge[i].w-=incf[t];edge[i^1].w+=incf[t];node.back()+=edge[i].cost;x=edge[i^1].to;}return d[t]*incf[t];
}void init()
{node.clear();sum.clear();memset(head,-1,sizeof(head));cnt=0;
}int solve(int st,int ed)
{int ans=0;while(spfa(st,ed)){node.push_back(0);update(st,ed);ans++;}return ans;
}int main()
{
#ifndef ONLINE_JUDGE
// freopen("input.txt","r",stdin);
// freopen("output.txt","w",stdout);
#endif
// ios::sync_with_stdio(false);int n,m,st=N-1,ed=st-1;while(scanf("%d%d",&n,&m)!=EOF){init();while(m--){int u,v,c;scanf("%d%d%d",&u,&v,&c);addedge(u,v,1,c);}addedge(st,1,inf,0);addedge(n,ed,inf,0);int max_flow=solve(st,ed);sum.push_back(0);for(int i=0;i<node.size();i++)sum.push_back(sum.back()+node[i]);int q;scanf("%d",&q);while(q--){LL u,v;scanf("%lld%lld",&u,&v);if(u*max_flow<v){puts("NaN");continue;}int j=v/u;LL ans1=sum[j]*u+(sum[j+1]-sum[j])*(v-j*u);LL ans2=v;LL gcd=__gcd(ans1,ans2);ans1/=gcd;ans2/=gcd;printf("%lld/%lld\n",ans1,ans2);}}return 0;
}
?
超強干貨來襲 云風專訪:近40年碼齡,通宵達旦的技術人生
總結
以上是生活随笔為你收集整理的牛客多校 - Minimum-cost Flow(最小费用最大流+贪心)的全部內容,希望文章能夠幫你解決所遇到的問題。
如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。