【bzoj3130】 SDOI2013费用流 最大流

kac出的好题，第一问裸最大流肯定没有问题，第二问我们首先可以考虑一个结论这些费用一定是只放在一条边上的，那么对于一个最大流方案我们肯定会选择放到流量最大的那条边上，那么小A就会想办法让自己流量最大的那条边最小，于是我们就可以二分一下答案，把所有边的流量全部替换为二分的答案，如果最大流发生变化就说明一定存在流量比二分的答案大的边，于是问题就解决了。不过这道题需要实数二分，不能整数二分，为什么呢？

（1,3,3）

（1,2,3）

（2,3,3）

（3,4,2）

（3,5,2）

（3,6,1）

（4,7,2）

（5,7,2）

（6,7,2）

最大边最小是2.5

kac出得一手神题！！！

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
#include<algorithm>
#define maxn 20010
#define inf 1000000000
#define eps 1e-6

using namespace std;

struct yts
{
	int x,y;
}a[maxn];

int to[maxn],next[maxn],head[maxn],d[maxn];
double c[maxn];
int q[maxn];
int n,m,p,ans,num,s,t;

bool bfs()
{
	int l=0,r=1;
	memset(d,-1,sizeof(d));
	q[1]=s;d[s]=0;
	while (l<r)
	{
		int x=q[++l];
		for (int p=head[x];p;p=next[p])
		  if (c[p] && d[to[p]]==-1)
		  {
		  	d[to[p]]=d[x]+1;
		  	q[++r]=to[p];
		  }
	}
	if (d[t]==-1) return 0; else return 1;
}

double find(int x,double low)
{
	if (x==t || low==0) return low;
	double totflow=0;
	for (int p=head[x];p;p=next[p])
	  if (d[to[p]]==d[x]+1 && c[p])
	  {
	  	double a=find(to[p],min(low,c[p]));
	  	low-=a;totflow+=a;
	  	c[p]-=a;c[p^1]+=a;
	  	if (low==0) return totflow;
	  }
	if (low>0) d[x]=-1;
	return totflow;
}

double Dinic()
{
	double ans=0;
	while (bfs()) ans+=find(s,inf);
	return ans;
}

void addedge(int x,int y,double z)
{
	num++;to[num]=y;c[num]=z;next[num]=head[x];head[x]=num;
	num++;to[num]=x;c[num]=0;next[num]=head[y];head[y]=num;
}

void build(double x)
{
	num=1;
	memset(head,0,sizeof(head));
	for (int i=1;i<=m;i++) addedge(a[i].x,a[i].y,x);
}

int main()
{
	scanf("%d%d%d",&n,&m,&p);num=1;s=1,t=n;
	for (int i=1;i<=m;i++)
	{
		int x,y,z;
		scanf("%d%d%d",&x,&y,&z);
		addedge(x,y,(double)z);
		a[i].x=x;a[i].y=y;
	}
	ans=(int)Dinic();
	printf("%d\n",ans);
	double l=0,r=50000,ans1;
	while (l+eps<r)
	{
		double mid=(l+r)/2;
		build(mid);
		if (Dinic()<(double)ans) ans1=mid,l=mid; else r=mid;
	}
	printf("%.4lf\n",(double)ans1*p);
	return 0;
}