|
MAPLE
|
# see http://www-maths.swan.ac.uk/pgrads/bb/project/node43.html perfp := [1, 2, 4, 6, 12, 16, 18, 30, 60, 88, 106, 126, 520, 606, 1278, 2202, 2280, 3216, 4252] : A000396 := proc(n) global perfp ; 2^op(n, perfp)*(2^(op(n, perfp)+1)-1) ; end: A134705 := proc(n) A000396(n)/2^n ; end: seq(A134705(n), n=1..12) ; - R. J. Mathar
|