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.

A107787 Golden semiprimes that are not brilliant numbers.

This page as a plain text file.
%I A107787 #11 Nov 29 2019 04:01:42
%S A107787 77,7303,10033,15229,644251,706609,836197,870071,936791,1027333,
%T A107787 1177993,1261807,1366553,1433143,1525441,1608161,62259511,62782091,
%U A107787 63462979,64441621,64973261,65427751,65599201,65837659,66043093,66834373,67966247,70635437,70882871
%N A107787 Golden semiprimes that are not brilliant numbers.
%H A107787 Amiram Eldar, <a href="/A107787/b107787.txt">Table of n, a(n) for n = 1..10000</a>
%e A107787 7303=67*109 is in the sequence because abs(67*phi-109) = 0.59172... < 1 and its prime factors have an unequal number of decimal digits.
%t A107787 f[p_] := Module[{x = GoldenRatio * p}, p1 = NextPrime[x, -1]; p2 = NextPrime[p1]; q = If[x - p1 < p2 - x, p1, p2]; If[Abs[q - x] < 1, q, 0]]; seq = {}; p = 1; Do[p = NextPrime[p]; q = f[p]; If[q > 0 && IntegerLength[p] != IntegerLength[q], AppendTo[seq, p*q]], {1000}]; seq (* _Amiram Eldar_, Nov 29 2019 *)
%Y A107787 Cf. A108540, A078972.
%K A107787 base,nonn
%O A107787 1,1
%A A107787 _Jason Earls_, Jun 14 2005
%E A107787 More terms from _Donovan Johnson_, Mar 06 2008