OFFSET

1,2

COMMENTS

Filter sequence for prime factorization patterns, including also information about gaps between prime factors. - Original name, gives the motivation for this sequence. Here the "gaps" refers to differences between the indices of primes present, not the prime gaps as usually understood.

For all i, j:

LINKS

EXAMPLE

a(10) = a(21) (= 6) because both have prime exponents [1, 1] and the difference between the prime indices is the same, as 10 = prime(1)*prime(3), and 21 = prime(2)*prime(4).

a(12) != a(18) because the prime exponents [2,1] and [1,2] do not occur in the same order.

a(140) = a(693) (= 71) because both numbers have prime exponents [2, 1, 1] (in this order) and the differences between the indices of the successive prime factors are same: 140 = prime(1)^2 * prime(3) * prime(4), 693 = prime(2)^2 * prime(4) * prime(5).

PROG

(PARI)

up_to = 65537;

rgs_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), u=1); for(i=1, length(invec), if(mapisdefined(om, invec[i]), my(pp = mapget(om, invec[i])); outvec[i] = outvec[pp] , mapput(om, invec[i], i); outvec[i] = u; u++ )); outvec; };

A246277(n) = if(1==n, 0, my(f = factor(n), k = primepi(f[1, 1])-1); for (i=1, #f~, f[i, 1] = prime(primepi(f[i, 1])-k)); factorback(f)/2);

v305897 = rgs_transform(vector(up_to, n, A246277(n)));

A305897(n) = v305897[n];

CROSSREFS

KEYWORD

nonn

AUTHOR

Antti Karttunen, Jun 14 2018

EXTENSIONS

Name changed by Antti Karttunen, Apr 30 2022

STATUS

approved