# This b-file was computed by Antti Karttunen, 2020-02-13 # 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); # v156552sigs = readvec("a156552.txt"); \\ Factorization file for A156552 prepared by Hans Havermann, available at https://oeis.org/A156552/a156552.txt # A323243(n) = if(n<=2,n-1,my(prsig=v156552sigs[n],ps=prsig[1],es=prsig[2]); prod(i=1,#ps,((ps[i]^(1+es[i]))-1)/(ps[i]-1))); # A048675(n) = { my(f = factor(n)); sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2; }; # A087808(n) = if(n<1, 0, if(n%2==0, 2*A087808(n/2), A087808((n-1)/2)+1)); # isA322225(n) = (A087808(A323243(n)) == 2*A048675(n)); # k=0; for(n=2,10000,if(isA322225(n),k++; write("b322225.txt", k, " ", n))); 1 4 2 9 3 12 4 20 5 44 6 52 7 60 8 108 9 124 10 125 11 132 12 140 13 156 14 172 15 188 16 204 17 236 18 300 19 308 20 396 21 412 22 436 23 476 24 492 25 612 26 644 27 700 28 836 29 876 30 884 31 891 32 924 33 972 34 980 35 1004 36 1044 37 1092 38 1100 39 1116 40 1148 41 1188 42 1196 43 1236 44 1260 45 1268 46 1292 47 1300 48 1308 49 1372 50 1380 51 1476 52 1620 53 1628 54 1724 55 1860 56 1900 57 2140 58 2244 59 2324 60 2356 61 2444 62 2460 63 2652 64 2660 65 2700 66 2820 67 2860 68 2892 69 2900 70 3108 71 3124 72 3300 73 3332 74 3380 75 3388 76 3556 77 3588 78 3596 79 3740 80 3796 81 4095 82 4100 83 4164 84 4284 85 4332 86 4380 87 4444 88 4484 89 4548 90 4588 91 4620 92 4708 93 4756 94 4836 95 4876 96 5060 97 5076 98 5148 99 5236 100 5292 101 5508 102 5548 103 5612 104 5652 105 5684 106 5732 107 5868 108 5900 109 5929 110 5964 111 6052 112 6132 113 6372 114 6444 115 6460 116 6612 117 6644 118 6860 119 6916 120 6956 121 7068 122 7084 123 7196 124 7260 125 7332 126 7372 127 7476 128 7500 129 7540 130 7700 131 7772 132 7828 133 8092 134 8428 135 8556 136 8748 137 8820 138 8900 139 9156 140 9180 141 9196 142 9348 143 9500 144 9548 145 9628 146 9660 147 9844 148 9916 149 9940 150 9996