OFFSET
1,2
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..16384
FORMULA
a(p) = 1 for p prime, p > 2.
a(2*p) = 2 for p prime, p > 3.
a(2*3^k) = k + 2, k >= 1;
a(2*p^k) = 2, k >= 1, p prime, p >= 5.
a(2^n) = n + 1.
a(n) = Sum_{d|n} A322860(d). - Antti Karttunen, Sep 11 2023
EXAMPLE
MATHEMATICA
f[p_, e_] := (p^(e + 1) - 1)/(p - 1); pracQ[n_] := (ind = Position[(fct = FactorInteger[n])[[;; , 1]]/(1 + FoldList[Times, 1, f @@@ Most@fct]), _?(# > 1 &)]) == {}; a[n_] := DivisorSum[n, 1 &, pracQ[#] &]; Array[a, 100] (* Amiram Eldar, Aug 21 2023 *)
PROG
(Magma) sk:=func<n, k|&+[Divisors(n)[i]:i in [1..k]]>; f:=func<n|forall{k: k in [2..#Divisors(n)]|sk(n, k-1) ge Divisors(n)[k]-1}>; [#[d:d in Divisors(n)|f(d)]:n in [1..100]];
a(n) = sumdiv(n, d, is_A005153(d)); \\ Michel Marcus, Sep 11 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Marius A. Burtea, Aug 18 2023
EXTENSIONS
Data section extended up to a(105) by Antti Karttunen, Jun 02 2024
STATUS
approved