A319681 Restricted growth sequence transform of A319680.

1, 2, 3, 4, 5, 6, 6, 7, 8, 9, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 20, 21, 22, 23, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 53, 54, 55, 49, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 75, 76, 77, 78, 79, 80, 81

Offset

1,2

Comments

For all i, j: a(i) = a(j) => A291784(i) = A291784(j).

Links

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

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; };
A019565(n) = {my(j, v); factorback(Mat(vector(if(n, #n=vecextract(binary(n), "-1..1")), j, [prime(j), n[j]])~))}; \\ From A019565
A319680(n) = { my(m=1); fordiv(n, d, if((1==moebius(n/d)), m *= A019565(d))); m; };
v319681 = rgs_transform(vector(up_to, n, A319680(n)));
A319681(n) = v319681[n];

Crossrefs

Cf. A319680, A319682.

Sequence in context: A176842 A278055 A291764 * A340767 A099848 A082287

Adjacent sequences: A319678 A319679 A319680 * A319682 A319683 A319684

Keyword

nonn

Author

Antti Karttunen, Sep 29 2018

Status

approved