結果
| 問題 |
No.223 1マス指定の魔方陣
|
| コンテスト | |
| ユーザー |
kei
|
| 提出日時 | 2018-09-24 16:11:06 |
| 言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 1,525 bytes |
| コンパイル時間 | 389 ms |
| コンパイル使用メモリ | 12,800 KB |
| 実行使用メモリ | 44,232 KB |
| 最終ジャッジ日時 | 2024-09-22 05:12:33 |
| 合計ジャッジ時間 | 28,646 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 1 |
| other | AC * 29 WA * 17 |
ソースコード
'''
<url:https://yukicoder.me/problems/no/217>
問題文============================================================
=================================================================
解説=============================================================
参考サイトや既存の魔法陣のテンプレ貼り付けてほい
================================================================
'''
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,x,y,z = input().strip().split(' ')
n,x,y,z =[int(n),int(x),int(y),int(z)]
x-=1
y-=1
ans = magic(n)
kx,ky = np.where(ans==z)
kx = int(kx)
ky = int(ky)
kx += n - y
ky += n - x
ans = np.concatenate([ans,ans,ans])
ans = np.concatenate([ans,ans,ans],axis=1)
for i in range(n):
for j in range(n):
print(ans[i+kx][j+ky],end=" ")
print("")
kei