reviewed
approved
reviewed
approved
proposed
reviewed
editing
proposed
1, 4, 15, 7, 209, 13, 2911, 97, 901, 181, 564719, 193, 7865521, 2521, 6989, 18817, 1525870529, 2701, 21252634831, 37441, 6779137, 489061, 4122901604639, 37633, 274758906449, 6811741, 6575588101, 1037623, 11140078609864049, 40321, 155161278879431551
ab(n) = my(v=divisors(n)); prod(i=1, #v, T(v[i])^moebius(n/v[i]))
ba(n) = if(isprime(n)&&!(12%n), ab(n), ab(n)/gcd(n, ab(n)))
(PARI) T(n) = ([4, -1; 1, 0]^n)[2, 1]
a(n) = my(v=divisors(n)); prod(i=1, #v, T(v[i])^moebius(n/v[i]))
b(n) = if(isprime(n)&&!(12%n), a(n), a(n)/gcd(n, a(n)))
1, 4, 15, 7, 209, 13
1,2
allocated
nonn
Jianing Song, Aug 06 2019
approved
editing
allocated for Jianing Song
recycled
allocated
recycled
approved