A308800 Primes p such that A001177(p) = (p-1)/7.
2731, 3739, 4831, 6091, 11159, 13679, 14771, 16871, 19559, 20399, 24179, 26111, 29191, 31039, 33811, 34511, 34679, 35911, 40111, 41651, 49211, 55259, 56099, 60859, 62819, 69539, 71191, 71359, 71471, 73291, 74831, 85751, 87991, 96979, 97231, 97931, 104959, 108179
Offset: 1
Keywords
Links
- Amiram Eldar, Table of n, a(n) for n = 1..10000
Crossrefs
Programs
-
Mathematica
pn[n_] := For[k = 1, True, k++, If[Mod[Fibonacci[k], n] == 0, Return[k]]]; Reap[For[p = 2, p < 50000, p = NextPrime[p], If[Mod[p, 7] == 1, If[pn[p] == (p - 1)/7, Print[p]; Sow[p]]]]][[2, 1]] (* Jean-François Alcover, Jul 05 2019 *)
-
PARI
Entry_for_decomposing_prime(p) = my(k=1, M=[k, 1; 1, 0]); if(isprime(p)&&kronecker(k^2+4,p)==1, my(v=divisors(p-1)); for(d=1, #v, if((Mod(M,p)^v[d])[2,1]==0, return(v[d])))) forprime(p=2, 109000, if(Entry_for_decomposing_prime(p)==(p-1)/7, print1(p, ", ")))
Comments