POJ - 3662

Telephone Lines

【Point】

1 ≤ N ≤ 1,000  1 ≤ P ≤ 10,000

很明显这是一个稀疏图求最短路的问题，紫书上有Dijkstra的plogn写法，但是我看了一下时间，发现大多100ms就可以完成，那么n^2也是可以接受的，就试着暴力了一发，结果782ms完成。

可能是比较老的题目，当时的评测机速度不行，划完这波水，我再去认真写写

【Algorithm】

要求行走路上第k+1大的值，用DFS或者DP强行算都是不行的，事实上，傻傻的我都写了一遍，成功超时了。

那么，就要用到神·作弊方法——二分，先把答案猜出来，再用答案去验证是否可行，对应的就是得出答案后，计算行走路上有多少个权值是大于我们所猜测的答案的。只要在最后的数目是满足要求的，那这个答案就是可行的。

【Code】

#define _CRT_SECURE_NO_WARNINGS

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>

using namespace std;
const int maxn = 1000 + 10;
const int inf = 1e6 + 10;
int d[maxn], v[maxn];
int w[maxn][maxn];
vector<int> G[maxn];
int n, p, k;

void dijkstra(int lim)
{
    memset(v, 0, sizeof(v));
    d[1] = 0;
    for (int i = 2; i <= n; i++)
        d[i] = inf;
    for (int i = 1; i <= n; i++)
    {
        int x, m = inf;
        for (int y = 1; y <= n; y++) if (!v[y] && d[y] <= m) m = d[x = y];
        v[x] = 1;

        for (int y = 1; y <= n; y++)
        {
            if (w[x][y] == inf)
                continue;
            int flag = 0;
            if (w[x][y] > lim)
                flag = 1;
            d[y] = min(d[y], d[x] + flag);
        }
    }
}

int main()
{
    scanf("%d%d%d", &n, &p, &k);
    for (int i = 1; i <= n; i++)
        for (int j = 1; j <= n; j++)
            w[i][j] = i == j ? 0 : inf;
    int a, b, ll;
    for (int i = 0; i < p; i++)
    {
        scanf("%d%d%d", &a, &b, &ll);
        G[a].push_back(b);
        G[b].push_back(a);
        w[a][b] = w[b][a] = ll;
    }
    int l = 0, r = inf;
    int ans = -1;
    while (l <= r)
    {
        int mid = l + (r - l) / 2;

        dijkstra(mid);
        if (d[n] <= k)
        {
            ans = mid;
            r = mid - 1;
        }
        else
            l = mid + 1;
    }
    printf("%d\n", ans);

    return 0;
}

【Note】

懒人也有懒解法，方法虽然笨，但是代码量确实少

Thanks For Watching

–END–