This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
A002110(n) = prod(i=1,n,prime(i)); A108951(n) = { my(f=factor(n)); prod(i=1, #f~, A002110(primepi(f[i, 1]))^f[i, 2]) }; A276085(n) = { my(f = factor(n)); sum(k=1, #f~, f[k, 2]*A002110(primepi(f[k, 1])-1)); }; A328571(n) = { my(m=1, p=2); while(n, m *= (p^!!(n%p)); n = n\p; p = nextprime(1+p)); (m); }; A346093(n) = A276085(A328571(A108951(n)));
up_to = 65537;
rgs_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), u=1); for(i=1, length(invec), if(mapisdefined(om,invec[i]), my(pp = mapget(om, invec[i])); outvec[i] = outvec[pp] , mapput(om,invec[i],i); outvec[i] = u; u++ )); outvec; };
A002110(n) = prod(i=1,n,prime(i));
A108951(n) = { my(f=factor(n)); prod(i=1, #f~, A002110(primepi(f[i, 1]))^f[i, 2]) };
A328571(n) = { my(m=1, p=2); while(n, m *= (p^!!(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A346091(n) = A328571(A108951(n));
v351955 = rgs_transform(vector(up_to, n, A346091(n)));
A351955(n) = v351955[n];
For n = 24, which has primorial base representation (see A049345) "400" as 24 = 4*A002110(2) + 0*A002110(1) + 0*A002110(0) = 4*6 + 0*2 + 0*1, thus a(24) = prime(3)^4 * prime(2)^0 * prime(1)^0 = 5^4 = 625. For n = 35 = "1021" as 35 = 1*A002110(3) + 0*A002110(2) + 2*A002110(1) + 1*A002110(0) = 1*30 + 0*6 + 2*2 + 1*1, thus a(35) = prime(4)^1 * prime(2)^2 * prime(1) = 7 * 3*3 * 2 = 126.
b = MixedRadix[Reverse@ Prime@ Range@ 12]; Table[Function[k, Times @@ Power @@@ # &@ Transpose@ {Prime@ Range@ Length@ k, Reverse@ k}]@ IntegerDigits[n, b], {n, 0, 51}] (* Michael De Vlieger, Aug 23 2016, Version 10.2 *)
f[n_] := Block[{a = {{0, n}}}, Do[AppendTo[a, {First@ #, Last@ #} &@ QuotientRemainder[a[[-1, -1]], Times @@ Prime@ Range[# - i]]], {i, 0, #}] &@ NestWhile[# + 1 &, 0, Times @@ Prime@ Range[# + 1] <= n &]; Rest[a][[All, 1]]]; Table[Times @@ Flatten@ MapIndexed[Prime[#2]^#1 &, Reverse@ f@ n], {n, 0, 73}] (* Michael De Vlieger, Aug 30 2016, Pre-Version 10 *)
a[n0_] := Module[{m = 1, i = 1, n = n0, p}, While[n > 0, p = Prime[i]; m *= p^Mod[n, p]; n = Quotient[n, p]; i++]; m];
Table[a[n], {n, 0, 100}] (* Jean-François Alcover, Dec 01 2021, after Antti Karttunen's Sage code *)
A276086(n) = { my(i=0,m=1,pr=1,nextpr); while((n>0),i=i+1; nextpr = prime(i)*pr; if((n%nextpr),m*=(prime(i)^((n%nextpr)/pr));n-=(n%nextpr));pr=nextpr); m; }; \\ Antti Karttunen, May 12 2017
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); }; \\ (Better than above one, avoids unnecessary construction of primorials). - Antti Karttunen, Oct 14 2019
from sympy import prime
def a(n):
i=0
m=pr=1
while n>0:
i+=1
N=prime(i)*pr
if n%N!=0:
m*=(prime(i)**((n%N)/pr))
n-=n%N
pr=N
return m # Indranil Ghosh, May 12 2017, after Antti Karttunen's PARI code
from sympy import nextprime
def a(n):
m, p = 1, 2
while n > 0:
n, r = divmod(n, p)
m *= p**r
p = nextprime(p)
return m
print([a(n) for n in range(74)]) # Peter Luschny, Apr 20 2024
def A276086(n): m=1 i=1 while n>0: p = sloane.A000040(i) m *= (p**(n%p)) n = floor(n/p) i += 1 return (m) # Antti Karttunen, Oct 14 2019, after Indranil Ghosh's Python code above, and my own leaner PARI code from Oct 14 2019. This avoids unnecessary construction of primorials.
(define (A276086 n) (let loop ((n n) (t 1) (i 1)) (if (zero? n) t (let* ((p (A000040 i)) (d (modulo n p))) (loop (/ (- n d) p) (* t (expt p d)) (+ 1 i))))))
(definec (A276086 n) (if (zero? n) 1 (* (expt (A053669 n) (A276088 n)) (A276086 (A276093 n))))) ;; Needs macro definec from http://oeis.org/wiki/Memoization#Scheme
(definec (A276086 n) (if (zero? n) 1 (* (A053669 n) (A276086 (- n (A002110 (A276084 n))))))) ;; Needs macro definec from http://oeis.org/wiki/Memoization#Scheme
nn = 60; b = MixedRadix[Reverse@ Prime@ Range@ PrimePi[nn + 1]]; Table[FromDigits[#, b] &@ Reverse@ If[n == 1, {0}, Function[k, ReplacePart[Table[0, {PrimePi[k[[-1, 1]]]}], #] &@ Map[PrimePi@ First@ # -> Last@ # &, k]]@ FactorInteger@ n], {n, nn}] (* Version 10.2, or *)
f[w_List] := Total[Times @@@ Transpose@ {Map[Times @@ # &, Prime@ Range@ Range[0, Length@ w - 1]], Reverse@ w}]; Table[f@ Reverse@ If[n == 1, {0}, Function[k, ReplacePart[Table[0, {PrimePi[k[[-1, 1]]]}], #] &@ Map[PrimePi@ First@ # -> Last@ # &, k]]@ FactorInteger@ n], {n, 60}] (* Michael De Vlieger, Aug 30 2016 *)
A276085(n) = { my(f = factor(n), pr=1, i=1, s=0); for(k=1, #f~, while(i <= primepi(f[k, 1])-1, pr *= prime(i); i++); s += f[k, 2]*pr); (s); }; \\ Antti Karttunen, Nov 11 2024
from sympy import primorial, primepi, factorint
def a002110(n):
return 1 if n<1 else primorial(n)
def a(n):
f=factorint(n)
return sum(f[i]*a002110(primepi(i) - 1) for i in f)
print([a(n) for n in range(1, 101)]) # Indranil Ghosh, Jun 22 2017
a(12) = a(2^2) * a(3) = (2#)^2 * (3#) = 2^2 * 6 = 24 a(45) = (3#)^2 * (5#) = (2*3)^2 * (2*3*5) = 1080 (as 45 = 3^2 * 5).
a[n_] := a[n] = Module[{f = FactorInteger[n], p, e}, If[Length[f]>1, Times @@ a /@ Power @@@ f, {{p, e}} = f; Times @@ (Prime[Range[PrimePi[p]]]^e)]]; a[1] = 1; Table[a[n], {n, 1, 42}] (* Jean-François Alcover, Feb 24 2015 *)
Table[Times @@ Map[#1^#2 & @@ # &, FactorInteger[n] /. {p_, e_} /; e > 0 :> {Times @@ Prime@ Range@ PrimePi@ p, e}], {n, 42}] (* Michael De Vlieger, Mar 18 2017 *)
primorial(n)=prod(i=1,primepi(n),prime(i)) a(n)=my(f=factor(n)); prod(i=1,#f~, primorial(f[i,1])^f[i,2]) \\ Charles R Greathouse IV, Jun 28 2015
from sympy import primerange, factorint
from operator import mul
def P(n): return reduce(mul, [i for i in primerange(2, n + 1)])
def a(n):
f = factorint(n)
return 1 if n==1 else reduce(mul, [P(i)**f[i] for i in f])
print([a(n) for n in range(1, 101)]) # Indranil Ghosh, May 14 2017
def sharp_primorial(n): return sloane.A002110(prime_pi(n)) def p(f): return sharp_primorial(f[0])^f[1] [prod(p(f) for f in factor(n)) for n in range (1,51)] # Giuseppe Coppoletta, Feb 07 2015
2556 has primorial base expansion [1,1,1,1,0,0] as 1*A002110(5) + 1*A002110(4) + 1*A002110(3) + 1*A002110(2) = 2310 + 210 + 30 + 6 = 2556. That in turn is converted by A276086 to 13^1 * 11^1 * 7^1 * 5^1 = 5005, whose arithmetic derivative is 5' * 1001 + 1001' * 5 = 1*1001 + 311*5 = 2556, thus 2556 is one of the rare fixed points (A328110) of this sequence.
Block[{b = MixedRadix[Reverse@ Prime@ Range@ 12]}, Array[Function[k, If[# < 2, 0, # Total[#2/#1 & @@@ FactorInteger[#]] ] &@ Abs[Times @@ Power @@@ # &@ Transpose@{Prime@ Range@ Length@ k, Reverse@ k}]]@ IntegerDigits[#, b] &, 65, 0]] (* Michael De Vlieger, Mar 12 2021 *)
A003415(n) = {my(fac); if(n<1, 0, fac=factor(n); sum(i=1, matsize(fac)[1], n*fac[i, 2]/fac[i, 1]))}; \\ From A003415 A276086(n) = { my(i=0,m=1,pr=1,nextpr); while((n>0),i=i+1; nextpr = prime(i)*pr; if((n%nextpr),m*=(prime(i)^((n%nextpr)/pr));n-=(n%nextpr));pr=nextpr); m; }; A327860(n) = A003415(A276086(n));
A327860(n) = { my(s=0, m=1, p=2, e); while(n, e = (n%p); m *= (p^e); s += (e/p); n = n\p; p = nextprime(1+p)); (s*m); }; \\ (Standalone version) - Antti Karttunen, Nov 07 2019
a(12) = 10 because the prime divisors of 12 are 2 and 3 so we have: 12/2 + 12/3 = 6 + 4 = 10. - _Geoffrey Critzer_, Mar 17 2015
[0] cat [n*&+[1/p: p in PrimeDivisors(n)]:n in [2..80]]; // Marius A. Burtea, Jan 21 2020
A069359 := n -> add(n/d, d = select(isprime, numtheory[divisors](n))): seq(A069359(i), i = 1..20); # Peter Luschny, Jan 31 2012 # second Maple program: a:= n-> n*add(1/i[1], i=ifactors(n)[2]): seq(a(n), n=1..100); # Alois P. Heinz, Oct 23 2019
f[list_, i_] := list[[i]]; nn = 100; a = Table[n, {n, 1, nn}]; b =
Table[If[PrimeQ[n], 1, 0], {n, 1, nn}]; Table[DirichletConvolve[f[a, n], f[b, n], n, m], {m, 1, nn}] (* Geoffrey Critzer, Mar 17 2015 *)
a(n) = n*sumdiv(n, d, isprime(d)/d); \\ Michel Marcus, Mar 18 2015
a(n) = my(ps=factor(n)[,1]~);sum(k=1,#ps,n\ps[k]) \\ Franklin T. Adams-Watters, Apr 09 2015
from sympy import primefactors def A069359(n): return sum(n//p for p in primefactors(n)) # Chai Wah Wu, Feb 05 2022
def A069359(n) : D = filter(is_prime, divisors(n)) return add(n/d for d in D) print([A069359(i) for i in (1..20)]) # Peter Luschny, Jan 31 2012
nn = 65; b = MixedRadix[Reverse@ Prime@ Range[IntegerLength[nn, 2] - 1]]; Table[FromDigits[IntegerDigits[n, 2], b], {n, 0, 65}] (* Version 10.2, or *)
Table[Total[Times @@@ Transpose@ {Map[Times @@ # &, Prime@ Range@ Range[0, Length@ # - 1]], Reverse@ #}] &@ IntegerDigits[n, 2], {n, 0, 65}] (* Michael De Vlieger, Aug 26 2016 *)
A276156(n) = { my(s=0, p=1, r=1); while(n, if(n%2, s += r); n>>=1; p = nextprime(1+p); r *= p); (s); }; \\ Antti Karttunen, Feb 03 2022
from sympy import prime, primorial, primepi, factorint
from operator import mul
def a002110(n): return 1 if n<1 else primorial(n)
def a276085(n):
f=factorint(n)
return sum([f[i]*a002110(primepi(i) - 1) for i in f])
def a019565(n): return reduce(mul, (prime(i+1) for i, v in enumerate(bin(n)[:1:-1]) if v == '1')) # after Chai Wah Wu
def a(n): return 0 if n==0 else a276085(a019565(n))
print([a(n) for n in range(101)]) # Indranil Ghosh, Jun 23 2017
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); }; A342002(n) = A342001(A276086(n)); \\ Uses also code from A342001.
A342002(n) = { my(s=0, m=1, p=2, e); while(n, e = (n%p); m *= p^(e>0); s += (e/p); n = n\p; p = nextprime(1+p)); (s*m); }; \\ Antti Karttunen, Mar 12 2021
A342002(n) = { my(s=0, p=2, e); while(n, e = (n%p); s += (e/p); n = n\p; p = nextprime(1+p)); numerator(s); }; \\ (Taking denominator instead would give A328571) - Antti Karttunen, Mar 12 2021
Comments