BZOJ3597: [Scoi2014]方伯伯运椰子

来源：互联网发布：龙江网络宽带客服电话编辑：程序博客网时间：2024/04/28 18:16

Description

Input

第一行包含二个整数N，M

接下来M行代表M条边，表示这个交通网络

每行六个整数，表示Ui,Vi,Ai,Bi,Ci,Di

接下来一行包含一条边，表示连接起点的边

Output

一个浮点数，保留二位小数。表示答案，数据保证答案大于0

Sample Input

5 10
1 5 13 13 0 412
2 5 30 18 396 148
1 5 33 31 0 39
4 5 22 4 0 786
4 5 13 32 0 561
4 5 3 48 0 460
2 5 32 47 604 258
5 7 44 37 75 164
5 7 34 50 925 441
6 2 26 38 1000 22

Sample Output

103.00

HINT

1<=N<=5000

0<=M<=3000

1<=Ui,Vi<=N+2

0<=Ai,Bi<=500

0<=Ci<=10000

0<=Di<=1000

Source

显然答案满足二分性，思考怎么check

如果存在一种费用更小的流法，那么一定会形成一个环

例如

+1->+1

↗ ↘

↘ ↗

-1->-1

这样就流量平衡，那么我们把增容的边弄成正向，减容的边弄成反向，那么一种最优的做法一定是一个负权环

所有边减去二分的答案后就是喜闻乐见的找负环了

#include <bits/stdc++.h>using namespace std;const int MAXN = 5200;const int MAXM = 6666;const double eps = 1e-3;double l, r, ans;int n, m, cnt, head[MAXN];struct edge{int to, nxt;double val;}e[MAXM];inline void addedge(int x, int y, double w) { e[ ++cnt ].to = y; e[ cnt ].nxt = head[ x ]; head[ x ] = cnt; e[ cnt ].val = w; }double dis[MAXN];bool vis[MAXN];inline bool dfs(int x){vis[ x ] = 1;for( int i = head[ x ] ; i ; i = e[ i ].nxt )if( dis[ e[ i ].to ] > dis[ x ] + e[ i ].val ){dis[ e[ i ].to ] = dis[ x ] + e[ i ].val;if( vis[ e[ i ].to ] ) return true;if( dfs( e[ i ].to ) ) return true;}vis[ x ] = 0;return false;}inline bool chk(){for( int i = 1 ; i <= n ; i++ ) dis[ i ] = 0, vis[ i ] = 0;for( int i = 1 ; i <= n ; i++ ) if( dfs( i ) ) return true;return false;}int main(){scanf( "%d%d", &n, &m );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 );addedge( u, v, b + d );if( c ) addedge( v, u, a - d );}l = 0, r = 1e9;while( r - l > eps ){double mid = ( l + r ) / 2.0;for( int i = 1 ; i <= cnt ; i++ ) e[ i ].val += mid;if( chk() ) l = mid;else r = mid;for( int i = 1 ; i <= cnt ; i++ ) e[ i ].val -= mid;}return printf( "%.2lf\n", r ), 0;}

0 0