# Computed 2018-06-11 in GP/PARI CALCULATOR Version 2.9.4 (released) amd64 running mingw (x86-64/GMP-6.1.2 kernel) 64-bit version, compiled: Dec 20 2017, gcc version 4.9.1 (GCC), threading engine: single # with the following program (by Antti Karttunen): # # allocatemem(2^30); # up_to = 65537; # rgs_transform(invec) = { my(occurrences = Map(), outvec = vector(length(invec)), u=1); for(i=1, length(invec), if(mapisdefined(occurrences,invec[i]), my(pp = mapget(occurrences, invec[i])); outvec[i] = outvec[pp] , mapput(occurrences,invec[i],i); outvec[i] = u; u++ )); outvec; }; # A278233(n) = { my(p=0, f=vecsort((factor(Pol(binary(n))*Mod(1, 2))[, 2]), , 4)); prod(i=1, #f, (p=nextprime(p+1))^f[i]); }; # v305788 = rgs_transform(vector(up_to,n,A278233(n))); # A305788(n) = v305788[n]; # A305812(n) = if(1==n,0, my(m=1); fordiv(n,d,if((d>1)&&(d