# This b-file was computed by Antti Karttunen, 2019-12-28 # 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^32); # upto_e = 128; \\ 128 --> 285 terms # A018819(n) = if( n<1, n==0, if( n%2, A018819(n-1), A018819(n/2)+A018819(n-1))); \\ From A018819 # v018819 = vector(upto_e,n,A018819(n)); # A050377(n) = factorback(apply(e -> v018819[e], factor(n)[, 2])); # A283980(n) = {my(f=factor(n)); prod(i=1, #f~, my(p=f[i, 1], e=f[i, 2]); if(p==2, 6, nextprime(p+1))^e)}; \\ From A283980 # A330688list(e) = { my(lista = List([1, 2]), i=2, u = 2^e, t, m=0, v025487); while(lista[i] != u, if(2*lista[i] <= u, listput(lista,2*lista[i]); t = A283980(lista[i]); if(t <= u, listput(lista,t))); i++); v025487 = vecsort(Vec(lista)); lista = List([]); for(i=1,#v025487,if((t=A050377(v025487[i]))>m, listput(lista,t); m=t)); Vec(lista); }; # v330688 = A330688list(upto_e); # A330688(n) = v330688[n]; # for(n=1,#v330688,write("b330688.txt", n, " ", A330688(n))); 1 1 2 2 3 4 4 6 5 8 6 10 7 12 8 14 9 16 10 20 11 24 12 28 13 40 14 56 15 60 16 80 17 84 18 104 19 112 20 120 21 144 22 160 23 168 24 184 25 200 26 208 27 216 28 224 29 240 30 260 31 288 32 320 33 360 34 368 35 400 36 416 37 432 38 460 39 480 40 520 41 576 42 600 43 624 44 640 45 720 46 736 47 800 48 864 49 920 50 960 51 1040 52 1104 53 1120 54 1152 55 1200 56 1440 57 1456 58 1472 59 1480 60 1600 61 1840 62 2016 63 2080 64 2400 65 2576 66 2880 67 2960 68 3360 69 3680 70 3760 71 4032 72 4144 73 4160 74 4800 75 5152 76 5264 77 5760 78 5920 79 6720 80 7360 81 7520 82 8064 83 8288 84 9600 85 10304 86 10528 87 11840 88 13440 89 15040 90 16576 91 18240 92 19200 93 21056 94 22400 95 23680 96 25536 97 26560 98 26880 99 30080 100 31360 101 31584 102 32320 103 33152 104 36480 105 37184 106 38400 107 39104 108 42112 109 44800 110 45120 111 45248 112 47360 113 47424 114 48000 115 51072 116 53120 117 54144 118 54720 119 60160 120 62720 121 64640 122 65664 123 67200 124 72960 125 75200 126 76800 127 78208 128 78400 129 80640 130 84224 131 89600 132 91200 133 94720 134 94848 135 96000 136 97760 137 102144 138 106240 139 108288 140 112000 141 120320 142 125440 143 129280 144 131328 145 132800 146 134400 147 135360 148 145920 149 150400 150 152320 151 156416 152 156800 153 161280 154 161600 155 164160 156 179200 157 182400 158 189696 159 192000 160 195520 161 201600 162 212480 163 216576 164 224000 165 232960 166 236800 167 240640 168 258560 169 262656 170 265600 171 268800 172 270720 173 276224 174 300800 175 304640 176 313600 177 322560 178 323200 179 328320 180 336128 181 364800 182 371840 183 382464 184 391040 185 403200 186 448000 187 452480 188 465408 189 474240 190 478080 191 481280 192 483840 193 531200 194 541440 195 548352 196 582400 197 601600 198 646400 199 656640 200 690560 201 729600 202 761600 203 806400 204 840320 205 896000 206 908800 207 956160 208 967680 209 990080 210 1062400 211 1082880 212 1163520 213 1164800 214 1203200 215 1292800 216 1313280 217 1370880 218 1381120 219 1459200 220 1523200 221 1612800 222 1635840 223 1680640 224 1792000 225 1817600 226 1912320 227 1935360 228 1980160 229 2124800 230 2165760 231 2327040 232 2329600 233 2362880 234 2585600 235 2626560 236 2741760 237 2762240 238 2767360 239 2792448 240 3046400 241 3225600 242 3271680 243 3361280 244 3635200 245 3824640 246 3960320 247 4121600 248 4224000 249 4249600 250 4654080 251 4725760 252 4887040 253 4992000 254 5171200 255 5483520 256 5491200 257 5524480 258 5946880 259 6092800 260 6543360 261 6722560 262 7006720 263 7270400 264 7603200 265 7649280 266 7756800 267 7920640 268 8360960 269 8448000 270 8985600 271 9308160 272 9451520 273 9715200 274 9774080 275 9815040 276 9984000 277 10368000 278 10967040 279 10982400 280 11481600 281 11635200 282 11893760 283 12072960 284 12185600 285 13086720