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.

A094889 Odd nonunitary abundant numbers.

This page as a plain text file.
%I A094889 #10 May 12 2019 06:51:04
%S A094889 33075,99225,165375,231525,259875,297675,363825,429975,467775,496125,
%T A094889 552825,562275,571725,606375,628425,675675,694575,716625,760725,
%U A094889 779625,798525,826875,848925,883575,893025,921375,937125,959175,1003275
%N A094889 Odd nonunitary abundant numbers.
%H A094889 Amiram Eldar, <a href="/A094889/b094889.txt">Table of n, a(n) for n = 1..10000</a>
%H A094889 J. A. Sellers, <a href="http://www.math.psu.edu/sellersj/mcedonan.pdf">A Note On Infinitely Many Odd Nonunitary Abundant Numbers</a>
%t A094889 f[p_, e_] := (p^(e + 1) - 1)/(p - 1); nuabQ[n_] := Times @@ (f @@@ FactorInteger[n]) - Times @@ (1 + Power @@@ FactorInteger[n]) > n; Select[Range[3, 10^6, 2], nuabQ] (* _Amiram Eldar_, May 12 2019 *)
%o A094889 (PARI) { usigma(n)=local(s=1,fac,i); fac=factor(n); for(i=1,matsize(fac)[1], s=s*(1+fac[i,1]^fac[i,2])); return(s); } nusigma(n)=sigma(n)-usigma(n); forstep(n=1,2^24,2,if(nusigma(n)>n,print1(n",")))
%Y A094889 Cf. A064597.
%K A094889 nonn
%O A094889 1,1
%A A094889 _Lekraj Beedassy_, Jun 16 2004
%E A094889 Corrected and extended by _Jason Earls_, Jun 18 2004