A045966 -id:A045966 - cp's OEIS frontend

A357852 Replace prime(k) with prime(k+2) in the prime factorization of n.

1, 5, 7, 25, 11, 35, 13, 125, 49, 55, 17, 175, 19, 65, 77, 625, 23, 245, 29, 275, 91, 85, 31, 875, 121, 95, 343, 325, 37, 385, 41, 3125, 119, 115, 143, 1225, 43, 145, 133, 1375, 47, 455, 53, 425, 539, 155, 59, 4375, 169, 605, 161, 475, 61, 1715, 187, 1625, 203

Offset: 1

Gus Wiseman, Oct 28 2022

This is the same as A045966 except the first term is 1 instead of 3.

			The terms together with their prime indices begin:
    1: {}
    5: {3}
    7: {4}
   25: {3,3}
   11: {5}
   35: {3,4}
   13: {6}
  125: {3,3,3}
   49: {4,4}
   55: {3,5}
   17: {7}
  175: {3,3,4}
   19: {8}
   65: {3,6}
   77: {4,5}
  625: {3,3,3,3}

Applying the transformation only once gives A003961.
A permutation of A007310.
Other multiplicative sequences: A064988, A064989, A357977, A357980, A357983.
A000040 lists the primes.
A056239 adds up prime indices, row-sums of A112798.
Cf. A000720, A003964, A066207, A076610, A215366, A296150, A299201, A357979.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Product[Prime[i+2],{i,primeMS[n]}],{n,30}]

a(n) = my(f=factor(n)); for (k=1, #f~, f[k,1] = nextprime(nextprime(f[k,1]+1)+1)); factorback(f); \\ Michel Marcus, Oct 28 2022

from math import prod
from sympy import nextprime, factorint
def A357852(n): return prod(nextprime(p,ith=2)**e for p, e in factorint(n).items()) # Chai Wah Wu, Oct 29 2022

a(n) = A003961(A003961(n)).

A045967 a(1)=4; if n = Product p_i^e_i, n > 1, then a(n) = Product p_{i+1}^{e_i+1}.

Original entry on oeis.org

4, 9, 25, 27, 49, 225, 121, 81, 125, 441, 169, 675, 289, 1089, 1225, 243, 361, 1125, 529, 1323, 3025, 1521, 841, 2025, 343, 2601, 625, 3267, 961, 11025, 1369, 729, 4225, 3249, 5929, 3375, 1681, 4761, 7225, 3969, 1849, 27225, 2209, 4563, 6125, 7569, 2809, 6075

Offset: 1

N. J. A. Sloane

If we had a(1) = 1 (instead of 4), then this would be multiplicative and a permutation of the odd powerful numbers (A062739). - Amiram Eldar, Aug 11 2022

From a puzzle proposed by Marc LeBrun.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000

Cf. A045966, A027748, A062739, A124010, A000040.

Cf. A151800.

a045967 1 = 4
a045967 n = product $ zipWith (^)
            (map a151800 $ a027748_row n) (map (+ 1) $ a124010_row n)
-- Reinhard Zumkeller, Jun 03 2013, Dec 23 2011

a[1]=4; a[n_] := Thread[f = FactorInteger[n]; Times @@ Power[f[[All, 1]] // NextPrime , f[[All, 2]] + 1]]; Array[a, 50] (* Jean-François Alcover, Feb 03 2015 *)

Sum_{n>=1} 1/a(n) = 2*zeta(2)*zeta(3)/(3*zeta(6)) - 3/4. - Amiram Eldar, Aug 11 2022

More terms from David W. Wilson

A306697 Square array T(n, k) read by antidiagonals, n > 0 and k > 0: T(n, k) is obtained by applying a Minkowski sum to sets related to the Fermi-Dirac factorizations of n and of k (see Comments for precise definition).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 4, 5, 4, 1, 1, 5, 9, 9, 5, 1, 1, 6, 7, 16, 7, 6, 1, 1, 7, 15, 25, 25, 15, 7, 1, 1, 8, 11, 36, 11, 36, 11, 8, 1, 1, 9, 27, 49, 35, 35, 49, 27, 9, 1, 1, 10, 25, 64, 13, 30, 13, 64, 25, 10, 1, 1, 11, 21, 81, 125, 77, 77, 125, 81

Offset: 1

Rémy Sigrist, Mar 05 2019

For any m > 0:
- let F(m) be the set of distinct Fermi-Dirac primes (A050376) with product m,
- for any i >=0 0 and j >= 0, let f(prime(i+1)^(2^i)) be the lattice point with coordinates X=i and Y=j (where prime(k) denotes the k-th prime number),
- f establishes a bijection from the Fermi-Dirac primes to the lattice points with nonnegative coordinates,
- let P(m) = { f(p) | p in F(m) },
- P establishes a bijection from the nonnegative integers to the set, say L, of finite sets of lattice points with nonnegative coordinates,
- let Q be the inverse of P,
- for any n > 0 and k > 0:
    T(n, k) = Q(P(n) + P(k))
              where "+" denotes the Minkowski addition on L.
This sequence has similarities with A297845, and their data sections almost match; T(6, 6) = 30, however A297845(6, 6) = 90.
This sequence has similarities with A067138; here we work on dimension 2, there in dimension 1.
This sequence as a binary operation distributes over A059896, whereas A297845 distributes over multiplication (A003991) and A329329 distributes over A059897. See the comment in A329329 for further description of the relationship between these sequences. - Peter Munn, Dec 19 2019

			Array T(n, k) begins:
  n\k|  1   2   3    4    5    6    7     8     9    10    11    12
  ---+-------------------------------------------------------------
    1|  1   1   1    1    1    1    1     1     1     1     1     1
    2|  1   2   3    4    5    6    7     8     9    10    11    12
    3|  1   3   5    9    7   15   11    27    25    21    13    45
    4|  1   4   9   16   25   36   49    64    81   100   121   144
    5|  1   5   7   25   11   35   13   125    49    55    17   175
    6|  1   6  15   36   35   30   77   216   225   210   143   540
    7|  1   7  11   49   13   77   17   343   121    91    19   539
    8|  1   8  27   64  125  216  343   128   729  1000  1331  1728
    9|  1   9  25   81   49  225  121   729   625   441   169  2025
   10|  1  10  21  100   55  210   91  1000   441   110   187  2100
   11|  1  11  13  121   17  143   19  1331   169   187    23  1573
   12|  1  12  45  144  175  540  539  1728  2025  2100  1573   720

Rémy Sigrist, PARI program for A306697
OEIS Wiki, "Fermi-Dirac representation" of n
Eric Weisstein's World of Mathematics, Distributive
Wikipedia, Minkowski addition

Columns (some differing for term 1) and equivalently rows: A003961(3), A000290(4), A045966(5), A045968(7), A045970(11).
Cf. A000040, A001221, A019565, A064547, A225546, A329050.
Related binary operations: A067138, A059896, A297845/A003991, A329329/A059897.

PARI
```
\\ See Links section.
```

For any m > 0, n > 0, k > 0, i >= 0, j >= 0:
- T(n, k) = T(k, n) (T is commutative),
- T(m, T(n, k)) = T(T(m, n), k) (T is associative),
- T(n, 1) = 1 (1 is an absorbing element for T),
- T(n, 2) = n (2 is an identity element for T),
- T(n, 3) = A003961(n),
- T(n, 4) = n^2 (A000290),
- T(n, 5) = A357852(n),
- T(n, 7) = A045968(n) (when n > 1),
- T(n, 11) = A045970(n) (when n > 1),
- T(n, 2^(2^i)) = n^(2^i),
- T(2^i, 2^j) = 2^A067138(i, j),
- T(A019565(i), A019565(j)) = A019565(A067138(i, j)),
- T(A000040(n), A000040(k)) = A000040(n + k - 1),
- T(2^(2^i), 2^(2^j)) = 2^(2^(i + j)),
- A001221(T(n, k)) <= A001221(n) * A001221(k),
- A064547(T(n, k)) <= A064547(n) * A064547(k).
From Peter Munn, Dec 05 2019:(Start)
T(A329050(i_1, j_1), A329050(i_2, j_2)) = A329050(i_1+i_2, j_1+j_2).
Equivalently, T(prime(i_1 - 1)^(2^(j_1)), prime(i_2 - 1)^(2^(j_2))) = prime(i_1+i_2 - 1)^(2^(j_1+j_2)), where prime(i) = A000040(i).
T(A059896(i,j), k) = A059896(T(i,k), T(j,k)) (T distributes over A059896).
T(A019565(i), 2^j) = A019565(i)^j.
T(A225546(i), A225546(j)) = A225546(T(i,j)).
(End)

A045968 a(1)=5; for n >= 2, if n = Product p_i^e_i, then a(n) = Product p_{i+3}^e_i.

Original entry on oeis.org

5, 7, 11, 49, 13, 77, 17, 343, 121, 91, 19, 539, 23, 119, 143, 2401, 29, 847, 31, 637, 187, 133, 37, 3773, 169, 161, 1331, 833, 41, 1001, 43, 16807, 209, 203, 221, 5929, 47, 217, 253, 4459, 53, 1309, 59, 931, 1573, 259, 61, 26411, 289, 1183, 319, 1127, 67, 9317, 247

Offset: 1

N. J. A. Sloane

			If n = 9 = 3^2, then a(n) = 11^2 = 121 (since 11 is the third prime after 3).

From a puzzle proposed by Marc LeBrun.

Cf. A045965, A045966, A045967, A045969, A045970, A045971, A045972, A045973.

f[p_, e_] := NextPrime[p, 3]^e; a[1] = 5; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 19 2023 *)

More terms from David W. Wilson

Erroneous linear recurrence deleted by Harvey P. Dale, May 07 2018

A045970 a(1)=7; if n = Product p_i^e_i, n > 1, then a(n) = Product p_{i+4}^e_i.

Original entry on oeis.org

7, 11, 13, 121, 17, 143, 19, 1331, 169, 187, 23, 1573, 29, 209, 221, 14641, 31, 1859, 37, 2057, 247, 253, 41, 17303, 289, 319, 2197, 2299, 43, 2431, 47, 161051, 299, 341, 323, 20449, 53, 407, 377, 22627, 59, 2717, 61, 2783, 2873, 451, 67, 190333, 361, 3179, 403

Offset: 1

N. J. A. Sloane

From a puzzle proposed by Marc LeBrun.

Cf. A045965, A045966, A045967, A045968, A045969, A045971, A045972, A045973.

f[p_, e_] := NextPrime[p, 4]^e; a[1] = 7; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 19 2023 *)

More terms from David W. Wilson

A045973 a(1)=10; if n = Product p_i^e_i, n > 1, then a(n) = Product p_{i+1}^e_i * Product p_{i+3}^e_i.

Original entry on oeis.org

10, 21, 55, 441, 91, 1155, 187, 9261, 3025, 1911, 247, 24255, 391, 3927, 5005, 194481, 551, 63525, 713, 40131, 10285, 5187, 1073, 509355, 8281, 8211, 166375, 82467, 1271, 105105, 1591, 4084101, 13585, 11571, 17017, 1334025, 1927, 14973, 21505, 842751

Offset: 1

N. J. A. Sloane

From a puzzle proposed by Marc LeBrun.

Cf. A045965, A045966, A045967, A045968, A045969, A045970, A045971, A045972.

f[p_, e_] := NextPrime[p]^e * NextPrime[p, 3]^e; a[1] = 10; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 19 2023 *)

Sum_{n>=1} 1/a(n) = -9/10 + Product_{k>=1} (1+1/(prime(k)*prime(k+4)-1)) = 0.2602421684... . - Amiram Eldar, Sep 19 2023

More terms from David W. Wilson

A045969 a(1)=6; if n = Product p_i^e_i, n>1, then a(n) = Product p_{i+1}^e_i * Product p_{i+2}^e_i.

Original entry on oeis.org

6, 15, 35, 225, 77, 525, 143, 3375, 1225, 1155, 221, 7875, 323, 2145, 2695, 50625, 437, 18375, 667, 17325, 5005, 3315, 899, 118125, 5929, 4845, 42875, 32175, 1147, 40425, 1517, 759375, 7735, 6555, 11011, 275625, 1763, 10005, 11305, 259875, 2021, 75075

Offset: 1

N. J. A. Sloane

From a puzzle proposed by Marc LeBrun.

Cf. A045965, A045966, A045967, A045968, A045970, A045971, A045972, A045973.

f[p_, e_] := (NextPrime[p] * NextPrime[p, 2])^e; a[1] = 6; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 19 2023 *)

Sum_{n>=1} 1/a(n) = -5/6 + Product_{k>=2} (1+1/(prime(k)*prime(k+1)-1)) = 0.31383788... . - Amiram Eldar, Sep 19 2023

More terms from David W. Wilson

A045971 a(1)=8; if n = Product p_i^e_i, n > 1, then a(n) = Product p_{i+1}^{e_i+2}.

Original entry on oeis.org

8, 27, 125, 81, 343, 3375, 1331, 243, 625, 9261, 2197, 10125, 4913, 35937, 42875, 729, 6859, 16875, 12167, 27783, 166375, 59319, 24389, 30375, 2401, 132651, 3125, 107811, 29791, 1157625, 50653, 2187, 274625, 185193, 456533, 50625, 68921, 328509, 614125

Offset: 1

N. J. A. Sloane

From a puzzle proposed by Marc LeBrun.

Cf. A045965, A045966, A045967, A045968, A045969, A045970, A045972, A045973, A065483.

f[p_, e_] := NextPrime[p]^(e + 2); a[1] = 8; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 19 2023 *)

Sum_{n>=1} 1/a(n) = (4/5) * A065483 - 7/8 = 0.196827322859... . - Amiram Eldar, Sep 19 2023

More terms from David W. Wilson

A045972 a(1)=9; if n = Product p_i^e_i, n > 1, then a(n) = Product p_{i+2}^{e_i+1}.

Original entry on oeis.org

9, 25, 49, 125, 121, 1225, 169, 625, 343, 3025, 289, 6125, 361, 4225, 5929, 3125, 529, 8575, 841, 15125, 8281, 7225, 961, 30625, 1331, 9025, 2401, 21125, 1369, 148225, 1681, 15625, 14161, 13225, 20449, 42875, 1849, 21025, 17689, 75625, 2209, 207025

Offset: 1

N. J. A. Sloane

From a puzzle proposed by Marc LeBrun.

Cf. A045965, A045966, A045967, A045968, A045969, A045970, A045971, A045973, A082695.

f[p_, e_] := NextPrime[p, 2]^(e + 1); a[1] = 9; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 19 2023 *)

Sum_{n>=1} 1/a(n) = (4/7) * (zeta(2)*zeta(3)/zeta(6)) - 8/9 = 0.221737646437... . - Amiram Eldar, Sep 19 2023

More terms from David W. Wilson

A357981 Numbers whose prime indices have only prime numbers as their own prime indices.

Original entry on oeis.org

1, 2, 4, 5, 8, 10, 11, 16, 20, 22, 23, 25, 31, 32, 40, 44, 46, 47, 50, 55, 59, 62, 64, 80, 88, 92, 94, 97, 100, 103, 110, 115, 118, 121, 124, 125, 127, 128, 137, 155, 160, 176, 179, 184, 188, 194, 197, 200, 206, 220, 230, 233, 235, 236, 242, 248, 250, 253, 254

Offset: 1

Gus Wiseman, Oct 23 2022

nonn

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.

Numbers whose prime indices are prime numbers are listed by A076610.

			The terms together with their prime indices begin:
     1: {}
     2: {1}
     4: {1,1}
     5: {3}
     8: {1,1,1}
    10: {1,3}
    11: {5}
    16: {1,1,1,1}
    20: {1,1,3}
    22: {1,5}
    23: {9}
    25: {3,3}
    31: {11}
    32: {1,1,1,1,1}

Contains all elements of A000079.
Contains all primes indexed by elements of A076610.
A000040 lists the prime numbers.
A056239 adds up prime indices, row-sums of A112798.
Cf. A003961, A045966, A064988, A066207, A215366, A357977, A357980, A357983.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[100],And@@PrimeQ/@Join@@primeMS/@primeMS[#]&]

cp's OEIS Frontend

A357852 Replace prime(k) with prime(k+2) in the prime factorization of n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A045967 a(1)=4; if n = Product p_i^e_i, n > 1, then a(n) = Product p_{i+1}^{e_i+1}.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Mathematica

Formula

Extensions

A306697 Square array T(n, k) read by antidiagonals, n > 0 and k > 0: T(n, k) is obtained by applying a Minkowski sum to sets related to the Fermi-Dirac factorizations of n and of k (see Comments for precise definition).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

Formula

A045968 a(1)=5; for n >= 2, if n = Product p_i^e_i, then a(n) = Product p_{i+3}^e_i.

Views

Author

Keywords

Examples

References

Crossrefs

Programs

Mathematica

Extensions

A045970 a(1)=7; if n = Product p_i^e_i, n > 1, then a(n) = Product p_{i+4}^e_i.

Views

Author

Keywords

References

Crossrefs

Programs

Mathematica

Extensions

A045973 a(1)=10; if n = Product p_i^e_i, n > 1, then a(n) = Product p_{i+1}^e_i * Product p_{i+3}^e_i.

Views

Author

Keywords

References

Crossrefs

Programs

Mathematica

Formula

Extensions

A045969 a(1)=6; if n = Product p_i^e_i, n>1, then a(n) = Product p_{i+1}^e_i * Product p_{i+2}^e_i.

Views

Author

Keywords

References

Crossrefs

Programs

Mathematica

Formula

Extensions

A045971 a(1)=8; if n = Product p_i^e_i, n > 1, then a(n) = Product p_{i+1}^{e_i+2}.

Views

Author

Keywords