A119628 Inverse permutation to sequence A083872.

1, 2, 3, 5, 9, 4, 31, 6, 17, 10, 271, 7, 793, 32, 11, 18, 5377, 19, 14689, 12, 33, 272, 96001, 8, 161, 794, 97, 34, 917281, 13, 2332801, 61, 273, 5378, 35, 20, 19595521, 14690, 795, 14, 102435841, 36, 258048001, 274, 21, 96002, 1258709761, 22, 1585, 162, 5379

Offset

1,2

Links

Rémy Sigrist, Table of n, a(n) for n = 1..10000

Formula

a(p) = A027423(p-1) + 1 for any prime number p. - Rémy Sigrist, Nov 08 2018

Prog

(PARI) {m=31; v=[]; for(n=1, m, v=concat(v, vecsort(eval(setminus(Set(divisors(n!)), Set(v)))))); z=100; u=vector(z); for(j=1, #v, if(v[j]<=z, u[v[j]]=j)); i=1; while(u[i]>0, print1(u[i], ", "); i++)} \\ Klaus Brockhaus, Aug 16 2006
(PARI) a(n) = my (ff=1, f=1); for (i=1, oo, f*=i; if (f%n==0, return (numdiv(ff) + sum(j=1, n, ff%j!=0 && f%j == 0))); ff=f) \\ Rémy Sigrist, Nov 08 2018

Crossrefs

Cf. A027423, A083872.

Sequence in context: A232562 A064358 A109736 * A021093 A342444 A011026

Adjacent sequences: A119625 A119626 A119627 * A119629 A119630 A119631

Keyword

nonn

Author

Leroy Quet, Jun 08 2006

Extensions

a(11) to a(36) from Klaus Brockhaus, Aug 16 2006

More terms from Rémy Sigrist, Nov 08 2018

Status

approved