login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A326227 Indices of nonsquarefree numerators (A001008) of harmonic numbers H(n) = Sum_{k=1..n} 1/k. 0
4, 6, 7, 10, 12, 16, 18, 22, 28, 29, 30, 36, 40, 42, 46, 52, 58, 60, 66, 70, 72, 78, 82, 88, 96, 100, 102, 106, 108, 112, 126, 130, 136, 138, 148, 150, 156, 162, 166, 172, 178, 180, 190, 192, 196, 198, 210, 222, 226, 228, 232, 238, 240, 250, 256, 262, 268, 270, 276 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

COMMENTS

It appears that the first term of A001008 having a cubic factor is A001008(848) = 11^3 * 1871 * C359.

By Wolstenholme's Theorem, p^2 divides A001008(p-1) whenever p >= 5 is prime (cf. A076637); see A308968 for illustration. Therefore, A006093 \ {1, 2} (primes - 1) is a subsequence. (Thanks to Bernard Schott.)

LINKS

Table of n, a(n) for n=1..59.

PROG

(PARI) is_A326227(n)={n>3&&vecmax(factor(A001008(n))[, 2])>1} \\ Add , 0 in factor() for much faster but possibly incorrect results [false negative].

for(n=1, oo, is_A326227(n) && print1(n", "))

CROSSREFS

Cf. A001008, A076637, A322434.

Cf. A308967 (number of prime factors), A308968 (table of factorization), A308969 (table of prime divisors), A308970 & A308971 (smallest & largest prime factor) of A001008(n).

Sequence in context: A047234 A089532 A285254 * A103401 A285091 A286485

Adjacent sequences:  A326224 A326225 A326226 * A326228 A326229 A326230

KEYWORD

nonn,more,hard

AUTHOR

M. F. Hasler, Jul 03 2019

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified November 21 04:33 EST 2019. Contains 329350 sequences. (Running on oeis4.)