OFFSET
1,1
COMMENTS
a(n)==1 (mod 4).
LINKS
Jinyuan Wang, Table of n, a(n) for n = 1..5000
PROG
(PARI) rad(n) = factorback(factorint(n)[, 1]); \\ A007947
lista(nn) = {for (n=1, nn, my(p=rad(n)); if (isprime(p) && isprime(p+6) && (p+1==rad(n+1)), print1(p, ", ")); ); } \\ Michel Marcus, Aug 22 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Giovanni Teofilatto, Aug 04 2008
EXTENSIONS
Definition clarified, more terms appended by R. J. Mathar, Aug 23 2010
A-number typo in definition corrected by R. J. Mathar, Aug 30 2010
More terms from Michel Marcus, Aug 22 2019
STATUS
approved