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.

A378057 Composite numbers k such that A378056(k) = gcd(lcm{d+1 : d|k}, lcm{d-1 : d > 1 and d|k}) = 2.

This page as a plain text file.
%I A378057 #7 Nov 15 2024 12:56:40
%S A378057 6,481,793,949,1417,2041,2257,2509,2701,2977,3133,3589,3601,4033,4069,
%T A378057 4453,4849,5161,5317,5809,5917,5941,6697,7033,7081,7141,7501,7957,
%U A378057 7969,8593,8917,9217,9529,9577,10249,10573,10777,11041,11401,11461,11581,11773,12469,12913,12961
%N A378057 Composite numbers k such that A378056(k) = gcd(lcm{d+1 : d|k}, lcm{d-1 : d > 1 and d|k}) = 2.
%C A378057 A378056(p) = 2 for all odd primes p.
%C A378057 6 is the only even term.
%C A378057 The least term that is not a semiprime is a(114) = 29341 = 13 * 37 * 61, and the least term that has more than 3 distinct prime factors is a(4087545) = 1038565321 = 37 * 61 * 421 * 1093.
%H A378057 Amiram Eldar, <a href="/A378057/b378057.txt">Table of n, a(n) for n = 1..10000</a>
%t A378057 s[n_] := Module[{d = Divisors[n]}, GCD[LCM @@ (d + 1), LCM @@ (Rest @ d - 1)]]; s[1] = 1; Select[Range[13000], CompositeQ[#] && s[#] == 2 &]
%o A378057 (PARI) is(k) = if(isprime(k), 0, my(d = divisors(k)); gcd(lcm(apply(x->x+1, d)), lcm(apply(x -> if(x > 1, x-1, x), d))) == 2);
%Y A378057 Cf. A378056.
%K A378057 nonn
%O A378057 1,1
%A A378057 _Amiram Eldar_, Nov 15 2024