A265750 -id:A265750 - cp's OEIS frontend

A003961 Completely multiplicative with a(prime(k)) = prime(k+1).

1, 3, 5, 9, 7, 15, 11, 27, 25, 21, 13, 45, 17, 33, 35, 81, 19, 75, 23, 63, 55, 39, 29, 135, 49, 51, 125, 99, 31, 105, 37, 243, 65, 57, 77, 225, 41, 69, 85, 189, 43, 165, 47, 117, 175, 87, 53, 405, 121, 147, 95, 153, 59, 375, 91, 297, 115, 93, 61, 315, 67, 111, 275, 729, 119

Offset: 1

Marc LeBrun

Meyers (see Guy reference) conjectures that for all r >= 1, the least odd number not in the set {a(i): i < prime(r)} is prime(r+1). - N. J. A. Sloane, Jan 08 2021
Meyers' conjecture would be refuted if and only if for some r there were such a large gap between prime(r) and prime(r+1) that there existed a composite c for which prime(r) < c < a(c) < prime(r+1), in which case (by Bertrand's postulate) c would necessarily be a term of A246281. - Antti Karttunen, Mar 29 2021
a(n) is odd for all n and for each odd m there exists a k with a(k) = m (see A064216). a(n) > n for n > 1: bijection between the odd and all numbers. - Reinhard Zumkeller, Sep 26 2001
a(n) and n have the same number of distinct primes with (A001222) and without multiplicity (A001221). - Michel Marcus, Jun 13 2014
From Antti Karttunen, Nov 01 2019: (Start)
More generally, a(n) has the same prime signature as n, A046523(a(n)) = A046523(n). Also A246277(a(n)) = A246277(n) and A287170(a(n)) = A287170(n).
Many permutations and other sequences that employ prime factorization of n to encode either polynomials, partitions (via Heinz numbers) or multisets in general can be easily defined by using this sequence as one of their constituent functions. See the last line in the Crossrefs section for examples.
(End)

			a(12) = a(2^2 * 3) = a(prime(1)^2 * prime(2)) = prime(2)^2 * prime(3) = 3^2 * 5 = 45.
a(A002110(n)) = A002110(n + 1) / 2.

Richard K. Guy, editor, Problems From Western Number Theory Conferences, Labor Day, 1983, Problem 367 (Proposed by Leroy F. Meyers, The Ohio State U.).

Indranil Ghosh, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
Index entries for sequences computed from indices in prime factorization.
Index entries for sequences related to Heinz numbers.

See A045965 for another version.
Row 1 of table A242378 (which gives the "k-th powers" of this sequence), row 3 of A297845 and of A306697. See also arrays A066117, A246278, A255483, A308503, A329050.
Cf. A064989 (a left inverse), A064216, A000040, A002110, A000265, A027746, A046523, A048673 (= (a(n)+1)/2), A108228 (= (a(n)-1)/2), A191002 (= a(n)*n), A252748 (= a(n)-2n), A286385 (= a(n)-sigma(n)), A283980 (= a(n)*A006519(n)), A341529 (= a(n)*sigma(n)), A326042, A049084, A001221, A001222, A122111, A225546, A260443, A245606, A244319, A246269 (= A065338(a(n))), A322361 (= gcd(n, a(n))), A305293.
Cf. A191555, A252738.
Cf. A249734, A249735 (bisections).
Cf. A246261 (a(n) is of the form 4k+1), A246263 (of the form 4k+3), A246271, A246272, A246259, A246281 (n such that a(n) < 2n), A246282 (n such that a(n) > 2n), A252742.
Cf. A275717 (a(n) > a(n-1)), A275718 (a(n) < a(n-1)).
Cf. A003972 (Möbius transform), A003973 (Inverse Möbius transform), A318321.
Cf. A300841, A305421, A322991, A250469, A269379 for analogous shift-operators in other factorization and quasi-factorization systems.
Cf. also following permutations and other sequences that can be defined with the help of this sequence: A005940, A163511, A122111, A260443, A206296, A265408, A265750, A275733, A275735, A297845, A091202 & A091203, A250245 & A250246, A302023 & A302024, A302025 & A302026.
A version for partition numbers is A003964, strict A357853.
A permutation of A005408.
Applying the same transformation again gives A357852.
Other multiplicative sequences: A064988, A357977, A357978, A357980, A357983.
A056239 adds up prime indices, row-sums of A112798.
Cf. A000720, A076610, A296150.

a003961 1 = 1
a003961 n = product $ map (a000040 . (+ 1) . a049084) $ a027746_row n
-- Reinhard Zumkeller, Apr 09 2012, Oct 09 2011
(MIT/GNU Scheme, with Aubrey Jaffer's SLIB Scheme library)
(require 'factor)
(define (A003961 n) (apply * (map A000040 (map 1+ (map A049084 (factor n))))))
;; Antti Karttunen, May 20 2014

a:= n-> mul(nextprime(i[1])^i[2], i=ifactors(n)[2]):
seq(a(n), n=1..80);  # Alois P. Heinz, Sep 13 2017

a[p_?PrimeQ] := a[p] = Prime[ PrimePi[p] + 1]; a[1] = 1; a[n_] := a[n] = Times @@ (a[#1]^#2& @@@ FactorInteger[n]); Table[a[n], {n, 1, 65}] (* Jean-François Alcover, Dec 01 2011, updated Sep 20 2019 *)
Table[Times @@ Map[#1^#2 & @@ # &, FactorInteger[n] /. {p_, e_} /; e > 0 :> {Prime[PrimePi@ p + 1], e}] - Boole[n == 1], {n, 65}] (* Michael De Vlieger, Mar 24 2017 *)

a(n)=local(f); if(n<1,0,f=factor(n); prod(k=1,matsize(f)[1],nextprime(1+f[k,1])^f[k,2]))

a(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ Michel Marcus, May 17 2014

use ntheory ":all";  sub a003961 { vecprod(map { next_prime($) } factor(shift)); }  # _Dana Jacobsen, Mar 06 2016

from sympy import factorint, prime, primepi, prod
def a(n):
    f=factorint(n)
    return 1 if n==1 else prod(prime(primepi(i) + 1)**f[i] for i in f)
[a(n) for n in range(1, 11)] # Indranil Ghosh, May 13 2017

If n = Product p(k)^e(k) then a(n) = Product p(k+1)^e(k).
Multiplicative with a(p^e) = A000040(A000720(p)+1)^e. - David W. Wilson, Aug 01 2001
a(n) = Product_{k=1..A001221(n)} A000040(A049084(A027748(n,k))+1)^A124010(n,k). - Reinhard Zumkeller, Oct 09 2011 [Corrected by Peter Munn, Nov 11 2019]
A064989(a(n)) = n and a(A064989(n)) = A000265(n). - Antti Karttunen, May 20 2014 & Nov 01 2019
A001221(a(n)) = A001221(n) and A001222(a(n)) = A001222(n). - Michel Marcus, Jun 13 2014
From Peter Munn, Oct 31 2019: (Start)
a(n) = A225546((A225546(n))^2).
a(A225546(n)) = A225546(n^2).
(End)
Sum_{k=1..n} a(k) ~ c * n^2, where c = (1/2) * Product_{p prime} ((p^2-p)/(p^2-nextprime(p))) = 2.06399637... . - Amiram Eldar, Nov 18 2022

A265752 a(n) = A007814(A265399(n)).

Original entry on oeis.org

0, 1, 0, 2, 1, 1, 1, 3, 0, 2, 2, 2, 3, 2, 1, 4, 5, 1, 8, 3, 1, 3, 13, 3, 2, 4, 0, 3, 21, 2, 34, 5, 2, 6, 2, 2, 55, 9, 3, 4, 89, 2, 144, 4, 1, 14, 233, 4, 2, 3, 5, 5, 377, 1, 3, 4, 8, 22, 610, 3, 987, 35, 1, 6, 4, 3, 1597, 7, 13, 3, 2584, 3, 4181, 56, 2, 10, 3, 4, 6765, 5, 0, 90, 10946, 3, 6, 145, 21, 5, 17711

Offset: 1

Antti Karttunen, Dec 15 2015

nonn

a(n) is the constant term of the reduction by x^2->x+1 of the polynomial encoded in the prime factorization of n. (Assuming here only polynomials with nonnegative integer coefficients, see e.g. A206296 for the details of the encoding).

Completely additive with a(prime(k)) = F(k-2), where F(k) denotes the k-th Fibonacci number, A000045(k) for k >= 0, or A039834(-k) for k <= 0. - Peter Munn, Apr 05 2021, incorporating comment by Antti Karttunen, Dec 15 2015

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

Cf. A007814, A192232, A192750, A206296, A265399, A265750, A265753.

Cf. also A000040, A000045/A039834.

\\ Needs also code from A265398 and A265399.
A265752 = n -> valuation(A265399(n),2);
for(n=1, 100, write("b265752.txt", n, " ", A265752(n)));

(define (A265752 n) (A007814 (A265399 n)))

a(n) = A007814(A265399(n)).
Other identities. For all n >= 1:
a(A000040(n+1)) = A000045(n-1). [Generalized by Peter Munn, Apr 05 2021]
a(A206296(n)) = A192232(n).
a(A265750(n)) = A192750(n).

A265753 a(n) = A007949(A265399(n)).

Original entry on oeis.org

0, 0, 1, 0, 1, 1, 2, 0, 2, 1, 3, 1, 5, 2, 2, 0, 8, 2, 13, 1, 3, 3, 21, 1, 2, 5, 3, 2, 34, 2, 55, 0, 4, 8, 3, 2, 89, 13, 6, 1, 144, 3, 233, 3, 3, 21, 377, 1, 4, 2, 9, 5, 610, 3, 4, 2, 14, 34, 987, 2, 1597, 55, 4, 0, 6, 4, 2584, 8, 22, 3, 4181, 2, 6765, 89, 3, 13, 5, 6, 10946, 1, 4, 144, 17711, 3, 9, 233, 35, 3, 28657, 3, 7, 21

Offset: 1

Antti Karttunen, Dec 15 2015

nonn

a(n) = Coefficient of x in the reduction under x^2->x+1 of the polynomial encoded in the prime factorization of n. (Assuming here only polynomials with nonnegative integer coefficients, see e.g. A206296 for the details).

Completely additive with a(prime(k)) = F(k-1), where F(k) denotes the k-th Fibonacci number, A000045(k). - Peter Munn, Mar 29 2021, incorporating comment by Antti Karttunen, Dec 15 2015

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

Cf. A007949, A112576, A192751, A206296, A265399, A265750, A265752.

Cf. also A000040, A000045.

\\ Needs also code from A265398 and A265399.
A265753 = n -> valuation(A265399(n),3);
for(n=1, 100, write("b265753.txt", n, " ", A265753(n)));

(define (A265753 n) (A007949 (A265399 n)))

a(n) = A007949(A265399(n)).
Other identities. For all n >= 1:
a(A000040(n)) = A000045(n-1). [Generalized by Peter Munn, Mar 29 2021]
a(A206296(n)) = A112576(n).
a(A265750(n)) = A192751(n).

A192750 Define a pair of sequences c_n, d_n by c_0=0, d_0=1 and thereafter c_n = c_{n-1}+d_{n-1}, d_n = c_{n-1}+4*n+2; sequence here is d_n.

Original entry on oeis.org

1, 6, 11, 21, 36, 61, 101, 166, 271, 441, 716, 1161, 1881, 3046, 4931, 7981, 12916, 20901, 33821, 54726, 88551, 143281, 231836, 375121, 606961, 982086, 1589051, 2571141, 4160196, 6731341, 10891541, 17622886, 28514431, 46137321, 74651756

Offset: 0

Clark Kimberling, Jul 09 2011

nonn

Old definition was: constant term of the reduction by x^2->x+1 of the polynomial p(n,x) defined recursively by p(n,x) = x*p(n-1,x) + 4n+2 for n>0, with p(0,x)=1.

For discussions of polynomial reduction, see A192232 and A192744.

Harvey P. Dale, Table of n, a(n) for n = 0..1000
Feng-Zhen Zhao, The log-behavior of some sequences related to the generalized Leonardo numbers, Integers (2024) Vol. 24, Art. No. A82.
Index entries for linear recurrences with constant coefficients, signature (2, 0, -1).

See A192751 for c_n.

Cf. A000045, A192744, A192232, A022088, A265750, A265752.

q = x^2; s = x + 1; z = 40;
p[0, n_] := 1; p[n_, x_] := x*p[n - 1, x] + 4 n + 2;
Table[Expand[p[n, x]], {n, 0, 7}]
reduce[{p1_, q_, s_, x_}] :=
FixedPoint[(s PolynomialQuotient @@ #1 +
       PolynomialRemainder @@ #1 &)[{#1, q, x}] &, p1]
t = Table[reduce[{p[n, x], q, s, x}], {n, 0, z}];
u1 = Table[Coefficient[Part[t, n], x, 0], {n, 1, z}]
  (* A192750 *)
u2 = Table[Coefficient[Part[t, n], x, 1], {n, 1, z}]
  (* A192751 *)
LinearRecurrence[{2,0,-1},{1,6,11},40] (* Harvey P. Dale, Dec 03 2023 *)

G.f.: ( 1+4*x-x^2 ) / ( (x-1)*(x^2+x-1) ). The first differences are in A022088. - R. J. Mathar, May 04 2014
a(n) = 5*Fibonacci(n+2)-4. - Gerry Martens, Jul 04 2015
a(n) = A265752(A265750(n)). - Antti Karttunen, Dec 15 2015

Entry revised by N. J. A. Sloane, Dec 15 2015

A192751 Define a pair of sequences c_n, d_n by c_0=0, d_0=1 and thereafter c_n = c_{n-1}+d_{n-1}, d_n = c_{n-1}+4*n+2; sequence here is c_n.

Original entry on oeis.org

0, 1, 7, 18, 39, 75, 136, 237, 403, 674, 1115, 1831, 2992, 4873, 7919, 12850, 20831, 33747, 54648, 88469, 143195, 231746, 375027, 606863, 981984, 1588945, 2571031, 4160082, 6731223, 10891419, 17622760, 28514301, 46137187, 74651618, 120788939

Offset: 0

Clark Kimberling, Jul 09 2011

nonn

Old definition was: coefficient of x in the reduction under x^2->x+1 of the polynomial p(n,x) defined recursively by p(n,x) = x*p(n-1,x) + 4n+2 for n>0, with p(0,x)=1.

For discussions of polynomial reduction, see A192232 and A192744.

Harvey P. Dale, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (3,-2,-1,1).

See A192750 for d_n.

Cf. A000045, A192744, A192232, A022088, A265750, A265752.

(See A192750.)
CoefficientList[Series[x (x^2-4x-1)/((x-1)^2(x^2+x-1)),{x,0,40}],x] (* or *) LinearRecurrence[{3,-2,-1,1},{0,1,7,18},40] (* Harvey P. Dale, Feb 23 2022 *)

G.f.: x*(x^2-4*x-1)/((x-1)^2*(x^2+x-1)). First differences are in A192750. [Colin Barker, Nov 13 2012]
a(n) = 5*Fibonacci(n+3) - (4*n+10). - N. J. A. Sloane, Dec 15 2015
a(n) = A265753(A265750(n)). - Antti Karttunen, Dec 15 2015

Description corrected by Antti Karttunen, Dec 15 2015

Entry revised by N. J. A. Sloane, Dec 15 2015

cp's OEIS Frontend

A003961 Completely multiplicative with a(prime(k)) = prime(k+1).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

Perl

Python

Formula

A265752 a(n) = A007814(A265399(n)).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Scheme

Formula

A265753 a(n) = A007949(A265399(n)).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Scheme

Formula

A192750 Define a pair of sequences c_n, d_n by c_0=0, d_0=1 and thereafter c_n = c_{n-1}+d_{n-1}, d_n = c_{n-1}+4*n+2; sequence here is d_n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A192751 Define a pair of sequences c_n, d_n by c_0=0, d_0=1 and thereafter c_n = c_{n-1}+d_{n-1}, d_n = c_{n-1}+4*n+2; sequence here is c_n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

Extensions