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.
%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