OFFSET
1,2
COMMENTS
LINKS
Jianing Song, Table of n, a(n) for n = 1..1000
FORMULA
EXAMPLE
12 = 3 * 2 * 2, so a(12) = 5^(3-1) * 13^(2-1) * 17^(2-1) = 5525.
15 = 5 * 3, so a(15) = 5^(5-1) * 13^(3-1) = 105625.
PROG
(PARI) a(n) = my(f=factor(n), w=omega(n), p=1, product=1); forstep(i=w, 1, -1, for(j=1, f[i, 2], p=nextprime(p+1); while(!(p%4==1), p=nextprime(p+1)); product *= p^(f[i, 1]-1))); product
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Jianing Song, Apr 24 2021
STATUS
approved