OFFSET
1,3
COMMENTS
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
FORMULA
EXAMPLE
a(8) = 3 because A025487(8) = 24 and 2^3 divides 24.
MATHEMATICA
max = 40000; A025487 = {1}; lpe = {}; Do[ pe = Sort[ FactorInteger[n][[All, 2]]]; If[FreeQ[lpe, pe], AppendTo[lpe, pe]; AppendTo[A025487, n]], {n, 2, max}]; a[n_] := FactorInteger[ A025487[[n]] ][[1, 2]]; a[1] = 0; Table[a[n], {n, 1, Length[A025487]}] (* Jean-François Alcover, Jun 14 2012, after Robert G. Wilson v *)
PROG
(Haskell)
a051282 = a007814 . a025487 -- Reinhard Zumkeller, Apr 06 2013
(PARI) isA025487(n)=my(k=valuation(n, 2), t); n>>=k; forprime(p=3, default(primelimit), t=valuation(n, p); if(t>k, return(0), k=t); if(k, n/=p^k, return(n==1)))
[valuation(n, 2) | n <- [1..1000], isA025487(n)]
\\ Or, for older versions:
apply(n->valuation(n, 2), select(isA025487, [1..1000])) \\ Charles R Greathouse IV, Nov 07 2014
CROSSREFS
KEYWORD
nice,nonn
AUTHOR
EXTENSIONS
More terms from Naohiro Nomoto, Mar 11 2001
STATUS
approved