|
|
A028948
|
|
An "extremely strange sequence": a(n+1) = [ A*a(n)+B ]/p^r, where p^r is the highest power of p dividing [ A*a(n)+B ] and p=2, A=4.001, B=1.2.
|
|
3
|
|
|
0, 1, 5, 21, 85, 341, 1365, 2731, 10927, 5465, 10933, 1367, 2735, 10943, 5473, 10949, 1369, 2739, 10959, 5481, 10965, 1371, 2743, 10975, 5489, 10981, 1373, 2747, 10991, 5497, 10997, 1375, 2751, 11007, 5505, 11013, 1377, 2755, 11023, 5513, 11029, 1379, 2759
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
LINKS
|
|
|
FORMULA
|
|
|
PROG
|
(PARI) vector(50, i, if(i>1, t=(t*4001+1200)\1000; t>>=valuation(t, 2), t=0)) \\ M. F. Hasler, Aug 07 2018
(Python)
from itertools import islice
def A028948_gen(): # generator of terms
x = 0
while True:
yield x
x = (y:=(x*4001+1200)//1000)>>(~y&y-1).bit_length()
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|