|
%I #20 Feb 25 2024 01:44:25
%S 1,1,1,1,1,1,1,2,1,1,1,1,1,1,1,2,2,1,1,1,1,1,2,1,1,1,1,1,1,1,2,1,1,2,
%T 2,1,2,1,1,1,1,1,1,1,3,1,1,1,2,1,2,1,1,1,1,1,1,1,2,1,1,1,1,1,1,3,1,2,
%U 1,1,1,1,2,1,1,1,1,1,1,1,3,1,1,1,1,1,1,2,2,1,1,1,1,1,2,1,2,1,1,1,1,2,1,1,1
%N First differences of indices of squarefree Jacobsthal numbers.
%C First differences of A107037.
%H Amiram Eldar, <a href="/A107039/b107039.txt">Table of n, a(n) for n = 1..898</a>
%p A001045 := proc(n) option remember : if n <= 1 then n; else A001045(n-1)+2*A001045(n-2) ; fi ; end: A107037 := proc(nmax) local a,n; a := [] ; n := 1: while nops(a) < nmax do if numtheory[issqrfree]( A001045(n) ) then a := [op(a),n] ; fi ; n := n+1: od : RETURN(a) ; end: A107039 := proc(nmax) local a107037 ; a107037 := A107037(nmax+1) : [seq( a107037[i]-a107037[i-1],i=2..nmax)] ; end: op(A107039(120)) ; # _R. J. Mathar_, Aug 12 2007
%t Position[LinearRecurrence[{1, 2}, {0, 1}, 131], _?SquareFreeQ]-1 // Rest // Flatten // Differences (* _Jean-François Alcover_, Dec 09 2023 *)
%Y Cf. A001045, A005117, A107037.
%K nonn
%O 1,8
%A _Paul Barry_, May 09 2005
%E More terms from _R. J. Mathar_, Aug 12 2007
%E Offset corrected by _Amiram Eldar_, Feb 25 2024
|
