|
MAPLE
|
read transforms : A000043 := proc(n) op(n, [2, 3, 5, 7, 13, 17, 19, 31, 61, 89, 107, 127, 521, 607, 1279, 2203, 2281, 3217, 4253, 4423, 9689, 9941, 11213, 19937]) ; end: A000668 := proc(n) 2^A000043(n)-1 ; end: A135623 := proc(n) digrev(A000668(n)) ; end: A135625 := proc(n) A000668(n)*A135623(n) ; end: for n from 1 to 13 do printf("%d, ", A135625 (n) ) ; od: - R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Feb 28 2008
|