【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();
}
上一篇:Charles代理配置截取https iOS11


下一篇:Java语言定义的线程状态分析