zoj 2526(一道很好的最短路应用题）

题目链接：http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=1538

题意：最短路问题，但是要求出最短路的条数，同时要求出所有可能的最短路选择中javabean最多的情况。

思路：求到终点的最短路径用Dijkstra，其路径更新条件（如果到某个点有多个路径长度一样的最短路径，则选择豆子总数最多的）就是直接加个else if条件就搞定了，最后就是dfs搜一下最短路的条数了，这个我就记忆化了一下。

然后我一开始使用spfa写的，wa了好多次，orz...改成dijkstra就过了。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<vector>
using namespace std;
#define MAXN 555
#define inf 1<<30
/*
struct Node{
    int v,w;
};
vector<Node>vet[MAXN];
*/
int map[MAXN][MAXN];
int path[MAXN];
int dist[MAXN];
int value[MAXN];
int sum_value[MAXN];
int dp[MAXN];
int mark[MAXN];
int n,m,st,ed;

/*
void SPFA(){
    memset(sum_value,0,sizeof(sum_value));
    memset(path,-1,sizeof(path));
    for(int i=0;i<n;i++)dist[i]=inf;
    dist[st]=0,sum_value[st]=value[st];
    queue<int>Q;
    Q.push(st);
    while(!Q.empty()){
        int u=Q.front();
        Q.pop();
        for(int i=0;i<vet[u].size();i++){
            int v=vet[u][i].v;
            int w=vet[u][i].w;
            if(dist[u]+w<=dist[v]){
                path[v]=u;
                Q.push(v);
                if(dist[u]+w<dist[v]){
                    dist[v]=dist[u]+w;
                    sum_value[v]=sum_value[u]+value[v];
                }else if(sum_value[v]<sum_value[u]+value[v]){
                    sum_value[v]=sum_value[u]+value[v];
                }
            }
        }
    }
}
*/

void Dijkstra(){
    memset(mark,false,sizeof(mark));
    memset(sum_value,0,sizeof(sum_value));
    memset(path,-1,sizeof(path));
    for(int i=0;i<n;i++)dist[i]=map[st][i];
    dist[st]=0,sum_value[st]=value[st],mark[st]=true;
    int u=st;
    for(int i=0;i<n-1;i++){
        int min=inf;
        for(int j=0;j<n;j++){
            if(!mark[j]&&map[u][j]<inf&&dist[u]+map[u][j]<=dist[j]){
                if(dist[u]+map[u][j]<dist[j]){
                    dist[j]=dist[u]+map[u][j];
                    sum_value[j]=sum_value[u]+value[j];
                    path[j]=u;
                }else if(sum_value[j]<sum_value[u]+value[j]){
                    sum_value[j]=sum_value[u]+value[j];
                    path[j]=u;
                }
            }
        }
        for(int j=0;j<n;j++){
            if(!mark[j]&&min>dist[j]){
                min=dist[j],u=j;
            }
        }
        mark[u]=true;
    }
}
    


/*
int dfs(int u){
    if(u==ed)return 1;
    if(dp[u])return dp[u];
    for(int i=0;i<vet[u].size();i++){
        int v=vet[u][i].v;
        int w=vet[u][i].w;
        if(dist[v]==dist[u]+w){
            dp[u]+=dfs(v);
        }
    }
    return dp[u];
}
*/

int dfs(int u){
    if(u==ed)return 1;
    if(dp[u])return dp[u];
    for(int i=0;i<n;i++)if(u!=i){
        if(dist[i]==dist[u]+map[u][i]){
            dp[u]+=dfs(i);
        }
    }
    return dp[u];
}


void Print(int u){
    if(path[u]==-1){
        printf("%d",u);
        return ;
    }
    Print(path[u]);
    printf(" %d",u);
}

int main(){
    int u,v,w;
    while(~scanf("%d%d%d%d",&n,&m,&st,&ed)){
    //    for(int i=0;i<n;i++)vet[i].clear();
        for(int i=0;i<n;i++){
            map[i][i]=0;
            for(int j=i+1;j<n;j++){
                map[i][j]=map[j][i]=inf;
            }
        }
        for(int i=0;i<n;i++)
            scanf("%d",&value[i]);
        for(int i=1;i<=m;i++){
            scanf("%d%d%d",&u,&v,&w);
            /*
            Node p1,p2;
            p1.v=v,p2.v=u;
            p1.w=p2.w=w;
            vet[u].push_back(p1);
            vet[v].push_back(p2);
            */
            map[u][v]=map[v][u]=w;
        }
    //    SPFA();
        Dijkstra();
        memset(dp,0,sizeof(dp));
        int ans=dfs(st);
        printf("%d %d\n",ans,sum_value[ed]);
        Print(ed);
        puts("");
    }
    return 0;
}