cp's OEIS Frontend

a(n) is the largest number m such that the arithmetic mean of {1, p, p^2, ..., p^k} is an integer for all k in 1..m.

Apparently, all the terms are of the form prime(k)-2 (A040976). Conjecture: The asymptotic density of the occurrences of prime(k)-2 is (1/s(k-1)-1/s(k)), where s(k) = A005867(k) = phi(prime(k)#), and prime(k)# is the k-th primorial number (A002110).

The sums of the first 10^k terms, for k = 1, 2, ..., are 15, 221, 2291, 23287, 233641, 2337007, 23379901, 233814475, 2338211029, 23382168187, ... . If the mentioned above conjecture is correct, then the asymptotic mean of this sequence is Sum_{k>=1} (prime(k)-2)*(1/s(k-1)-1/s(k)) = 2.33821872365981424748... .

Apparently, the indices of records after n = 1 occur at A000720(A073917(n)) (verified for the first 12 terms of A073917) with record values a(A000720(A073917(n))) = prime(n+1) - 2 (verified for the first 150 terms of A073917).

a(n) is the least number m whose number of divisors is A000005(m) = n such that the arithmetic mean of the first k divisors of m is an integer for all k in 1..n.

a(17) = 4084081^16 = 5.991...*10^105 is too large to include in the data section.

a(n) exists for all n >= 1. For n > 1, consider a prime p of the form m*lcm(1,2,...n-1) + 1, with m >= 1. Such a prime exists by Dirichlet's theorem on arithmetic progressions. Then, p^(n-1) has n divisors, and p^k == 1 (mod lcm(1..n-1)) for k = 0..(n-1). Therefore, Sum_{k=0..n-1} p^k == k (mod lcm(1,2,...n-1)), or equivalently, Sum_{k=0..n-1} p^k is divisible by k for k = 0..(n-1). Thus, p^(n-1) is in A359260.