# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/ Search: id:a346099 Showing 1-1 of 1 %I A346099 #8 Jul 19 2021 01:25:55 %S A346099 1,2,3,4,5,3,7,1,1,5,11,3,13,7,5,1,17,1,19,5,7,11,23,3,1,13,3,7,29,5, %T A346099 31,2,11,17,7,1,37,19,13,1,41,1,43,11,1,23,47,1,1,1,17,13,53,9,11,7, %U A346099 19,29,59,1,61,31,1,4,13,11,67,17,23,1,71,3,73,37,25,19,1,13,79,1,1,41,83,1,17,43,29,11,89 %N A346099 a(n) = gcd(n, A346098(n)). %C A346099 Only powers of primes (A000961) occur as terms. A346100 lists the exponents. %H A346099 Index entries for sequences computed from indices in prime factorization %H A346099 Index entries for sequences related to primorial base %H A346099 Index entries for sequences related to primorial numbers %F A346099 a(n) = gcd(n, A346098(n)) = gcd(n, A064989(A319626(A324886(n)))). %o A346099 (PARI) %o A346099 A064989(n) = { my(f = factor(n)); if((n>1 && f[1,1]==2), f[1,2] = 0); for (i=1, #f~, f[i,1] = precprime(f[i,1]-1)); factorback(f); }; %o A346099 A319626(n) = (n / gcd(n, A064989(n))); %o A346099 A346099(n) = gcd(n, A064989(A319626(A324886(n)))); \\ Rest of program in A324886. %Y A346099 Cf. A000961, A064989, A319626, A324886, A346095, A346096, A346097, A346098, A346100. %Y A346099 Cf. A346090 (positions of ones). %K A346099 nonn %O A346099 1,2 %A A346099 _Antti Karttunen_, Jul 07 2021 # Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE