Heavy Transportation POJ - 1797 && Frogger POJ - 2253 Dijstra 算法变形

这两道题都对 Dijstra 进行了改造，继承了 Dijstra 的精髓，也就是贪心和松弛。我们先找贪心的从未访问过的点中找到最符合条件的点，因为通过这个点确定的边进行松弛之后有可能改善之前的情况，比如最短路中先找到当前最短的边，然后利用这条边进行松弛，就可能通过这条边找到到其他点的最短路径，然后重复此过程即可。

POJ 1797 题目链接

题意：

在一个图中，每条边都有一定的承载力，一辆车从源点往终点运送货物，货物重量不能超过他经过的边的承载力，问从源点到目标点的路线中，它能运送货物的最大重量是多少。

思路：

我们用数组 d[N] 表示起点到各个点的最大承载力的值，0 代表起点，起初我们置 d[0] 为无穷大。就可以认为 0 号点是当前已知最大承载力的点，也就是最符合条件的，然后以此点开始进入循环，并进行松弛操作，通过 v 松弛后的最大承载力是 min(d[v],limit[v][i])，而我们希望通过松弛得到更大的承载力，因此我们应做如下更新 d[i] = max(d[i],min(d[v],limit[v][i]))。松弛操作完成之后重复选出未访问过的点钟最符合条件的点进行松弛操作即可。

代码如下：

#include <cstdio>
#include <iostream>
#include <cstring>
#define MAX_V 1005
#define fp(p1,p2,p3) for(int p1 = p2;p1 < p3;p1 ++)
using namespace std;
int r,round,n,x,y,w,m,used[MAX_V],d[MAX_V],limit[MAX_V][MAX_V];
int dijstra(){
    fill(used,used + MAX_V,0);
    fill(d,d + MAX_V,0);
    d[0] = (1 << 28);
    while(1){
        int v = -1;
        fp(i,0,n){
            if(!used[i]&&(v==-1||d[v]<d[i])){
                v=i;
            }
        }
        if(v==-1) break;
        used[v]=1;
        fp(i,0,n){
            if(!used[i]){
                d[i] = max(d[i],min(d[v],limit[v][i]));
            }
        }
    }
    return d[n - 1];
}
int main(void)
{
    cin >> round;
    while(r ++ < round){
        memset(limit,0,sizeof(limit));
        scanf("%d%d",&n,&m);
        fp(i,0,m){
            scanf("%d%d%d",&x,&y,&w);
            x --,y --;
            limit[x][y] = limit[y][x] = w;
        }
        printf("Scenario #%d:\n%d\n\n",r,dijstra());
    }
    return 0;
}

POJ 2253 题目链接

题意：

一直青蛙想通过在池塘中石头上跳来跳去的方式去拜访一个朋友，问为了能拜访他朋友，他的最小所需跳跃能力为多少。

思路：

首先我们应将题目中给出的坐标转化为邻接矩阵。然后我们用数组 d[N] 来表示从起点 0 到达某点所需要的最小跳跃能力，显然 d[0] = 0,然后从未访问过的各点找到所需最小所需跳跃能力最小的点，置访问位为真（因为这个点不可能通过其他点松弛进而得到一个更小的最小所需跳跃能力），然后通过此点进行松弛，更新其他点的最小所需跳跃能力，通过松弛后的最小所需跳跃能力为 max(d[v],costs[v][i])（因为他必须有能力跳过路径中最长的边才行），我们期望得到的最小所需跳跃能力尽量小，所以做如下更新 d[i]=min(d[i],max(d[v],costs[v][i])) ,然后重复这个操作即可。
代码：

#include <iostream>
#include <cstdio>
#include <cmath>
#define N 205
#define INF (1 <<28)
#define fp(p,q,r) for(int p = q;p<r;p++)
using namespace std;
typedef long long ll;
double co[N][2];
ll ro;
bool used[N];
double costs[N][N],d[N],n;
double dijstra(){
    fill(used,used+N,0);
    fill(d,d+N,INF);
    d[0]=0;
    while(1){
        int v=-1;
        fp(i,0,n){
            if(!used[i]&&(v==-1||d[v]>d[i])) v=i;
        }
        if(v==-1) break;
        used[v]=1;
        fp(i,0,n){
            d[i]=min(d[i],max(d[v],costs[v][i]));
        }
    }
    return d[1];
}
int main(void)
{
    while(cin >> n && n){
        fill(costs[0],costs[0] + N * N,INF);
        fp(i,0,n){
            ll x,y;
            cin >> x >> y;
            co[i][0] = x,co[i][1]=y;
            fp(j,0,i){
                costs[i][j]=costs[j][i]=sqrt((co[i][0]-co[j][0])*(co[i][0]-co[j][0])+(co[i][1]-co[j][1])*(co[i][1]-co[j][1]));
            }
        }
        printf("Scenario #%lld\nFrog Distance = %.3lf\n\n",++ ro,dijstra());
    }
    return 0;
}