HPU1413- StarFarming [最短路]

HPU1413: StarFarming [最短路]

时间限制: 1 Sec 内存限制: 128 MB

提交: 359 解决: 46 统计

题目描述

星农（StarFarming）公司计划要给员工发路费津贴，发放的规则是这样的：1到n-1代表各个员工家的序号，n代表公司。路费津贴只发给上班的最短路与回家的最短路的总路程最长的人。该市的路建造的有些奇怪，修路只修单行道，即只允许往某一个方向通行。

现在给你城市的有向图的地图，TLG请你帮忙计算谁能得到津贴，以及他上班和回家的总路程是多少。

输入
有多组测试数据。

每组第一行输入两个整数N，M。表示点的个数，与单行道的数量（可能有重复）

接下来m行，每行输入三个整数x,y,z。表示从x到y城市有一条单行道，距离为z。

题目保证至少一人存在来回的路径。不存在的不发津贴（班都没法好好上还想要钱？！）

1≤N≤10001≤N≤1000
1≤M≤1000001≤M≤100000
1≤x,y≤N1≤x,y≤N
1≤z≤2001≤z≤200

输出

对于每组数据，输出两个整数，分别表示获得津贴的人的序号以及总路程。（如果有多个人路程相同，取序号最小的）

样例输入

4 7

1 2 2

2 3 2

1 3 4

4 1 2

4 2 2

3 4 1

4 3 5

样例输出

1 7

提示

对于样例，

1来回需要的最短路程是7：1->2->3->4->1

2来回需要的最短路程是5：2->3->4->2

3来回需要的最短路程是5：3->4->2->3

所以输出1 7

传送门

分析：对于这道题目，比赛的时候一直卡,先是因为没有while内没有取反超时了，一直没有发现，后来各种优化，什么链式前向星，Dijkstra堆优化，全用上了，还是T了，最后终于发现了,改了过来,之后开开心心去交了。mmp的，wrong answer, 比赛之后才发现 原来最后筛选最大值的时候忘记在小于inf的范围内取了，因为有的N到其他点到不了,dis的值为inf,在筛选最大值的时候加个判定条件就好了。唉，现在想想自己还是太菜，路还长着呢。

说下思路,题目需要求得是所有1到n-1的所有点到n点距离加上n点到其他n-1点的距离和最大的值和序号。跑两遍Dijkstra,注意正向建边和反向建边。（石乐志的少年！）

ac代码如下(有点乱)：

#include<iostream>
#include<algorithm>
#include<queue>
#include<stdio.h>
#include<cstring>
#define inf 1<<30
using namespace std;
typedef long long ll;
struct Edge{
    int nex,to,w;
}edge[1000005][2];
int head[1000005][2];
int cnt[3];
ll dis[1000005];
ll dis1[1000005];
ll dis2[1000005];
int n;
void add(int u,int v,int w,int p){
    edge[cnt[p]][p].to=v;
    edge[cnt[p]][p].w=w;
    edge[cnt[p]][p].nex=head[u][p];
    head[u][p]=cnt[p]++;
}
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
struct mycmp{
public:
    bool operator() (const int& lhs, const int&rhs) const  {
        return ( dis[lhs] > dis[rhs]);
    }
};
void dij(int s,int p){
    priority_queue<int ,vector<int>,mycmp>q;
    bool vis[1000005];
    memset(vis,0,sizeof vis);
    for(int i=1;i<=n;i++)dis[i]=inf;
    dis[s]=0;
    q.push(s);
    while(!q.empty()){
        int f=q.top();
        q.pop();
        if(vis[f])continue;
        vis[f]=1;
        for(int i=head[f][p];~i;i=edge[i][p].nex){
            int v=edge[i][p].to;
            int w=edge[i][p].w;
            if(dis[v]>dis[f]+w&&!vis[v]){
                dis[v]=dis[f]+w;
                q.push(v);
            }
        }
    }
}
int main()
{
    int m;
    while(~scanf("%d%d",&n,&m)){
        cnt[1]=0;
        cnt[2]=0;
        memset(head,-1,sizeof head);
        int u,v,w;
        for(int i=1;i<=m;i++){
            u=read(),v=read(),w=read();
            add(u,v,w,1);
            add(v,u,w,2);
        }
        memset(dis2,0,sizeof dis2);
        dij(n,1);
        for(int i=1;i<=n;i++)dis1[i] = dis[i];
        dij(n,2);
        //for(int i=1;i<=n;i++) dis2[i]=dis[i];
        /*加和*/
        for(int i = 1 ; i <= n ;i++){
            dis[i] += dis1[i];
        }
        ll index=0,num = -1;
        for(int i = 1; i<= n;i++){
            if(dis[i]>num&&dis[i]<inf){//一定要注意这里！！！
                index = i;
                num = dis[i];
            }
        }
        printf("%lld %lld\n",index,num);
    }
    return 0;
}