|
|
A205947
|
|
Carmichael numbers not congruent to 1 modulo 6.
|
|
2
|
|
|
561, 2465, 62745, 162401, 656601, 1909001, 5444489, 11921001, 19384289, 26719701, 45318561, 84350561, 151530401, 174352641, 221884001, 230996949, 275283401, 434932961, 662086041, 684106401, 689880801, 710382401
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Now known to be infinite, see Wright.
No member of this sequence is divisible by any prime of the form 6k+1, hence all prime factors for this sequence are members of A045410. (End)
|
|
LINKS
|
|
|
FORMULA
|
Wright shows that there are at least x^(K/(log log log x)^2) terms up to x, for an explicitly computable (though not computed) constant K. - Charles R Greathouse IV, Jul 20 2015
|
|
MAPLE
|
korselt:= proc(n) uses numtheory; local p;
if isprime(n) or not issqrfree(n) then return false fi;
for p in factorset(n) do
if n-1 mod (p-1) <> 0 then return false fi
od;
true
end proc:
select(korselt, [seq(seq(6*i+j, j=[3, 5]), i=1..10^5)]); # Robert Israel, Jul 20 2015
|
|
MATHEMATICA
|
Select[Range[100000], !PrimeQ[#] && IntegerQ[(#-1)/CarmichaelLambda[#]] && !Mod[#, 6]==1&]
|
|
PROG
|
(PARI) Korselt(n, f=factor(n))=for(i=1, #f[, 1], if(f[i, 2]>1||(n-1)%(f[i, 1]-1), return(0))); 1
list(lim)={
my(v=List(), p=2);
forstep(n=561, lim, [12, 6],
if(Korselt(n), listput(v, n))
);
forprime(q=3, lim,
forstep(n=p+if(p%6<5, 4, 6), q-2, 6,
if(Korselt(n), listput(v, n))
);
p=q
);
vecsort(Vec(v))
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|