%I#60 2022年10月22日08:05:54

%S 1,1,1,1,5,1,7,1,1,5,11,13,7,5,1,17,1,19,5,7,11,23,1,25,13,7,29,5，

%电话31,1,17,35,1,37,19,13,5,41,7,43,11,5,23,47,1,49,25,17,13,53,1,55，

%U 7,19,29,59,5,61,31,7,165,11,67,17,23,35,71,1,73,37,25,19,77,13,79,5,1

%N a（N）=最大值{k|gcd（N，k）=k，gcd（k，6）=1}。

%C Bennett、Filaseta和Trifonov证明，如果n>8，则a（n^2+n）>n^0.285_Charles R Greathouse IV，2014年5月21日

%H Reinhard Zumkeller，n的表，n=1..10000的a（n）</a>

%H M.A.Bennett、M.Filaseta和O.Trifonov，<A href=“http://www.math.ubc.ca网站/~bennett/BFTpaper0207.pdf“>关于连续整数的因式分解，J.Reine Angew.Math.629（2009），第171-200页。

%F a（n）*A065331（n）=n。

%F与a（2^e）=1，a（3^e）=1，a（p^e）=p^e，p>3_Vladeta Jovovic_，2001年11月2日

%F A106799（n）=A001222（a（n））。-_Reinhard Zumkeller_，2005年5月19日

%F a（1）=1；则a（2n）=a（n），a（2n+1）=a_Benoit Cloitre_，2007年6月4日

%F Dirichlet g.F.zeta（s-1）*（1-2^（1-s））*（1-3 ^（1-s））/（1-2 ^（-s））*_R.J.Mathar，2011年7月4日

%F a（n）=A038502（A000265（n））_Reinhard Zumkeller，2011年7月6日

%F a（n）=n/GCD（n，6^n）_Stanislav Sykora，2016年2月8日

%F和{k=1..n}a（k）~（1/4）*n^2.-_Amiram Eldar，2022年10月22日

%e a（30）=5。

%p A065330:=程序（n）

%p局部a，f，p，e；

%pa:=1；

%ifactors（n）[2]do中f的p

%p p：=op（1，f）；

%p e：=op（2，f）；

%如果p>3，则为p

%p a：=a*p ^e；

%p end if；

%p端do：

%p a；

%p end程序：#_R.J.Mathar_，2012年7月12日

%p与（padic）：a:=n->n/（2^ordp（n，2）*3^ordp；

%p序列（a（n），n=1..81）；#_Peter Luschny_，2014年3月25日

%t f[n_]：=时间@@（第一个@#^最后一个@#&/@选择[因子整数@n，第一个@#！=2&&First@#！=3 &]); 数组[f，81]（*_Robert G.Wilson v_，2006年8月18日*）

%t f[n_]：=分母[6^n/n]；阵列[f，100]（*_Vladimir Joseph Stephan Orlovsky_，2011年2月16日*）

%t表[n/GCD[n，6^n]，{n，100}]（*_文森佐图书馆，2016年2月9日*）

%o（PARI）a（n）=如果（n<2,1，如果（n%2，如果（n%3，n，a（n/3）），a（n/2）））

%o（PARI）a（n）=n\gcd（n，6^n）\\效率不高，但很简单_Stanislav Sykora，2016年2月8日

%o（PARI）a（n）=n>>估价（n，2）/3^估价

%o（哈斯克尔）

%o a065330=a038502。a000265---Reinhard Zumkeller，2011年7月6日

%o（岩浆）[n div Gcd（n，6^n）：n in[1..100]]；//_文森佐·利班迪（Vincenzo Librandi），2016年2月9日

%Y参见A065331、A000265、A038502、A165725。

%K mult，nonn（非n）

%O 1,5型

%A _Reinhard Zumkeller，2001年10月29日