%I#17 2021年11月30日09:54:54

%S 1,2,4,5,9,7,15,12,18,15,28,16,36,23,31,30,51,26,59,34,50,43,75,37,77，

%电话：52,72,55102,42112,69,90,73106,61141,84109,80159,66169,97119，

%U 108187,84185103155121218,97193126179142248,95262152185单位

%N a（N）=和{1<=k<=N，gcd（k，N）=1}层（N/k）。

%C和{k|n}a（k）=sum{k=1到n}d（k），其中d（k）是k的正除数。

%C等于A054525*A006218（A006218Mobius变换）_Gary W.Adamson_，2008年8月7日

%H Charles R Greathouse IV，n表，n=1..10000的a（n）</a>

%F a（n）也是Sum_{k|n}mu（n/k）（Sum_{j=1..k}d（j））和Sum__{k=1..n}phi（n，n/k。

%e a（6）=7，因为相对于6的素数和不超过6的数字是1和5，从而得出floor（6/1）+floor（5/5）=7。

%p a：=进程（n）局部s，j：s：=0：对于从1到n的j do，如果gcd（j，n）=1，则s：=s+楼层（n/j），否则s：=s:fiod：s：结束：seq（a（n），n=1..75）；

%t表[a:=选择[Range[n]，GCD[n，#]==1&]；总和[Floor[n/a[i]]]，{i，1，Length[a]}]，{n，1，60}]

%o（PARI）A116477（n）=总和（k=1，n，（gcd（k，n）==1）*楼层（n/k））\\迈克尔·波特，2010年3月1日

%o（PARI）A006218（n）=总和（k=1，平方（n），n\k）*2-平方

%o a（n，f=系数（n））=我的（K=1，n）；对于（i=1，#f~，如果（f[i，2]>1，K*=f[i，1]^（f[i，2]-1）；f[i，2]=1））；sumdiv（N=N/K，K，moebius（N/K）*A006218（K*K））\\查尔斯·格里塔斯四世，2021年11月30日

%Y参考A006218。A122191的行总和。

%Y参考A054525.-_Gary W.Adamson_，2008年8月7日

%K容易，不是

%O 1,2号机组

%A _罗伊查询，2006年3月18日

%E更多术语摘自德国电气公司和斯特凡·斯坦伯格公司，2006年4月1日