Morris' Blog 我在摸索我存在的意義、生命的意義、內涵及價值。 沒有真正神手般的厲害，套用模式是大部分的情況。 沒那麼厲害的我存在的意義為何，襯托出神的存在？ 有時候不能怪別人不給我機會，我已經開始省思， 為什麼我沒有給別人機會的原因，而是全要靠別人。 有時候我認為我與眾不同，那是因為我做不到很多人能辦到的事情。 我聽不到、玩不起來、說得不夠明瞭、理解得不夠多， 我缺少的語文能力，就是一切一切無能的判定。 為此，我不曉得能做些對別人有幫助的事情。 我知道我自己是個不服輸的性格，但現在我不得不認輸。 我沒有能力，因此我必須認輸。 我到底來這世界作什麼？處處充滿不確定性， 而我總是無法出類拔萃，想當個一般人似乎也回不了頭了。 越強大的挑戰需要越強力的支持，而我的支持是什麼？ 最近嘗試與人談話，但都是失敗的結局，沒人能接受， 不需要言語，不需要溝通，一切所想都表達不出來。

49愛的鼓勵 6訂閱站台

a259. Q10917: A Walk Through the Forest

內容 ：

輸入說明 ：

輸出說明 ：

範例輸入 ：

範例輸出 ：

提示 ：

出處 ：

---

作法 : SPFA + DP
搜尋所有點到 2 的距離，接著使用 DP，計算到點 1 的方法數

/**********************************************************************************/
/*  Problem: a259 "Q10917: A Walk Through the Forest" from UVa/ACM , Lucky cat 翻譯*/
/*  Language: C (1888 Bytes)                                                      */
/*  Result: AC(440ms, 608KB) judge by this@ZeroJudge                              */
/*  Author: morris1028 at 2011-10-28 14:26:37                                     */
/**********************************************************************************/


#include<stdio.h>
#include<stdlib.h>
#define oo 1000000000
typedef struct {
    int x, y, v;
}Pair;
Pair Path[1000000];
int cmp(const void *i, const void *j) {
    Pair *x, *y;
    x = (Pair *)i, y = (Pair *)j;
    if(x->x < y->x)        return -1;
    else if(x->x == y->x)
            return x->v - y->v;
    else    return 1;
}
int Min(int x, int y) {
    return x < y ? x : y;
}
int NodeT[1001], NodeL[1001], Q[1000001];
int SPFA(int st, int ed, int N) {
    int i, j, k;
    static int Used[1001], V[1001], QIdx;
    static int DP[1001], tNode;
    for(i = 1; i <= N; i++)
        Used[i] = 0, V[i] = oo, DP[i] = 0;
    Used[st] = 1, V[st] = 0, QIdx = 0, Q[QIdx] = st;
    for(i = 0; i <= QIdx; i++) {
        tNode = Q[i], Used[tNode] = 0;
        for(j = 0, k = NodeL[tNode]; j < NodeT[tNode]; j++, k++) {
            if(V[tNode] + Path[k].v < V[Path[k].y]) {
                V[Path[k].y] = V[tNode] + Path[k].v;
                if(!Used[Path[k].y]) {
                    Used[Path[k].y] = 1, Q[++QIdx] = Path[k].y;
                }
            }
        }
    }
    static Pair Tmp[1001];
    for(i = 1; i <= N; i++)
        Tmp[i-1].x = V[i], Tmp[i-1].y = i;;
    DP[st] = 1;
    qsort(Tmp, N, sizeof(Pair), cmp);
    for(i = 0; i < N; i++) {
        tNode = Tmp[i].y;
        for(j = 0, k = NodeL[tNode]; j < NodeT[tNode]; j++, k++) {
            if(V[tNode] < V[Path[k].y] && V[Path[k].y] != oo) {
                DP[Path[k].y] += DP[tNode];
            }
        }
        if(tNode == ed)    break;
    }
    return DP[ed];
}
int main() {
    int N, M, i;
    while(scanf("%d", &N) == 1 && N) {
        scanf("%d", &M);
        for(i = 1; i <= N; i++)
            NodeL[i] = oo, NodeT[i] = 0;
        for(i = 0; i < M; i++) {
            scanf("%d %d %d", &Path[i].x, &Path[i].y, &Path[i].v);
            Path[M+i].x = Path[i].y;
            Path[M+i].y = Path[i].x;
            Path[M+i].v = Path[i].v;
        }
        qsort(Path, 2*M, sizeof(Pair), cmp);
        for(i = 0, M *= 2; i < M; i++) {
            NodeL[Path[i].x] = Min(i, NodeL[Path[i].x]);
            NodeT[Path[i].x] ++;
        }
        printf("%d\n", SPFA(2, 1, N));
    }
    return 0;
}