timeit (N=10000) で計測すると、mjyさんの回答を、以下のようにnp.zeros
からnp.empty
に変更するのが一番早くなりました。260msぐらいで、他の回答より2倍以上早いです。今回のケースでは結果が全部置き換えられるのでnp.empty
を使って初期化しなくても問題ありません。
def by_np_slice2(N):
nums = np.arange(0, N, dtype=np.int)
result = np.empty((N, N-1), dtype=np.int)
for i, a in enumerate(result):
a[:i] = nums[:i]
a[i:] = nums[i+1:]
return result
また、コメントにあるnp.arange
を使ってfor
で繰り返し処理をする方法も360msぐらいで悪くはありません。for
の繰り返しでresult[i]
とnumpyのarray
をインデックスで扱う処理が遅いので、numba
を入れてJITコンパイルすると速くなります。230msぐらいで処理できるようになり上の処理よりも少し早くなりました。
@numba.jit
def by_np_arange(N):
result = np.empty((N, N-1),dtype=np.int)
for i in range(N):
result[i] = np.concatenate((np.arange(0,i), np.arange(i+1,N)))
return result
for文を使わない方法だと以下の単位3角行列を使うのが単純で速いのですが390msぐらいです。for文を使うなというのが常識ですが、場合によるようです。
Z = -np.tri(N, N-1, -1, dtype=np.int)
Z += np.arange(1,N)
Z = [range(0,i) + range(i+1,N) for i in range(N)]
とした場合、N=10000
で 3.79秒 -> 1.14秒になりました(timeit
で計測)。