A277886 -id:A277886 - cp's OEIS frontend

A097248 a(n) is the eventual stable point reached when iterating k -> A097246(k), starting from k = n.

1, 2, 3, 3, 5, 6, 7, 6, 5, 10, 11, 5, 13, 14, 15, 5, 17, 10, 19, 15, 21, 22, 23, 10, 7, 26, 15, 21, 29, 30, 31, 10, 33, 34, 35, 15, 37, 38, 39, 30, 41, 42, 43, 33, 7, 46, 47, 15, 11, 14, 51, 39, 53, 30, 55, 42, 57, 58, 59, 7, 61, 62, 35, 15, 65, 66, 67, 51, 69, 70, 71, 30, 73, 74, 21

Offset: 1

Reinhard Zumkeller, Aug 03 2004

nonn

a(n) = r(n,m) with m such that r(n,m)=r(n,m+1), where r(n,k) = A097246(r(n,k-1)), r(n,0)=n. (The original definition.)
A097248(n) = r(n,a(n)).
From Antti Karttunen, Nov 15 2016: (Start)
The above remark could be interpreted to mean that A097249(n) <= a(n).
All terms are squarefree, and the squarefree numbers are the fixed points.
These are also fixed points eventually reached when iterating A277886.
(End)

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

Range of values is A005117.
Cf. A008683, A019565, A048675, A097246, A097247, A097249, A277905, A332824.
A003961, A225546, A277885, A277886, A331590 are used to express relationship between terms of this sequence.
The formula section also details how the sequence maps the terms of A007913, A260443, A329050, A329332.
See comments/formulas in A283475, A283478, A331751 giving their relationship to this sequence.

Table[FixedPoint[Times @@ Map[#1^#2 & @@ # &, Partition[#, 2, 2] &@ Flatten[FactorInteger[#] /. {p_, e_} /; e >= 2 :> {If[OddQ@ e, {p, 1}, {1, 1}], {NextPrime@ p, Floor[e/2]}}]] &, n], {n, 75}] (* Michael De Vlieger, Mar 18 2017 *)

A097246(n) = { my(f=factor(n)); prod(i=1, #f~, (nextprime(f[i,1]+1)^(f[i,2]\2))*((f[i,1])^(f[i,2]%2))); };
A097248(n) = { my(k=A097246(n)); while(k<>n, n = k; k = A097246(k)); k; };
\\ Antti Karttunen, Mar 18 2017

from sympy import factorint, nextprime
from operator import mul
def a097246(n):
    f=factorint(n)
    return 1 if n==1 else reduce(mul, [(nextprime(i)**int(f[i]/2))*(i**(f[i]%2)) for i in f])
def a(n):
    k=a097246(n)
    while k!=n:
        n=k
        k=a097246(k)
    return k # Indranil Ghosh, May 15 2017

;; with memoization-macro definec
;; Two implementations:
(definec (A097248 n) (if (not (zero? (A008683 n))) n (A097248 (A097246 n))))
(definec (A097248 n) (if (zero? (A277885 n)) n (A097248 (A277886 n))))
;; Antti Karttunen, Nov 15 2016

a(A005117(n)) = A005117(n).
From Antti Karttunen, Nov 15 2016: (Start)
If A008683(n) <> 0 [when n is squarefree], a(n) = n, otherwise a(n) = a(A097246(n)).
If A277885(n) = 0, a(n) = n, otherwise a(n) = a(A277886(n)).
A007913(a(n)) = a(n).
a(A007913(n)) = A007913(n).
A048675(a(n)) = A048675(n).
a(A260443(n)) = A019565(n).
(End)
From Peter Munn, Feb 06 2020: (Start)
a(1) = 1; a(p) = p, for prime p; a(m*k) = A331590(a(m), a(k)).
a(A331590(m,k)) = A331590(a(m), a(k)).
a(n^2) = a(A003961(n)) = A003961(a(n)).
a(A225546(n)) = a(n).
a(n) = A225546(2^A048675(n)) = A019565(A048675(n)).
a(A329050(n,k)) = prime(n+k-1) = A000040(n+k-1).
a(A329332(n,k)) = A019565(n * k).
Equivalently, a(A019565(n)^k) = A019565(n * k).
(End)
From Antti Karttunen, Feb 22-25 & Mar 01 2020: (Start)
a(A019565(x)*A019565(y)) = A019565(x+y).
a(A332461(n)) = A332462(n).
a(A332824(n)) = A019565(n).
a(A277905(n,k)) = A277905(n,1) = A019565(n), for all n >= 1, and 1 <= k <= A018819(n).
(End)

Name changed and the original definition moved to the Comments section by Antti Karttunen, Nov 15 2016

A277905 Irregular table: Each row n (n >= 0) lists in ascending order all A018819(n) numbers k for which A048675(k) = n.

Original entry on oeis.org

1, 2, 3, 4, 6, 8, 5, 9, 12, 16, 10, 18, 24, 32, 15, 20, 27, 36, 48, 64, 30, 40, 54, 72, 96, 128, 7, 25, 45, 60, 80, 81, 108, 144, 192, 256, 14, 50, 90, 120, 160, 162, 216, 288, 384, 512, 21, 28, 75, 100, 135, 180, 240, 243, 320, 324, 432, 576, 768, 1024, 42, 56, 150, 200, 270, 360, 480, 486, 640, 648, 864, 1152, 1536, 2048, 35, 63, 84, 112, 125, 225, 300, 400

Offset: 1

Antti Karttunen, Nov 14 2016

Each row beginning with an odd number (rows with even index) is followed by a row of the same length, with the same terms, but multiplied by 2. See also comments in the Formula section of A018819.
Note that although the indexing of rows start from zero, the indexing of this sequence starts from 1, with a(1) = 1.
Also Heinz numbers of integer partitions whose binary rank is n, where the binary rank of a partition y is given by Sum_i 2^(y_i-1). For example, row n = 6 is 15, 20, 27, 36, 48, 64, corresponding to the partitions (3,2), (3,1,1), (2,2,2), (2,2,1,1), (2,1,1,1,1), (1,1,1,1,1,1). - Gus Wiseman, May 25 2024
Also, row n lists in ascending order all A018819(n) numbers k for which A097248(k) = A019565(n). - Flávio V. Fernandes, Jul 19 2025

			The irregular table begins as:
  row terms
   0   1;
   1   2;
   2   3,  4;
   3   6,  8;
   4   5,  9,  12,  16;
   5  10, 18,  24,  32;
   6  15, 20,  27,  36,  48,  64;
   7  30, 40,  54,  72,  96, 128;
   8   7, 25,  45,  60,  80,  81, 108, 144, 192, 256;
   9  14, 50,  90, 120, 160, 162, 216, 288, 384, 512;
  10  21, 28,  75, 100, 135, 180, 240, 243, 320, 324, 432,  576,  768, 1024;
  11  42, 56, 150, 200, 270, 360, 480, 486, 640, 648, 864, 1152, 1536, 2048;
...

Antti Karttunen, Table of n, a(n) for n = 1..166
Index entries for sequences that are permutations of the natural numbers

Cf. A019565 (the left edge, the only terms that are squarefree).
Cf. A000079 (the trailing edge).
Cf. A048675, A260443, A277886, A277896, A277903, A277904.
Row lengths are A018819 (number of partitions of binary rank n).
A000009 counts strict partitions, ranks A005117.
A029837 stc_sum or A070939 bin_len, opposite A070940 binexp_lastpos_1.
A048675 gives binary rank of prime indices, distinct A087207.
A048793 lists binary indices, product A096111, reverse A272020.
A061395 gives greatest prime index, least A055396.
A112798 lists prime indices, cf. A001222, A003963, A056239, A296150.
A372890 adds up binary ranks of partitions, strict A372888.
Cf. A000040, A000041, A001511, A005940, A118462, A215366, A225620, A246868.
Cf. A097248, A048675.

prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n], {p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Select[Range[0,2^k],Total[2^(prix[#]-1)]==k&],{k,0,10}] (* Gus Wiseman, May 25 2024 *)

(definec (A277905 n) (A277905bi (A277903 n) (A277904 n)))
(define (A277905bi row col) (let outloop ((k (A019565 row)) (col col)) (if (zero? col) k (let inloop ((j (+ 1 k))) (if (= (A048675 j) row) (outloop j (- col 1)) (inloop (+ 1 j))))))) ;; Very slow implementation.
;; Implementation based on a naive recurrence:
(definec (A277905 n) (if (= 1 n) n (let ((maybe_next (A277896 (A277905 (- n 1))))) (if (not (zero? maybe_next)) maybe_next (A019565 (A277903 n))))))

a(1) = 1; for n > 1, if A277896(a(n-1)) > 0, then a(n) = A277896(a(n-1)), otherwise a(n) = A019565(A277903(n)). [A naive recurrence for a one-dimensional version.]
Other identities. For all n >= 1:
A048675(a(n)) = A277903(n).

A097246 Replace factors of n that are squares of a prime with the prime succeeding this prime.

Original entry on oeis.org

1, 2, 3, 3, 5, 6, 7, 6, 5, 10, 11, 9, 13, 14, 15, 9, 17, 10, 19, 15, 21, 22, 23, 18, 7, 26, 15, 21, 29, 30, 31, 18, 33, 34, 35, 15, 37, 38, 39, 30, 41, 42, 43, 33, 25, 46, 47, 27, 11, 14, 51, 39, 53, 30, 55, 42, 57, 58, 59, 45, 61, 62, 35, 27, 65, 66, 67, 51, 69, 70, 71, 30, 73

Offset: 1

Reinhard Zumkeller, Aug 03 2004

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

Cf. A000035, A000040, A002110, A049084, A097250, A000079, A000188, A000244, A003961, A004526, A005117, A007814, A007913, A048675, A064989, A151800.

Cf. A097247, A097248 (fixed points of iteration), A097249 (number of iterations needed to reach them for each n), A277886, A277899.

Table[Times @@ Map[#1^#2 & @@ # &, Partition[#, 2, 2] &@ Flatten[ FactorInteger[n] /. {p_, e_} /; e >= 2 :> {If[OddQ@ e, {p, 1}, {1, 1}], {NextPrime@ p, Floor[e/2]}}]], {n, 73}] (* Michael De Vlieger, Mar 18 2017 *)

A097246(n) = { my(f=factor(n)); prod(i=1, #f~, (nextprime(f[i, 1]+1)^(f[i,2]\2))*((f[i,1])^(f[i,2]%2))); }; \\ Antti Karttunen, Mar 18 2017

from sympy import factorint, nextprime
from operator import mul
def a(n):
    f=factorint(n)
    return 1 if n==1 else reduce(mul, [(nextprime(i)**int(f[i]/2))*(i**(f[i]%2)) for i in f]) # Indranil Ghosh, May 15 2017

(definec (A097246 n) (if (= 1 n) 1 (* (A000244 (A004526 (A007814 n))) (A000079 (A000035 (A007814 n))) (A003961 (A097246 (A064989 n))))))
(define (A097246 n) (* (A003961 (A000188 n)) (A007913 n)))
;; Antti Karttunen, Nov 15 2016

Multiplicative with p^e -> NextPrime(p)^floor(e/2) * p^(e mod 2), where p prime and NextPrime(p)=A000040(A049084(p)+1).
a(n) <= n; a(n) = n iff n is squarefree: a(A005117(n)) = A005117(n);
a(m*n) <= a(m)*a(n); a(m*n) = a(m)*a(n) iff m and n are coprime;
a(A000040(k)^n) = A000040(k+1)^floor(n/2)*A000040(k)^(n mod 2); a(2^n) = 3^floor(n/2) * (1 + n mod 2);
a(A000040(k)*A002110(n)/A002110(k-1)) = A000040(k+1)*A002110(n)/A002110(k) for k <= n, see also A097250.
From Antti Karttunen, Nov 15 2016: (Start)
a(1) = 1; for n > 1, a(n) = 2^A000035(A007814(n)) * 3^A004526(A007814(n)) * A003961(a(A064989(n))).
a(n) = A003961(A000188(n)) * A007913(n).
A048675(a(n)) = A048675(n).
(End)
Sum_{k=1..n} a(k) ~ c * n^2, where c = (1/2) * Product_{p prime} (p^4-p^2)/(p^4-nextprime(p)) = 0.4059779303..., where nextprime is A151800. - Amiram Eldar, Nov 29 2022

A277896 a(n) = the least k > n for which A048675(k) = A048675(n), 0 if no such number exists (when n is a power of 2).

Original entry on oeis.org

0, 0, 4, 0, 9, 8, 25, 0, 12, 18, 49, 16, 121, 50, 20, 0, 169, 24, 289, 27, 28, 98, 361, 32, 45, 242, 36, 75, 529, 40, 841, 0, 44, 338, 63, 48, 961, 578, 52, 54, 1369, 56, 1681, 147, 60, 722, 1849, 64, 175, 90, 68, 363, 2209, 72, 99, 150, 76, 1058, 2809, 80, 3481, 1682, 84, 0, 117, 88, 3721, 507, 92, 126, 4489, 96, 5041, 1922, 100, 867, 275

Offset: 1

Antti Karttunen, Nov 15 2016

nonn

Apart from zeros, a permutation of A013929.

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

Cf. A000079, A013929, A048675, A209229, A277886, A277893, A277905.

Numbers not in this sequence: A005117 (A019565).

(define (A277896 n) (if (= 1 (A209229 n)) 0 (let ((v (A048675 n))) (let loop ((k (+ 1 n))) (if (= (A048675 k) v) k (loop (+ 1 k)))))))

a(A000079(n)) = 0.

For all n, except powers of two, a(n) >= A277893(n).

cp's OEIS Frontend

A097248 a(n) is the eventual stable point reached when iterating k -> A097246(k), starting from k = n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Python

Scheme

Formula

Extensions

A277905 Irregular table: Each row n (n >= 0) lists in ascending order all A018819(n) numbers k for which A048675(k) = n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Scheme

Formula

A097246 Replace factors of n that are squares of a prime with the prime succeeding this prime.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Python

Scheme

Formula

A277896 a(n) = the least k > n for which A048675(k) = A048675(n), 0 if no such number exists (when n is a power of 2).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Scheme

Formula