A226025 Odd composite numbers that are not squares of primes.
15, 21, 27, 33, 35, 39, 45, 51, 55, 57, 63, 65, 69, 75, 77, 81, 85, 87, 91, 93, 95, 99, 105, 111, 115, 117, 119, 123, 125, 129, 133, 135, 141, 143, 145, 147, 153, 155, 159, 161, 165, 171, 175, 177, 183, 185, 187, 189, 195, 201, 203, 205, 207, 209, 213, 215, 217
Offset: 1
Keywords
Links
- Arkadiusz Wesolowski, Table of n, a(n) for n = 1..10000
Programs
-
Haskell
a226025 n = a226025_list !! (n-1) a226025_list = filter ((/= 2) . a100995) a071904_list -- Reinhard Zumkeller, Jun 15 2013
-
Magma
[n: n in [3..217 by 2] | not IsPrime(n) and not IsSquare(n) or IsSquare(n) and not IsPrime(Floor(n^(1/2)))];
-
Maple
select(n -> not(isprime(n)) and (not(issqr(n)) or not(isprime(sqrt(n)))), [seq(2*i+1,i=1..1000)]); # Robert Israel, Sep 08 2014
-
Mathematica
Select[Range[3, 217, 2], ! PrimeQ[#] && ! PrimeQ@Sqrt[#] &] r = Prime@Range[2, 6]^2; Complement[Select[Range[3, Last[r] - 2, 2], ! PrimeQ[#] &], Most[r]] Select[Range[3,251,2],NoneTrue[{#,Sqrt[#]},PrimeQ]&] (* Harvey P. Dale, Sep 06 2021 *)
-
PARI
is_A226025(n)={bittest(n,0)&&!isprime(n,0)&&!(issquare(n)&&isprime(sqrtint(n)))&&n>1} \\ - M. F. Hasler, Oct 20 2013
Comments