A330751 Number of values of k, 1 <= k <= n, with A319626(k) = A319626(n), where A319626(n) gives the numerator of rational valued primorial deflation of n.

1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 2, 3, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 1, 2, 1, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 3

Offset

1,6

Comments

Ordinal transform of A319626.

Links

Antti Karttunen, Table of n, a(n) for n = 1..65537

Prog

(PARI)
up_to = 65537;
ordinal_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), pt); for(i=1, length(invec), if(mapisdefined(om, invec[i]), pt = mapget(om, invec[i]), pt = 0); outvec[i] = (1+pt); mapput(om, invec[i], (1+pt))); outvec; };
A064989(n) = {my(f); 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)};
A319626(n) = (n / gcd(n, A064989(n)));
v330751 = ordinal_transform(vector(up_to, n, A319626(n)));
A330751(n) = v330751[n];

Crossrefs

Cf. A064989, A319626, A330750, A330752.

Sequence in context: A349340 A326297 A060128 * A327407 A290104 A308479

Adjacent sequences: A330748 A330749 A330750 * A330752 A330753 A330754

Keyword

nonn

Author

Antti Karttunen, Dec 30 2019

Status

approved