160426

Ruby

Wilson's theorem

オンライン整数列大辞典の
A061006（http://oeis.org/A061006/list）
と比較し、答え合わせしてみる。

def A061006(n)
  ary = [0]
  (2..n).each{|i|
    j = 1
    (1..i - 1).each{|k|
      j *= k
      j %= i
    }
    ary << j
  }
  ary
end
ary = A061006(94)

# OEIS A061006のデータ
ary0 =
[0,1,2,2,4,0,6,0,0,0,10,0,12,0,0,0,16,0,18,0,0,0,
 22,0,0,0,0,0,28,0,30,0,0,0,0,0,36,0,0,0,40,0,42,0,
 0,0,46,0,0,0,0,0,52,0,0,0,0,0,58,0,60,0,0,0,0,0,
 66,0,0,0,70,0,72,0,0,0,0,0,78,0,0,0,82,0,0,0,0,0,
 88,0,0,0,0,0]
# 一致の確認
p ary == ary0