A290082 - OEIS

The OEIS is supported by the many generous donors to the OEIS Foundation.

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”).

Other ways to Give

A290082

Restricted growth sequence transform of A289625(A003961(n)).

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

(list; graph; refs; listen; history; text; internal format)

OFFSET

1,2

COMMENTS

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

LINKS

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

PROG

(PARI)

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

write_to_bfile(start_offset, vec, bfilename) = { for(n=1, length(vec), write(bfilename, (n+start_offset)-1, " ", vec[n])); }

A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); }; \\ This function from Michel Marcus

A289625(n) = { my(m=1, p=2, v=znstar(n)[2]); for(i=1, length(v), m *= p^v[i]; p = nextprime(p+1)); (m); };

write_to_bfile(1, rgs_transform(vector(16384, n, A289625(A003961(n)))), "b290082_upto16384.txt");

CROSSREFS

Cf. A003961, A003972, A289625, A290076.

Sequence in context: A356991 A083058 A374843 * A127035 A325106 A309687

Adjacent sequences: A290079 A290080 A290081 * A290083 A290084 A290085

KEYWORD

nonn

AUTHOR

Antti Karttunen, Jul 20 2017

STATUS

approved