%I#31 2022年10月2日13:38:54

%S 1,3,5,9,11,15,17,23,27,31,33,41,43,47,51,60,62,70,72,80,84,88,90102，

%电话：106110116124126134136148152156160176178182186198200，

%电话：20821021822623023225025426266276288292308312314330

%N gcd（x，y）与x*y的除数之和。

%C A124315的部分金额。

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

%H Masum Billal，<a href=“https://arxiv.org/abs/2206.05023“>GCD-和推广的渐近结果，arXiv:2206.05023[math.NT]，2022。

%H Adrian W.Dudek，<a href=“http://arxiv.org/abs/1602.03555“>关于错误处理欧几里德引理的成功，arXiv:1602.03555[math.HO]，2016。见第3页备注1。

%H Adrian W.Dudek，<a href=“https://doi.org/10.4169/amer.math.monthly.123.9.924“>关于错误处理欧几里德引理的成功，《美国数学月刊》，第123卷，第9期（2016年），924-927。

%H Randell Heyman，<a href=“https://arxiv.org/abs/2003.13937“>涉及除数函数和GCD函数的求和，arXiv:2003.13937[math.NT]，2020。

%F a（n）=总和{k=1..层（平方（n））}_Daniel Suteu，2019年1月8日

%F a（n）=n*zeta（2）*（log（n）+2*gamma-1+2*zeta’（2）/zeta（二））+O（sqrt（n）*log（n）），其中gamma是Euler-Mascheroni常数A001620_Daniel Suteu，2019年1月11日

%t表[Total@Flatten@Map[Function[k，DivisorSigma[0，GCD[#，k]]&/@Select[Range@n，#k<=n&]]，Range@n]，{n，60}]（*_Michael De Vlieger_，2016年2月12日*）

%o（PARI）a（n）=总和（k=1，n，sumdiv（k，d，numdiv（gcd（d，k/d）））；

%o（PARI）a（n）=总和（k=1，平方（n），2*总和_Daniel Suteu，2019年1月8日

%o（PARI）a（n）=总和（k=1，n，总和（j=1，sqrt（n/k），楼层（n/k/j^2））；\\_Benoit Cloitre_，2022年10月2日

%Y参考A000005，A124315。

%Y参考A001620、A013661、A306016。

%K nonn公司

%O 1,2号机组

%2016年2月12日，马库斯