# This b-file was computed by Antti Karttunen, 2020-02-07 # in GP/PARI CALCULATOR Version 2.9.4 (released) amd64 running linux (x86-64/GMP-6.1.2 kernel) 64-bit version, compiled: Dec 19 2017, gcc version 7.3.0 (Ubuntu 7.3.0-1ubuntu1), threading engine: pthread # with the following program: # # default(parisizemax,2^31); # A097246(n) = { my(f=factor(n)); prod(i=1, #f~, (nextprime(f[i, 1]+1)^(f[i, 2]\2))*((f[i, 1])^(f[i, 2]%2))); }; # A097248(n) = { my(k=A097246(n)); while(k<>n, n = k; k = A097246(k)); k; }; # isA331751(n) = (A097248(2*n)==A097248(sigma(n))); # k=0; for(n=1,2^25,if(isA331751(n),k++; write("b331751.txt", k, " ", n); print(n " -> ", factor(n)))); 1 2 2 6 3 27 4 28 5 84 6 270 7 496 8 1053 9 1120 10 1488 11 1625 12 1638 13 3360 14 3780 15 4875 16 8128 17 10530 18 24384 19 66960 20 147420 21 167400 22 406224 23 611226 24 775000 25 872960 26 943250 27 1097280 28 1245699 29 1255338 30 1303533 31 1464320 32 1686400 33 1740024 34 1922375 35 1952500 36 2011625 37 2193408 38 2325000 39 2611440 40 2618880 41 2829750 42 2941029 43 4392960 44 4909212 45 5059200 46 5767125 47 5857500 48 6034875 49 6528600 50 6656832 51 12456990 52 12733875 53 13035330 54 14100480 55 17137575 56 26316675 57 27220480 58 29362905 59 29410290 60 33550336