login
A366807
a(n) = A020639(A120944(n))*A120944(n).
2
12, 20, 28, 45, 63, 44, 52, 60, 99, 68, 175, 76, 117, 84, 92, 153, 275, 171, 116, 124, 325, 132, 207, 140, 148, 539, 156, 164, 425, 172, 261, 637, 279, 188, 475, 204, 315, 212, 220, 333, 228, 575, 236, 833, 244, 369, 387, 260, 931, 268, 276, 423, 284, 1573, 725
OFFSET
1,1
COMMENTS
Define f(x) to be lpf(k)*k, where lpf(k) = A020639(k). This sequence contains the mappings of f(x) across composite squarefree numbers A120944.
Define sequence R_k = { m : rad(m) | k }, where rad(n) = A007947(n) and k is squarefree. Then the sequence k*R_k contains all numbers divisible by squarefree k that are also not divisible by any prime q coprime to k. It is plain to see that k is the first term in the sequence k*R_k. This sequence gives the second term in k*R_k since lpf(k) is the second term in R_k.
Permutation of A366825. Contains numbers whose prime signature has at least 2 terms, of which is 2, the rest of which are 1s.
Proper subset of A364996, which itself is contained in A126706.
LINKS
FORMULA
a(n) = A065642(A120944(n)), n > 1.
a(n) = A285109(A120944(n)).
EXAMPLE
Let b(n) = A120944(n).
a(1) = 12 = 2^2*3^1 = b(1)*lpf(b(1)) = 6*lpf(6) = 6*2. In {6*A003586}, 12 is the second term.
a(2) = 20 = 2^2*5^1 = b(2)*lpf(b(2)) = 10*lpf(10) = 10*2. In {10*A003592}, 20 is the second term.
a(4) = 45 = 3^2*5^1 = b(4)*lpf(b(4)) = 15*lpf(15) = 15*3. In {15*A003593}, 45 is the second term, etc.
MATHEMATICA
nn = 150; s = Select[Range[nn], And[SquareFreeQ[#], CompositeQ[#]] &];
Array[#*FactorInteger[#][[1, 1]] &[s[[#]]] &, Length[s]]
PROG
(Python)
from math import isqrt
from sympy import primepi, mobius, primefactors
def A366807(n):
def f(x): return n+1+primepi(x)+x-sum(mobius(k)*(x//k**2) for k in range(1, isqrt(x)+1))
m, k = n+1, f(n+1)
while m != k:
m, k = k, f(k)
return m*min(primefactors(m)) # Chai Wah Wu, Aug 02 2024
KEYWORD
nonn,easy
AUTHOR
Michael De Vlieger, Dec 16 2023
STATUS
approved