|
%I #34 Apr 21 2021 07:14:08
%S 1,1,2,1,1,2,3,3,6,1,1,2,1,1,2,3,3,6,4,5,10,4,5,10,12,15,30,1,1,2,1,1,
%T 2,3,3,6,1,1,2,1,1,2,3,3,6,4,5,10,4,5,10,12,15,30,5,7,14,5,7,14,15,21,
%U 42,5,7,14,5,7,14,15,21,42,20,35,70,20,35,70
%N The second of a pair of coprime numbers whose factorizations depend on the ternary representation of n (See Comments for precise definition).
%C For n >= 0, with ternary representation Sum_{i=1..k} t_i * 3^e_i (all t_i in {1, 2} and all e_i distinct and in increasing order):
%C - let S(0) = A000961 \ { 1 },
%C - and S(i) = S(i-1) \ { p^(f + j), with p^f = the (e_i+1)-th term of S(i-1) and j > 0 } for any i=1..k,
%C - then a(n) = Product_{i=1..k such that t_i=2} "the (e_i+1)-th term of S(k)".
%C See A289815 for the first coprime number and additional comments.
%C The number of distinct prime factors of a(n) equals the number of twos in the ternary representation of n.
%H Rémy Sigrist, <a href="/A289816/b289816.txt">Table of n, a(n) for n = 0..10000</a>
%F a(n) = A289815(A004488(n)) for any n >= 0.
%F a(A005836(n)) = 1 for any n > 0.
%F a(2 * A005836(n)) = A289272(n-1) for any n > 0.
%e For n=42:
%e - 42 = 2*3^1 + 1*3^2 + 1*3^3,
%e - S(0) = { 2, 3, 4, 5, 7, 8, 9, 11, 13, 16, 17, 19, 23, 25, 27, 29, ... },
%e - S(1) = S(0) \ { 3^(1+j) with j > 0 }
%e = { 2, 3, 4, 5, 7, 8, 11, 13, 16, 17, 19, 23, 25, 29, ... },
%e - S(2) = S(1) \ { 2^(2+j) with j > 0 }
%e = { 2, 3, 4, 5, 7, 11, 13, 17, 19, 23, 25, 29, ... },
%e - S(3) = S(2) \ { 5^(1+j) with j > 0 }
%e = { 2, 3, 4, 5, 7, 11, 13, 17, 19, 23, 29, ... },
%e - a(42) = 3.
%o (PARI) a(n) = my (v=1, x=1); \
%o for (o=2, oo, \
%o if (n==0, return (v)); \
%o if (gcd(x,o)==1 && omega(o)==1, \
%o if (n % 3, x *= o); \
%o if (n % 3==2, v *= o); \
%o n \= 3; \
%o ); \
%o );
%o (Python)
%o from sympy import gcd, primefactors
%o def omega(n): return 0 if n==1 else len(primefactors(n))
%o def a(n):
%o v, x, o = 1, 1, 2
%o while True:
%o if n==0: return v
%o if gcd(x, o)==1 and omega(o)==1:
%o if n%3: x*=o
%o if n%3==2:v*=o
%o n //= 3
%o o+=1
%o print([a(n) for n in range(101)]) # _Indranil Ghosh_, Aug 02 2017
%Y Cf. A000961, A004488, A289815.
%K nonn,base,look
%O 0,3
%A _Rémy Sigrist_, Jul 12 2017
|
