login
Record high values in A078142.
1

%I #13 Sep 07 2022 17:25:13

%S 0,2,4,6,8,10,12,14,17,20,23,28,29,32,33,40,41,42,48,52,54,57,58,65,

%T 67,72,74,77,80,89,91,92,98,102,108,112,113,114,117,122,126,127,132,

%U 138,140,143,148,150,152,153,161,168,171,173,180,182,188,191,197,203,209

%N Record high values in A078142.

%H Amiram Eldar, <a href="/A076828/b076828.txt">Table of n, a(n) for n = 1..10000</a>

%H Amiram Eldar, <a href="/A076828/a076828.txt">Table of n, k, a(n) = A078142(k) for n=1..10000</a>

%t s[n_] := Total[Ceiling[Sqrt[(p = FactorInteger[n][[;; , 1]])]]^2 - p]; seq={}; sm = -1; Do[s1 = s[n]; If[s1 > sm, sm = s1; AppendTo[seq, s1]], {n, 1, 10^4}]; seq (* _Amiram Eldar_, Dec 08 2019 *)

%t f[n_]:=Module[{difs=Transpose[FactorInteger[n]][[1]]},Total[Ceiling[Sqrt[difs]]^2-difs]];DeleteDuplicates[Array[f,12000],GreaterEqual] (* _Harvey P. Dale_, Sep 07 2022 *)

%Y Cf. A078142.

%K easy,nonn

%O 1,2

%A _Jason Earls_, Nov 21 2002