A083356 -id:A083356 - cp's OEIS frontend

A090252 The Two-Up sequence: a(n) is the least positive number not already used that is coprime to the previous floor(n/2) terms.

1, 2, 3, 5, 4, 7, 9, 11, 13, 17, 8, 19, 23, 25, 21, 29, 31, 37, 41, 43, 47, 53, 16, 59, 61, 67, 71, 73, 55, 79, 27, 49, 83, 89, 97, 101, 103, 107, 109, 113, 127, 131, 137, 139, 149, 151, 26, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199, 211, 85, 121, 223, 227, 57, 229

Offset: 1

Amarnath Murthy, Nov 27 2003

nonn

a(n) is coprime to the next n terms. - David Wasserman, Oct 24 2005
All values up to a(1000000) are either prime powers or semiprimes; this suggests the sequence is unlikely to be a permutation of the integers.
It appears that a(n) is even iff n = 3*2^k-1 for some k (A083356). - N. J. A. Sloane, Nov 01 2014
The even terms in the present sequence are listed in A354255.
We have a(1) = 1 and a(2) = 2. At step k >= 2, the sequence is extended by adding two terms: a(2*k-1) = smallest unused number which is relatively prime to a(k), a(k+1), ..., a(2*k-2), and a(2*k) = smallest unused number which is relatively prime to a(k), a(k+1), ..., a(2*k-1). So at step k=2 we add a(3)=3, a(4)=5; at step k=3 we add a(5)=4, a(6)=7; and so on. - N. J. A. Sloane, May 21 2022
Comments from N. J. A. Sloane, May 23 2022: (Start)
Conjecture 1. A090252 is a subsequence of A354144 (prime powers and semiprimes).
Conjecture 2. The terms of A354144 that are missing from A090252 are 6, 10, 14, 15, 22, 33, 34, 35, 38, 39, 46, 51, 58, 62, 65, 69, 74, 77, 82, 86, 87, 91, 93, 94, 95, 106, 111, 115, 118, 119, 122, 123, 129, 133, 134, 141, 142, 143, 145, 146, 155, 158, 166, 177, 178, 183, 185, 187, 194, 201, 202, 203, 209, 213, 214, 215, 218, 219, 221, ...
But since there is no proof that any one of these numbers is really missing, this list cannot yet have an entry in the OEIS.
Let S_p = list of indices of terms in A090252 that are divisible by the prime p.
Conjecture 3. For a prime p, there are constants v_1, v_2, ..., v_K and c such that
   S_p = { v_1, v_2, ..., v_k, lambda*2^i - 1, i >= c}.
For example, from Michael S. Branicky's 10000-term b-file, it appears that:
   S_2 = { 3*2^k-1, k >= 0 }  cf. A083329
   S_3 = { 2^k-1, k >= 2 }  cf. A000225
   S_5 = { 4 then 15*2^k-1 k >= 0 }  cf. A196305
   S_7 = { 6, 15, then 33*2^k-1, k >= 0  }
   S_11 = { 8, 29, then 61*2^k-1, k >= 0 }
   S_13 = { 9, 47, 97*2^n-1, n >= 0 }
   S_17 = { 10, 59, 121*2^n-1, n >= 0 }
   S_19 = { 12, 63, 129*2^n-1, n >= 0 }
   S_23 = { 13, 65, 133*2^n-1, n >= 0 }
   S_29 = { 16, 121, 245*2^n-1, n >= 0 }
   S_31 = { 17, 131, 265*2^n-1, n >= 0 }
The initial primes p and the corresponding values of lambda are:
    p:   2   3   5   7   11   13   17   19   23   29   31
lambda:..3...1..15..33...61...97..121..129..133..245..265
(This sequence of lambdas does not seem to have any simpler explanation, is not in the OEIS, and cannot be since the terms shown are all conjectural.)
Conjecture 2 is a consequence of Conjecture 3. For example, 6 does not appear in A090252, since the sets S_2 and S_3 are disjoint.
Also 10 does not appear, since S_2 and S_5 are disjoint.
In fact 2*p for 3 <= p <= 11 does not appear, but 26 = 2*13 does appear since S_2 and S_13 have 47 in common.
Assuming the numbers that appear to be missing (see Conjecture 2) really are missing, the numbers that take a record number of steps to appear are 1, 2, 3, 4, 7, 8, 16, 26, 32, 64, 128, 206, 256, 478, 512, 933, ..., and the indices where they appear are 1, 2, 3, 5, 6, 11, 23, 47, 95, 191, 383, 767, 1535, 3071, 6143, 8191, .... These two sequences are not yet in the OEIS, and cannot be added since the terms are all conjectural.
(End)
From N. J. A. Sloane, Jun 06 2022 (Start)
Theorem: (a) a(n) <= prime(n-1) for all n >= 2 (cf. A354154).
  (b) A stronger upper bound is the following. Let c(n) = A354166(n) denote the number of nonprime terms among a(1) .. a(n). Note c(1)=1. Then a(n) <= prime(n-c(n)) for n <> 7 and 14.
It appears that a(n) = prime(n-c(n)) for almost all n. That is, this is the equation to the line in the graph that contains most of the terms.
For example, a(34886) = 408710 (see the b-file) = prime(34886 - A354166(34886)) = prime(34886 - 374) = prime(34512) = 408710.
Another example: Consider Russ Cox's table of the first N = 5764982 terms. We see that a(5764982) = 99999989 = prime(5761455) = prime(N - 3527) which agrees with c(N) = 3527 (from the first Russ Cox link).
(End)
If we consider the May 23 2022 comment, note the conjectured indices show near complete overlap with terms of A081026: 1, 2, 3, 5, 6, 11, 23, 47, 95, 191, 383, 767, 1535, 3071, 6143, 8191. - Bill McEachen, Aug 09 2024

Michael S. Branicky, Table of n, a(n) for n = 1..34886
Russ Cox, Table of nonprime entries in A090252: n, A090252(n), # of prime factors, n = 1..3527.
Russ Cox, Table of n, a(n) for n = 1..5764982, up to the first term that is greater than 10^8 [gzipped file]
Michael De Vlieger, Thomas Scheuerle, Rémy Sigrist, N. J. A. Sloane, and Walter Trump, The Binary Two-Up Sequence, arXiv:2209.04108 [math.CO], Sep 11 2022.
N. J. A. Sloane, Blog post about the Two-Up sequence, June 13 2022.
Hugo van der Sanden, Perl program to calculate this sequence and A249064 (requires Math::Pari)
Hugo van der Sanden, Faster Perl program on github, used to compute 10^9 terms. [Link changed by _N. J. A. Sloane_, Jun 19 2022]
Hugo van der Sanden, Table of nonprime entries in the first 10^9 terms of A090252 [See beginning of the file for description. The blog in the above link has comments from _Hugo van der Sanden_ describing the algorithm used to generate this table.]

Cf. A083329, A084937, A196305, A249064, A354146, A354148-A354151, A354154, A354159-A354167, A354255 (even terms), A355893.
See also A354169, A354764, A354765, A355057.
See A247665 for the case when the numbers are required to be at least 2. A353730 is another version.
For a squarefree analog, see A354790, A354791, A354792.

nn = 120; c[] = 0; a[1] = c[1] = 1; u = 2; Do[k = u; While[Nand[c[k] == 0, AllTrue[Array[a[i - #] &, Floor[i/2]], CoprimeQ[#, k] &]], k++]; Set[{a[i], c[k]}, {k, i}]; If[k == u, While[c[u] > 0, u++]], {i, 2, nn}]; Array[a, nn] (* _Michael De Vlieger, May 21 2022 *)

A090252_first(N, U=[0], L=List())=vector(N, i, for(k=U[1]+1, oo, setsearch(U, k) && next; foreach(L,m, gcd(k,m)>1 && next(2)); bitand(i,1) || listpop(L,1); listput(L,k); if( k>U[1]+1, U=setunion(U,[k]), U[1]++; while(#U>1 && U[2]==U[1]+1, U=U[^1]));break); L[#L]) \\ M. F. Hasler, Jun 14 2022

from math import gcd, prod
from itertools import count, islice
def agen(): # generator of terms
    alst = [1]; aset = {1}; yield 1
    mink = 2
    for n in count(2):
        k, prodall = mink, prod(alst[n-n//2-1:n-1])
        while k in aset or gcd(prodall, k) != 1: k += 1
        alst.append(k); aset.add(k); yield k
        while mink in aset: mink += 1
print(list(islice(agen(), 64))) # Michael S. Branicky, May 21 2022

More terms from David Wasserman, Oct 24 2005

A143127 a(n) = Sum_{k=1..n} k*d(k) where d(k) is the number of divisors of k.

Original entry on oeis.org

1, 5, 11, 23, 33, 57, 71, 103, 130, 170, 192, 264, 290, 346, 406, 486, 520, 628, 666, 786, 870, 958, 1004, 1196, 1271, 1375, 1483, 1651, 1709, 1949, 2011, 2203, 2335, 2471, 2611, 2935, 3009, 3161, 3317, 3637, 3719, 4055, 4141, 4405, 4675, 4859, 4953, 5433

Offset: 1

Gary W. Adamson, Jul 26 2008

a(n) is also the sum of all parts of all partitions of all positive integers <= n into equal parts. - Omar E. Pol, May 29 2017

a(n) is also the sum of the multiples of k, not exceeding n, for k = 1, 2, ..., n. See a formula and an example below. - Wolfdieter Lang, Oct 18 2021

			a(3) = 11 = (1 + 4 + 6), where n*d(n) = (1, 4, 6, 12, 10, 24, ...).
a(4) = 23 = (8 + 7 + 5 + 3), where (8, 7, 5, 3) = row 4 of triangle A110661.
a(4) = 23 is the sum of [1 2 3 4|2 4|3|4] (multiples of k=1..4, not exceeding n). - _Wolfdieter Lang_, Oct 18 2021
a(4) = [1] + [2 + 1 + 1] + [3 + 1 + 1 + 1] + [4 + 2 + 2 + 1 + 1 + 1 + 1] = 23. - _Omar E. Pol_, Oct 18 2021

Enrique Pérez Herrero, Table of n, a(n) for n = 1..1000
Jacob F. F. Bulmer, Javier Martínez-Cifuentes, Bryn A. Bell, and Nicolás Quesada, Simulating lossy and partially distinguishable quantum optical circuits: theory, algorithms and applications to experiment validation and state preparation, arXiv:2412.17742 [quant-ph], 2024. See p. 29.
Vaclav Kotesovec, Graph - The asymptotic ratio (1000000 terms)

Partial sums of A038040.
Cf. A000005, A083356.
Row sums of triangle A110661.
Row sums of triangle A143310. - Gary W. Adamson, Aug 06 2008
Cf. A018804.

a143127 n = a143127_list !! (n-1)
a143127_list = scanl1 (+) a038040_list
-- Reinhard Zumkeller, Jan 21 2014

Accumulate[DivisorSigma[0, Range[48]] Range[48]] (* Giovanni Resta, May 29 2018 *)

a(n) = sum(k=1, n, k*numdiv(k)); \\ Michel Marcus, May 29 2018

from math import isqrt
def A143127(n): return -((k:=isqrt(n))*(k+1)>>1)**2+sum(i*(m:=n//i)*(1+m) for i in range(1,k+1)) # Chai Wah Wu, Jul 11 2023

a(n) = Sum_{k=1..n} A038040(k).
a(n) = Sum_{m=1..floor(sqrt(n))} m*(m+floor(n/m))*(floor(n/m)+1-m) - A000330(floor(sqrt(n))) = 2*A083356(n) - A000330(floor(sqrt(n))). - Max Alekseyev, Jan 31 2012
G.f.: x*f'(x)/(1 - x), where f(x) = Sum_{k>=1} x^k/(1 - x^k). - Ilya Gutkovskiy, Apr 13 2017 [Sum_{k>=1} k*x^k/((1-x)*(1-x^k)^2), see A038040. - Wolfdieter Lang, Oct 18 2021]
a(n) = Sum_{k=1..n} k/2 * floor(n/k) * floor(1 + n/k). - Daniel Suteu, May 28 2018
a(n) ~ log(n) * n^2 / 2 + (gamma - 1/4)*n^2, where gamma is the Euler-Mascheroni constant A001620. - Vaclav Kotesovec, Sep 08 2018
From Daniel Hoying, May 21 2020: (Start)
a(n) = (Sum_{i=1..floor(sqrt(n))} i*floor(n/i)*(1+floor(n/i))) - (floor(sqrt(n))*(1+floor(sqrt(n)))/2)^2;
     = (Sum_{i=1..floor(sqrt(n))} i*floor(n/i)*(1+floor(n/i))) - A000537(floor(sqrt(n))).
a(n) = A000537(floor(sqrt(n))) ; n=1;
     = A000537(floor(sqrt(n))) + n*(n+1) - floor(n/2)*(floor(n/2)+1) ; 1     = A000537(floor(sqrt(n))) + n*(n+1) - floor(n/2)*(floor(n/2)+1) + Sum_{i=floor(sqrt(n))+1..floor(n/2)} i*floor(n/i)*(1+floor(n/i)) ; n>=6. (End)
a(n) = Sum_{i=1..n} A018804(i)*floor(n/i). - Ridouane Oudra, Mar 15 2021
a(n) = Sum_{k=1..n} b(n,k), with b(n, k) = Sum_{j=1..floor(n/k)} j*k = k * floor(n/k) * (floor(n/k) + 1)/2. See the formula by Daniel Suteu above. - Wolfdieter Lang, Oct 18 2021

More terms from Carl Najafi, Dec 24 2011

Edited by Max Alekseyev, Jan 31 2012

A060872 Sum of dd' over all unordered pairs (d,d') with dd' = n.

Original entry on oeis.org

1, 2, 3, 8, 5, 12, 7, 16, 18, 20, 11, 36, 13, 28, 30, 48, 17, 54, 19, 60, 42, 44, 23, 96, 50, 52, 54, 84, 29, 120, 31, 96, 66, 68, 70, 180, 37, 76, 78, 160, 41, 168, 43, 132, 135, 92, 47, 240, 98, 150, 102, 156, 53, 216, 110, 224, 114, 116, 59, 360, 61, 124, 189, 256

Offset: 1

N. J. A. Sloane, May 04 2001

a(n) is also the sum of all parts of all partitions of n into consecutive parts that differ by 2. - Omar E. Pol, May 05 2020

			a(4)=8 because pairs of factors are 1*4 and 2*2 and 1*4 + 2*2 = 8.
For n = 16 there are three partitions of 16 into consecutive parts that differ by 2 (including 16 as a valid partition). They are [16], [9, 7] and [7, 5, 3, 1]. The sum of the parts is [16] + [9 + 7] + [7 + 5 + 3 + 1] = 48, so a(16) = 48. - _Omar E. Pol_, May 05 2020

Vincenzo Librandi, Table of n, a(n) for n = 1..5000

Cf. A038548, A060866.

First differences of A083356.

[n*Ceiling(DivisorSigma(0, n)/2): n in [1..70]]; // Vincenzo Librandi, Apr 12 2017

Table[ n * Ceiling[ DivisorSigma[0, n] /2 ], {n, 1, 73} ]

a(n) = n*ceil(numdiv(n)/2); \\ Michel Marcus, Jul 12 2023

from sympy import divisor_count
def A060872(n): return n*(divisor_count(n)+1>>1) # Chai Wah Wu, Jul 11 2023

a(n) = n * ceiling( d(n)/2) where d is the number of divisors function.
G.f.: x*f'(x), where f(x) = Sum_{k>=1} x^k^2/(1 - x^k). - Ilya Gutkovskiy, Apr 10 2017
a(n) = n*A038548(n). - Omar E. Pol, May 05 2020
Dirichlet g.f.: (zeta(2*s-2) + zeta(s-1)^2)/2. - Vaclav Kotesovec, Oct 21 2024

More terms from Robert G. Wilson v, Jun 23 2001

A354146 Even numbers in A353730 in order of appearance.

Original entry on oeis.org

2, 4, 8, 16, 26, 32, 64, 128, 206, 256, 454, 446, 512, 1024, 2048, 3142

Offset: 1

N. J. A. Sloane, May 21 2022

A090252(1535) = 256 and A090252(3071) = 478 are also even terms in A090252; the latter breaks the correspondence with this sequence. - Michael S. Branicky, May 21 2022

Cf. A083356, A090252, A153894, A353730, A354255.

cp's OEIS Frontend

A083357 Numbers n such that A083356(n) (the total area of all incongruent integer-sided rectangles of area <= n) is a square.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Extensions

A083358 Square root of n-th perfect square in A083356.

Views

Author

Keywords

Comments

Examples

Crossrefs

Formula

Extensions

A090252 The Two-Up sequence: a(n) is the least positive number not already used that is coprime to the previous floor(n/2) terms.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Python

Extensions

A143127 a(n) = Sum_{k=1..n} k*d(k) where d(k) is the number of divisors of k.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Python

Formula

Extensions

A060872 Sum of d*d' over all unordered pairs (d,d') with d*d' = n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Python

Formula

Extensions

A354146 Even numbers in A353730 in order of appearance.

Views

Author

Keywords

Comments

Crossrefs

Extensions

A060872 Sum of dd' over all unordered pairs (d,d') with dd' = n.