【LOJ】 #2308. 「APIO2017」商旅

题解

分数题可以想到分数规划，我们预处理出从i到j卖什么货物赚的最多，然后把每条边的边权改成“利润 - 效率 × 时间”
用spfa找正环即可

代码

#include <bits/stdc++.h>
#define fi first
#define se second
#define pii pair<int,int>
#define mp make_pair
#define pb push_back
#define enter putchar('\n')
#define space putchar(' ')
#define MAXN 300005
#define mo 994711
#define eps 1e-8
//#define ivorysi
using namespace std;
typedef long long int64;
typedef long double db;
typedef unsigned int u32;
template<class T>
void read(T &res) {
    res = 0;char c = getchar();T f = 1;
    while(c < '0' || c > '9') {
        if(c == '-') f = -1;
        c = getchar();
    }
    while(c >= '0' && c <= '9') {
        res = res * 10 + c - '0';
        c = getchar();
    }
    res *= f;
}
template<class T>
void out(T x) {
    if(x < 0) {putchar('-');x = -x;}
    if(x >= 10) out(x / 10);
    putchar('0' + x % 10);
}
int N,M,K;
int64 g[105][105],h[105][105],inf = 1e14;
int64 B[105][1005],S[105][1005],dis[105];
bool vis[105],mark[105];
struct node {
    int to,next;int64 val;
}E[20005];
int head[105],sumE;
void add(int u,int v,int64 c) {
    E[++sumE].to = v;
    E[sumE].next = head[u];
    E[sumE].val = c;
    head[u] = sumE;
}
bool SPFA_dfs(int u) {
    vis[u] = 1;
    mark[u] = 1;
    for(int i = head[u] ; i ; i = E[i].next) {
    int v = E[i].to;
    if(dis[v] <= dis[u] + E[i].val) {
        dis[v] = dis[u] + E[i].val;
        if(vis[v]) return false;
        if(!SPFA_dfs(v)) return false;
    }
    }
    vis[u] = 0;
    return true;
}
bool check(int64 mid) {
    memset(head,0,sizeof(head));sumE = 0;
    for(int i = 1 ; i <= N ; ++i) {
    for(int j = 1 ; j <= N ; ++j) {
        if(g[i][j] < inf)  add(i,j,h[i][j] - g[i][j] * mid);
    }
    }
    memset(vis,0,sizeof(vis));memset(mark,0,sizeof(mark));
    for(int i = 1 ; i <= N ; ++i) dis[i] = 0;
    for(int i = 1 ; i <= N ; ++i) {
    if(!mark[i]) {
        if(!SPFA_dfs(i)) return true;
    }
    }
    return false;
}
void Solve() {
    read(N);read(M);read(K);
    for(int i = 1 ; i <= N ; ++i) {
    for(int j = 1 ; j <= N ; ++j) {
        g[i][j] = inf;
    }
    }
    for(int i = 1 ; i <= N ; ++i) {
    for(int j = 1 ; j <= K ; ++j) {
        read(B[i][j]);read(S[i][j]);
    }
    }
    int v,w;int64 t;
    for(int i = 1 ; i <= M ; ++i) {
    read(v);read(w);read(t);
    g[v][w] = t;
    }
    for(int k = 1 ; k <= N ; ++k) {
    for(int i = 1 ; i <= N ; ++i) {
        for(int j = 1 ; j <= N ; ++j) {
        g[i][j] = min(g[i][j],g[i][k] + g[k][j]);
        }
    }
    }
    for(int i = 1 ; i <= N ; ++i) {
    for(int j = 1 ; j <= N ; ++j) {
        for(int k = 1 ; k <= K ; ++k) {
        if(B[i][k] != -1 && S[j][k] != -1) h[i][j] = max(S[j][k] - B[i][k],h[i][j]);
        }
    }
    }
    int64 L = 0,R = 1e9;
    while(L < R) {
    int64 mid = (L + R + 1) >> 1;
    if(check(mid)) L = mid;
    else R = mid - 1;
    }
    out(L);enter;
}

int main() {
#ifdef ivorysi
    freopen("f1.in","r",stdin);
#endif
    Solve();
}

转载于:https://www.cnblogs.com/ivorysi/p/9766128.html