Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
A356161
a(n) = 1 if the smallest positive k such that n divides k*A276086(k) is equal to n, otherwise 0, where A276086 is primorial base exp-function.
4
1, 0, 0, 1, 1, 0, 1, 1, 0, 0, 1, 0, 1, 0, 0, 1, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1
OFFSET
1
FORMULA
a(n) = [A355944(n) == n] = [A355945(n) == 0], where [ ] is the Iverson bracket.
PROG
(PARI)
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A356161(n) = for(k=1, oo, if((k*A276086(k))%n==0, return(n==k)));
CROSSREFS
Characteristic function of A356160.
Sequence in context: A089012 A083035 A359422 * A187074 A188398 A288929
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jul 28 2022
STATUS
approved