[SDOI2013]费用流

传送门
这个题好神奇呀
我一开始看错题了，然后就一直在想怎么做，然后看了题解发现自己的思路全错了。。
题解是二分+最大流
其实我们可以先从第二个人开始考虑
假如已经确定了一个最大流方案，那么他要怎么分配费用呢？
由于总费用固定，那么我们就要让这些费用影响到的流量尽可能地多
所以，显然的结论就是：最优策略就是将所有的流量都放在最大流量的边上
所以第一个人要让费用最小，那么其实就是要让流量最大的边流量最小
这样就是一个很经典的二分的问题了
我们只需要二分答案，然后重置cap和flow，直接跑一边最大流，判断和原来的最大流是否相同即可
话说我还是没有卡到洛谷rank1，空间没rank1小。。
直接上代码：

#prag\
ma \
GCC \
optimize("O3")
#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstdlib>
#define R register
using namespace std;
char buf[200000],*pt=buf,*o=buf;
inline void read(int &s){
    s=0; while(*pt<'0'||*pt>'9')pt++;
    while(*pt>='0'&&*pt<='9')s=(s<<3)+(s<<1)+*pt++-48;
}
const int N=101,M=2001;
const double eps=1e-5;
int n,m,s,t,p,tot=-1;
int head[N],cur[N],d[N],q[N],to[M],Next[M],Flow[M];
double flow[M];
inline void addedge(int x,int y,int l){
    to[++tot]=y;Next[tot]=head[x];Flow[tot]=l;flow[tot]=l;head[x]=tot;
    to[++tot]=x;Next[tot]=head[y];Flow[tot]=0;flow[tot]=0;head[y]=tot;
}
inline bool bfs(){
    for(R int i=1;i<=n;i++)d[i]=0x3f3f3f3f;
    int l=1,r=1;d[s]=0;q[l]=s;
    while(l<=r){
        int x=q[l++];
        for(R int i=head[x];i!=-1;i=Next[i]){
            int u=to[i];
            if(fabs(flow[i])>eps&&d[u]>d[x]+1){
                d[u]=d[x]+1;
                q[++r]=u;
            }
        }
    }
    return d[t]!=0x3f3f3f3f;
}
inline double dfs(int x,double a){
    if(fabs(a)<eps)return 0;
    if(x==t)return a;
    double F=0,f;
    for(R int &i=cur[x];i!=-1;i=Next[i]){
        int u=to[i];
        if(fabs(flow[i])>eps&&d[u]==d[x]+1&&fabs(f=dfs(u,min(a,flow[i])))>0){
            flow[i]-=f;
            flow[i^1]+=f;
            F+=f;
            a-=f;
            if(fabs(a)<eps)return F;
        }
    }
    return F;
}
inline double dinic(){
    double F=0;
    while(bfs()){
        for(R int i=1;i<=n;i++)cur[i]=head[i];
        F+=dfs(s,0x3f3f3f3f);
    }
    return F;
}
int main(){
    fread(buf,1,200000,stdin);
    memset(head,-1,sizeof(head));
    read(n);read(m);read(p);s=1;t=n;
    int x,y,c;
    for(R int i=1;i<=m;i++){
        read(x);read(y);read(c);
        addedge(x,y,c);
    }
    double ans=dinic();
    printf("%d\n",(int)ans);
    double l=0,r=50001,mid;
    while(fabs(l-r)>eps){
        mid=(l+r)/2;
        for(R int i=0;i<=tot;i++)flow[i]=min((double)Flow[i],mid);
        double num=dinic();
        if(fabs(num-ans)>eps)l=mid;
        else r=mid;
    }
    l*=(double)p;
    printf("%.6lf",l);
    return 0;
}