OFFSET
1,2
COMMENTS
LINKS
Sean A. Irvine, Table of n, a(n) for n = 1..59
EXAMPLE
Initial solutions:
n a(n) palindrome
--------------------------------------------
1 1 [1]
2 9 [5, 5]
3 3 [5, 1, 5]
4 8 [1, 5, 5, 1]
5 3 [5, 1, 5, 1, 5]
6 30 [5, 1, 5, 5, 1, 5]
7 3 [5, 1, 5, 1, 5, 1, 5]
8 29 [1, 5, 1, 5, 5, 1, 5, 1]
9 16 [5, 5, 1, 1, 5, 1, 1, 5, 5]
10 87 [5, 1, 5, 1, 5, 5, 1, 5, 1, 5]
11 23 [5, 5, 1, 5, 1, 5, 1, 5, 1, 5, 5]
12 86 [5, 5, 1, 5, 1, 5, 5, 1, 5, 1, 5, 5]
MAPLE
with(numtheory):
Nmax := 50: Pmax := 500000:
b := Array(1..Pmax):
for k from 1 to Pmax do b[k] := irem(ithprime(k), 6) end do:
a := Array(1..Nmax, fill=infinity):
for k from 1 to Pmax do
r := 0:
while k-r >=1 and k+r <= Pmax and b[k-r]=b[k+r] do
n := 2*r+1: start := k-r:
if n <= Nmax and start < a[n] then a[n] := start: end if:
r := r+1:
end do:
r := 0:
while k-r >=1 and k+r+1 <= Pmax and b[k-r]=b[k+r+1] do
n := 2*r+2: start := k-r:
if n <= Nmax and start < a[n] then a[n] := start: end if:
r := r+1:
end do:
end do:
seq(a[n], n=1..Nmax);
PROG
(PARI) a(n, v) = for (i=1, #v-n, my(w = vector(n, k, v[i+k-1])); if (Vecrev(w) == w, return(i)); );
lista(nn) = my(v = primes(nn) % 6); for (n=1, 50, my(x = a(n, v)); if (x, print1(x, ", "), return));
lista(10^7) \\ Michel Marcus, Mar 10 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Mar 09 2026
EXTENSIONS
More terms from Michel Marcus, Mar 10 2026
STATUS
approved
