A250474 -id:A250474 - cp's OEIS frontend

A249822 Square array of permutations: A(row,col) = A078898(A246278(row,col)), read by antidiagonals A(1,1), A(1,2), A(2,1), A(1,3), A(2,2), A(3,1), ...

1, 2, 1, 3, 2, 1, 4, 3, 2, 1, 5, 5, 3, 2, 1, 6, 4, 9, 3, 2, 1, 7, 8, 4, 14, 3, 2, 1, 8, 6, 12, 4, 28, 3, 2, 1, 9, 14, 5, 21, 4, 36, 3, 2, 1, 10, 13, 42, 5, 33, 4, 57, 3, 2, 1, 11, 11, 17, 92, 5, 45, 4, 67, 3, 2, 1, 12, 7, 19, 33, 305, 5, 63, 4, 93, 3, 2, 1, 13, 23, 6, 25, 39, 455, 5, 80, 4, 139, 3, 2, 1, 14, 9, 59, 6, 43, 61, 944, 5, 116, 4, 154, 3, 2, 1, 15, 17, 7, 144, 6, 52, 70, 1238, 5, 148, 4, 210, 3, 2, 1

Offset: 1

Antti Karttunen, Nov 06 2014

			The top left corner of the array:
1, 2, 3,  4,  5,   6,   7,    8,    9,   10,  11,   12,  13,   14,   15, ...
1, 2, 3,  5,  4,   8,   6,   14,   13,   11,   7,   23,   9,   17,   18, ...
1, 2, 3,  9,  4,  12,   5,   42,   17,   19,   6,   59,   7,   22,   26, ...
1, 2, 3, 14,  4,  21,   5,   92,   33,   25,   6,  144,   7,   32,   39, ...
1, 2, 3, 28,  4,  33,   5,  305,   39,   43,   6,  360,   7,   48,   50, ...
1, 2, 3, 36,  4,  45,   5,  455,   61,   52,   6,  597,   7,   63,   68, ...
1, 2, 3, 57,  4,  63,   5,  944,   70,   76,   6, 1053,   7,   95,   84, ...
1, 2, 3, 67,  4,  80,   5, 1238,   96,   99,   6, 1502,   7,  106,  121, ...
...

Inverse permutations can be found from table A249821.
Row k+1 is a right-to-left composition of the first k rows of A251722.
Row 1: A000027 (an identity permutation), Row 2: A048673, Row 3: A249824, Row 4: A249826.
Column 4: A250474, Column 6: A250477, Column 8: A250478.
Cf. A002260, A004736, A078898, A246278, A249818.

(define (A249822 n) (A249822bi (A002260 n) (A004736 n)))
(define (A249822bi row col) (A078898 (A246278bi row col))) ;; Code for A246278bi given in A246278.

A251721 Square array of permutations: A(row,col) = A249822(row, A249821(row+1, col)), read by antidiagonals A(1,1), A(1,2), A(2,1), A(1,3), A(2,2), A(3,1), ...

Original entry on oeis.org

1, 2, 1, 3, 2, 1, 5, 3, 2, 1, 4, 4, 3, 2, 1, 7, 6, 4, 3, 2, 1, 11, 7, 5, 4, 3, 2, 1, 6, 9, 6, 5, 4, 3, 2, 1, 13, 10, 7, 6, 5, 4, 3, 2, 1, 17, 5, 8, 7, 6, 5, 4, 3, 2, 1, 10, 12, 10, 8, 7, 6, 5, 4, 3, 2, 1, 19, 15, 11, 9, 8, 7, 6, 5, 4, 3, 2, 1, 9, 8, 13, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 8, 16, 14, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 23, 19, 15, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1

Offset: 1

Antti Karttunen, Dec 07 2014

These are the "first differences" between permutations of array A249821, in a sense that by composing the first k rows of this array [from left to right, as in a(n) = row_1(row_2(...(row_k(n))))], one obtains row k+1 of A249821.

On row n, the first A250473(n) terms are fixed, and the first non-fixed term comes at A250474(n).

			The top left corner of the array:
1, 2, 3, 5, 4, 7, 11, 6, 13, 17, 10, 19, 9, 8, 23, 29, 14, 15, 31, 22, ...
1, 2, 3, 4, 6, 7, 9, 10, 5, 12, 15, 8, 16, 19, 21, 22, 13, 24, 11, 27, ...
1, 2, 3, 4, 5, 6, 7, 8, 10, 11, 13, 14, 15, 9, 16, 18, 20, 21, 23, 24, ...
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 15, 16, 17, 18, 19, 20, 22, ...
...

Inverse permutations can be found from array A251722.
Row 1: A064216, Row 2: A249745, Row 3: A250475.
Cf. A000027, A002260, A004736, A078898, A083221, A246277, A246278, A249821, A249822, A250473, A250474.

(define (A251721bi row col) (A249822bi row (A249821bi (+ row 1) col)))
(define (A251721 n) (A251721bi (A002260 n) (A004736 n)))
;; Code for A249821bi and A249822bi given in A249821, A249822.

A(row,col) = A249822(row, A249821(row+1, col)).

A(row,col) = A078898(A246278(row, A246277(A083221(row+1, col)))).

A251722 Square array of permutations: A(row,col) = A249822(row+1, A249821(row, col)), read by antidiagonals A(1,1), A(1,2), A(2,1), A(1,3), A(2,2), A(3,1), ...

Original entry on oeis.org

1, 2, 1, 3, 2, 1, 5, 3, 2, 1, 4, 4, 3, 2, 1, 8, 9, 4, 3, 2, 1, 6, 5, 5, 4, 3, 2, 1, 14, 6, 6, 5, 4, 3, 2, 1, 13, 12, 7, 6, 5, 4, 3, 2, 1, 11, 7, 8, 7, 6, 5, 4, 3, 2, 1, 7, 8, 14, 8, 7, 6, 5, 4, 3, 2, 1, 23, 19, 9, 9, 8, 7, 6, 5, 4, 3, 2, 1, 9, 10, 10, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 17, 17, 21, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 18, 42, 11, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1

Offset: 1

Antti Karttunen, Dec 07 2014

These are the "first differences" between permutations of array A249822, in a sense that by composing the first k rows of this array [from right to left, as in a(n) = row_k(...(row_2(row_1(n))))], one obtains row k+1 of A249822.

On row n the first non-fixed term is A250474(n+1) at position A250474(n), i.e., on row 1 it is 5 at n=4, on row 2 it is 9 at n=5, on row 3 it is 14 at n=9, etc. All the previous A250473(n) terms are fixed.

			The top left corner of the array:
1, 2, 3, 5, 4, 8, 6, 14, 13, 11, 7, 23, 9, 17, 18, 41, 10, 38, 12, 32, ...
1, 2, 3, 4, 9, 5, 6, 12, 7, 8, 19, 10, 17, 42, 11, 13, 22, 26, 14, 29, ...
1, 2, 3, 4, 5, 6, 7, 8, 14, 9, 10, 21, 11, 12, 13, 15, 33, 16, 25, 17, ...
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 28, 14, 15, 16, 17, 18, 19, ...
...

Inverse permutations can be found from array A251721.
Row 1: A048673, Row 2: A249746, Row 3: A250476.
Cf. A000027, A002260, A004736, A078898, A083221, A246277, A246278, A249821, A249822, A250473, A250474.

(define (A251722bi row col) (A249822bi (+ row 1) (A249821bi row col)))
(define (A251722 n) (A251722bi (A002260 n) (A004736 n)))
;; Code for A249821bi and A249822bi given in A249821, A249822.

A(row,col) = A249822(row+1, A249821(row, col)).

A(row,col) = A078898(A246278(row+1, A246277(A083221(row, col)))).

A053726 "Flag numbers": number of dots that can be arranged in successive rows of K, K-1, K, K-1, K, ..., K-1, K (assuming there is a total of L > 1 rows of size K > 1).

Original entry on oeis.org

5, 8, 11, 13, 14, 17, 18, 20, 23, 25, 26, 28, 29, 32, 33, 35, 38, 39, 41, 43, 44, 46, 47, 48, 50, 53, 56, 58, 59, 60, 61, 62, 63, 65, 67, 68, 71, 72, 73, 74, 77, 78, 80, 81, 83, 85, 86, 88, 89, 92, 93, 94, 95, 98, 101, 102, 103, 104, 105, 107, 108, 109, 110, 111, 113, 116

Offset: 1

Dan Asimov, asimovd(AT)aol.com, Apr 09 2003

Numbers of the form F(K, L) = KL+(K-1)(L-1), K, L > 1, i.e. 2KL - (K+L) + 1, sorted and duplicates removed.
If K=1, L=1 were allowed, this would contain all positive integers.
Positive numbers > 1 but not of the form (odd primes plus one)/2. - Douglas Winston (douglas.winston(AT)srupc.com), Sep 11 2003
In other words, numbers n such that 2n-1, or equally, A064216(n) is a composite number. - Antti Karttunen, Apr 17 2015
Note: the following comment was originally applied in error to the numerically similar A246371. - Allan C. Wechsler, Aug 01 2022
From Matthijs Coster, Dec 22 2014: (Start)
Also area of (over 45 degree) rotated rectangles with sides > 1. The area of such rectangles is 2ab - a - b + 1 = 1/2((2a-1)(2b-1)+1).
Example: Here a = 3 and b = 5. The area = 23.
           *
          ***
         *****
          *****
           *****
            ***
             *
(End)
The smallest integer > k/2 and coprime to k, where k is the n-th odd composite number. - Mike Jones, Jul 22 2024
Numbers k such that A193773(k-1) > 1. - Allan C. Wechsler, Oct 22 2024

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

Essentially same as A104275, but without the initial one.
A144650 sorted into ascending order, with duplicates removes.
Cf. A006254 (complement, apart from 1, which is in neither sequence).
Cf. also A000720, A008508, A064216, A071904, A199593, A250474.
Differs from its subsequence A246371 for the first time at a(8) = 20, which is missing from A246371.

select( {is_A053726(n)=n>4 && !isprime(n*2-1)}, [1..115]) \\ M. F. Hasler, Aug 02 2022

from sympy import isprime
def ok(n): return n > 1 and not isprime(2*n-1)
print(list(filter(ok, range(1, 117)))) # Michael S. Branicky, May 08 2021

from sympy import primepi
def A053726(n):
    if n == 1: return 5
    m, k = n, (r:=primepi(n)) + n + (n>>1)
    while m != k:
        m, k = k, (r:=primepi(k)) + n + (k>>1)
    return r+n # Chai Wah Wu, Aug 02 2024

;; with Antti Karttunen's IntSeq-library.
(define A053726 (MATCHING-POS 1 1 (lambda (n) (and (> n 1) (not (prime? (+ n n -1)))))))
;; Antti Karttunen, Apr 17 2015

;; with Antti Karttunen's IntSeq-library.
(define (A053726 n) (+ n (A000720 (A071904 n))))
;; Antti Karttunen, Apr 17 2015

a(n) = A008508(n) + n + 1.
From Antti Karttunen, Apr 17 2015: (Start)
a(n) = n + A000720(A071904(n)). [The above formula reduces to this. A000720(k) gives number of primes <= k, and A071904 gives the n-th odd composite number.]
a(n) = A104275(n+1). (End)
a(n) = A116922(A071904(n)). - Mike Jones, Jul 22 2024
a(n) = A047845(n+1)+1. - Amiram Eldar, Jul 30 2024

More terms from Douglas Winston (douglas.winston(AT)srupc.com), Sep 11 2003

A250477 Number of times prime(n) (the n-th prime) occurs as the least prime factor among numbers 1 .. (prime(n)^2 * prime(n+1)): a(n) = A078898(A251720(n)).

Original entry on oeis.org

6, 8, 12, 21, 33, 45, 63, 80, 116, 148, 182, 232, 265, 296, 356, 433, 490, 548, 625, 674, 740, 829, 919, 1055, 1187, 1252, 1313, 1376, 1446, 1657, 1897, 2029, 2134, 2301, 2484, 2605, 2785, 2946, 3110, 3301, 3439, 3654, 3869, 3978, 4086, 4349, 4811, 5147, 5273, 5395, 5604, 5787, 6049, 6403, 6684, 6954, 7153

Offset: 1

Antti Karttunen, Dec 14 2014

nonn

a(n) = Position of 6 on row n of array A249821. This is always larger than A250474(n), the position of 4 on row n, as 4 is guaranteed to be the first composite term on each row of A249821.
From Antti Karttunen, Mar 29 2015: (Start)
a(n) = 1 + number of positive integers <= (prime(n)*prime(n+1)) whose smallest prime factor is at least prime(n).
That a(n) >  A250474(n) can also be seen by realizing that prime(n) must occur at least as many times as the smallest prime factor for the numbers in range 1 .. (prime(n)^2 * prime(n+1)) than for numbers in (smaller) range 1 .. (prime(n)^3), and also by realizing that a(n) cannot be equal to A250474(n) because each row of A249822 is a permutation of natural numbers.
Or more simply, by considering the comment given in A256447 which follows from the new interpretation given above.
(End)

Antti Karttunen, Table of n, a(n) for n = 1..564
Antti Karttunen, Ratio a(n)/A250474(n+1), plotted up to n=65 with OEIS Plot2-utility

Column 6 of A249822. Cf. also A250474 (column 4), A250478 (column 8).
First differences: A256446. Cf. also A256447, A256448.
Cf. A002110, A006094, A020639, A078898, A249821, A251720, A281890.

f[n_] := Count[Range[Prime[n]^2*Prime[n + 1]], x_ /; Min[First /@ FactorInteger[x]] == Prime@ n]; Array[f, 20] (* Michael De Vlieger, Mar 30 2015 *)

allocatemem(234567890);
A002110(n) = prod(i=1, n, prime(i));
A250477(n) = { my(m); m = (prime(n) * prime(n+1)); sumdiv(A002110(n-1), d, (moebius(d)*(m\d))); };
for(n=1, 23, print1(A250477(n),", "));
\\ A more practical program:

allocatemem(234567890);
vecsize = (2^24)-4;
v020639 = vector(vecsize);
v020639[1] = 1; for(n=2,vecsize, v020639[n] = vecmin(factor(n)[, 1]));
A020639(n) = v020639[n];
A250477(n) = { my(p=prime(n),q=prime(n+1),u=p*q,k=1,s=1); while(k <= u, if(A020639(k) >= p, s++); k++); s; };
for(n=1, 564, write("b250477.txt", n, " ", A250477(n)));
\\ Antti Karttunen, Mar 29 2015

a(n) = A078898(A251720(n)).
a(1) = 1, a(n) = Sum_{d | A002110(n-1)} moebius(d) * floor(A006094(n) / d). [Follows when A251720, (p_n)^2 * p_{n+1} is substituted to the similar formula given for A078898. Here p_n is the n-th prime (A000040(n)), A006094(n) gives the product p_n * p{n+1} and A002110(n) gives the product of the first n primes. Because the latter is always squarefree, one could use here also Liouville's lambda (A008836) instead of Moebius mu (A008683)].
a(n) = A250474(n) + A256447(n).

A250473 Length of the maximal prefix of noncomposite numbers on row n of A249821.

Original entry on oeis.org

3, 4, 8, 13, 27, 35, 56, 66, 92, 138, 153, 209, 252, 271, 316, 395, 472, 503, 592, 657, 686, 791, 865, 978, 1140, 1228, 1269, 1355, 1396, 1495, 1848, 1946, 2110, 2158, 2456, 2513, 2694, 2879, 3006, 3203, 3397, 3472, 3827, 3903, 4046, 4120, 4582, 5060, 5227, 5308, 5473, 5742, 5831, 6268, 6542, 6815, 7106, 7196, 7487, 7685, 7783, 8294, 9028, 9247, 9353, 9567, 10350

Offset: 1

Antti Karttunen, Nov 23 2014

nonn

a(n) = Length of the initial fixed portion of the permutation at row n of A251721 & A251722.

Also, for n > 1, the length of the maximal increasing prefix on row n of A249821.

One more than A054272.
One less than A250474.
Cf. A000879, A249821, A251721, A251722.

a(n) = A054272(n) + 1.
a(n) = A250474(n) - 1.
a(n) = 2 + A000879(n) - n.

A250478 Number of times prime(n) occurs as the least prime factor among numbers 1 .. prime(n)^4: a(n) = A078898(A030514(n)).

Original entry on oeis.org

8, 14, 42, 92, 305, 455, 944, 1238, 2085, 3995, 4710, 7757, 10273, 11558, 14742, 20701, 28019, 30444, 39680, 46534, 49856, 62350, 71394, 86977, 111352, 124421, 130649, 145076, 151939, 167759, 236113, 257098, 291830, 302611, 370060, 382610, 427214, 475078

Offset: 1

Antti Karttunen, Dec 14 2014

Jon E. Schoenfield, Table of n, a(n) for n = 1..100

Column 8 of A249822.
Cf. also A250474 (column 4), A250477 (column 6).
Cf. A002110, A030514, A030078, A078898.

allocatemem(234567890);
A002110(n) = prod(i=1, n, prime(i));
A250478(n) = { my(p3); p3 = (prime(n)^3); sumdiv(A002110(n-1), d, (moebius(d)*(p3\d))); };
for(n=1, 23, print1(A250478(n),", "));

(define (A250478 n) (A078898 (A030514 n)))

a(n) = A078898(A030514(n)).

a(1) = 1, a(n) = sum_{d | A002110(n-1)} moebius(d) * floor(prime(n)^3 / d). [Follows when A030514, prime(n)^4 is substituted to the similar formula given for A078898. Here A002110(n) gives the product of the first n primes. Because the latter is always squarefree, one could use here also Liouville's lambda (A008836) instead of Moebius mu (A008683).]

More terms from Jon E. Schoenfield, Dec 14 2014

A256449 a(n) = A256447(n) - A256448(n).

Original entry on oeis.org

3, 2, 1, 0, 2, -3, 2, 0, 0, 3, 0, 1, 5, 3, -1, -3, 3, -1, -1, 3, 1, 0, -7, -10, 4, 5, 0, -1, -1, -31, -2, 0, -2, -14, -3, 1, -5, 5, 9, 0, 7, 7, 6, 5, 4, -6, -22, 5, 9, 7, -9, -1, -3, -6, 9, -15, 2, 5, 14, -4, 11, -24, 13, 0, 4, -9, -8, -10, 6, -2, 0, -2, 16, 11, -7, -13, 7, -11, 0

Offset: 1

Antti Karttunen, Mar 29 2015

sign

Positions of zeros: 4, 8, 9, 11, 22, 27, 32, 40, 64, 71, 79, 104, 113, 126, 140, 201, 225, 332, 333, 394, 451, ...

Corresponding primes: 7, 19, 23, 31, 79, 103, 131, 173, 311, 353, 401, 569, 617, 701, 809, 1229, 1427, 2237, 2239, 2707, 3187, ...

Antti Karttunen, Table of n, a(n) for n = 1..564
A. Karttunen, Sequences A256449 and A251723 compared with OEIS Plot2-script (See also the ratio-plots linked in A256447.)

Cf. A250474, A250477, A251723, A256446, A256447, A256448, A256470.

(define (A256449 n) (- (A256447 n) (A256448 n)))
(define (A256449 n) (- (* 2 (A250477 n)) (A250474 n) (A250474 (+ n 1))))

a(n) = A256447(n) - A256448(n).
a(n) = 2*A250477(n) - A250474(n) - A250474(n+1).
a(n) = 3 - A256470(n).

cp's OEIS Frontend

A249822 Square array of permutations: A(row,col) = A078898(A246278(row,col)), read by antidiagonals A(1,1), A(1,2), A(2,1), A(1,3), A(2,2), A(3,1), ...

Views

Author

Keywords

Examples

Crossrefs

Programs

Scheme

A251721 Square array of permutations: A(row,col) = A249822(row, A249821(row+1, col)), read by antidiagonals A(1,1), A(1,2), A(2,1), A(1,3), A(2,2), A(3,1), ...

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Scheme

Formula

A251722 Square array of permutations: A(row,col) = A249822(row+1, A249821(row, col)), read by antidiagonals A(1,1), A(1,2), A(2,1), A(1,3), A(2,2), A(3,1), ...

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Scheme

Formula

A053726 "Flag numbers": number of dots that can be arranged in successive rows of K, K-1, K, K-1, K, ..., K-1, K (assuming there is a total of L > 1 rows of size K > 1).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Python

Python

Scheme

Scheme

Formula

Extensions

A250477 Number of times prime(n) (the n-th prime) occurs as the least prime factor among numbers 1 .. (prime(n)^2 * prime(n+1)): a(n) = A078898(A251720(n)).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

A250473 Length of the maximal prefix of noncomposite numbers on row n of A249821.

Views

Author

Keywords

Comments

Crossrefs

Formula

A250478 Number of times prime(n) occurs as the least prime factor among numbers 1 .. prime(n)^4: a(n) = A078898(A030514(n)).

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Scheme

Formula

Extensions

A256449 a(n) = A256447(n) - A256448(n).

Views

Author

Keywords

Comments

Links