收费站
【题目描述】
在某个遥远的国家里,有n个城市。编号为1,2,3,……,n。
这个国家的政府修建了m条双向的公路。每条公路连接着两个城市。沿着某条公路,开车从一个城市到另一个城市,需要花费一定的汽油。
开车每经过一个城市,都会被收取一定的费用(包括起点和终点城市)。所有的收费站都在城市中,在城市间的公路上没有任何的收费站。
小红现在要开车从城市u到城市v(1<=u,v<=n)。她的车最多可以装下s升的汽油。在出发的时候,车的油箱是满的,并且她在路上不想加油。
在路上,每经过一个城市,她要交一定的费用。如果她某次交的费用比较多,她的心情就会变得很糟。所以她想知道,在她能到达目的地的前提下,她交的费用中最多的一次最少是多少。这个问题对于她来说太难了,于是她找到了聪明的你,你能帮帮她吗?
【输入】
第一行5个正整数,n,m,u,v,s。分别表示有n个城市,m条公路,从城市u到城市v,车的油箱的容量为s升。
接下来有n行,每行1个正整数,fi。表示经过城市i,需要交费fi元。
再接下来有m行,每行3个正整数,ai,bi,ci(1<=ai,bi<=n)。表示城市ai和城市bi之间有一条公路,如果从城市ai到城市bi,或者从城市bi到城市ai,需要用ci升汽油。
【输出】
仅一个整数,表示小红交费最多的一次的最小值。
如果她无法到达城市v,输出-1。
【样例输入输出1】
cost.in | cost.out |
4 4 2 3 8 8 5 6 10 2 1 2 2 4 1 1 3 4 3 4 3 | 8 |
【样例输入输出2】
cost.in | cost.out |
4 4 2 3 3 8 5 6 10 2 1 2 2 4 1 1 3 4 3 4 3 | -1 |
【数据规模】
对于60%的数据,满足n<=200,m<=10000,s<=200
对于100%的数据,满足n<=10000,m<=50000,s<=1000000000
对于100%的数据,满足ci<=1000000000,fi<=1000000000,可能有两条边连接着相同的城市。
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include <algorithm>
using namespace std;
FILE *in, *out;
typedef long long int ll;
ll q[3000001];
ll f[10001];
ll a[10001];
ll d[10001];
ll v[10001];
ll b[10001][251];
ll c[10001][251];
ll n, m, st, ed, tt, limit, tot;
void find_it(ll s, ll t)
{
ll i, j, x, tmp;
i = s;
j = t;
x = f[(i + j) / 2];
while (i <= j) {
while (f[i] < x) { ++i; }
while (f[j] > x) { --j; }
if (i <= j) {
tmp = f[i];
f[i] = f[j];
f[j] = tmp;
++i;
--j;
}
}
if (s < j) find_it(s, j);
if (i < t) find_it(i, t);
}
void spfa(ll s)
{
ll i, head, tail, now;
sizeof(v, 0, sizeof(v));
for (i = 1; i <= n; ++i) d[i] = 0x7fffffff;
d[s] = 0;
q[1] = s;
head = 1;
tail = 1;
v[s] = 1;
while (head <= tail) {
now = q[head];
for (i = 1; i <= b[now][0]; ++i) {
if (a[b[now][i]]<=limit) {
if (d[b[now][i]]>d[now]+c[now][i]) {
d[b[now][i]] = d[now]+c[now][i];
if (v[b[now][i]] == 0) {
v[b[now][i]] = 1;
++tail;
q[tail] = b[now][i];
}
}
}
}
++head;
v[now] = 0;
}
}
void process(ll s, ll t)
{
ll mid;
if (s == t) { fprintf(out, "%d\n", f[s]); fclose(in);fclose(out);return ;}
mid = (s + t) / 2;
limit = f[mid];
spfa(st);
if (d[ed] > tt) process(mid + 1, t);
else process(s, mid);
}
int main()
{
in = fopen("cost.in", "rt");
out = fopen("cost.out", "wt");
ll s, t, l, i, lim;
fscanf(in, "%lld%lld%lld%lld%lld", &n, &m, &st, &ed, &tt);
for (i = 1; i <= n; ++i) fscanf(in, "%lld", &a[i]);
if (a[st] > a[ed]) lim = a[st];
else lim = a[ed];
for (i = 1; i <= n; ++i) {
if (a[i] >= lim) {
++tot;
f[tot] = a[i];
}
}
find_it(1, tot);
for (i = 1; i <= m; ++i) {
fscanf(in, "%lld%lld%lld", &s, &t, &l);
++b[s][0];
b[s][b[s][0]] = t;
c[s][b[s][0]] = l;
++b[t][0];
b[t][b[t][0]] = s;
c[t][b[t][0]] = l;
}
limit = 1000000001;
spfa(st);
if (d[ed] > tt) { fprintf(out, "-1\n"); system("pause");fclose(in);fclose(out);return 0;}
process(1, tot);
system("pause");
return 0;
}
/*
4 4 2 3 8
8 5 6 10
2 1 2
2 4 1
1 3 4
3 4 3
4 4 2 3 3
8 5 6 10
2 1 2
2 4 1
1 3 4
3 4 3
*/