BZOJ 3130: [Sdoi2013]费用流 网络流 二分 最大流

https://www.lydsy.com/JudgeOnline/problem.php?id=3130

本来找费用流的题，权当复习一下网络流好了。

有点麻烦的是double，干脆判断大小或者二分增加下限都用eps=1e-8操作好了（毕竟只要求精确到4位）。

我普通最大流都快忘了，板子写错了一次超时了。

网络流板子的细节要记清楚：1.增广的时候访问完哪个点就把dep标记改为-1防止同一次增广再次访问。2.往下层传递的时候用val-cnt而不是val。3.cnt=val时及时返回。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
using namespace std;
const int maxn=1010;
const double eps=1e-8;
int n,m,p;double mx;
struct nod{
    int y,rev,next;double v,v1;
}e[maxn*2];
int head[maxn]={},tot=0;double Ans=0;
int num[maxn]={},vis[maxn*2]={};
queue<int>q;
void init(int x,int y,double v){
    e[++tot].y=y;e[tot].v=v;e[tot].rev=tot+1;e[tot].next=head[x];head[x]=tot;
    e[++tot].y=x;e[tot].v=0;e[tot].rev=tot-1;e[tot].next=head[y];head[y]=tot;
}
bool fir(){
    memset(num,-1,sizeof(num));
    q.push(1);num[1]=0;
    while(!q.empty()){
        int x=q.front();q.pop();
        for(int i=head[x];i;i=e[i].next){
            if(e[i].v1<eps)continue;
            if(num[e[i].y]==-1){
                num[e[i].y]=num[x]+1;
                q.push(e[i].y);
            }
        }
    }
    return num[n]!=-1;
}
double dfs(int x,double val){
    if(x==n)return val;
    double cnt=0,z;
    for(int i=head[x];i;i=e[i].next){
        if((e[i].v1<eps)||num[e[i].y]!=num[x]+1)continue;
        z=dfs(e[i].y,min(val-cnt,e[i].v1));
        cnt+=z; e[i].v1-=z; e[e[i].rev].v1+=z;
        if(val-cnt<eps)return cnt;
    }num[x]=-1;
    return cnt;
}
bool Check(double shu){
    mx=shu;double cnt=0;
    for(int i=1;i<tot;i+=2)e[i].v1=min(mx,e[i].v);
    for(int i=2;i<=tot;i+=2)e[i].v1=0;
    while(fir()){
        cnt+=dfs(1,50000000);
    }
    return Ans-cnt<eps;
}
double erfen(double mx){
    double l=0,r=mx;
    while(r-l>eps){
        double mid=(l+r)/2;
        if(Check(mid))r=mid;
        else l=mid+eps;
    }
    return l;
}
int main(){
    scanf("%d%d%d",&n,&m,&p);
    int x,y,z;
    for(int i=1;i<=m;i++){scanf("%d%d%d",&x,&y,&z);init(x,y,(double)z);mx=max(mx,(double)z);}
    for(int i=1;i<=tot;i++)e[i].v1=e[i].v;
    while(fir())Ans+=dfs(1,mx);
    if(Ans==0){
        printf("0
0.0000
");
    }
    else{
        printf("%.0f
",Ans);
        printf("%.4lf
",erfen(mx)*(double)p);
    }
    return 0;
}