最短路Dijksta算法温故

POJ2387

最基本得最短路问题，双向图，无负权回路，可以用来熟悉一下自己得思路，和预处理得过程

#include <iostream>
#include <string.h>
#define inf 0x3f3f3f3f
using namespace std;
const int maxn = 1010;
int mp[maxn][maxn];
int dis[maxn];
int vis[maxn];
int N,T;

void init(int n)
{
    //memset(vis,0,sizeof(vis));
    //memset(dis,0,sizeof(dis));
    for(int i = 1;i <= n;i++)
    {
        for(int j = 1;j <= n;j++)
        {
            if(i == j)mp[i][j] = 0;
            else mp[i][j] = inf;
        }
        vis[i] = 0;
    }
}
void Dijkstra(int s,int e)
{
    vis[s] = 1;
    for(int i = 1;i <= N;i++)dis[i] = mp[s][i];//根据中央节点初始化dis
    int net;
    for(int i = 1;i <= N;i++)//i其实就是一个松弛得次数,每次松弛肯定会优化一个d松弛N次就足够了
    {
        int minlen = inf;
        for(int j = 1;j <= N;j++)
        {
            if(!vis[j] && dis[j] < minlen)
            {
                minlen = dis[j];
                net = j;
            }
        }
        if(minlen == inf)break;
        vis[net] = 1;
        for(int j = 1;j <= N;j++)
        {
            if(!vis[j] && dis[j] > mp[net][j] + dis[net])
            {
                dis[j] = mp[net][j] + dis[net];
            }
        }
    }
}
int main()
{
    int a,b,x;
    while(~scanf("%d%d",&T,&N))
    {
        init(N);
        while(T--)
        {
            scanf("%d%d%d",&a,&b,&x);
            if(mp[a][b] > x)mp[a][b] = mp[b][a] = x;
        }
        Dijkstra(1,N);
        printf("%d\n",dis[N]);
    }
    return 0;
}

本来习惯用memset，其实在预处理mp得时候，就可以顺带着处理vis清零得问题

在dijksta算法中，以中央节点1为起点预处理dis数组，含义是j点到中央节点的最短距离

第二次i得for循环是循环的松弛次数，所以可能已经松弛完毕，你需要判断一下可以用minlen == inf 来判断，也可以预处理net = 0，在进行一次筛选点的时候，出来发现net任然是0也可以判断是松弛完毕

POJ2253

这个题题意我就理解了半天是真的半天，青蛙要从石头1跳到石头2，需要借助其他的石头，那么我们肯定有不止一条路，面对这个多的路，每一条路都有一个最大的跳跃距离比如1跳到3距离为99，3跳到2距离为66，那么这条路径得最大跳跃距离是99，问你所有路径中最大跳跃路径得最小值~~

这里求得并不是dis[2]最短路了，而是dis[2]所有路径中最大路径中的最小值，所以我dis得意义也被改变了，dis【i】表示得是从1到i节点的所有路径中最短的跳跃路径

也就是普通的Dijkstra算法求得net点时，一层for循环寻找可以优化得j点

‘那么dis[j] = min(dis[j],max(dis[net],mp[net][j]))

右边的max里dis【net】表示得所有得到net点的路径方案中最小的跳跃距离，那么这些路径再通一个j点加上了mp【net】【j】时，dis【j】得结果可能时dis【j】也有可能会因此变小

可能还时有点不理解为什么内部时max而不是min，那就来分类讨论一下，如果在原来的所有的路径上，我加上了mp[net][j]那么如果mp[net][j]就会时原来所有路径得一部分，一个结尾，如果这个结尾比dis[net]大那就代表由net到j最大跳跃距离要跟新了，因为原来的那个值（dis[net]）并不是那一条路径上的最大跳跃距离了，所以更新了，如果比dis[net]小，那就不必更新直接去比较就好~~

#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#define inf 0x3f3f3f3f
using namespace std;
const int maxn = 210;
double mp[maxn][maxn];
double stone[maxn][2];
double dis[maxn];
int vis[maxn];
int n;
double far(double ax,double ay,double bx,double by)
{
    return sqrt((ax - bx) * (ax - bx) + (ay - by) * (ay - by));
}
void Dijkstra(int s)
{
    vis[s] = 1;
    for(int i = 1;i <= n;i++)dis[i] = mp[s][i];

    for(int i = 1;i <= n;i++)
    {
        int minlen = inf;
        int net = 0;
        for(int j = 1;j <= n;j++)
        {
            if(!vis[j] && dis[j] < minlen)
            {
                minlen = dis[j];
                net = j;
            }
        }
        if(net == 0)break;
        vis[net] = 1;
        for(int j = 1;j <= n;j++)
        {
            if(!vis[j])
            {
                dis[j] = min(dis[j],max(dis[net],mp[net][j]));
            }
        }
    }
    /*现在要求求出所有通路的最大距离，
    并把这些最大距离作比较，
    把最小的一个最大距离作为青蛙的最小跳远距离。*/
}



int main()
{
    int cas = 1;
    while(~scanf("%d",&n),n)
    {
        for(int i = 1;i <= n;i++)
        {
            scanf("%lf %lf",&stone[i][0],&stone[i][1]);
        }
        for(int i = 1;i <= n;i++)
        {
            for(int j = 1;j <= n;j++)
            {
                if(i == j)mp[i][j] = 0.0;
                else mp[i][j] = mp[j][i] = far(stone[i][0],stone[i][1],stone[j][0],stone[j][1]);
            }
            vis[i] = 0;
            dis[i] = inf;
        }
        Dijkstra(1);
        printf("Scenario #%d\nFrog Distance = %.3lf\n\n",cas++,dis[2]);
    }
    return 0;
}