%I #24 Feb 11 2021 22:53:54
%S 1,5,79,445,16055,27067045,1208198549,1208198609,1208198615,
%T 2636913002983,885992692751479,885992692751483,885992692751833,
%U 43525513764815035,43525513764815057,43525513764815173,43525513764816367
%N Numbers k such that prime(k) == 6 (mod k).
%C a(18) > 2.8*10^20. - _Giovanni Resta_, Feb 23 2020
%t NextPrim[n_] := Block[{k = n + 1}, While[ !PrimeQ[k], k++ ]; k]; p = 1; Do[ If[ Mod[p = NextPrim[p], n] == 6, Print[n]], {n, 1, 10^9}] (* _Robert G. Wilson v_, Feb 18 2004 *)
%o (Sage)
%o def A023148(max) :
%o terms = []
%o p = 2
%o for n in range(1, max+1) :
%o if (p - 6) % n == 0 : terms.append(n)
%o p = next_prime(p)
%o return terms
%o # _Eric M. Schmidt_, Feb 05 2013
%Y Cf. A092048, A023143, A023144, A023145, A023146, A023147, A023149, A023150, A023151, A023152.
%K nonn,more
%O 1,2
%A _David W. Wilson_
%E Extended by _Robert G. Wilson v_, Feb 18 2004
%E a(7)-a(9) from _Robert G. Wilson v_, Feb 22 2006
%E First two terms inserted by _Eric M. Schmidt_, Feb 05 2013
%E a(10)-a(17) from _Giovanni Resta_, Feb 23 2020
|