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.

Showing 1-4 of 4 results.

A323525 Number of ways to arrange the parts of a multiset whose multiplicities are the prime indices of n into a square matrix.

Original entry on oeis.org

1, 1, 0, 0, 0, 0, 1, 0, 6, 4, 0, 12, 0, 0, 0, 24, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 9, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 36, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 84, 0, 0, 72, 0, 0, 0, 0, 0, 0, 0, 0, 126, 252, 0, 0, 0, 0, 0, 0
Offset: 1

Views

Author

Gus Wiseman, Jan 17 2019

Keywords

Comments

This multiset (row n of A305936) is generally not the same as the multiset of prime indices of n. For example, the prime indices of 12 are {1,1,2}, while a multiset whose multiplicities are {1,1,2} is {1,1,2,3}.

Examples

			The a(9) = 6 matrices:
  [1 1] [1 2] [1 2] [2 1] [2 1] [2 2]
  [2 2] [1 2] [2 1] [1 2] [2 1] [1 1]
The a(38) = 9 matrices:
  [1 1 1] [1 1 1] [1 1 1] [1 1 1] [1 1 1] [1 1 1] [1 1 2] [1 2 1] [2 1 1]
  [1 1 1] [1 1 1] [1 1 1] [1 1 2] [1 2 1] [2 1 1] [1 1 1] [1 1 1] [1 1 1]
  [1 1 2] [1 2 1] [2 1 1] [1 1 1] [1 1 1] [1 1 1] [1 1 1] [1 1 1] [1 1 1]
		

Crossrefs

The positions of 0's are numbers whose sum of prime indices is not a perfect square (A323527).
The positions of 1's are primes indexed by squares (A323526).

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,Reverse[primeMS[n]]];
    Table[If[IntegerQ[Sqrt[Total[primeMS[n]]]],Length[Permutations[nrmptn[n]]],0],{n,100}]

Formula

If A056239(n) is a perfect square, a(n) = A318762(n). Otherwise, a(n) = 0.

A323528 Numbers whose sum of prime indices is a perfect square.

Original entry on oeis.org

1, 2, 7, 9, 10, 12, 16, 23, 38, 51, 53, 65, 68, 77, 78, 94, 97, 98, 99, 104, 105, 110, 125, 126, 129, 132, 135, 140, 150, 151, 162, 168, 172, 176, 178, 180, 200, 205, 216, 224, 227, 240, 246, 249, 259, 288, 298, 311, 320, 328, 332, 333, 341, 361, 370, 377, 384
Offset: 1

Views

Author

Gus Wiseman, Jan 17 2019

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798. The sum of prime indices of n is A056239(n).
Also Heinz numbers of integer partitions of perfect squares, where the Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).

Examples

			10 is in the sequence because 10 = 2*5 = prime(1)*prime(3) and 1 + 3 = 4 is a square.
		

Crossrefs

Programs

  • Maple
    select(k-> issqr(add(numtheory[pi](i[1])*i[2], i=ifactors(k)[2])), [$1..400])[]; # Alois P. Heinz, Jan 22 2019
  • Mathematica
    Select[Range[100],IntegerQ[Sqrt[Sum[PrimePi[f[[1]]]*f[[2]],{f,FactorInteger[#]}]]]&]
  • PARI
    isok(n) = {my(f=factor(n)); issquare(sum(k=1, #f~, primepi(f[k, 1])*f[k,2]));} \\ Michel Marcus, Jan 18 2019

A323521 Numbers whose number of prime factors counted with multiplicity (A001222) is not a perfect square.

Original entry on oeis.org

4, 6, 8, 9, 10, 12, 14, 15, 18, 20, 21, 22, 25, 26, 27, 28, 30, 32, 33, 34, 35, 38, 39, 42, 44, 45, 46, 48, 49, 50, 51, 52, 55, 57, 58, 62, 63, 64, 65, 66, 68, 69, 70, 72, 74, 75, 76, 77, 78, 80, 82, 85, 86, 87, 91, 92, 93, 94, 95, 96, 98, 99, 102, 105, 106
Offset: 1

Views

Author

Gus Wiseman, Jan 17 2019

Keywords

Crossrefs

Programs

  • Mathematica
    Select[Range[100],!IntegerQ[Sqrt[PrimeOmega[#]]]&]

A323526 One and prime numbers whose prime index is a perfect square.

Original entry on oeis.org

1, 2, 7, 23, 53, 97, 151, 227, 311, 419, 541, 661, 827, 1009, 1193, 1427, 1619, 1879, 2143, 2437, 2741, 3083, 3461, 3803, 4211, 4637, 5051, 5519, 6007, 6481, 6997, 7573, 8161, 8737, 9341, 9931, 10627, 11321, 12049, 12743, 13499, 14327, 15101, 15877, 16747, 17609
Offset: 1

Views

Author

Gus Wiseman, Jan 17 2019

Keywords

Crossrefs

Programs

  • Mathematica
    Array[Prime[#^2]&,20]
  • PARI
    vector(50, n, if (n==1, 1, prime((n-1)^2))) \\ Michel Marcus, Feb 15 2019

Formula

a(n) = A011757(n-1) for n > 1. - Alois P. Heinz, Jan 17 2019
Showing 1-4 of 4 results.