結果

問題 No.217 魔方陣を作ろう
ユーザー keikei
提出日時 2018-09-24 15:17:56
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 509 ms / 5,000 ms
コード長 884 bytes
コンパイル時間 90 ms
コンパイル使用メモリ 12,800 KB
実行使用メモリ 44,480 KB
最終ジャッジ日時 2024-09-22 03:54:29
合計ジャッジ時間 10,746 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 507 ms
44,140 KB
testcase_01 AC 506 ms
44,012 KB
testcase_02 AC 502 ms
44,260 KB
testcase_03 AC 503 ms
43,756 KB
testcase_04 AC 501 ms
44,264 KB
testcase_05 AC 507 ms
43,884 KB
testcase_06 AC 503 ms
43,884 KB
testcase_07 AC 507 ms
43,888 KB
testcase_08 AC 504 ms
44,140 KB
testcase_09 AC 508 ms
44,392 KB
testcase_10 AC 509 ms
44,480 KB
testcase_11 AC 509 ms
43,760 KB
testcase_12 AC 504 ms
43,888 KB
testcase_13 AC 501 ms
44,136 KB
testcase_14 AC 500 ms
43,880 KB
testcase_15 AC 500 ms
43,760 KB
testcase_16 AC 498 ms
44,136 KB
testcase_17 AC 502 ms
44,428 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import numpy as np
def magic(n):
    n = int(n)
    if n % 2 == 1:
        p = np.arange(1, n+1)
        return n*np.mod(p[:, None] + p - (n+3)//2, n) + np.mod(p[:, None] + 2*p-2, n) + 1
    elif n % 4 == 0:
        J = np.mod(np.arange(1, n+1), 4) // 2
        K = J[:, None] == J
        M = np.arange(1, n*n+1, n)[:, None] + np.arange(n)
        M[K] = n*n + 1 - M[K]
    else:
        p = n//2
        M = magic(p)
        M = np.block([[M, M+2*p*p], [M+3*p*p, M+p*p]])
        i = np.arange(p)
        k = (n-2)//4
        j = np.concatenate((np.arange(k), np.arange(n-k+1, n)))
        M[np.ix_(np.concatenate((i, i+p)), j)] = M[np.ix_(np.concatenate((i+p, i)), j)]
        M[np.ix_([k, k+p], [0, k])] = M[np.ix_([k+p, k], [0, k])]
    return M

n = int(input())
ret = magic(n);

for i in range(n):
    for j in range(n):
        print(ret[i][j],end=" ")
    print("")
0