Alois P. Heinz, <a href="/A115103/b115103_1.txt">Table of n, a(n) for n = 1..10000</a> (first 1000 terms from Vincenzo Librandi)
Alois P. Heinz, <a href="/A115103/b115103_1.txt">Table of n, a(n) for n = 1..10000</a> (first 1000 terms from Vincenzo Librandi)
editing
approved
Vincenzo Librandi, Alois P. Heinz, <a href="/A115103/b115103_1.txt">Table of n, a(n) for n = 1..10000</a> (first 1000</a> terms from Vincenzo Librandi)
approved
editing
editing
approved
# second Maple program:
q:= p-> isprime(p) and (f-> f(p+1)=f(p-1))(numtheory[bigomega]):
select(q, [$1..3000])[]; # Alois P. Heinz, May 08 2022
approved
editing
editing
approved
isA115103 := proc(n)
if not type(n, prime) then
return false;
end if;
if numtheory[bigomega](n-1) <> numtheory[bigomega](n+1) then
false;
else
true ;
end if ;
end proc:
for n from 2 to 3000 do
if isA115103(n) then
printf("%d, ", n) ;
end if;
end do: # R. J. Mathar, Feb 13 2019
approved
editing
proposed
approved
editing
proposed
Cf. A067386 (without multiplicity).
approved
editing