OFFSET
1,1
COMMENTS
Differs from A095990 starting with the number 2250.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
Canonical factorization of 2250 = 2^1 * 3^2 * 5^3 or 2,3,5 raised to powers 1,2,3, an arithmetic progression.
MAPLE
filter:= proc(n) local F, m, i;
F:= map(t->t[2], sort(ifactors(n)[2], (a, b) -> a[1]<b[1]));
nops(F) > 1 and F = [$1..nops(F)]
end proc:
select(filter, [$2..10^4]); # Robert Israel, Nov 23 2016
MATHEMATICA
filterQ[n_] := CompositeQ[n] && With[{f = FactorInteger[n][[All, 2]]}, f == Range[Length[f]]];
Select[Range[10000], filterQ] (* Jean-François Alcover, Aug 28 2020 *)
PROG
(PARI) omnipprog(n, m) = { local(a, x, j, nf, fl=0); for(x=1, n, a=factor(x); nf=omega(x); for(j=1, nf, if(a[j, 2]==j, fl=1, fl=0; break); ); if(fl&nf>1, print1(x", ")) ) }
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Cino Hilliard, Feb 27 2005
EXTENSIONS
Name changed by Robert Israel, Nov 23 2016
STATUS
approved