; A003415 o=0:a(n)=n’=n的算术导数:a(0)=a(1)=0,a(素数)=1,a(mn)=m*a(n)+n*a(m)。; 由Antti Karttunen手动编码2021-02-28,https://github.com/karttu; 对于Christian Krause的LODA系统(一种汇编语言、计算模型和挖掘整数序列的工具); 可在获取https://github.com/ckrause/loda; 这是在dif-instruction或任何后来的增强之前,用该语言的最早版本之一编写的。;; 基本上,我们实现了以下PARI脚本的算法:;; A003415(n)={my(s=0,m=1,spf);而(n>1,spf=A020639(n);n/=spf;s+=m*n;m*=spf)(s);};; 这是根据Zumkeller的非分支递归公式得出的迭代解:; a(n)=a(A032742(n))*A020639(n)+A032741(n),对于n>1。;; 移动$1.0;初始化result-register,结果(这是一个和)就是这样构造的(LODA中实际上不需要这样的零位!)mov 2.2美元;这是迄今为止我们遇到的最小的主护目镜。移动5.1美元;上述PARI脚本中的当前产品m。lpb$0,1;启动主回路。当0美元中没有剩余的东西时,我们就停止了。保证每次迭代都会减少。移动$3,$0;剩下的0美元是求其最小素因子的安全上限。lpb$3.1;在这个子循环中,找到除$0的下一个素数>=$2,即=A020639($0)。移动$4,$0模块$4,$2cmp$4.0cmp$4.03美元以下,4美元以下;如果$2没有除n,则减去1(继续搜索),否则减去0,退出循环。增加2.1美元;注意:当我们退出时,这个递增的最后一个实例被丢弃,并且在$2中仍然有spf-divisor液化石油气; 现在,对于$0>0,我们在$2中有lpf=A020639(n)。分区$0,$2;从0美元中除以该(当前)最小质数因子的一个实例。移动$4,$0mul 4美元,5美元加上1美元、4美元;s+=m*n;(如PARI脚本)。5美元,2美元;m*=spf;(同上)液化石油气; 结果现在是1美元。