|
%I #8 Dec 23 2018 16:13:08
%S 1,1,1,1,1,2,1,1,2,1,1,2,1,1,2,1,1,2,1,1,2,1,1,2,1,1,2,1,1,2,1,1,2,1,
%T 1,2,1,1,2,1,1,2,1,1,2,1,1,2,1,1,2,1,1,2,1,1,2,1,1,2,1,1,2,1,1,2,1,1,
%U 2,1,1,2,1,1,2,1,1,2,1,1,2,1,1,2,1,1,2,1,1,2,1,1,2,1,1,2,1,1,2,1,1,2,1,1,2
%N a(n) = A107750(n+1) - A107750(n).
%H Antti Karttunen, <a href="/A107751/b107751.txt">Table of n, a(n) for n = 0..65537</a>
%F a(n) = if n<=4 then 1 else 0^(a(n-1)-1) + 0^(a(n-2)-1).
%t A107751[ nmax_ ] := Length /@ Split[ (Total[ 1 - IntegerDigits[ #1, 2 ] ] &) /@ Range[ 0, nmax ] ]; A107751[ 200 ] - Peter Pein (petsie(AT)dordos.net), Oct 12 2007
%o (PARI)
%o up_to = 65537;
%o A107751list(up_to) = { my(v=vector(up_to)); v[1]=v[2]=v[3]=v[4]=1; for(n=5,up_to,v[n] = (0^(v[n-1]-1) + 0^(v[n-2]-1))); (v); };
%o v107751 = A107751list(up_to);
%o A107751(n) = if(!n,1,v107751[n]); \\ _Antti Karttunen_, Dec 23 2018
%Y Cf. A100063, A099837, A100051, A107688.
%K nonn
%O 0,6
%A _Reinhard Zumkeller_, May 23 2005
|
