OFFSET
1,1
COMMENTS
Contains all terms of A004611 except 1. - Robert Israel, Aug 03 2016
LINKS
Clark Kimberling, Table of n, a(n) for n = 1..10000
EXAMPLE
39 = 3^1 13^1, so that the number of distinct primes 3*k+1 is 1 and the number of distinct primes 3*k + 2 is 0.
MAPLE
filter:= proc(n) local P1, P2;
P1, P2:= selectremove(t -> t mod 3 = 1, numtheory:-factorset(n));
nops(P1) > nops(P2 minus {3})
end proc:
select(filter, [$1..1000]); # Robert Israel, Aug 03 2016
MATHEMATICA
g[n_] := Map[First, FactorInteger[n]] ; z = 5000;
p1 = Select[Prime[Range[z]], Mod[#, 3] == 1 &];
p2 = Select[Prime[Range[z]], Mod[#, 3] == 2 &];
q1[n_] := Length[Intersection[g[n], p1]]
q2[n_] := Length[Intersection[g[n], p2]]
Select[Range[z], q1[#] == q2[#] &]; (* A274435 *)
Select[Range[z], q1[#] < q2[#] &]; (* A274436 *)
Select[Range[z], q1[#] > q2[#] &]; (* A274437 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Clark Kimberling, Jul 19 2016
STATUS
approved