A085450 a(n) is the smallest m > 1 such that m divides Sum_{k=1..m} prime(k)^n.

23, 19, 25, 2951, 25, 4357, 25, 43, 281525, 269, 25, 37, 23, 295, 17, 3131, 395191, 37, 25, 19, 139, 1981, 23, 37, 25, 455, 25, 41, 124403, 61, 17, 511, 193, 535, 23, 5209, 1951, 19, 25, 301, 891, 9805, 25, 527, 23, 83, 17, 37, 131, 43, 25, 193, 53, 37, 25, 19

Offset

1,1

Comments

By definition a(1) is A045345(2).

This sequence has a very interesting behavior. If Mod(n, 2)(Mod(n, 20)-1)(Mod(n, 20)-9)(Mod(n, 20)-13)(Mod(n, 20)-17)!=0, a(n)=17, 23 or 25; in other cases a(n) may be too large. If Mod[n, 16] = 15, a(n) = 17. For example, a(n) = 17 for n = 15, 31, 47, 63, 79, 95, 111, 127, 143, 159, 175, 191, ...; also, a(n) = 23 for n = 1, 13, 23, 35, 45, 57, 67, 89, 101, 123, 133, 145, 155, 167, 177, 189, 199, ...; a(n) = 25 for n = 3, 5, 7, 11, 19, 25, 27, 39, 43, 51, 55, 59, 65, 71, 75, ..., . For a(n) = 19 for n = 2, 20, 38, 56, 74, 92, 110, 128, 146, 164, 182, 200, 218, ..., == 2 (mod 18).

From Alexander Adamchuk, Jul 20 2008: (Start)

Conjectures:

a(n) exists for all n; a(n) >= 17.

a(325)-a(575) = {25,19,25,5851,1843,61,23,821,89,301,17,37,131,455,25,1607,297,37,23,19,25,

325,25,37,353,47,17,1663,23,691,25,691,509269,155,25,269,105893,19,25,3971,

23,213215,17,26021,327,79,25,37,151,83,23,161,101,37,25,19,327,265,17,37,25,

43,23,41,169,61,25,113,21761,6289,25,47,23,19,17,4073,1137,565,25,527,25,

325,25,37,23,455,25,431,13195,37,17,19,53,155,23,37,89,455,25,18839,25,6221,

25,41,18597,229,17,811,623173,19,25,193,2079,673,25,881,23,47,25,37,25,97,

17,79,131,37,25,19,23,56501,25,37,299,455,25,167,2707,446963,17,157,25,325,

25,41,53,19,25,5917,103,1051,23,607,101,155,17,37,6233,455,25,9049,23,37,25,

19,327,5359,25,37,43,455,17,9187,23,193,25,1861,7923,301,25,113,25,19,23,41,

89,61,17,43,1785,131,25,37,1417,455,23,151,53,37,25,19,25,79,17,37,23,455,

25,289,59,47,25,511,47,83,25,739,23,19,17,301,25,269,25,41,707,2735,23,37,

299,43,25,283,69723,37,17,19,1785,479,23,37,25,455,25,1867,131,61,25,31799,

23,161,17}.

a(n) is currently unknown and a(n)>10^7 for n = {324, 576, ...}. (End)

All but one of the terms up to n=1000 are known and they are less than 10^8. Currently the only unknown term for n<=1000 is a(656)>10^8. - Alexander Adamchuk, May 24 2009

a(656) > 23,491,000,000. - Robert Price, Apr 22 2014

Links

Alexander Adamchuk and Robert Price, Table of n, a(n) for n = 1..655 (first 323 terms from Alexander Adamchuk)

Formula

For[m=2, Mod[Sum[Prime[k]^n, {k, m}], m]!=0, m++ ]; m

Example

a(3) = 25 because 2^3+3^3+5^3...+prime(25)^3 == 0 (mod 25) and for 1 < n < 25 2^3+3^3+...+prime(n)^3 is not congruent to zero (mod n).

Mathematica

a[n_] := Block[{m = 2, s = 2^n}, While[s = s + Prime[m]^n; Mod[s, m] != 0, m++ ]; m]; Table[ a[n], {n, 1, 56}]
a[n_] := Block[{m = 2, s = 2^n}, While[s = s + Prime[m]^n; Mod[s, m] != 0&& m<10000000, m++ ]; m]; Table[ a[n], {n, 1, 1000}] (* Alexander Adamchuk, Jul 20 2008 *)

Prog

(PARI) a(n)=my(s=2^n, m=1); forprime(p=3, , if((s+=p^n)%m++==0, return(m))) \\ Charles R Greathouse IV, Feb 06 2015

Crossrefs

Cf. A045345, A111441, A122140, A125907, A122142, A125825, A125826, A125828, A131263, A131264, A125827, A131272, A131273, A131274, A131275, A131276, A131277, A131278, A131279.

Sequence in context: A204598 A160436 A105818 * A077146 A077576 A004512

Adjacent sequences: A085447 A085448 A085449 * A085451 A085452 A085453

Keyword

nonn

Author

Farideh Firoozbakht, Aug 14 2003

Extensions

Edited and extended by Robert G. Wilson v, Aug 14 2003

More terms: a(324) = 18642551, a(576) = 12824827. Alexander Adamchuk, May 24 2009

Status

approved