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

49愛的鼓勵 6訂閱站台

內容 ：

輸入說明 ：

輸出說明 ：

範例輸入 ：

範例輸出 ：

提示 ：

出處 ：

擷取至 Matrix67 文章點此

脑海中出现的另一个牛B算法题则是POJ3318： 给你三个n*n的矩阵A、B、C，问你A*B是不是等于C。数据保证O(n^3)铁定超时，因此你需要想一个不用把A和B乘起来就可以验证的算法。一个基 于概率的算法是随机生成一个n乘1的矩阵R，然后判断A*B*R是否等于C*R，而前者相当于A*(B*R)，与后者一样都可以在O(n^2)的时间里算 出来。如果算出来的结果相等，几乎可以肯定A*B和C也是相等的。

/**********************************************************************************/
/*  Problem: a473 "外援計畫" from                                             */
/*  Language: C (1112 Bytes)                                                      */
/*  Result: AC(312ms, 6.2MB) judge by this@ZeroJudge                              */
/*  Author: morris1028 at 2012-10-08 14:39:11                                     */
/**********************************************************************************/


#include <stdio.h>
#define FOR(i,n) for(i=0;i<n;i++)
#define LL long long
#define N 505
LL A[N][N], B[N][N], C[N][N], R[N], D[N], E[N];
int n, i, j;
int solve() {
    FOR(i,n)
        R[i] = i+1, D[i] = E[i] = 0;
    FOR(i,n)FOR(j,n)
        E[i] += C[i][j]*R[j];
    FOR(i,n)FOR(j,n)
        D[i] += B[i][j]*R[j];
    FOR(i,n)    R[i] = 0;
    FOR(i,n)FOR(j,n)
        R[i] += A[i][j]*D[j];
    FOR(i,n)
        if(R[i] != E[i])
            return 0;
    return 1;
}
LL ReadInt(LL *x) {
    static char c, neg;
    while((c = getchar()) < '-')    {if(c == EOF) return EOF;}
    neg = (c == '-') ? -1 : 1;
    *x = (neg == 1) ? c-'0' : 0;
    while((c = getchar()) >= '0')
        *x = (*x << 3) + (*x << 1) + c-'0';
    *x *= neg;
    return 1;
}
int main() {
    while(scanf("%d", &n) == 1) {
        FOR(i,n)FOR(j,n)
            ReadInt(&A[i][j]);
        FOR(i,n)FOR(j,n)
            ReadInt(&B[i][j]);
        FOR(i,n)FOR(j,n)
            ReadInt(&C[i][j]);
        if(solve())
            puts("YES");
        else
            puts("NO");
    }
    return 0;
}