luogu3305/bzoj3130 费用流 (二分答案+dinic)

Bob肯定想挑一个流量最大的边，然后把所有的费用都加给它呗

那Alice就让流量最大的边尽量小呗

那就二分一下答案再dinic呗

 #include<bits/stdc++.h>
 #define pa pair<int,int>
 #define CLR(a,x) memset(a,x,sizeof(a))
 using namespace std;
 typedef long long ll;
 const int maxn=,maxm=;
 const double inf=1e9;

 inline ll rd(){
     ll x=;char c=getchar();int neg=;
     while(c<''||c>''){if(c=='-') neg=-;c=getchar();}
     while(c>=''&&c<='') x=x*+c-'',c=getchar();
     return x*neg;
 }

 struct Edge{
     int a,b,ne;
     double l;
 }eg[maxm];
 int egh[maxn],ect=;
 double oril[maxm];
 int N,M,P,S,T;
 int dep[maxn],cur[maxn];
 double maxl;
 queue<int> q;

 inline void adeg(int a,int b,int c){
     eg[++ect].a=a,eg[ect].b=b;oril[ect]=eg[ect].l=c,eg[ect].ne=egh[a];egh[a]=ect;
     eg[++ect].a=b,eg[ect].b=a;oril[ect]=eg[ect].l=,eg[ect].ne=egh[b];egh[b]=ect;
 }

 inline bool bfs(){
     CLR(dep,);CLR(cur,-);
     dep[S]=;q.push(S);
     while(!q.empty()){
         int p=q.front();q.pop();
         for(int i=egh[p];i;i=eg[i].ne){
             int b=eg[i].b;
             if(!dep[b]&&eg[i].l){
                 dep[b]=dep[p]+;
                 q.push(b);
             }
         }
     }
     return dep[T]!=;
 }

 double dinic(int x,double y){
     if(x==T) return y;
     double tmp=y;
     if(cur[x]==-) cur[x]=egh[x];
     for(int &i=cur[x];i;i=eg[i].ne){
         int b=eg[i].b;
         if(dep[b]!=dep[x]+||!eg[i].l) continue;
         double re=dinic(b,min(eg[i].l,tmp));
         tmp-=re,eg[i].l-=re,eg[i^].l+=re;
         if(tmp<=1e-) break;
     }return y-tmp;
 }

 inline bool judge(double m){
     for(int i=;i<=ect;i++)
         eg[i].l=min(oril[i],m);
     double l=;
     while(bfs()) l+=dinic(S,inf);
     return fabs(maxl-l)<=1e-;
 }

 int main(){
     //freopen("","r",stdin);
     int i,j,k;
     N=rd(),M=rd(),P=rd();
     S=,T=N;
     for(i=;i<=M;i++){
         int a=rd(),b=rd(),c=rd();
         adeg(a,b,c);
     }
     while(bfs()) maxl+=dinic(S,inf);
     double l=,r=5e4,ans=5e4;
     while(r-l>=1e-){
         double m=(l+r)/;
         if(judge(m)) ans=m,r=m-1e-;
         else l=m+1e-;
     }
     printf("%d\n",(int)maxl);
     printf("%lf\n",ans*P);
     return ;
 }