OFFSET
1,2
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
expQ[e1_, e2_] := Module[{m = Length[e2], ans = 1}, Do[If[BitAnd[e1[[i]], e2[[i]]] < e2[[i]], ans = 0; Break[]], {i, 1, m}]; ans];
e[n_] := e[n] = FactorInteger[n!][[;; , 2]]; q[n_] := Sum[expQ[e[n], e[m]], {m, 2, n}] <= 1; Select[Range[240], q]
PROG
(PARI) isexp(e1, e2) = {my(m = #e2, ans = 1); for(i=1, m, if(bitand(e1[i], e2[i]) < e2[i], ans = 0; break)); ans; }
e(n) = factor(n!)[, 2];
is(n) = sum(m = 2, n, isexp(e(n), e(m))) <= 1;
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Aug 22 2024
STATUS
approved