|
（PARI）usigma（n）={局部（f，s=1）；f=因子（n）；对于（i=1，矩阵大小（f）[1]，s*=1+f[i，1]^f[i、2]）；返回}
{a=0；对于（n=11000，a+=usigma（n）；写入（“b064609.txt”，n，“”，a））}\\哈里·史密斯2009年9月20日
（Python）
从sympy.theory.factor导入审核sigma
定义a（n）：返回和（范围（1，n+1）中j的udivisor_sigma（j，1））
打印（[a（n）代表范围（1101）中的n）]#因德拉尼尔·戈什2017年3月18日
| 