OFFSET
1,1
COMMENTS
The probability that a random integer is cubefree is 1/zeta(3) (see A088453).
The asymptotic density of this sequence is 1 - 15/(7*zeta(3)) + (7/8)*Product_{p odd prime} (1-2/p^3) = 1 - (15/7)*A088453 + (7/6)*A340153 = 0.00704953764810822631953... . - Amiram Eldar, Aug 20 2026
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..700 from R. J. Mathar)
EXAMPLE
312 is a term because 2^3 divides 312 and 5^3 divides 312*2+1.
MATHEMATICA
noncfQ[n_] := AnyTrue[FactorInteger[n][[;; , 2]], # > 2 &]; Select[Range[7000], noncfQ[#] && noncfQ[2*# + 1] &] (* Amiram Eldar, May 25 2025 *)
PROG
(PARI) iscubefree(k) = if(k == 1, 1, vecmax(factor(k)[, 2]) < 3);
isok(k) = !iscubefree(k) && !iscubefree(2*k+1); \\ Amiram Eldar, Aug 20 2026
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Jonathan Vos Post, Mar 03 2006
EXTENSIONS
Edited by Don Reble, Mar 05 2006
2875 inserted by R. J. Mathar, Dec 08 2015
STATUS
approved