|
%I #20 Dec 09 2022 07:07:25
%S 6,10,12,14,15,18,20,21,22,24,26,28,30,33,34,35,38,39,40,42,44,45,46,
%T 48,50,51,52,54,55,56,57,58,60,62,63,65,66,68,69,70,72,74,75,76,77,78,
%U 80,82,84,85,86,87,88,90,91,92,93,94,95,96,98,99,102,104,105,106
%N Composite numbers that are not perfect powers.
%H Amiram Eldar, <a href="/A106543/b106543.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) = n + O(n/log n). - _Charles R Greathouse IV_, Oct 03 2011
%t perfPQ[n_]:=GCD@@FactorInteger[n][[All,2]]>1; Select[Range[110], CompositeQ[ #] && !perfPQ[#]&] (* _Harvey P. Dale_, Oct 10 2017 *)
%o (PARI) lista(nn)=forcomposite(i=1, nn, if (! ispower(i), print1(i, ", "));); \\ _Michel Marcus_, Jun 27 2013
%o (PARI) is(n)=!isprime(n) && !ispower(n) && n>1 \\ _Charles R Greathouse IV_, Oct 19 2015
%o (CASE)
%o def A106543_list(n) : return [k for k in (2..n) if not k.is_prime() and not k.is_perfect_power()]
%o A106543_list(106) # _Terry D. Grant_, Jul 17 2016
%Y Intersection of A002808 and A007916.
%Y Cf. A105642, A089229.
%K easy,nonn
%O 1,1
%A _Alexandre Wajnberg_, May 08 2005
|
