传送门
这个题好神奇呀
我一开始看错题了,然后就一直在想怎么做,然后看了题解发现自己的思路全错了。。
题解是二分+最大流
其实我们可以先从第二个人开始考虑
假如已经确定了一个最大流方案,那么他要怎么分配费用呢?
由于总费用固定,那么我们就要让这些费用影响到的流量尽可能地多
所以,显然的结论就是:最优策略就是将所有的流量都放在最大流量的边上
所以第一个人要让费用最小,那么其实就是要让流量最大的边流量最小
这样就是一个很经典的二分的问题了
我们只需要二分答案,然后重置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;
}