login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A121707 Numbers n > 1 such that n^3 divides Sum_{k=1..n-1} k^n = A121706(n). 24
35, 55, 77, 95, 115, 119, 143, 155, 161, 187, 203, 209, 215, 221, 235, 247, 253, 275, 287, 295, 299, 319, 323, 329, 335, 355, 371, 377, 391, 395, 403, 407, 413, 415, 437, 455, 473, 475, 493, 497, 515, 517, 527, 533, 535, 539, 551, 559, 575, 581, 583, 589, 611 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

COMMENTS

All terms belong to A038509 (Composite numbers with smallest prime factor >= 5). Many but not all terms belong to A060976 (Odd nonprimes, c, which divide Bernoulli(2*c)).

Many terms are semiprimes:

- the non-semiprimes are {275, 455, 475, 539, 575, 715, 775, 875, 935, ...}: see A321487;

- semiprime terms that are multiples of 5 have indices {7, 11, 19, 23, 31, 43, 47, 59, 67, 71, 79, 83, ...} = A002145 (Primes of form 4*k + 3, except 3, or k > 0; or Primes which are also Gaussian primes);

- semiprime terms that are multiples of 7 have indices {5, 11, 17, 23, 29, 41, 47, 53, 59, 71, 83, 89, ...} = A003627 (Primes of form 3*k - 1, except 2, or k > 1);

- semiprime terms that are multiples of 11 have indices {5, 7, 13, 17, 19, 23, 37, 43, 47, 53, 59, 67, 73, 79, 83, ...} = Primes of the form 4*k + 1 and 4*k - 1. [Edited by Michel Marcus, Jul 21 2018, M. F. Hasler, Nov 09 2018]

Conjecture: odd numbers n > 1 such that n divides Sum_{k=1..n-1} k^(n-1). - Thomas Ordowski and Robert Israel, Oct 09 2015. Professor Andrzej Schinzel (in a letter to me, dated Dec 29 2015) proved this conjecture. - Thomas Ordowski, Jul 20 2018

Note that n^2 divides Sum_{k=1..n-1} k^n for every odd number n > 1. - Thomas Ordowski, Oct 30 2015

Conjecture: these are "anti-Carmichael numbers" defined; n > 1 such that p - 1 does not divide n - 1 for every prime p dividing n. Equivalently, odd numbers n > 1 such that n is coprime to A027642(n-1). A number n > 1 is an "anti-Carmichael" if and only if gcd(n, b^n - b) = 1 for some integer b. - Thomas Ordowski, Jul 20 2018

It seems that these numbers are all composite terms of A317358. - Thomas Ordowski, Jul 30 2018

a(62) = 697 is the first term not in A267999: see A306097 for all these terms. - M. F. Hasler, Nov 09 2018

If the conjecture from Thomas Ordowski is true, then no term is a multiple of 2 or 3. - Jianing Song, Jan 27 2019

Conjecture: an odd number n > 1 is a term iff gcd(n, A027642(n-1)) = 1. - Thomas Ordowski, Jul 19 2019

LINKS

Giovanni Resta, Table of n, a(n) for n = 1..10000 (first 1371 terms from Robert Israel)

Don Reble, Comments on A121707

MAPLE

filter:= n -> add(k &^ n mod n^3, k=1..n-1) mod n^3 = 0:

select(filter, [$2..1000]); # Robert Israel, Oct 08 2015

MATHEMATICA

fQ[n_] := Mod[Sum[PowerMod[k, n, n^3], {k, n - 1}], n^3] == 0; Select[

Range[2, 611], fQ] (* Robert G. Wilson v, Apr 04 2011 and slightly modified Aug 02 2018 *)

PROG

(PARI) is(n)=my(n3=n^3); sum(k=1, n-1, Mod(k, n3)^n)==0 \\ Charles R Greathouse IV, May 09 2013

(PARI) for(n=2, 1000, if(sum(k=1, n-1, k^n) % n^3 == 0, print1(n", "))) \\ Altug Alkan, Oct 15 2015

(Sage) # after Andrzej Schinzel

def isA121707(n):

    if n == 1 or is_even(n): return False

    return n.divides(sum(k^(n-1) for k in (1..n-1)))

[n for n in (1..611) if isA121707(n)] # Peter Luschny, Jul 18 2019

CROSSREFS

Cf. A000312, A002145, A002997, A027642, A031971, A038509, A060976, A121706, A267999 (probably a subsequence).

Cf. A306097 for terms of this sequence A121707 not in sequence A267999, A321487 for terms which are not semiprimes.

Cf. A191677 (n divides Sum_{k<n} k^(n-1)).

Sequence in context: A043915 A318572 A171082 * A267999 A319386 A157352

Adjacent sequences:  A121704 A121705 A121706 * A121708 A121709 A121710

KEYWORD

nonn

AUTHOR

Alexander Adamchuk, Aug 16 2006

EXTENSIONS

Sequence corrected by Robert G. Wilson v, Apr 04 2011

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 March 28 17:51 EDT 2020. Contains 333103 sequences. (Running on oeis4.)