A324187 a(n) = A106315(A163511(n)).

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

Links

Antti Karttunen, Table of n, a(n) for n = 0..16384

Antti Karttunen, Data supplement: n, a(n) computed for n = 0..65537

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));
A324187(n) = A106315(A163511(n));
(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. A054429, A106315, A163511, A324057, A324183, A324184, A324188, A324189.

Cf. A324199 (positions of zeros).

Sequence in context: A285641 A257856 A021661 * A219120 A324996 A359426

Adjacent sequences: A324184 A324185 A324186 * A324188 A324189 A324190

Keyword

nonn

Author

Antti Karttunen, Feb 17 2019

Status

approved