|
%I #21 Mar 07 2020 11:17:46
%S 1,0,1,-1,2,-1,1,0,0,1,-2,3,-2,2,-2,3,-3,4,-3,3,-2,1,1,-2,4,-4,5,-4,4,
%T -3,4,-5,6,-5,5,-5,6,-6,7,-6,6,-5,4,-2,1,0,-1,2,-1,0,2,-3,2,0,-1,3,-4,
%U 5,-4,3,-1,0,1,-1,2,-3,5,-6,7,-7,8,-7,7,-6,5,-3
%N a(1) = 1 and for n > 0, a(n+1) = a(k) - a(n) where k is the number of terms equal to a(n) among the first n terms.
%C In other words, for n > 0, a(n+1) = a(o(n)) - a(n) where o is the ordinal transform of the sequence.
%C The sequence has interesting graphical features (see plot in Links section).
%H Rémy Sigrist, <a href="/A329985/b329985.txt">Table of n, a(n) for n = 1..25000</a>
%H Rémy Sigrist, <a href="/A329985/a329985.png">Density plot of the first 2^22 terms</a>
%H N. J. A. Sloane (in collaboration with Scott R. Shannon), <a href="/A331452/a331452.pdf">Art and Sequences</a>, Slides of guest lecture in Math 640, Rutgers Univ., Feb 8, 2020. Mentions this sequence.
%e The first terms, alongside their ordinal transform, are:
%e n a(n) o(n)
%e -- ---- ----
%e 1 1 1
%e 2 0 1
%e 3 1 2
%e 4 -1 1
%e 5 2 1
%e 6 -1 2
%e 7 1 3
%e 8 0 2
%e 9 0 3
%e 10 1 4
%t A={1};For[n=2,n<=76,n++,A=Append[A,Part[A,Count[Table[Part[A,i],{i,1,n-1}],Part[A,n-1]]]-Part[A,n-1]]];A (* _Joshua Oliver_, Nov 26 2019 *)
%t Nest[Append[#, #[[Count[#, #[[-1]] ] ]] - #[[-1]]] &, {1}, 75] (* _Michael De Vlieger_, Dec 01 2019 *)
%o (PARI) for (n=1, #(a=vector(76)), print1 (a[n]=if (n==1, 1, a[sum(k=1, n-1, a[k]==a[n-1])]-a[n-1])", "))
%Y o(n) is A330334.
%Y See A329981 for similar sequences.
%K sign,look
%O 1,5
%A _Rémy Sigrist_, Nov 26 2019
|
