【*】第三题：收费站（提高组第一试2011年10月21日）（2011年NOIP冲刺模拟试题）（

 收费站

【题目描述】

    在某个遥远的国家里，有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
*/