cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A332622 Numbers where records occur for the product of exponential divisors function (A157488).

This page as a plain text file.
%I A332622 #40 Jun 06 2020 03:19:13
%S A332622 1,2,3,4,8,9,12,16,20,28,36,72,100,144,324,400,576,900,1764,2700,3528,
%T A332622 3600,7056,10800,14400,28224,32400,44100,88200,108900,129600,176400,
%U A332622 396900,435600,608400,705600,1587600,3920400,5336100,5475600,6350400,14288400,15681600
%N A332622 Numbers where records occur for the product of exponential divisors function (A157488).
%C A332622 The corresponding record values are 1, 2, 3, 8, 16, 27, 72, 128, 200, 392, 46656, 186624, 1000000, ...
%H A332622 Amiram Eldar, <a href="/A332622/b332622.txt">Table of n, a(n) for n = 1..68</a> (terms below 10^10)
%F A332622 The first 8 terms of A157488 are 1, 2, 3, 8, 5, 6, 7 and 16. The record values occur at 1, 2, 3, 4 and 8 - the first 5 terms of this sequence.
%t A332622 f[p_, e_] := p^(DivisorSigma[1, e]/DivisorSigma[0, e]); exprod[n_] := (Times @@ (f @@@ (fct = FactorInteger[n])))^(Times @@ DivisorSigma[0, Last /@ fct]); em = 0; s = {}; Do[If[(e = exprod[n]) > em, em = e; AppendTo[s, n]], {n, 1, 10^6}]; s
%Y A332622 Cf. A034287, A061742, A157488, A328134, A318278, A322791.
%K A332622 nonn
%O A332622 1,2
%A A332622 _Amiram Eldar_, Jun 05 2020