The OEIS is supported by the many generous donors to the OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A356866 Smallest Carmichael number (A002997) with n prime factors that is also a strong pseudoprime to base 2 (A001262). 0
 15841, 5310721, 440707345, 10761055201, 5478598723585, 713808066913201, 1022751992545146865, 5993318051893040401, 120459489697022624089201, 27146803388402594456683201, 14889929431153115006659489681 (list; graph; refs; listen; history; text; internal format)
 OFFSET 3,1 LINKS Table of n, a(n) for n=3..13. FORMULA a(n) >= max(A180065(n), A006931(n)). PROG (PARI) carmichael_strong_psp(A, B, k, base) = A=max(A, vecprod(primes(k+1))\2); (f(m, l, p, k, k_exp, congr, u=0, v=0) = my(list=List()); if(k==1, forprime(q=u, v, my(t=m*q); if((t-1)%l == 0 && (t-1)%(q-1) == 0, my(tv=valuation(q-1, 2)); if(tv > k_exp && Mod(base, q)^(((q-1)>>tv)< k_exp && Mod(base, q)^(((q-1)>>tv)<u, u=r); list=concat(list, f(t, L, r, k-1, k_exp, congr, u, v)))))))); list); my(res=f(1, 1, 3, k, 0, 1)); for(v=0, logint(B, 2), res=concat(res, f(1, 1, 3, k, v, -1))); vecsort(Vec(res)); a(n, base=2) = if(n < 3, return()); my(x=vecprod(primes(n+1))\2, y=2*x); while(1, my(v=carmichael_strong_psp(x, y, n, base)); if(#v >= 1, return(v[1])); x=y+1; y=2*x); CROSSREFS Cf. A001262, A002997, A006931, A063847, A180065. Sequence in context: A063847 A367231 A265669 * A184612 A277350 A101320 Adjacent sequences: A356863 A356864 A356865 * A356867 A356868 A356869 KEYWORD nonn,more AUTHOR Daniel Suteu, Oct 01 2022 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified June 23 07:16 EDT 2024. Contains 373629 sequences. (Running on oeis4.)