OFFSET
1,2
COMMENTS
The relationship between a practical number and the largest primitive practical number that divides it such that the radical of the quotient is a divisor of same primitive provides an equivalence relation. Practical numbers that have the same progenitive primitive. This property arises from the characteristic of practical numbers and, in particular, primitives that says a practical number multiplied by power combinations of any of its divisors is also practical. This sequence identifies the primitive progenitor of each practical number A005153(n).
Note that this sequence and A378202 are similar but the first difference is at a(63) as explained in the example.
LINKS
Frank M Jackson, Table of n, a(n) for n = 1..10000
EXAMPLE
a(63) = 66. A005153(63) = 264 and the largest primitive practical number that divides the practical number 264 is 88. However the radical of the quotient 264/88 is 3 and 3 is not a divisor of 88. The next greatest primitive divisor of 264 is 66 and the radical of the quotient 264/66 is 2 and 2 is a divisor of 66.
a(131) = 306. A005153(131) = 612 and it is divisible by two primitive practical numbers 204 and 306 with their quotient a divisor of their primitive in both cases but 306 is chosen as the larger primitive.
MATHEMATICA
plst=Last/@ReadList["https://oeis.org/A005153/b005153.txt", {Number, Number}]; pplst=Last/@ReadList["https://oeis.org/A267124/b267124.txt", {Number, Number}]; Rad[n_] := Times @@ First /@ FactorInteger[n]; getpplst[n_] := Module[{}, Select[pplst, #<=n &]]; lst1={}; Do[lst=getpplst[plst[[n]]]; lnh=Length@lst; m=0; While[Mod[j=plst[[n]], k=lst[[lnh-m]]]!=0||Mod[k, Rad[j/k]]!=0, m++]; AppendTo[lst1, {j, k}], {n, 1, 100}]; Last/@lst1
CROSSREFS
KEYWORD
nonn
AUTHOR
Frank M Jackson, Dec 20 2024
STATUS
approved