【BZOJ4898】[Apio2017]商旅 分数规划+SPFA

【BZOJ4898】[Apio2017]商旅

Description

在广阔的澳大利亚内陆地区长途跋涉后，你孤身一人带着一个背包来到了科巴。你被这个城市发达而美丽的市场所深深吸引，决定定居于此，做一个商人。科巴有个集市，集市用从1到N的整数编号，集市之间通过M条单向道路连接，通过每条道路都需要消耗一定的时间。在科巴的集市上，有K种不同的商品，商品用从1到K的整数编号。每个集市对每种商品都有自己的定价，买入和卖出商品的价格可以是不同的。并非每个集市都可以买卖所有的商品：一个集市可能只提供部分商品的双向交易服务；对于一种商品，一个集市也可能只收购而不卖出该商品或只卖出而不收购该商品。如果一个集市收购一种商品，它收购这种商品的数量是不限的，同样，一个集市如果卖出一种商品，则它卖出这种商品的数量也是不限的。为了更快地获得收益，你决定寻找一条盈利效率最高的环路。环路是指带着空的背包从一个集市出发，沿着道路前进，经过若干个市场并最终回到出发点。在环路中，允许多次经过同一个集市或同一条道路。在经过集市时，你可以购买或者卖出商品，一旦你购买了一个商品，你需要把它装在背包里带走。由于你的背包非常小，任何时候你最多只能持有一个商品。在购买一个商品时，你不需要考虑你是否有足够的金钱，但在卖出时，需要注意只能卖出你拥有的商品。从环路中得到的收益为在环路中卖出商品得到的金钱减去购买商品花费的金钱，而一条环路上消耗的时间则是依次通过环路上所有道路所需要花费的时间的总和。环路的盈利效率是指从环路中得到的收益除以花费的时间。需要注意的是，一条没有任何交易的环路的盈利效率为0。你需要求出所有消耗时间为正数的环路中，盈利效率最高的环路的盈利效率。答案向下取整保留到整数。如果没有任何一条环路可以盈利，则输出0。

Input

第一行包含3个正整数N,M和K,分别表示集市数量、道路数量和商品种类数量。

接下来的N行，第行中包含2K个整数描述一个集市Bi,1 Si,1 Bi,2 Si,2...Bik Si,k。

对于任意的1<=j<=k,整数和分别表示在编号为的集市上购买、卖出编号为的商品时的交易价格。

如果一个交易价格为-1，则表示这个商品在这个集市上不能进行这种交易。

接下来M行，第行包含3个整数Vp,Wp和Tp,表示存在一条从编号为Vp的市场出发前往编号为Wp的市场的路径花费Tp分钟。

1<=N<=100,1<=M<=9900

如果在编号为的集市i中，编号为j的商品既可以购买又可以卖出则0<Si,j<=Bi,j<=10^9

对于编号为P(1<=P<=M)的道路，保证Vp<>Wp且1<=Tp<=10^7

不存在满足1<=P<Q<=M的P,Q,使得(Vp,Wp)=(Vq,Wq) 。

Output

输出包含一个整数，表示盈利效率最高的环路盈利效率，答案向下取整保留到整数。如果没有任何一条环路可以盈利，则输出0。

 

Sample Input

4 5 2
10 9 5 2
6 4 20 15
9 7 10 9
-1 -1 16 11
1 2 3
2 3 3
1 4 1
4 3 1
3 1 1

Sample Output

2
在样例中，我们考虑下面两条环路，“1 - 2 - 3 - 1” 和 “1 - 4 - 3 - 1”。
考虑环路 “1 - 2 - 3 - 1” ：这条环路消耗的总时间是 分钟。在这条环路中，最佳的交易方式是：在编号为 1 的集市中购买编号为 2 的商品（花费的金钱为 5 ）；在编号为 2 的集市中卖出编号为 2 的商品（得到的金钱为 15 ），然后立即购买编号为 1 的商品（花费的金钱为 6 ）；带着编号为 1 的商品经过编号为 3 的集市，在回到编号为 1 的城市后卖出（得到的金钱为 9 ）。在这个环路中，总盈利为13。 这个环路的盈利效率为13/7 ，向下取整后为 1 。
考虑环路 “1 - 4 - 3 - 1” ：这条环路消耗的总时间是 分钟。在这条环路中，最佳的交易方式是：在编号为 1 的集市中购买编号为 2 的商品（花费的金钱为 5 ）；在编号为 4 的集市中卖出编号为 2 的商品（得到的金钱为 11 ）；然后经过编号为 3 的集市回到编号为 1 的城市。在这个环路中，总盈利为 6。 这个环路的盈利效率为6/3 ，向下取整后为 2 。
综上所述，盈利效率最高的环路的盈利效率为 2 。

题解：考场上最简单的题，嗯，当我把正解想出来时，距离考试结束还有不到20分钟~

先预处理出任意两个点之间的最短路径以及最优购买策略。那么二分答案mid，从i到j连（收益-长度*mid）的边，（如果i能到j的话），然后用SPFA判正环即可。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <queue>
#include <cmath>
using namespace std;
typedef long long ll;
int n,m,K,cnt;
int S[110][1010],B[110][1010],to[20010],next[20010],head[110],inq[110],len[110];
double val[20010],dis[110];
int map[110][110],td[110][110];
queue<int> q;
inline void add(int a,int b,double c)
{
	to[cnt]=b,val[cnt]=c,next[cnt]=head[a],head[a]=cnt++;
}
bool check(double x)
{
	int i,j,u;
	memset(head,-1,sizeof(head)),memset(dis,0,sizeof(dis)),cnt=0;
	for(i=1;i<=n;i++)	for(j=1;j<=n;j++)	if(i!=j&&td[i][j]>=0)	add(i,j,td[i][j]-x*map[i][j]);
	for(i=1;i<=n;i++)	q.push(i),inq[i]=len[i]=1;
	while(!q.empty())
	{
		u=q.front(),q.pop(),inq[u]=0;
		for(i=head[u];i!=-1;i=next[i])	if(dis[to[i]]<dis[u]+val[i])
		{
			dis[to[i]]=dis[u]+val[i],len[to[i]]=len[u]+1;
			if(len[to[i]]>n)	return 1;
			if(!inq[to[i]])	inq[to[i]]=1,q.push(to[i]);
		}
	}
	return 0;
}
inline int rd()
{
	int ret=0,f=1;	char gc=getchar();
	while(gc<'0'||gc>'9')	{if(gc=='-')f=-f;	gc=getchar();}
	while(gc>='0'&&gc<='9')	ret=ret*10+gc-'0',gc=getchar();
	return ret*f;
}
int main()
{
	n=rd(),m=rd(),K=rd();
	int i,j,k,a,b,c;
	memset(map,0x3f,sizeof(map)),memset(td,0xc0,sizeof(td));
	for(i=1;i<=n;i++)
	{
		map[i][i]=0;
		for(j=1;j<=K;j++)	B[i][j]=rd(),S[i][j]=rd();
	}
	for(i=1;i<=m;i++)	a=rd(),b=rd(),c=rd(),map[a][b]=min(map[a][b],c);
	for(k=1;k<=n;k++)	for(i=1;i<=n;i++)	for(j=1;j<=n;j++)	map[i][j]=min(map[i][j],map[i][k]+map[k][j]);
	double l=0,r=0,mid;
	for(i=1;i<=n;i++)	for(j=1;j<=n;j++)	if(i!=j&&map[i][j]<0x3f3f3f3f)
	{
		td[i][j]=0;
		for(k=1;k<=K;k++)	if(B[i][k]!=-1&&S[j][k]!=-1&&B[i][k]<S[j][k])	td[i][j]=max(td[i][j],S[j][k]-B[i][k]);
		r=(r>td[i][j])?r:td[i][j];
	}
	for(i=1;i<=30;i++)
	{
		mid=(l+r)/2;
		if(check(mid))	l=mid;
		else	r=mid;
	}
	printf("%d",int(floor(l+1e-9)));
	return 0;
}