P2868 [USACO07DEC]Sightseeing Cows G(最优比率环)

luogu传送门

奶牛从任意点出发,经过一些路回到起点

价值是 ∑ f i ∑ t i \frac{\sum f_i}{\sum t_i} tifi

那么按照 01 01 01分数的基本套路来就是

f ( r ) = r ∗ ∑ t i − ∑ f i f(r)=r*\sum t_i-\sum f_i f(r)=rtifi

显然对于不同的环对应不同的直线

二分答案 m i d mid mid,作直线 x = m i d x=mid x=mid相交上述直线若干点

若存在交点的 y y y坐标小于零说明最大值还在 m i d mid mid右边

否则还在左边,所以这题的关键就是去找负环

这样把每个点拆点,彼此连一条权值为 f i f_i fi的边就成功转化为全是边的模型了

然后原图的所有边扩大为原来的 m i d mid mid

但是这样有一个问题,环内有个点的权值会算两次

进去一次出来一次

#include <bits/stdc++.h>
using namespace std;
#define int long long
const int maxn=4e5+10;
const int inf=1e16;
const double eps=1e-6;
struct edge{
	int to,nxt; double w;
}d[maxn]; int head[maxn],cnt=1;
double p[maxn];
void add(int u,int v,double w)
{
	d[++cnt]=(edge){v,head[u],w},head[u]=cnt;
}
double dis[maxn];
int num[maxn],n,m,vis[maxn];
bool spfa(double mid)
{
	for(int i=1;i<=n;i++)	dis[i]=inf,num[i]=0,vis[i]=0;
	num[0]=1; dis[0]=0;
	queue<int>q; q.push(0);
	while( !q.empty() )
	{
		int u=q.front(); q.pop();
		vis[u]=0;
		for(int i=head[u];i;i=d[i].nxt )
		{
			int v=d[i].to;
			double val = d[i].w*mid-p[v];
			if( dis[v]>dis[u]+eps+val )
			{
				dis[v]=dis[u]+val;
				if( !vis[v] )
				{
					vis[v]=1,q.push(v);
					num[v]++;
					if( num[v]>=n )	return true;
				}
			}
		}
	}
	return false;
}
signed main()
{
	ios::sync_with_stdio(false);
	cin.tie(0); cout.tie(0);
	cin >> n >> m;
	for(int i=1;i<=n;i++)	cin >> p[i];
	for(int i=1;i<=m;i++)
	{
		int l,r; double w;
		cin >> l >> r >> w;
		add(l,r,w); 
	}
	for(int i=1;i<=n;i++)	add(0,i,0); 
	double l=0,r=20Q0000000009,ans=-1;
	while( r>=l+eps )
	{
		double mid = (l+r)/2.0;
		if( spfa(mid) )	ans=mid,l=mid+eps;
		else	r=mid-eps;
	}
	printf("%.2lf",ans);
}
©️2020 CSDN 皮肤主题: 1024 设计师:上身试试 返回首页