算法提高 道路和航路

时间:2022-03-11 11:16:46
 算法提高 道路和航路  
时间限制:1.0s   内存限制:256.0MB
提交此题   锦囊1   锦囊2
问题描述
农夫约翰正在针对一个新区域的牛奶配送合同进行研究。他打算分发牛奶到T个城镇(标号为1..T),这些城镇通过R条标号为(1..R)的道路和P条标号为(1..P)的航路相连。

每一条公路i或者航路i表示成连接城镇Ai(1<=A_i<=T)和Bi(1<=Bi<=T)代价为Ci。每一条公路,Ci的范围为0<=Ci<=10,000;由于奇怪的运营策略,每一条航路的Ci可能为负的,也就是-10,000<=Ci<=10,000。

每一条公路都是双向的,正向和反向的花费是一样的,都是非负的。

每一条航路都根据输入的Ai和Bi进行从Ai->Bi的单向通行。实际上,如果现在有一条航路是从Ai到Bi的话,那么意味着肯定没有通行方案从Bi回到Ai。

农夫约翰想把他那优良的牛奶从配送中心送到各个城镇,当然希望代价越小越好,你可以帮助他嘛?配送中心位于城镇S中(1<=S<=T)。

输入格式
输入的第一行包含四个用空格隔开的整数T,R,P,S。

接下来R行,描述公路信息,每行包含三个整数,分别表示Ai,Bi和Ci。

接下来P行,描述航路信息,每行包含三个整数,分别表示Ai,Bi和Ci。

输出格式
输出T行,分别表示从城镇S到每个城市的最小花费,如果到不了的话输出NO PATH。
样例输入
6 3 3 4
1 2 5
3 4 5
5 6 10
3 5 -100
4 6 -100
1 3 -10
样例输出
NO PATH
NO PATH
5
0
-95
-100
数据规模与约定
对于20%的数据,T<=100,R<=500,P<=500;

对于30%的数据,R<=1000,R<=10000,P<=3000;

对于100%的数据,1<=T<=250001<=R<=500001<=P<=50000



只得了85分,尽管刚学习了如何用输入输出外挂,但还是不行,有三组数据超时,我用的bellman方法,看来还是得学spfa算法。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <vector>
#include <set>
#include <queue>
#define max_v 650005
#define max_e 1000005
#define INF 0x3f3f3f3f

using namespace std;

struct edge{int from,to,cost;};
edge ee[max_e];//存储每条边的信息
long long int d[max_v];//表示s到每个顶点的最短距离
int v,e;//顶点数和边数
int r,p,s;

template <class T>
inline bool scan_d(T &ret) {
    char c; int sgn;
    if(c=getchar(),c==EOF) return 0; //EOF
    while(c!='-'&&(c<'0'||c>'9')) c=getchar();
    sgn=(c=='-')?-1:1;
    ret=(c=='-')?0:(c-'0');
    while(c=getchar(),c>='0'&&c<='9') ret=ret*10+(c-'0');
    ret*=sgn;
    return 1;
}
inline void out(long long int x) {
    if(x>9) out(x/10);
    putchar(x%10+'0');
}

bool bellman(int s){
    memset(d,INF,sizeof(d));//赋最大值的技巧
    d[s]=0;
    int j=0;//添加的计数器
    while(true){
        bool update=false;
        for(int i=0;i<e+r;i++){
            edge et=ee[i];
            if(d[et.from]!=INF&&d[et.to]>d[et.from]+et.cost){
                d[et.to]=d[et.from]+et.cost;
                update=true;
            }
        }
        j++;
        if(!update) return true;
        if(j==v) return false;//当进行第v次循环时,说明存在负圈
    }
}

int main()
{
    //freopen("data.txt", "r", stdin);
    //freopen("data.out", "w", stdout);
    //scanf("%d %d %d %d",&v,&r,&p,&s);
    scan_d(v);
    scan_d(r);
    scan_d(p);
    scan_d(s);
    e=r+p;
    int ii=0;
    int t1,t2,t3;
    for(int i=0;i<e;i++){
        //scanf("%d %d %d",&t1,&t2,&t3);
        scan_d(t1);
        scan_d(t2);
        scan_d(t3);
        t1-=1;
        t2-=1;
        if(i<r){
            ee[ii].from=t1;
            ee[ii].to=t2;
            ee[ii].cost=t3;
            ee[ii+e].from=t2;
            ee[ii+e].to=t1;
            ee[ii+e].cost=t3;
            ii++;
        }else{
            ee[ii].from=t1;
            ee[ii].to=t2;
            ee[ii].cost=t3;
            ii++;
        }
    }
    if(bellman(s-1)){
        for(int i=0;i<v;i++){
            if(d[i]<249990009){
                //printf("%d\n",d[i]);
                if(d[i]<0){
                    printf("-");
                    out(d[i]*-1);
                    printf("\n");
                }else{
                    out(d[i]);
                    printf("\n");
                }
            }else{
                printf("NO PATH\n");
            }

        }
    }
    //fclose(stdin);
    //fclose(stdout);
    return 0;
}
/*

*/