OFFSET
1,2
LINKS
FORMULA
b(1)=0, b(n)=(b(n-1) * 214013 + 2531011) mod 2^32. a(n) = (floor(b(n)/65536) mod 32768). The sequence is periodic with period length 2^32.
a(n) = floor(A096557(n)/2^16) mod 2^15 = floor((2531011*(214013^(n-1)-1)/214012 mod 2^32)/2^16) mod 2^15. - M. F. Hasler, May 14 2015
MAPLE
b:= proc(n) option remember; `if`(n<2, 0,
irem(214013 *b(n-1) +2531011, 4294967296))
end:
a:= n-> irem(iquo(b(n), 65536), 32768):
seq(a(n), n=1..50); # Alois P. Heinz, Jun 10 2014
MATHEMATICA
A096557 = NestList[Mod[#*214013 + 2531011, 2^32] &, 0, 50];
Mod[BitShiftRight[A096557, 16], 2^15] (* Paolo Xausa, Aug 29 2024 *)
PROG
(PARI) a(n)=A096557(n)>>16%2^15 \\ M. F. Hasler, May 14 2015
(PARI) A096558(n)=lift((Mod(214013, 2^34)^(n-1)-1)*13129821757)>>18%32768 \\ M. F. Hasler, May 14 2015
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Hugo Pfoertner, Jul 21 2004
STATUS
approved