bzoj3597 [Scoi2014]方伯伯运椰子 01分数规划

来源:互联网 发布:打谱软件 mac 编辑:程序博客网 时间:2024/05/16 01:19

题意:给你一个满流的残量网络,让你通过调整以后,使得调整前后的平均费用之差(调整后比调整前小)最大,同时要求调整后仍然满流。
其实这题不难,只不过网上的题解大多数都模棱两可,写法也参差不一,搞得我很蛋疼。。
事实上我们可以发现,当前给我们的图肯定不是最小流,否则无法增广。
然后,根据费用流的消圈定理,我们只用增广一个负环就好了,所以用01分数规划一下,然后spfa判断是否有负环。

根据01分数规划,本题目标是:最大化λ=(X-Y)/k
则有f(λ)=λk+Y-X
扩边一次的费用为b[i]+d[i]+λ
缩边一次的费用为a[i]-d[i]+λ
每一次二分的时候重新连一次就好了。

#include<cstdio>#include<algorithm>#include<cstring>#define fo(i,a,b) for(int i=a;i<=b;i++)#define fd(i,a,b) for(int i=a;i>=b;i--)using namespace std;const int N=3e5+5;const double eps=1e-3;int n,m;int head[N],go[N],tot,next[N];double dis[N],val[N];int vis[N]; struct node{    int a,b,c,d,u,v;}e[N];inline void add(int x,int y,double z){    go[++tot]=y;    next[tot]=head[x];    val[tot]=z;    head[x]=tot;}inline bool dfs(int x){    vis[x]=1;    for(int i=head[x];i;i=next[i])    {        int v=go[i];        if (dis[x]+val[i]<dis[v])        {            if(vis[v])return 1;            dis[v]=dis[x]+val[i];            if (dfs(v))return 1;        }    }    vis[x]=0;    return 0;}bool pd(double x){    memset(head,0,sizeof(head));    tot=0;    fo(i,1,m)    {        add(e[i].u,e[i].v,e[i].b+e[i].d+x);        if (e[i].c)        add(e[i].v,e[i].u,e[i].a-e[i].d+x);    }    fo(i,1,n+2)dis[i]=vis[i]=0;    fo(i,1,n+2)if (dfs(i))return 1;    return 0;}int main(){    scanf("%d%d",&n,&m);    fo(i,1,m)    {        int u,v,a,b,c,d;        scanf("%d%d%d%d%d%d",&u,&v,&a,&b,&c,&d);        e[i].a=a;        e[i].b=b;        e[i].c=c;        e[i].d=d;        e[i].u=u;        e[i].v=v;    }    double l=0.0,r=1e9;    while (r-l>eps)    {        double mid=(l+r)/2;        if (pd(mid))l=mid;        else r=mid;    }    printf("%.2lf\n",l);} 
阅读全文
0 0
原创粉丝点击