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.

A181623 Sequence starting with 1 such that the sum of any two distinct elements has four distinct prime factors.

This page as a plain text file.
%I A181623 #12 Sep 16 2017 00:33:27
%S A181623 1,209,1121,2989,11381,34889,47701,62453,188785,878185,1761737,
%T A181623 3931385,5630905,7990481,32892077,204570037,253223785,1353794333,
%U A181623 2877954833
%N A181623 Sequence starting with 1 such that the sum of any two distinct elements has four distinct prime factors.
%C A181623 Choose the first number not leading to a contradiction.
%e A181623 Each of the three pairwise sums of the subset {1, 209, 1121} is the product of four distinct prime factors: {2*3*5*7,  2*3*11*17,  2*3*5*137}.
%p A181623 with(numtheory):nn:=100000:T:=array(1..nn): U:=array(1..nn): for p from 1 to
%p A181623   nn do: T[p]:=p:U[p]:=1:od:for u from 1 to 30 do: k:=1+u:for n from u+1 to nn
%p A181623   do:s:=T[n]+T[u]:s1:=nops(factorset(s)):s2:=bigomega(s):if s1=4 and s2=4 then
%p A181623   U[k]:=T[n]:k:=k+1:else fi:od:for i from 1 to nn do:T[i]:=U[i]:od:od:for j from
%p A181623   1 to 30 do:printf(`%d, `, T[j]):od:
%Y A181623 Cf. A180514, A180565, A180615, A181620, A181622.
%K A181623 nonn
%O A181623 1,2
%A A181623 _Michel Lagneau_, Jan 31 2011
%E A181623 a(9)-a(19) from _Donovan Johnson_, Feb 14 2011