%I#16 2023年12月8日12:31:31

%S 1,1,2,3,4,0,6,7,8,2,10,2,12,4,6,15,16,3,18,8,10,8,22,6,24,10,26,14，

%电话28，-12,30,31,18,14,22,17,36,16,22,20,40，-12,42,26,27,20,46,14,48,17，

%U 30,32,52,12,38,34,26,58，-18,60,28,43,63,46，-12,66,44,42，-4,70,45,72,34,41,50,58，-12,78,44,80,38,82，-14,62

%N a（N）=N-A308135（N），其中A308135。

%H Antti Karttunen，n的表，n=1..20000的a（n）</a>

%F a（n）=n-A308135（n）=n-（σ（n）-A057723（n））。

%F a（n）=A336563（n）+A033879（n）。【由_Georg Fischer修订，2022年12月13日】

%F和{k=1..n}a（k）~c*n^2/2，其中c=A065487-A013661+1=0.586357_Amiram Eldar，2023年12月8日

%tf[p，e]：=（p^（e+1）-1）/（p-1）；fc[p，e_]：=f[p，e]-1；a[1]=1；a[n_]：=n-次数@@f@@@（fct=FactorInteger[n]）+次数@@fc@@fct；阵列[a，100]（*_Amiram Eldar_，2023年12月8日*）

%o（PARI）

%o A007947（n）=因子回收（因子（n）[，1]）；

%o A057723（n）={my（r=A007947（n））；（r*西格玛（n/r））；}；

%o A308135（n）=（西格玛（n）-A057723（n））；

%o A336564（n）=（n-A308135（n））；

%Y参见A000203、A033879、A057723、A308135、A336563、A336565、A336566。

%Y参考A013661、A065487。

%K符号

%氧1,3

%2020年7月27日，安蒂·卡图内