OFFSET
1,1
COMMENTS
a(n) up to a(7) = 23 coincides with A007665[n+1] = Tower of Hanoi with 5 pegs. It appears that a(n) includes all A007665[n] = {1, 3, 5, 7, 11, 15, 19, 23, 27, 31, 39, 47, 55, 63, 71, 79, 87, 95, 103, 111, 127, 143, 159, 175, 191, 207, 223, 239, 255, 271, 287, 303, 319, 335, 351, 383, 415, 447, 479, 511, 543, 575, 607, 639, 671, 703, 735, 767, 799, ...} except A007665[1] = 1.
EXAMPLE
There are 25 primes p < 100, p(n) = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97}.
MATHEMATICA
Select[Range[300], IntegerQ[Sum[ Prime[k]^#1, {k, 1, 25}]/25]&]
PROG
(PARI) for(n=1, 100, if(sum(k=1, 25, prime(k)^n)%25==0, print1(n, ", ")));
print; print("Alternative method not using primes:");
for(n=1, 100, m=(n-1)%6; print1((n-m)*3+(n-m+if(m>1, (m-1)*12-1, m*6-1))/3, ", ")) \\ K. Spage, Oct 23 2009
CROSSREFS
KEYWORD
nonn
AUTHOR
Alexander Adamchuk, Aug 21 2006, Sep 18 2006, Sep 21 2006
STATUS
approved