%I #13 Aug 07 2015 03:00:37
%S 1,0,2,5,4,5,7,7,2,2,2,4,4,4,6,5,6,7,7,2,2,3,2,6,4,4,6,6,7,6,4,7,7,4,
%T 5,3,4,6,5,6,7,7,2,2,2,3,2,5,3,3,2,7,4,2,3,6,5,2,4,4,5,4,7,6,3,4,8,5,
%U 5,7,3,4,6,5,7,5,2,6,7,3,4,3,3,6,4,5,7,7,6,2,2,2,2,3,2,7,7,6,2,5,2,2,3,4,3
%N a(n) = 3 + (H(n) mod 6) + floor(r) where H()=A005185() and r = (H(n) - 2*H(n+1) + H(n+2) - 4) / H(n).
%p A005185 := proc(n)
%p option remember;
%p if n<=2 then
%p 1
%p elif n > procname(n-1) and n > procname(n-2) then
%p procname(n-procname(n-1))+procname(n-procname(n-2));
%p end if;
%p end proc:
%p A086267 := proc(n)
%p local H ;
%p H := A005185(n) ;
%p H-2*A005185(n+1)+A005185(n+2)-4;
%p %/H ;
%p 3+ floor(%)+ (H mod 6) ;
%p end proc:
%p seq(A086267(n),n=1..50) ; # _R. J. Mathar_, Oct 10 2011
%t Hofstadter[n_Integer?Positive] := Hofstadter[n] = Hofstadter[n - Hofstadter[n-1]] + Hofstadter[n - Hofstadter[n-2]] Hofstadter[1] = Hofstadter[2] = 1 Digits=502 a=Table[Hofstadter[n], {n, 1, Digits}]; b=Table[Floor[(a[[n]]-2*a[[n+1]]+a[[n+2]]-4)/a[[n]]]+Mod[a[[n]], 6]+3, {n, 1, Digits-2}] ListPlot[b]
%K nonn,obsc
%O 1,3
%A _Roger L. Bagula_, Aug 28 2003