「SCOI2014」方伯伯运椰子 解题报告

「SCOI2014」方伯伯运椰子

可以看出是分数规划

然后我们可以看出其实只需要改变1的流量就可以了，因为每次改变要保证流量守恒，必须流成一个环，在正负性确定的情况下，变几次是无所谓的。

然后按照套路，设

\[ans=\frac{X-Y}{k}\\
ans\times k =X-Y\\
ans\times k=-\sum w_i\\
\sum ans-w_i=0
\]

从第二部到第三步是把X和Y中的共同边都减掉了

\(w\)是根据扩容或者缩容建的边权为\(b+d,a-d\)的边权集合

注意一点，缩小容量必须\(c_i>0\)

然后发现环的边数就是\(k\)，减过去就可以二分ans了

---

Code:

#include <cstdio>
#include <cctype>
#include <cstring>
template <class T>
void read(T &x)
{
	x=0;char c=getchar();
	while(!isdigit(c)) c=getchar();
	while(isdigit(c)) x=x*10+c-'0',c=getchar();
}
const int N=5010;
const int M=6010;
int n,m,u[N],v[N],a[N],b[N],c[N],d[N];
int head[N],to[M],Next[M],cnt;double edge[M];
void add(int u,int v,double w)
{
	to[++cnt]=v,edge[cnt]=w,Next[cnt]=head[u],head[u]=cnt;
}
int used[N],vis[N],q[N*N],l,r;double dis[N];
bool spfa()
{
	for(int i=1;i<=n+2;i++) dis[i]=1e12;
	memset(vis,0,sizeof vis);
	memset(used,0,sizeof used);
	dis[q[l=r=1]=n+1]=0;
	while(l<=r)
	{
		int now=q[l++];
		used[now]=0;
		for(int v,i=head[now];i;i=Next[i])
			if(dis[v=to[i]]>dis[now]+edge[i])
			{
				dis[v]=dis[now]+edge[i];
				if((++vis[v])==n+2) return true;
				if(!used[v]) used[q[++r]=v]=1;
			}
	}
	return false;
}
bool check(double x)
{
	memset(head,0,sizeof head),cnt=0;
	for(int i=1;i<=m;i++)
	{
	    if(u[i]!=n+1)
		{
		    if(c[i]!=0) add(v[i],u[i],x+a[i]-d[i]);
		    add(u[i],v[i],x+b[i]+d[i]);
		}
        else
            add(u[i],v[i],0);
	}
	return spfa();
}
int main()
{
	read(n),read(m);
	double l=0,r=0;
	for(int i=1;i<=m;i++)
	{
		read(u[i]),read(v[i]),read(a[i]);
		read(b[i]),read(c[i]),read(d[i]);
		r+=1.0*c[i]*d[i];
	}
	while(l+1e-6<r)
	{
		double mid=(l+r)/2;
		if(check(mid)) l=mid;
		else r=mid;
	}
	printf("%.2lf\n",l);
	return 0;
}

---

2019.2.24