|
|
|
|
1, 2, 2, 1, 2, 9, 1, 5, 2, 3, 9, 25, 1, 15, 5, 7, 2, 81, 3, 125, 9, 25, 25, 49, 1, 9, 15, 35, 5, 21, 7, 11, 2, 81, 81, 125, 3, 375, 125, 343, 9, 225, 25, 245, 25, 49, 49, 121, 1, 135, 9, 175, 15, 105, 35, 7, 5, 21, 21, 55, 7, 33, 11, 13, 2, 729, 81, 3125, 81, 625, 125, 2401, 3, 1125, 375, 343, 125, 147, 343, 1331
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
COMMENTS
|
|
|
LINKS
|
|
|
FORMULA
|
|
|
PROG
|
(PARI)
A005940(n) = { my(p=2, t=1); n--; until(!n\=2, if((n%2), (t*=p), p=nextprime(p+1))); t };
A054429(n) = ((3<<#binary(n\2))-n-1);
(Python)
from math import gcd
from sympy import nextprime
c, p, k = 1, 1, n
while k:
c *= (p:=nextprime(p))**(s:=(~k&k-1).bit_length())
k >>= s+1
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,frac
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|