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.

A265285 Carmichael numbers (A002997) k such that k-1 is a square.

This page as a plain text file.
%I A265285 #54 Feb 16 2025 08:33:27
%S A265285 46657,2433601,67371265,351596817937,422240040001,18677955240001,
%T A265285 458631349862401,286245437364810001,20717489165917230086401
%N A265285 Carmichael numbers (A002997) k such that k-1 is a square.
%C A265285 This sequence contains all Carmichael numbers n such that for all primes p dividing n, p-1 divides n-1 and furthermore, n-1 is a square.
%C A265285 Numbers sqrt(a(n)-1) form a subsequence of A135590. - _Max Alekseyev_, Apr 25 2024
%H A265285 G. Tarry, I. Franel, A. Korselt, and G. Vacca, <a href="https://oeis.org/wiki/File:Probl%C3%A8me_chinois.pdf">Problème chinois</a>, L'intermédiaire des mathématiciens 6 (1899), pp. 142-144.
%H A265285 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/CarmichaelNumber.html">Carmichael Number</a>.
%H A265285 <a href="/index/Ca#Carmichael">Index entries for sequences related to Carmichael numbers</a>.
%e A265285 46657 is a term because 46657 - 1 = 46656 = 216^2.
%e A265285 2433601 is a term because 2433601 - 1 = 2433600 = 1560^2.
%p A265285 isA002997:= proc(n) local F,p;
%p A265285          if n::even or isprime(n)  then return false fi;
%p A265285          F:= ifactors(n)[2];
%p A265285          if max(seq(f[2],f=F)) > 1 then return false fi;
%p A265285          andmap(f -> (n-1) mod (f[1]-1) = 0,  F)
%p A265285 end proc:
%p A265285 select(isA002997, [seq(4*i^2+1,i=1..10^6)]); # _Robert Israel_, Dec 08 2015
%o A265285 (PARI) is_c(n) = { my(f); bittest(n, 0) && !for(i=1, #f=factor(n)~, (f[2, i]==1 && n%(f[1, i]-1)==1)||return) && #f>1 }
%o A265285 for(n=1, 1e10, if(is_c(n) && issquare(n-1), print1(n, ", ")))
%o A265285 (PARI) lista(kmax) = {my(m); for(k = 2, kmax, m = k^2 + 1; if(!isprime(m), f = factor(k); for(i = 1, #f~, f[i, 2] *= 2); fordiv(f, d, if(!(m % (d+1)) && isprime(d+1), m /= (d+1))); if(m == 1, print1(k^2 + 1, ", ")))); } \\ _Amiram Eldar_, May 02 2024
%Y A265285 Subsequence of A265237 and of A265328.
%Y A265285 Cf. A002997, A135590, A265237, A303791.
%K A265285 nonn,hard,more
%O A265285 1,1
%A A265285 _Altug Alkan_, Dec 06 2015
%E A265285 a(4)-a(5), using A002997 b-file, from _Michel Marcus_, Dec 07 2015
%E A265285 a(6) and a(7) from _Robert Israel_, Dec 08 2015
%E A265285 a(8) from _Max Alekseyev_, Apr 30 2018
%E A265285 a(9) from _Daniel Suteu_ confirmed by _Max Alekseyev_, Apr 25 2024