OFFSET
1,1
COMMENTS
Taking m=1 in the definition we get the primes 2, 3, 5.
If n is in A226960, then n is a product of terms of this sequence.
If k is only allowed to be 0 or 1, we get 2, 3, 7, 43 and no more. - Jianing Song, Feb 21 2021
Also prime factors of terms in A341858. It is conjectured that this sequence is infinite. - Jianing Song, Feb 22 2021
LINKS
Ray Chandler, Table of n, a(n) for n = 1..1000
MATHEMATICA
fa = FactorInteger; free[n_] := n == Product[fa[n][[i, 1]], {i, Length[fa[n]]}] ; Os[b_, 1] = True; Os[b_, b_] = True; Os[b_, n_] := Os[b, n] = PrimeQ[n] && free[(n - 1)/b^IntegerExponent[n - 1, b]] &&IntegerExponent[n - 1, b] < 3 && Union@Table[Os[b, fa[n - 1][[i, 1]]], {i, Length[fa[n - 1]]}] == {True}; G[b_] := Select[Prime[Range[1000]], Os[b, #] &]; G[2]
PROG
(PARI) is(n)=if(!isprime(n), return(0)); if(n<13, return(1)); my(k=valuation(n-1, 2), m=n>>k, f); if(k>2, return(0)); f=factor(m); if(lcm(f[, 2])>1, return(0)); for(i=1, #f~, if(!is(f[i, 1]), return(0))); 1 \\ Charles R Greathouse IV, Oct 28 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
José María Grau Ribas, Oct 05 2013
EXTENSIONS
Revised definition from Charles R Greathouse IV, Nov 13 2013
Terms corrected by José María Grau Ribas, Nov 14 2013
STATUS
approved