ad[n_]:=使用[{f=FactorInteger[n]},n*Total[f[[All,2]]/f[[Al,1]]];
okQ[n_,k_]:=如果[Not[CompositeQ[k-n]&&CompositeQ[k+n]],False,ad[k-n]==ad[k+n]];
a[n_]:=For[k=1,True,k++,If[okQ[n,k],Print[“a(”,n,“)=“,k];返回[k]]];
数组[a,46](*Jean-François Alcover公司2017年12月20日*)