%I #25 Feb 09 2021 02:40:27
%S 49,99,125,243,244,343,351,424,476,549,604,725,775,845,846,847,1025,
%T 1251,1275,1324,1376,1421,1449,1520,1665,1675,1681,1682,1683,1849,
%U 1863,1925,2008,2024,2057,2107,2151,2224,2276,2349,2367,2524,2528,2575,2645,2825
%N Squareful numbers with squareful neighbors.
%H Amiram Eldar, <a href="/A235578/b235578.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) = A070258(n) + 1. - _Amiram Eldar_, Feb 09 2021
%e 243 is in the sequence because 242 = 2 * 11^2, 243 = 3^5, 244 = 2^2 * 61.
%e 244 is also in the sequence because 245 = 5 * 7^2.
%e But 245 is not in the sequence because 246 = 2 * 3 * 41.
%t Select[Range@3051, And @@ Not /@ SquareFreeQ /@ {#-1, #, #+1} &] (* _Giovanni Resta_, Jan 12 2014 *)
%t Select[Range[2000], Union[MoebiusMu[{# - 1, #, # + 1}]] == {0} &] (* _Alonso del Arte_, Jan 12 2014 *)
%o (PARI) isok(n) = !issquarefree(n) && !issquarefree(n-1) && !issquarefree(n+1); \\ _Michel Marcus_, Jan 12 2014
%Y Subsequence of A013929.
%Y Cf. A001694, A070258, A228649 (squarefree numbers with squarefree neighbors).
%K nonn
%O 1,1
%A _Juri-Stepan Gerasimov_, Jan 12 2014
%E Terms corrected by _Michel Marcus_, Jan 12 2014