OFFSET
1,2
LINKS
FORMULA
PROG
(PARI)
default(primelimit, 2^30);
A064989(n) = {my(f); f = factor(n); if((n>1 && f[1, 1]==2), f[1, 2] = 0); for (i=1, #f~, f[i, 1] = precprime(f[i, 1]-1)); factorback(f)};
for(n=1, 8191, write("b254118.txt", n, " ", A254118(n)));
(Scheme) ;; Two versions.
(Python)
from sympy import factorint, prevprime, floor
from operator import mul
from functools import reduce
def a064216(n):
f=factorint(2*n - 1)
return 1 if n==1 else reduce(mul, [prevprime(i)**f[i] for i in f])
def a254103(n):
if n==0: return 0
if n%2==0: return 3*a254103(n//2) - 1
else: return floor((3*(1 + a254103((n - 1)/2)))//2)
def a254116(n): return a064216(a254103(n))
def a(n): return (a254116(2*n + 1) - 1)//2
print([a(n) for n in range(1, 101)]) # Indranil Ghosh, Jun 06 2017
CROSSREFS
KEYWORD
nonn,look
AUTHOR
Antti Karttunen, Feb 05 2015
STATUS
approved
