【BZOJ】3130: [Sdoi2013]费用流

传送门：http://www.lydsy.com/JudgeOnline/problem.php?id=3130

费用明显就是堆到流最大那条边上即可。
【然后想都没想就二分最大边冲上去了】

后来感觉奇怪，才发现流不一定是整数

红色容量为2
蓝色容量为1
最大流量边就是1.5了

所以做法就是二分最大边，然后跑实数的流即可
【我输漏了一个数调了一天QwQ】

#include<stdio.h>
#include<cstring>
#include<algorithm>
#define M 2005
#define N 105
#define eps 0.000001
#define cint const int &
using namespace std;

int tot=1,n,m,p,s[N],tm[N],ts,Q[N],lev[N],cur[N];
double L,R;

struct edge{int v,n;double c;}E[M],e[M];

inline void push(cint u,cint v,cint c)
{
    E[++tot]=(edge){v,s[u],c};s[u]=tot;
    E[++tot]=(edge){u,s[v],0};s[v]=tot;
}

inline bool bfs()
{
    memcpy(cur,s,sizeof(cur));
    tm[Q[1]=1]=++ts;
    for (int l=1,r=1;l<=r;l++) for (int i=s[Q[l]];i;i=e[i].n) if (eps<e[i].c && tm[e[i].v]<ts)
        tm[Q[++r]=e[i].v]=ts,lev[Q[r]]=lev[Q[l]]+1;
    return tm[n]==ts;
}

double dfs(cint u,double c)
{
    if (u==n) return c;
    double g=0,tmp;
    for (int i=cur[u];eps<c && i;cur[u]=i,i=e[i].n) 
      if (tm[e[i].v]==ts && lev[e[i].v]==lev[u]+1 && eps<e[i].c && (tmp=dfs(e[i].v,min(c,e[i].c))))
        g+=tmp,c-=tmp,e[i].c-=tmp,e[i^1].c+=tmp;
    return g;
}

inline double Judge(const double &k)
{
    memcpy(e,E,sizeof(e));
    for (int i=2;i<=tot;i++) if (k<e[i].c) e[i].c=k;
    double ans=0;
    while (bfs()) for (double rt=dfs(1,233333333);eps<rt;rt=dfs(1,233333333)) ans+=rt;
    return ans;
}

int main()
{
    scanf("%d%d%d",&n,&m,&p);
    for (int u,v,c;m--;) scanf("%d%d%d",&u,&v,&c),push(u,v,c),R=max((double)c,R);
    double lim=Judge(R);
    for (double mid=(L+R)/2;L+eps<R;mid=(L+R)/2) if (lim<Judge(mid)+eps) R=mid;else L=mid;
    printf("%.0lf\n%.4lf",lim+eps,L*p);
}