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.
%I A183099 #15 Aug 29 2023 04:20:31 %S A183099 0,0,0,4,0,0,0,12,9,0,0,4,0,0,0,28,0,9,0,4,0,0,0,12,25,0,36,4,0,0,0, %T A183099 60,0,0,0,49,0,0,0,12,0,0,0,4,9,0,0,28,49,25,0,4,0,36,0,12,0,0,0,4,0, %U A183099 0,9,124,0,0,0,4,0,0,0,129,0,0,25,4,0,0,0,28,117,0,0,4,0,0,0,12,0,9,0,4,0,0,0,60,0,49,9,129 %N A183099 a(n) = sum of powerful divisors d (excluding 1) of n. %C A183099 a(n) = sum of divisors d of n from set A001694(m) - powerful numbers for m >=2. %H A183099 Antti Karttunen, <a href="/A183099/b183099.txt">Table of n, a(n) for n = 1..16385</a> %H A183099 <a href="/index/Su#sums_of_divisors">Index entries for sequences related to sums of divisors</a>. %F A183099 a(n) = A000203(n) - A183100(n) = A183097(n) - 1. %F A183099 a(1) = 0, a(p) = 0, a(pq) = 0, a(pq...z) = 0, a(p^k) = ((p^(k+1)-1) / (p-1))-p-1, for p, q = primes, k = natural numbers, pq...z = product of k (k > 2) distinct primes p, q, ..., z. %e A183099 For n = 12, set of such divisors is {4}; a(12) = 4. %t A183099 f[p_, e_] := (p^(e+1)-1)/(p-1) - p; a[1] = 0; a[n_] := Times @@ f @@@ FactorInteger[n] - 1; Array[a, 100] (* _Amiram Eldar_, Aug 29 2023 *) %o A183099 (PARI) A183099(n) = (sumdiv(n, d, ispowerful(d)*d) - 1); \\ _Antti Karttunen_, Oct 07 2017 %Y A183099 Cf. A000203, A001694, A183097, A183100. %K A183099 nonn,easy %O A183099 1,4 %A A183099 _Jaroslav Krizek_, Dec 25 2010