#2018-08-30在GP/PARI CALCULATOR版本2.9.4（发布）和运行linux（x86-64/GMP-6.1.2内核）64位版本的d64中计算，编译时间：2017年12月19日，gcc版本7.3.0（Ubuntu 7.3.0-1ubuntu1），线程引擎：pthread#使用以下程序：##up_to=2^14；#A064549（n）={my（f=因子（n））；对于（i=1，#f~，f[i，2]++）；因子回复（f）；}；#DirSqrt（v）={my（n=#v，u=向量（n））；u[1]=1；对于（n=2，n，u[n]=（v[n]/v[1]-sumdiv（n，d，如果（d>1&&d-->