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

49愛的鼓勵 6訂閱站台

這題不算難，難是難在卡時間，很清楚地我們知道要進行矩陣連乘積，
也就是 N*N 的矩陣連乘 E-S 次，最後乘上一個 N*1 的矩陣。
不能直接做矩陣乘法，因為會消耗 O(N*N*N) 計算兩個 N*N 相乘。

在此倒過來運算 (N*N)*(N*1) = (N*1) 消耗 O(N*N)，
逆序運算得 O(N*N*(S-E))。
直接算 O(N*N*N*(S-E))。

或許有人會想用線段樹降低 O(E-S)，在此先說不僅僅是 MLE 而且還會 TLE。

---

#include <stdio.h>
#include <string.h>
struct mx {
    unsigned char v[55][55];
};
mx mult(mx A, mx B, int a, int b, int c){
    mx C;
    int i, j, k;
    memset(&C, 0, sizeof(C));
    for(i = 0; i < a; i++) {
        for(j = 0; j < c; j++) {
            for(k = 0; k < b; k++)
                C.v[i][j] += A.v[i][k]*B.v[k][j];
        }
    }
    return C;
}
mx A[855], B, one;
int M, N, K, S, E;
/*
3 2 1
10 100 200 255
8 6 2 3
1 0 0 1
1 0 2 8
*/
int main() {
    int i, j, k, x;
    memset(&one, 0, sizeof(one));
    for(i = 0; i < 55; i++)
        one.v[i][i] = 1;
    while(scanf("%d %d %d", &M, &N, &K) == 3) {
        for(i = 0; i < M; i++) {
            for(j = 0; j < N; j++)
                for(k = 0; k < N; k++) {
                    scanf("%d", &x);
                    A[i].v[j][k] = (unsigned char)x;
                }
        }
        mx ans;
        while(K--) {
            scanf("%d %d", &S, &E);
            for(i = 0; i < N; i++) {
                scanf("%d", &x);
                B.v[i][0] = x;
            }
            ans = B;
            if(S <= E) {
                for(i = E; i >= S; i--)
                ans = mult(A[i], ans, N, N, 1);
            } else {
                for(i = E; i >= 0; i--)
                ans = mult(A[i], ans, N, N, 1);
                for(i = M-1; i >= S; i--)
                ans = mult(A[i], ans, N, N, 1);
            }
            for(i = 0; i < N; i++) {
                if(i)
                    printf(" ");
                printf("%d", ans.v[i][0]);
            }
            puts("");
        }
    }
    return 0;
}