#该b文件由Antti Karttunen计算,2020-02-07#在GP/PARI CALCULATOR版本2.9.4(发布)和运行linux(x86-64/GMP-6.1.2内核)64位版本的amd64中,编译时间:2017年12月19日,gcc版本7.3.0(Ubuntu 7.3.0-1ubuntu1),线程引擎:pthread#使用以下程序:##默认值(parisizemax,2^31);#A097246(n)={my(f=因子(n));prod(i=1,#f~,(下一素数(f[i,1]+1)^;#A097248(n)={my(k=A097246(n));而(k<>n,n=k;k=A097 246(k));k;};#isA331751(n)=(A097248(2*n)==A09724八(σ(n)));#k=0;对于(n=1,2^25,如果(isA331751(n),k++;写入(“b331751.txt”,k,“”,n);打印(n“->”,系数(n)));1 22 63 274 285 846 2707 4968 10539 112010 148811 162512 163813 336014 378015 487516 812817 1053018 2438419 6696020 14742021 16740022 40622423 61122624 77500025 87296026 94325027 109728028 124569929 125533830 130353331 146432032 168640033 174002434 192237535 195250036 201162537 219340838 232500039 261144040 261888041 282975042 294102943 439296044 490921245 505920046 576712547 585750048 603487549 652860050 665683251 1245699052 1273387553 1303533054 1410048055 1713757556 2631667557 2722048058 2936290559 2941029060 33550336