OFFSET
1,1
LINKS
Muniru A Asiru, Table of n, a(n) for n = 1..71
MATHEMATICA
a := Select[Range[1000000], IntegerQ[HarmonicMean[Divisors[#]]] &]; Table[Part[a, n + 1] - Part[a, n], {n, 1, 30}] (* G. C. Greubel, Aug 29 2016 *)
PROG
(GAP) H:=Concatenation([1], Filtered([2, 4..3*10^6], n->IsInt(n*Tau(n)/Sigma(n))));;
a:=List([1..Length(H)-1], i->H[i+1]-H[i]);; Print(a); # Muniru A Asiru, Dec 23 2018
(Python)
from numpy import diff
from sympy import divisor_sigma as sigma
def ore(n): return (n*sigma(n, 0))%sigma(n, 1) == 0
print(diff(list(filter(ore, range(1, 9000))))) # Michael S. Branicky, May 22 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Omar E. Pol, Jan 20 2009
EXTENSIONS
More terms from R. J. Mathar, Feb 19 2009
STATUS
approved