題目描述
四川的方伯伯為了致富,決定引進海南的椰子樹。方伯伯的椰子園十分現代化,椰子園中有一套獨特的交通系統。
現在用點來表示交通節點,邊來表示道路。這樣,方伯伯的椰子園就可以看作一個有 n + 2 個交通節點,m條邊的有向無環圖。n +1 號點為入口,n +2 號點為出口。每條道路都有 6 個參數,ui,vi,ai,bi,ci,di,分別表示,該道路從 ui 號點通向 vi 號點,將它的容量壓縮一次要 ai 的花費,容量擴大一次要 bi 的花費,該條道路當前的運輸容量上限為 ci,并且每單位運輸量通過該道路要 di 的費用。
在這個交通網絡中,只有一條道路與起點相連。因為弄壞了這條道路就會導致整個交通網絡癱瘓,聰明的方伯伯決定絕不對這條道路進行調整,也就是說,現在除了這條道路之外,對其余道路都可以進行調整。
有兩種調整方式:
選擇一條道路,將其進行一次壓縮,這條道路的容量會下降 1 單位。選擇一條道路,將其進行一次擴容,這條道路的容量會上升 1 單位。一條道路可以被多次調整。
由于很久以前,方伯伯就請過一個工程師,對這個交通網絡進行過一次大的優化調整。所以現在所有的道路都被完全的利用起來了,即每條道路的負荷都是滿的(每條道路的流量等于其容量)。
但方伯伯一想到自己的海南椰子會大豐收,就十分擔心巨大的運輸量下,會導致過多的花費。因此,方伯伯決定至少進行一次調整,調整之后,必須要保持每條道路滿負荷,且總交通量不會減少。
設調整后的總費用是 Y,調整之前的總費用是 X。現在方伯伯想知道,最優調整比率是多少,即假設他進行了 k 次調整,(X - Y)/k最大能是多少?
注:總費用 = 交通網絡的運輸花費 + 調整的花費 輸入輸出格式 輸入格式:
第一行包含二個整數N,M接下來M行代表M條邊,表示這個交通網絡每行六個整數,表示Ui,Vi,Ai,Bi,Ci,Di接下來一行包含一條邊,表示連接起點的邊
輸出格式:
一個浮點數,保留二位小數。表示答案,數據保證答案大于0
分析: 1.因為本題要求的是分數的最大值,所以要用到分數規劃 2.因為與源點相連的只有一條邊,所以這個圖的流量是守恒的(總量不會變化) 3.那么壓縮相當于退流,擴容相當于增廣,增廣相當于加了一條。 4.擴容1的費用為b+d,壓縮1的費用為a-d 5.化一下給的式子: x?yk>mid x?y?k?mid>0 y?x+k?mid<0 y-x即是擴容的費用加上壓縮的費用。由于最后要除上操作總數k,因此對于相同的邊操作多次是沒有意義的。又因為存在mid的影響,所以每次加上mid再判斷是否有負權環即可。
#include<cstdio>#include<algorithm>#include<cstring>#include<cmath>#include<queue>using namespace std;const int maxn=5010;const int maxm=3010;const int INF=1e9;int to[maxm*2],Next[maxm*2],w[maxm*2],Begin[maxn],e;int n,m;double dis[maxn];void add(int x,int y,int z){ to[++e]=y; Next[e]=Begin[x]; Begin[x]=e; w[e]=z;}int inq[maxn],cnt[maxn];int s,t;bool SPFA(double add){ memset(inq,0,sizeof(inq)); memset(cnt,0,sizeof(cnt)); for(int i=1;i<=n;i++) dis[i]=INF*1.0; queue<int>q; q.push(s);inq[s]=1;cnt[s]++; dis[s]=0; while(!q.empty()){ int u=q.front();q.pop(); for(int i=Begin[u];i;i=Next[i]){ int v=to[i]; if(cnt[v]>n) return true; double length=w[i]+add; if(dis[v]>dis[u]+length){ dis[v]=dis[u]+length; if(!inq[v]) q.push(v),cnt[v]++; } } } return false;}int main(){ scanf("%d%d",&n,&m); double L=0,R=0; s=n+1,t=n+2; n+=2; for(int i=1;i<=m;i++){ int u,v,a,b,c,d; scanf("%d%d%d%d%d%d",&u,&v,&a,&b,&c,&d); if(c) add(v,u,a-d); add(u,v,b+d); if(a-d<0) R+=d-a; } double ans=0; while(R-L>=1e-3){ double mid=(L+R)/2; if(SPFA(mid)){ ans=mid;L=mid; }else R=mid; } ^_^