OFFSET
1,1
COMMENTS
For a guide to related sequences, see A204892.
a(n) <= n+1. If n is divisible by p^2 then a(n) <= p+n/p. - Robert Israel, May 14 2021
LINKS
Robert Israel, Table of n, a(n) for n = 1..3000
EXAMPLE
1 divides 2^5-1^5 -> k=2, j=1
2 divides 3^5-1^5 -> k=3, j=1
3 divides 4^5-1^5 -> k=4, j=1
4 divides 4^5-2^5 -> k=4, j=2
5 divides 6^5-1^5 -> k=6, j=1
6 divides 7^5-1^5 -> k=7, j=1
MAPLE
N:= 100: # for a(1)..a(N)
V:= Vector(N):
count:= 0:
for k from 1 while count < N do
for j from 1 to k-1 while count < N do
Q:= select(t -> t <= N and V[t] = 0, numtheory:-divisors(k^5-j^5));
if Q <> {} then
newcount:= nops(Q);
count:= count + newcount;
V[convert(Q, list)]:= k;
fi
od od:
convert(V, list); # Robert Israel, May 14 2021
MATHEMATICA
s = Table[n^4, {n, 1, 120}] ;
lk = Table[
NestWhile[# + 1 &, 1,
Min[Table[Mod[s[[#]] - s[[j]], z], {j, 1, # - 1}]] =!= 0 &], {z, 1,
Length[s]}]
Table[NestWhile[# + 1 &, 1,
Mod[s[[lk[[j]]]] - s[[#]], j] =!= 0 &], {j, 1, Length[lk]}]
(* Peter J. C. Moses, Jan 27 2012 *)
Array[(k=1; While[FreeQ[Mod[Table[k^5-j^5, {j, k-1}], #], 0], k++]; k)&, 100] (* Giorgos Kalogeropoulos, May 14 2021 *)
CROSSREFS
KEYWORD
nonn,look
AUTHOR
Clark Kimberling, Feb 01 2012
STATUS
approved