login
a(n) is the smallest multiple of n that is a term of A072873.
2

%I #15 Feb 16 2024 05:27:26

%S 1,4,27,4,3125,108,823543,16,27,12500,285311670611,108,

%T 302875106592253,3294172,84375,16,827240261886336764177,108,

%U 1978419655660313589123979,12500,22235661,1141246682444,20880467999847912034355032910567,432,3125,1211500426369012,27,3294172

%N a(n) is the smallest multiple of n that is a term of A072873.

%H Amiram Eldar, <a href="/A365636/b365636.txt">Table of n, a(n) for n = 1..388</a>

%F Multiplicative with a(p^e) = p^(p*ceiling(e/p)).

%F a(n) = n * A365637(n).

%F a(n) >= n with equality if and only if n is in A072873.

%F Sum_{n>=1} 1/a(n) = Product_{p prime} (1 + p/(p^p-1)) = 1.86196549645040699446... .

%t f[p_, e_] := p^(Ceiling[e/p]*p); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 30]

%o (PARI) a(n) = {my(f = factor(n)); prod(i = 1, #f~, f[i,1]^(f[i,1] * ceil(f[i,2] / f[i,1])));}

%Y Cf. A072873, A365637.

%K nonn,easy,mult

%O 1,2

%A _Amiram Eldar_, Sep 14 2023

%E Data, formulas and codes corrected by _Amiram Eldar_, Feb 15 2024