[go: up one dir, main page]

login
A324187
a(n) = A106315(A163511(n)).
7
0, 1, 5, 2, 2, 1, 0, 4, 18, 28, 30, 13, 16, 12, 4, 6, 3, 42, 72, 32, 51, 78, 21, 33, 12, 36, 24, 44, 36, 20, 8, 10, 67, 2, 168, 1, 176, 504, 128, 172, 84, 10, 312, 102, 32, 198, 75, 97, 108, 120, 144, 58, 48, 72, 128, 20, 50, 66, 48, 4, 0, 36, 16, 12, 4, 731, 372, 3126, 625, 6, 785, 801, 456, 1332, 768, 1720, 540, 232, 688, 932, 145, 660
OFFSET
0,3
FORMULA
a(n) = A106315(A163511(n)) = (A163511(n)*A324183(n)) mod A324184(n).
For n > 0, a(n) = A324057(A054429(n)).
PROG
(PARI)
A106315(n) = (n*numdiv(n) % sigma(n));
A163511(n) = if(!n, 1, my(p=2, t=1); while(n>1, if(!(n%2), (t*=p), p=nextprime(1+p)); n >>= 1); (t*p));
(PARI)
A324183(n) = if(!n, 1, n = ((3<<#binary(n\2))-n-1); my(e=0, m=1); while(n>0, if(!(n%2), m *= (1+e); e=0, e++); n >>= 1); (m*(1+e)));
A324184(n) = if(!n, 1, my(p=2, mp=p*p, m=1); while(n>1, if(n%2, p=nextprime(1+p); mp = p*p, if((2==n)||!(n%4), mp *= p, m *= (mp-1)/(p-1))); n >>= 1); (m*(mp-1)/(p-1)));
A324187(n) = ((A163511(n)*A324183(n))%A324184(n));
CROSSREFS
Cf. A324199 (positions of zeros).
Sequence in context: A285641 A257856 A021661 * A219120 A324996 A359426
KEYWORD
nonn
AUTHOR
Antti Karttunen, Feb 17 2019
STATUS
approved