cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A245500 Concatenation of multiplicities of prime divisors of highly composite numbers A002182(n).

Original entry on oeis.org

0, 1, 2, 11, 21, 31, 22, 41, 211, 311, 221, 411, 321, 421, 3111, 2211, 4111, 3211, 4211, 3311, 5211, 4311, 6211, 4221, 32111, 4411, 5221, 42111, 33111, 52111, 43111, 62111, 42211, 53111, 44111, 52211, 63111, 421111, 331111, 521111, 431111, 621111, 422111
Offset: 1

Views

Author

Michael De Vlieger, Jul 24 2014

Keywords

Comments

For prime decomposition of A002182(n) = 2^a * 3^b * 5^c * ..., a(n) = "abc..." converted to a decimal number.
In other words, each "place" read from left to right represents the n-th prime, starting with 2 at left and increasing to the right. A number in the "place" represents the multiplicity of the corresponding prime in A002182(n).
This notation is corrupt when any multiplicity exceeds 9. The smallest instance of this is at n = 221.
Similar to A054841 but multiplicities are in reverse order.
Given that the exponents e(i) (a,b,c... in the above) of the prime factorization are (weakly) decreasing, their concatenation remains unambiguous way beyond n = 221 (first instance where e(1) >= 10) and even beyond n = 8869 (first instance where e(2) >= 10). Only when e2 >= 10 + e(3) for the first time, in principle the first digit of e(2) could be mistaken for the last digit of e(1); yet it is unlikely if not impossible that e(2) < 10 and e(1) > 100. So the first ambiguous decomposition would require concat(e(1),e(2),e(3)) = concat(a',b',c') with a' > b' >= e(3), thus e(2) significantly larger than 10 + e(3) and e(1) much larger than 100. - M. F. Hasler, Jan 03 2020

Examples

			A002182(4)   = 12 = 2^2 * 3^1, thus a(4) = 21;
A002182(17)  = 2520 = 2^3 * 3^2 * 5^1 * 7^1, thus a(17) = 3211;
A002182(220) = 2^10 * 3^4 * 5^3 * 7^2 * 11 * ... * 53 (skipping no primes), thus a(220) cannot be represented using a single decimal place for the multiplicity 10.
		

Crossrefs

Programs

  • Mathematica
    encodePrimeSignature[n_Integer] :=
      Catch[FromDigits[Reverse[IntegerDigits[Apply[Plus,
           Which[n == 0, Throw["undefined"],
              n == 1, 0,
              Max[Last /@ FactorInteger @ n ] > 9, Throw["overflow"],
              True, Power[10, PrimePi[Abs[#]] - 1]] & /@
           Flatten[ConstantArray @@@ FactorInteger[n]] ]]]]];
    lst = FoldList[Max, 1, Table[DivisorSigma[0, n], {n, 2, 100000000}]];
    Map[encodePrimeSignature, Flatten[Position[lst, #, 1, 1] & /@ Union[lst]]]
  • PARI
    apply( A245500(n)=fromdigits(factor(A002182(n))[,2]~), [1..99]) \\ For n >= 8869, fromdigits must be replaced by (s->if(s,eval(concat([Str(e)|e<-s])))). - M. F. Hasler, Jan 03 2020

Extensions

Offset corrected to 1 by M. F. Hasler, Jan 03 2020