login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A331183
Number of values of k, 1 <= k <= n, with A323898(k) = A323898(n).
2
1, 1, 1, 2, 1, 1, 1, 3, 2, 2, 1, 1, 1, 1, 1, 4, 1, 1, 1, 2, 1, 2, 1, 1, 1, 3, 1, 1, 1, 1, 1, 5, 1, 3, 1, 1, 1, 4, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 3, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 6, 1, 1, 1, 3, 1, 1, 1, 1, 1, 5, 1, 4, 1, 1, 1, 2, 1, 6, 1, 1, 2, 3, 2, 2, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 3, 1
OFFSET
1,4
COMMENTS
Ordinal transform of A323898, or equally, of the ordered pair [A000120(n), A083254(n)].
PROG
(PARI) up_to = 16384;
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; };
A083254(n) = (2*eulerphi(n)-n);
A323898aux(n) = [hammingweight(n), A083254(n)];
v331183 = ordinal_transform(vector(up_to, n, A323898aux(n)));
A331183(n) = v331183[n];
CROSSREFS
Cf. also A331182, A331184.
Sequence in context: A345272 A373381 A373370 * A284082 A047040 A047020
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jan 12 2020
STATUS
approved