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.

A339818 Carmichael numbers k for which the 2-adic valuation of phi(k) does not exceed the 2-adic valuation of k-1.

This page as a plain text file.
%I A339818 #12 Dec 26 2020 16:06:22
%S A339818 1729,15841,3057601,3828001,5310721,8355841,8830801,9439201,14676481,
%T A339818 15829633,17236801,40280065,78091201,83099521,84350561,92625121,
%U A339818 94536001,104852881,118901521,129762001,157731841,163954561,180115489,193708801,214852609,221884001,279377281,382304161,382536001,438359041,481239361,511338241
%N A339818 Carmichael numbers k for which the 2-adic valuation of phi(k) does not exceed the 2-adic valuation of k-1.
%H A339818 Amiram Eldar, <a href="/A339818/b339818.txt">Table of n, a(n) for n = 1..22784</a> (terms below 2^64)
%t A339818 carmichaels = Cases[Import["https://oeis.org/A002997/b002997.txt", "Table"], {_, _}][[;; , 2]]; q[n_] := IntegerExponent[EulerPhi[n], 2] <= IntegerExponent[n - 1, 2]; Select[carmichaels, q] (* _Amiram Eldar_, Dec 26 2020 *)
%o A339818 (PARI)
%o A339818 A002322(n) = lcm(znstar(n)[2]); \\ From A002322
%o A339818 isA339818(n) = ((n>1)&&issquarefree(n)&&!isprime(n)&&(valuation(eulerphi(n),2)<=valuation(n-1,2))&&(0==((n-1)%A002322(n))));
%Y A339818 Intersection of A002997 and A339817 (see comments in latter).
%Y A339818 Cf. also A339869, A339878, A339909.
%K A339818 nonn
%O A339818 1,1
%A A339818 _Antti Karttunen_, Dec 20 2020