%I#34 2023年9月23日03:11:12

%S 3,2,15,6,35,3,7,5,11,30143,42195,14255,2323,10399110483,33,23，

%电话：39,3182,87210899,151023,171155,341295,381443,9515991051763，

%电话：462215506235138,47,6,51,262703,782915,21333520333638703599

%N N^2-1除以其最大平方因子。

%换句话说，n^2-1的平方自由部分。

%C x^2-m*y^2=1有x=n解的最小m。

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

%F a（n）=A007913（n^2-1）。

%F a（n）=A005563（n-1）/A008833（n^2-1））。【Reinhard Zumkeller，2011年11月26日；Georg Fischer修正，2022年12月10日】

%e a（6）=35，因为6^2-1=35本身是平方自由的。

%e 7^2-1=48=A005563（6），其最大平方因子为A008833（48）=16，因此a（7）=48/16=3。

%t a[n_]：=次数@@（#[[1]]^Mod[#[2]]，2]和/@FactorInteger[n^2-1]）

%t表[（n^2-1）/Max[选择[Divisors[n^2-1]，IntegerQ[Sqrt[#]]&]]，{n，2,60}]（*_哈维P.Dale_，2019年12月8日*）

%o（PARI）a（n）=核心（n*n-1）；\\_David Wasserman，2005年3月7日

%o（哈斯克尔）

%o a068310 n=f 1$a027746_row（n^2-1），其中

%o f y[]=y

%o f y[p]=y*p

%o f y（p:ps'@（p'：ps））|p==p'=f y ps

%o|否则=f（y*p）ps'

%o--_Reinhard Zumkeller，2011年11月26日

%Y参见A002350、A007913、A067872、A033314、A027746、A175607。

%放松，好，不

%氧2,1

%A _Lekraj Beedassy，2002年2月25日

%E编辑：Dean Hickerson，2002年3月19日

%E条目由N.J.A.Sloane修订，2007年4月27日