cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A260524 Pseudoprimes to bases 2, 3, 5 and 7 that are congruent to 5 (modulo 6) but are not Carmichael numbers (A002997).

This page as a plain text file.
%I A260524 #16 Sep 07 2015 14:43:05
%S A260524 468950021,493108481,659846021,5936122901,8144063621,11408333333,
%T A260524 12601267541,14252656133,18074903681,27223783841,30633711701,
%U A260524 31093792133,31797754721,61426533761,65085388961,86610942881,91945013333,92380393121,102538073177
%N A260524 Pseudoprimes to bases 2, 3, 5 and 7 that are congruent to 5 (modulo 6) but are not Carmichael numbers (A002997).
%H A260524 Charles R Greathouse IV, <a href="/A260524/b260524.txt">Table of n, a(n) for n = 1..10000</a>
%t A260524 fQ[n_] := !PrimeQ[n] && PowerMod[2, n - 1, n] == 1 &&
%t A260524 > PowerMod[3, n - 1, n] == 1 && PowerMod[5, n - 1, n] == 1 && PowerMod[7, n - 1, n] == 1 && Mod[n, CarmichaelLambda[n]] != 1; k = 1; lst = {}; While[k < 25000000001, If[ fQ@ k, AppendTo[lst, k]; Print@ k]; k += 6]; lst
%o A260524 (PARI) Korselt(n)=my(f=factor(n)); for(i=1, #f[, 1], if(f[i, 2]>1||(n-1)%(f[i, 1]-1), return(0))); 1
%o A260524 is(n)=n%6==5 && Mod(2,n)^n==2 && Mod(3,n)^n==3 && Mod(5,n)^(n-1)==1 && Mod(7,n)^(n-1)==1 && !isprime(n) && !Korselt(n) \\ _Charles R Greathouse IV_, Jul 29 2015
%o A260524 (Perl) use ntheory ":all"; foroddcomposites { say if $_%6 == 5 && is_pseudoprime($_,2,3,5,7) && $_ % carmichael_lambda($_) != 1; } 1e9; # _Dana Jacobsen_, Sep 07 2015
%Y A260524 Cf. A153581.
%K A260524 nonn
%O A260524 1,1
%A A260524 _Ray Chandler_, _Artur Jasinski_ and _Robert G. Wilson v_, Jul 28 2015
%E A260524 a(9)-a(19) from _Charles R Greathouse IV_, Jul 29 2015