bzoj3597[Scoi2014]方伯伯运椰子 01分数规划+spfa判负环

3597: [Scoi2014]方伯伯运椰子

Time Limit: 30 Sec  Memory Limit: 64 MB
Submit: 594  Solved: 360
[Submit][Status][Discuss]

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

 

很容易就可以想到01分数规划，然后思考怎么判断是否有可行解
一次完整的修改应该是找两条路径，一条路径容量扩大1，流量扩大1，另一条流量-1，容量缩小1
如果第一条路径上的边为e1，第二条路径上的边为e2，代价就是 $\sum{a_{e1}}+\sum{d_{e1}}-\sum{d_{e2}}+\sum{b_{e2}}$

可以发现这两条路径除掉前面的公共路径之后可以形成一个无向环，并且一条边的扩容、缩容付出的代价是独立的
得到一个思路：建立双向边，正向边权是扩容代价，如果容量上限不为0，反向边权是缩容代价，否则不建反向边

由于题目的答案分式和变化的容量无直接关系，所以容量变化1和变化x是没有区别的，直接检查容量变化1是否可行
01规划判断的时候，把每个环的权值定义为 $边数*mid+\sum边权$
把每条边权值加上mid后找负环，如果存在负环就有可行解，否则没有

#include<bits/stdc++.h>
#define N 5010
using namespace std;
const double eps=1e-3;
int n,m,tot,fg,hd[N],vis[N];double d[N];
struct edge{int v,w,next;}e[N<<1];
void adde(int u,int v,int w){
    e[++tot].v=v;
    e[tot].next=hd[u];
    e[tot].w=w;
    hd[u]=tot;
}
void dfs(int u,double x){
    vis[u]=1;
    if(fg)return;
    for(int i=hd[u];i;i=e[i].next){
        int v=e[i].v;
        if(d[v]>d[u]+x+e[i].w){
            d[v]=d[u]+x+e[i].w;
            if(vis[v]){fg=1;return;}
            dfs(v,x);
        }
    }
    vis[u]=0;
}
bool check(double x){
    fg=0;memset(d,0,sizeof(d));
    memset(vis,0,sizeof(vis));
    for(int i=1;i<=n&&!fg;i++)dfs(i,x);
    return fg;
}
int main(){
    scanf("%d%d",&n,&m);
    n+=2;int u,v,a,b,c,d;
    for(int i=1;i<=m;i++){
        scanf("%d%d%d%d%d%d",&u,&v,&a,&b,&c,&d);
        adde(u,v,b+d);if(c)adde(v,u,a-d);
    }
    double l=0,r=1e6,mid,ans;
    while(l+eps<=r){
        mid=(l+r)/2;
        if(check(mid))l=ans=mid;
        else r=mid;
    }
    printf("%.2lf\n",ans);
    return 0;
}