|
%I #41 Apr 03 2023 10:36:09
%S 1,2,19,103,311,691,1321,2309,3671,5519,7919,10957,14753,19403,24809,
%T 31319,38873,47657,57559,69031,81799,96137,112291,130073,149717,
%U 171529,195043,220861,248851,279431,312583,347707,386093,427169,470933,517553
%N a(0)=1, a(n) = prime(n^3).
%C A sequence of increments for Shell sort that produces good results. A bit better than Sedgewick's A036562 and A003462.
%H David A. Corneth, <a href="/A055875/b055875.txt">Table of n, a(n) for n = 0..2200</a> (first 1001 terms from Ivan Panchenko)
%H <a href="/index/So#sorting">Index entries for sequences related to sorting</a>
%H Andrew Booker, <a href="https://t5k.org/nthprime/index.php#n-th">Nth Prime Page</a>.
%F a(n) = A000040(A000578(n)), n>0.
%t {1}~Join~Array[Prime[#^3] &, 35] (* _Michael De Vlieger_, Apr 13 2021 *)
%o (Magma) [NthPrime(n^3): n in [0..50] ]; // _Vincenzo Librandi_, Apr 22 2011
%o (PARI) first(n) = { my(res = vector(n), t = 0); forprime(p = 2, oo, t++; if(ispower(t, 3, &i), print1([i, p]", "); res[i] = p; if(i >= n, return(concat(1, res))))) } \\ _David A. Corneth_, Apr 13 2021
%Y Cf. A033844, A036562, A003462, A036567, A036569, A055876.
%Y Sequences used for Shell sort: A003462, A033622, A036562, A036564, A036569, A055875.
%Y Cf. A011757, A109724.
%K nonn
%O 0,2
%A Steven Pigeon (pigeon(AT)iro.umontreal.ca), Jul 14 2000
%E More terms from _Jonathan Vos Post_, Aug 13 2005
|
