[APIO2017]商旅

link

这题卡我精度，调了一晚上才调对，因为没有想到图还可以不连通

其实可以预处理出好多东西，距离（$dis(u,v)$），买卖物品(从$u$到$v$买卖物品的最大利润，例($max{S_{u,i}-B_{v,i}}$)，然后其实可以发现就是一个十分普通普遍的分数规划式子，就每次二分$k$,然后建边，$spfa$判正环就行

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<climits>
#include<queue>
#define int long long
using namespace std;
inline int read(){
    int f=,ans=;char c=getchar();
    while(c<''||c>''){if(c=='-')f=-;c=getchar();}
    while(c>=''&&c<=''){ans=ans*+c-'';c=getchar();}
    return f*ans;
}
const int MAXN=;
int n,m,k;
int dis[MAXN][MAXN],b[MAXN][MAXN],s[MAXN][MAXN],st[MAXN][MAXN],l,r,maxn,D[MAXN][MAXN],num[MAXN],Dis[MAXN],vis[MAXN];
queue<int> que;
bool check(int P){
    while(!que.empty()) que.pop();
    for(int i=;i<=n;i++)
        for(int j=;j<=n;j++) D[i][j]=st[i][j]-dis[i][j]*P;
    memset(Dis,,sizeof(Dis));
    memset(num,,sizeof(num));
    memset(vis,,sizeof(vis));
    for(int i=;i<=n;i++) num[i]=vis[i]=,que.push(i);
    while(!que.empty()){
        int xx=que.front();que.pop();vis[xx]=;
        for(int i=;i<=n;i++){
            int d=D[xx][i];
            if(dis[xx][i]>INT_MAX) continue;
            if(Dis[xx]+d>=Dis[i]){
                Dis[i]=d+Dis[xx];
                num[i]=num[xx]+;
                if(num[i]>n+) return ;
                if(!vis[i]){
                    vis[i]=;
                    que.push(i);
                }

            }
        }
    }return ;
}
signed main(){
    memset(dis,/,sizeof(dis));
    n=read(),m=read(),k=read();
    for(int i=;i<=n;i++)
        for(int j=;j<=k;j++) b[i][j]=read(),s[i][j]=read();
    for(int i=;i<=k;i++){
        for(int x=;x<=n;x++)
           for(int y=;y<=n;y++)
            if(b[x][i]!=-&&s[y][i]!=-) {
                st[x][y]=max(st[x][y],s[y][i]-b[x][i]),r=max(r,st[x][y]);
            }
    }
    for(int i=;i<=m;i++){
        int u=read(),v=read(),w=read();
        dis[u][v]=min(dis[u][v],w);
    }
    for(int ss=;ss<=n;ss++)
        for(int i=;i<=n;i++)
            for(int j=;j<=n;j++)
                dis[i][j]=min(dis[i][j],dis[i][ss]+dis[ss][j]);
    while(l<=r){
        int mid=l+r>>;
        if(check(mid)) maxn=max(maxn,mid),l=mid+;
        else r=mid-;
    }
    printf("%lld\n",maxn);
}