%I #19 May 07 2022 16:33:25
%S 1,1,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,
%T 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
%U 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,1
%N a(n) = 1 if A351546(n) is a unitary divisor of n, otherwise 0. Here A351546(n) is the largest unitary divisor of sigma(n) coprime with A003961(n).
%C For all known triperfect numbers, n = 1..6, a(A005820(n)) = 1.
%C For all known 5-multiperfect numbers, n = 1..65, a(A046060(n)) = 1.
%C For all known multiperfects m such that sigma(m) is also multiperfect, n = 1..23, a(A323653(n)) = 1.
%C Observation: Apparently, for no other odd terms than 1 of A006872, a(2*A006872(n)) = 1.
%H Antti Karttunen, <a href="/A353633/b353633.txt">Table of n, a(n) for n = 1..100000</a>
%H <a href="/index/Ch#char_fns">Index entries for characteristic functions</a>
%H <a href="/index/Pri#prime_indices">Index entries for sequences computed from indices in prime factorization</a>
%H <a href="/index/Si#SIGMAN">Index entries for sequences related to sigma(n)</a>
%F a(n)) = [1 == A353668(n)] * [1 == gcd(A351546(n),A353667(n))], where [ ] are the Iverson brackets.
%o (PARI)
%o A003961(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
%o A351546(n) = { my(f=factor(sigma(n)),u=A003961(n)); prod(k=1,#f~,f[k,1]^((0!=(u%f[k,1]))*f[k,2])); };
%o A353633(n) = { my(u=A351546(n)); (!(n%u) && 1==gcd(u,n/u)); };
%Y Characteristic function of A351551.
%Y Cf. A000203, A003961, A005820, A351546, A353667, A353668.
%Y Cf. also A005820, A006872, A046060, A323653.
%K nonn
%O 1
%A _Antti Karttunen_, May 04 2022