OFFSET
1,1
COMMENTS
LINKS
Robert G. Wilson v, Table of n, a(n) for n = 1..10000 (first 185 terms from Jaroslav Krizek)
EXAMPLE
Prime 433 is a term because 432 = 2^4 * 3^3 and 4/2 + 3/3 = 3 (integer).
PROG
(Magma) [n: n in [3..10^8] | IsPrime(n) and Denominator(&+[p[2]/p[1]: p in Factorization(n-1)]) eq 1]
(PARI) isA072873(n)=my(f=factor(n)); for(i=1, #f~, if(f[i, 2]%f[i, 1], return(0))); 1
lista(nn) = {forprime(p=2, nn, if (isA072873(p-1), print1(p, ", ")); ); } \\ Michel Marcus, Jan 21 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Jaroslav Krizek, Jan 11 2016
STATUS
approved