|
%I #17 Oct 01 2023 01:57:35
%S 15,21,30,33,35,39,42,51,55,57,60,65,66,69,70,77,78,84,85,87,91,93,95,
%T 102,105,110,111,114,115,119,120,123,129,130,132,133,135,138,140,141,
%U 143,145,154,155,156,159,161,165,168,170,174,177,182,183,185,186
%N Numbers having more than one odd prime factor to an odd power in their prime factorization.
%C The sequence is of asymptotic density one, a(n) ~ n.
%H N. J. A. Sloane, <a href="/A285800/b285800.txt">Table of n, a(n) for n = 1..19999</a>
%e 15 = 3*5, 21 = 3*7, 30 = 2*15, 33 = 3*11 are the smallest positive integers having at least two prime factors to an odd power in their factorization.
%e a(10) = 57, a(100) = 287, a(10^3) = 1950, a(10^4) = 15701, a(10^5) = 138540, a(10^6) = 1284998.
%p s800:=[]; s801:=[1];
%p for n from 2 to 1000 do
%p c:=0;
%p t2:=ifactors(n)[2];
%p for t3 in t2 do if t3[1]>2 and (t3[2] mod 2 = 1) then c:=c+1; fi; od:
%p if c <= 1 then s801:=[op(s801),n]; else s800:=[op(s800),n]; fi;
%p od:
%p s800; # A285800
%p s801; # A285801 - _N. J. A. Sloane_, Sep 30 2017
%o (PARI) is(n)=1<#select(t->bittest(t,0),factor(n>>valuation(n,2))[,2])
%Y Subsequence of A248150.
%Y Complement of A285801.
%K nonn,easy
%O 1,1
%A _M. F. Hasler_, Apr 26 2017
|
