OFFSET
1,1
LINKS
Donovan Johnson, Table of n, a(n) for n = 1..1000
EXAMPLE
n = pqrst, p<q<r<s<t, primes, p+q+r+s+t = kt; n = 8970 = 2*3*5*13*23, sum = 46 = 2*23.
MATHEMATICA
ffi[x_] := Flatten[FactorInteger[x]] lf[x_] := Length[FactorInteger[x]] ba[x_] := Table[Part[ffi[x], 2*w-1], {w, 1, lf[x]}] sb[x_] := Apply[Plus, ba[x]] ma[x_] := Part[Reverse[Flatten[FactorInteger[x]]], 2] amo[x_] := Abs[MoebiusMu[x]] Do[s=sb[n]/ma[n]; If[IntegerQ[s]&&Equal[lf[n], 5]&& !Equal[amo[n], 0], Print[{n, ba[n]}]], {n, 2, 1000000}]
sdpQ[n_]:=Module[{fi=FactorInteger[n][[;; , 1]]}, Length[fi]==5&&SquareFreeQ[n]&&Mod[Total[ fi], Max[fi]]==0]; Select[Range[150000], sdpQ] (* Harvey P. Dale, May 04 2023 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, May 13 2002
STATUS
approved