bzoj 4898: [Apio2017]商旅

题意：

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

题解：

分数规划+floyd+判负环。
分数规划就挺套路了，那么先做一次floyd，再算出两个点间最大获利，spfa判负环就可以了。
还有k应该是1000。
code:

#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<cstring>
#define LL long long
using namespace std;
LL n,m,k;
LL b[110][1010],s[110][1010],val[110][110];
LL f[110][110],dis[110];
bool vis[110],ty;
void dfs(LL x,LL ans)
{
    vis[x]=true;
    for(LL i=1;i<=n;i++)
    {
        if(ty) return;
        if(i==x||f[x][i]>999999999) continue;
        LL c=ans*f[x][i]-val[x][i];
        if(dis[i]>=dis[x]+c)
        {
            if(vis[i]) {ty=true;return;}
            dis[i]=dis[x]+c;
            dfs(i,ans);
        }
    }
    vis[x]=false;
}
bool check(LL x)
{
    memset(dis,0,sizeof(dis));
    memset(vis,false,sizeof(vis));
    ty=false;
    for(LL i=1;i<=n;i++)
    {
        dfs(i,x);
        if(ty) return true;
    }
    return false;
}
int main()
{
    scanf("%lld %lld %lld",&n,&m,&k);
    for(LL i=1;i<=n;i++)
        for(LL j=1;j<=k;j++)
            scanf("%lld %lld",&b[i][j],&s[i][j]);
    for(LL i=1;i<=n;i++)
        for(LL j=1;j<=n;j++)
        {
            val[i][j]=0;
            for(LL u=1;u<=k;u++)
                if(b[i][u]!=-1&&s[j][u]!=-1)
                    val[i][j]=max(val[i][j],s[j][u]-b[i][u]);
        }
    LL l=1,r=1000000000,ans=0;
    memset(f,63,sizeof(f));
    for(LL i=1;i<=m;i++)
    {
        LL x,y,c;scanf("%lld %lld %lld",&x,&y,&c);
        f[x][y]=c;
    }
    for(LL K=1;K<=n;K++)
        for(LL i=1;i<=n;i++)
            for(LL j=1;j<=n;j++)
                f[i][j]=min(f[i][j],f[i][K]+f[K][j]);
    while(l<=r)
    {
        LL mid=(l+r)/2;
        if(check(mid)) ans=mid,l=mid+1;
        else r=mid-1;
    }
    printf("%lld",ans);
}