[seqfan] observation about A003160 a(1)=a(2)=1, a(n)=n-a(a(n-1))-a(a(n-2))

Joerg Arndt arndt at jjj.de
Fri Apr 30 13:01:40 CEST 2010


Compute A003160:
? N=100; v=vector(N); v[1]=v[2]=1;
? for(n=3,N,v[n]=n-v[v[n-1]]-v[v[n-2]]); v
  [1, 1, 1, 2, 3, 4, 4, 4, 5, 5, 5, 6, 6, 6, 7, 8, ...]

Positions where seq increases:
? for(n=2,N,if(v[n]-v[n-1]>0,print1(n,", ")))
  4, 5, 6, 9, 12, 15, 16, 17, 20, 21, 22, 25, 26, ...

This appears to be A003156 (except initial term 1):
 a(n) is the number of ones before the n-th zero in
 the Feigenbaum sequence A035263. 

Initial term is obtained if we prepend a zero at the
beginning of A003160





More information about the SeqFan mailing list