%I #9 Mar 30 2012 18:35:54
%S 1,2,3,6,58,124,254,309,519,1029,1179,1569,1986,3795,10008,31133,
%T 39260,76772,126798,190293,613553,873413,1324947,16893137,23186977,
%U 65348522,91513433,168375480,836588442,844570409
%N Numbers starting with 1 such that the sum of any two distinct elements has an odd number of distinct prime factors.
%C Numbers starting with 4 :
%C 4, 5, 12, 25, 85, 126, 145, 186, 252, 1146, ...
%C Numbers starting with 5 :
%C 5, 6, 11, 26, 55, 424, 444, 589, 722, 1573, ...
%C Numbers starting with 7 :
%C 7, 9, 10, 20, 22, 118, 350, 1012, 1433, 2043, ...
%e The set {6, 58, 124} gives the number of distinct prime factors {1, 3, 3}.
%p with(numtheory):nn:=1000000:T:=array(1..nn): U:=array(1..nn): for p from 1
%p to 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
%p nn do:s:=T[n]+T[u]:s1:=nops(factorset(s)):z:=irem(s1, 2):if z=1 then U[k]:=T[n]:k:=k+1:else
%p fi:od:for i from 1 to nn do:T[i]:=U[i]:od:od:for j from 1 to 30 do:printf(`%d, `, T[j]):od:
%t t={1}; k=1; Do[k++; While[! And @@ OddQ[Length /@ FactorInteger[t+k]], k++]; AppendTo[t, k], {10}]; t
%Y Cf. A180514.
%K nonn
%O 1,2
%A _Michel Lagneau_, Jan 21 2011
%E a(23)-a(30) from _Donovan Johnson_, Jan 25 2011