# This b-file was computed by Antti Karttunen, 2020-04-06 # 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); # search_up_to = 2^23; # A333003list(up_to) = { my(u=vector(up_to), v=vector(up_to)); u[1] = v[1] = 1; for(n=2,up_to, my(ps=factor(n)[, 1]~); u[n] = vecsum(apply(p -> u[n-n/p], ps)); v[n] = (u[n]*n)+vecsum(apply(p -> v[n-n/p], ps))); vector(up_to, n, denominator(v[n]/u[n])); }; # v333003 = A333003list(search_up_to); # A333003(n) = v333003[n]; # k=0; for(n=1,search_up_to,if(1==A333003(n),k++;write("b333785.txt", k, " ", n))); 1 1 2 2 3 3 4 4 5 5 6 8 7 14 8 16 9 17 10 20 11 25 12 32 13 50 14 64 15 68 16 82 17 83 18 84 19 100 20 101 21 125 22 128 23 162 24 163 25 170 26 235 27 243 28 256 29 257 30 272 31 289 32 456 33 457 34 512 35 548 36 621 37 1024 38 1028 39 1040 40 2048 41 4096 42 4112 43 7225 44 8192 45 8738 46 9248 47 13058 48 16384 49 16480 50 16481 51 17476 52 17477 53 32128 54 32768 55 65536 56 65537 57 131072 58 132098 59 262144 60 262148 61 264196 62 278528 63 279616 64 297092 65 524288 66 933424 67 1048576 68 1048592 69 2097152 70 4194304 71 4210688 72 4456456 73 4473856 74 5242880 75 8388608