Нашел решение, может кому пригодится:
Код:
default(parisizemax,10^9)
n=10
l=vector(2^n,i,0)
for(i=2,2^n,l[i]=l[i\((1+sqrt(5))/2)]+1)
l1=vector(2^n,i,0)
l1[1]=1
for(i=1,2^n-1,l1[i+1]=2+l[i])
fib=vector(2^n,i,0)
fib[1]=1
fib[2]=1
for(i=2,2^n-1,fib[i+1]=fib[i]+fib[i-1])
l2=vector(2^n,i,0)
l2[1]=1
for(i=1,2^n-1,l2[i+1]=fib[l1[i+1]])
l3=vector(2^n,i,0)
for(i=0,2^n-1,l3[i+1]=i+1-l2[i+1])
l4=vector(2^n,i,0)
for(i=1,2^n-1,l4[i+1]=(l1[i+1]-l1[l3[i+1]])%2)
l5=vector(2^n,i,0)
l5[1]=2
for(i=2,2^n,l5[i]=l5[l3[i]]+l4[i])
l6=vector(2^n,i,0)
l6[1]=1
for(i=1,2^n-1,l6[i+1]=l5[i]*l6[l3[i]+1])
print(l5)
Последняя функция дает уникальные (после удаления небольшого числа дубликатов например в Excel) члены последовательности
вот отсюда в порядке их появления в ней. В исходной последовательности они ограничены
, в здешней -
(числа Фибоначчи). Вот их число для малых
:
Последовательности нет в OEIS. С ростом
число удаляемых дубликатов растет.