%I #54 Sep 08 2024 10:02:16
%S 1,2,3,4,1,6,1,8,9,2,1,12,1,2,3,16,1,18,1,4,3,2,1,24,1,2,27,4,1,6,1,
%T 32,3,2,1,36,1,2,3,8,1,6,1,4,9,2,1,48,1,2,3,4,1,54,1,8,3,2,1,12,1,2,9,
%U 64,1,6,1,4,3,2,1,72,1,2,3,4,1,6,1,16,81,2,1,12,1,2,3,8,1,18,1,4,3,2,1,96
%N Largest 3-smooth divisor of n.
%C Bennett, Filaseta, & Trifonov show that if n > 8 then a(n^2 + n) < n^0.715. - _Charles R Greathouse IV_, May 21 2014
%H Charles R Greathouse IV, <a href="/A065331/b065331.txt">Table of n, a(n) for n = 1..10000</a>
%H M. A. Bennett, M. Filaseta, and O. Trifonov, <a href="http://www.math.ubc.ca/~bennett/BFTpaper0207.pdf">On the factorization of consecutive integers</a>, J. Reine Angew. Math. 629 (2009), pp. 171-200.
%H Jon Maiga, <a href="http://sequencedb.net/s/A065331">Computer-generated formulas for A065331</a>, Sequence Machine.
%F a(n) = n / A065330(n).
%F a(n) = A006519(n) * A038500(n).
%F a(n) = (2^A007814 (n)) * (3^A007949(n)).
%F Multiplicative with a(2^e)=2^e, a(3^e)=3^e, a(p^e)=1, p>3. - _Vladeta Jovovic_, Nov 05 2001
%F Dirichlet g.f.: zeta(s)*(1-2^(-s))*(1-3^(-s))/ ( (1-2^(1-s))*(1-3^(1-s)) ). - _R. J. Mathar_, Jul 04 2011
%F a(n) = gcd(n,6^n). - _Stanislav Sykora_, Feb 08 2016
%F a(A225546(n)) = A225546(A053165(n)). - _Peter Munn_, Jan 17 2020
%F Sum_{k=1..n} a(k) ~ n*(log(n)^2 + (2*gamma + 3*log(2) + 2*log(3) - 2)*log(n) + (2 + log(2)^2/6 + 3*log(2)*(log(3) - 1) - 2*log(3) + log(3)^2/6 + gamma*(3*log(2) + 2*log(3) - 2) - 2*sg1)) / (6*log(2)*log(3)), where gamma is the Euler-Mascheroni constant A001620 and sg1 is the first Stieltjes constant (see A082633). - _Vaclav Kotesovec_, Sep 19 2020
%F a(n) = A003586(A322026(n)), A322026(n) = A071521(a(n)). - _Antti Karttunen_, Sep 08 2024
%p A065331 := proc(n) n/A065330(n) ; end: # _R. J. Mathar_, Jun 24 2009
%p seq(2^padic:-ordp(n,2)*3^padic:-ordp(n,3), n=1..100); # _Robert Israel_, Feb 08 2016
%t Table[GCD[n, 6^n], {n, 100}] (* _Vincenzo Librandi_, Feb 09 2016 *)
%t a[n_] := Times @@ ({2, 3}^IntegerExponent[n, {2, 3}]); Array[a, 100] (* _Amiram Eldar_, Sep 19 2020 *)
%o (PARI) a(n)=3^valuation(n,3)<<valuation(n,2) \\ _Charles R Greathouse IV_, Aug 21 2011
%o (PARI) a(n)=gcd(n,6^n) \\ Not very efficient, but simple. _Stanislav Sykora_, Feb 08 2016
%o (PARI) a(n)=gcd(6^logint(n,2),n) \\ 'optimized' version of Sykora's script; _Charles R Greathouse IV_, Jul 23 2024
%o (Haskell)
%o a065331 = f 2 1 where
%o f p y x | r == 0 = f p (y * p) x'
%o | otherwise = if p == 2 then f 3 y x else y
%o where (x', r) = divMod x p
%o -- _Reinhard Zumkeller_, Nov 19 2015
%o (Magma) [Gcd(n,6^n): n in [1..100]]; // _Vincenzo Librandi_, Feb 09 2016
%Y Cf. A003586, A065330, A006519, A038500, A007814, A007949, A007310, A047229, A071521, A322026.
%Y Related to A053165 via A225546.
%Y Cf. A126760 (ordinal transform of this sequence, from its term a(1) = 1 onward).
%K mult,nonn,easy
%O 1,2
%A _Reinhard Zumkeller_, Oct 29 2001