【BZOJ】3130: [Sdoi2013]费用流

来源：互联网发布：微山湖武工队知乎编辑：程序博客网时间：2024/05/29 17:15

传送门：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);}

0 0