|
%I #12 Apr 02 2020 03:11:29
%S 3,7,8,11,15,17,19,23,24,26,31,35,39,43,44,47,48,51,53,55,59,62,63,67,
%T 71,74,79,80,83,87,89,91,95,97,99,103,107,111,115,116,119,120,123,124,
%U 127,131,134,139,143,146,149,151,152,155,159,161,163,167,168,170,174
%N Numbers k such that the squarefree part of k is greater than the squarefree part of (k+1).
%H Amiram Eldar, <a href="/A069122/b069122.txt">Table of n, a(n) for n = 1..10000</a>
%t f[p_, e_] := If[OddQ[e], p, 1]; sqf[n_] := Times @@ (f @@@ FactorInteger[n]); Select[Range[200], sqf[#] > sqf[# + 1] &] (* _Amiram Eldar_, Apr 02 2020 *)
%o (PARI) for(n=1,220,if(core(n)>core(n+1),print1(n,",")))
%Y Cf. A007913.
%K easy,nonn
%O 1,1
%A _Benoit Cloitre_, Apr 07 2002
|
