login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A072777 Powers of squarefree numbers that are not squarefree. 11

%I #50 Jan 04 2024 14:13:34

%S 4,8,9,16,25,27,32,36,49,64,81,100,121,125,128,169,196,216,225,243,

%T 256,289,343,361,441,484,512,529,625,676,729,841,900,961,1000,1024,

%U 1089,1156,1225,1296,1331,1369,1444,1521,1681,1764,1849

%N Powers of squarefree numbers that are not squarefree.

%C For all n exists k: a(n) = A072774(k) and A072776(k) > 1.

%H Stanislav Sykora and Reinhard Zumkeller, <a href="/A072777/b072777.txt">Table of n, a(n) for n = 1..20000</a> (first 10000 terms from Reinhard Zumkeller)

%F Sum_{n>=1} 1/a(n) = Sum_{n>=2} mu(n)^2/(n*(n-1)) = Sum_{n>=2} (zeta(n)/zeta(2*n) - 1) = 0.8486338679... (A368250). - _Amiram Eldar_, Jul 22 2020

%e The number 144 = 12^2 is not a member because 12 is not squarefree.

%e 64 = 2^6 and 49 = 7^2 are members because, though not squarefree, they are powers of the squarefree numbers 2 and 7, respectively. Note that 64 is included even though it is also a square of a nonsquarefree number. - _Stanislav Sykora_, Jul 11 2014

%t Select[Range[2000], Length[u = Union[FactorInteger[#][[All, 2]]]] == 1 && u[[1]] > 1 &] (* _Jean-François Alcover_, Mar 27 2013 *)

%o (Haskell)

%o import Data.Map (singleton, findMin, deleteMin, insert)

%o a072777 n = a072777_list !! (n-1)

%o a072777_list = f 9 (drop 2 a005117_list) (singleton 4 (2, 2)) where

%o f vv vs'@(v:ws@(w:_)) m

%o | xx < vv = xx : f vv vs' (insert (bx*xx) (bx, ex+1) $ deleteMin m)

%o | xx > vv = vv : f (w*w) ws (insert (v^3) (v, 3) m)

%o where (xx, (bx, ex)) = findMin m

%o -- _Reinhard Zumkeller_, Apr 06 2014

%o (PARI) BelongsToA(n) = {my(f, k, e); if(n == 1, return(0));

%o f = factor(n); e = f[1, 2]; if(e == 1, return(0));

%o for(k = 2, #f[, 2], if(f[k, 2] != e, return(0))); return(1);}

%o Ntest(nmax, test) = {my(k = 1, n = 0, v); v = vector(nmax); while(1, n++; if(test(n), v[k] = n; k++; if(k > nmax, break)); ); return(v); }

%o a = Ntest(20000, BelongsToA) \\ Note: not very efficient. - _Stanislav Sykora_, Jul 11 2014

%o (PARI) is(n)=ispower(n,,&n) && issquarefree(n) \\ _Charles R Greathouse IV_, Oct 16 2015

%Y Cf. A013929, A368250.

%Y Cf. A005117, subsequence of A001597 and A072774.

%K nonn,nice

%O 1,1

%A _Reinhard Zumkeller_, Jul 10 2002

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 00:26 EDT 2024. Contains 371264 sequences. (Running on oeis4.)