cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 12 results. Next

A064216 Replace each p^e with prevprime(p)^e in the prime factorization of odd numbers; inverse of sequence A048673 considered as a permutation of the natural numbers.

Original entry on oeis.org

1, 2, 3, 5, 4, 7, 11, 6, 13, 17, 10, 19, 9, 8, 23, 29, 14, 15, 31, 22, 37, 41, 12, 43, 25, 26, 47, 21, 34, 53, 59, 20, 33, 61, 38, 67, 71, 18, 35, 73, 16, 79, 39, 46, 83, 55, 58, 51, 89, 28, 97, 101, 30, 103, 107, 62, 109, 57, 44, 65, 49, 74, 27, 113, 82, 127, 85, 24, 131
Offset: 1

Views

Author

Howard A. Landman, Sep 21 2001

Keywords

Comments

a((A003961(n) + 1) / 2) = n and A003961(a(n)) = 2*n - 1 for all n. If the sequence is indexed by odd numbers only, it becomes multiplicative. In this variant sequence, denoted b, even indices don't exist, and we get b(1) = a(1) = 1, b(3) = a(2) = 2, b(5) = 3, b(7) = 5, b(9) = 4 = b(3) * b(3), ... , b(15) = 6 = b(3) * b(5), and so on. This property can also be stated as: a(x) * a(y) = a(((2x - 1) * (2y - 1) + 1) / 2) for x, y > 0. - Reinhard Zumkeller [re-expressed by Peter Munn, May 23 2020]
Not multiplicative in usual sense - but letting m=2n-1=product_j (p_j)^(e_j) then a(n)=a((m+1)/2)=product_j (p_(j-1))^(e_j). - Henry Bottomley, Apr 15 2005
From Antti Karttunen, Jul 25 2016: (Start)
Several permutations that use prime shift operation A064989 in their definition yield a permutation obtained from their odd bisection when composed with this permutation from the right. For example, we have:
A243505(n) = A122111(a(n)).
A243065(n) = A241909(a(n)).
A244153(n) = A156552(a(n)).
A245611(n) = A243071(a(n)).
(End)

Examples

			For n=11, the 11th odd number is 2*11 - 1 = 21 = 3^1 * 7^1. Replacing the primes 3 and 7 with the previous primes 2 and 5 gives 2^1 * 5^1 = 10, so a(11) = 10. - _Michael B. Porter_, Jul 25 2016
		

Crossrefs

Odd bisection of A064989 and A252463.
Row 1 of A251721, Row 2 of A249821.
Cf. A048673 (inverse permutation), A048674 (fixed points).
Cf. A246361 (numbers n such that a(n) <= n.)
Cf. A246362 (numbers n such that a(n) > n.)
Cf. A246371 (numbers n such that a(n) < n.)
Cf. A246372 (numbers n such that a(n) >= n.)
Cf. A246373 (primes p such that a(p) >= p.)
Cf. A246374 (primes p such that a(p) < p.)
Cf. A246343 (iterates starting from n=12.)
Cf. A246345 (iterates starting from n=16.)
Cf. A245448 (this permutation "squared", a(a(n)).)
Cf. A253894, A254044, A254045 (binary width, weight and the number of nonleading zeros in base-2 representation of a(n), respectively).
Cf. A285702, A285703 (phi and sigma applied to a(n).)
Here obviously the variant 2, A151799(n) = A007917(n-1), of the prevprime function is used.
Cf. also A003961, A270430, A270431.

Programs

  • Mathematica
    Table[Times @@ Power[If[# == 1, 1, NextPrime[#, -1]] & /@ First@ #, Last@ #] &@ Transpose@ FactorInteger[2 n - 1], {n, 69}] (* Michael De Vlieger, Dec 18 2014, revised Mar 17 2016 *)
  • PARI
    a(n) = {my(f = factor(2*n-1)); for (k=1, #f~, f[k,1] = precprime(f[k,1]-1)); factorback(f);} \\ Michel Marcus, Mar 17 2016
    
  • Python
    from sympy import factorint, prevprime
    from operator import mul
    def a(n):
        f=factorint(2*n - 1)
        return 1 if n==1 else reduce(mul, [prevprime(i)**f[i] for i in f]) # Indranil Ghosh, May 13 2017
  • Scheme
    (define (A064216 n) (A064989 (- (+ n n) 1))) ;; Antti Karttunen, May 12 2014
    

Formula

a(n) = A064989(2n - 1). - Antti Karttunen, May 12 2014
Sum_{k=1..n} a(k) ~ c * n^2, where c = Product_{p prime > 2} ((p^2-p)/(p^2-q(p))) = 0.6621117868..., where q(p) = prevprime(p) (A151799). - Amiram Eldar, Jan 21 2023

Extensions

More terms from Reinhard Zumkeller, Sep 26 2001
Additional description added by Antti Karttunen, May 12 2014

A246277 Column index of n in A246278: a(1) = 0, a(2n) = n, a(2n+1) = a(A064989(2n+1)).

Original entry on oeis.org

0, 1, 1, 2, 1, 3, 1, 4, 2, 5, 1, 6, 1, 7, 3, 8, 1, 9, 1, 10, 5, 11, 1, 12, 2, 13, 4, 14, 1, 15, 1, 16, 7, 17, 3, 18, 1, 19, 11, 20, 1, 21, 1, 22, 6, 23, 1, 24, 2, 25, 13, 26, 1, 27, 5, 28, 17, 29, 1, 30, 1, 31, 10, 32, 7, 33, 1, 34, 19, 35, 1, 36, 1, 37, 9, 38, 3, 39, 1, 40, 8, 41, 1, 42
Offset: 1

Views

Author

Antti Karttunen, Aug 21 2014

Keywords

Comments

If n >= 2, n occurs in column a(n) of A246278.
By convention, a(1) = 0 because 1 does not occur in A246278.

Crossrefs

Terms of A348717 halved. A305897 is the restricted growth sequence transform.
Positions of terms 1 .. 8 in this sequence are given by the following sequences: A000040, A001248, A006094, A030078, A090076, A251720, A090090, A030514.
Cf. A078898 (has the same role with array A083221 as this sequence has with A246278).
This sequence is also used in the definition of the following permutations: A246274, A246276, A246675, A246677, A246683, A249815, A249817 (A249818), A249823, A249825, A250244, A250245, A250247, A250249.
Also in the definition of arrays A249821, A251721, A251722.
Sum of prime indices of a(n) is A359358(n) + A001222(n) - 1, cf. A326844.
A112798 lists prime indices, length A001222, sum A056239.

Programs

  • Mathematica
    a246277[n_Integer] := Module[{f, p, a064989, a},
      f[x_] := Transpose@FactorInteger[x];
      p[x_] := Which[
        x == 1, 1,
        x == 2, 1,
        True, NextPrime[x, -1]];
      a064989[x_] := Times @@ Power[p /@ First[f[x]], Last[f[x]]];
      a[1] = 0;
      a[x_] := If[EvenQ[x], x/2, NestWhile[a064989, x, OddQ]/2];
    a/@Range[n]]; a246277[84] (* Michael De Vlieger, Dec 19 2014 *)
  • PARI
    A064989(n) = {my(f); f = factor(n); if((n>1 && f[1,1]==2), f[1,2] = 0); for (i=1, #f~, f[i,1] = precprime(f[i,1]-1)); factorback(f)};
    A246277(n) = { if(1==n, 0, while((n%2), n = A064989(n)); (n/2)); };
    
  • PARI
    A246277(n) = if(1==n, 0, my(f = factor(n), k = primepi(f[1,1])-1); for (i=1, #f~, f[i,1] = prime(primepi(f[i,1])-k)); factorback(f)/2); \\ Antti Karttunen, Apr 30 2022
    
  • Python
    from sympy import factorint, prevprime
    from operator import mul
    from functools import reduce
    def a064989(n):
        f=factorint(n)
        return 1 if n==1 else reduce(mul, [1 if i==2 else prevprime(i)**f[i] for i in f])
    def a(n): return 0 if n==1 else n//2 if n%2==0 else a(a064989(n))
    print([a(n) for n in range(1, 101)]) # Indranil Ghosh, Jun 15 2017
  • Scheme
    ;; two different variants, the second one employing memoizing definec-macro)
    (define (A246277 n) (if (= 1 n) 0 (let loop ((n n)) (if (even? n) (/ n 2) (loop (A064989 n))))))
    (definec (A246277 n) (cond ((= 1 n) 0) ((even? n) (/ n 2)) (else (A246277 (A064989 n)))))
    

Formula

a(1) = 0, a(2n) = n, a(2n+1) = a(A064989(2n+1)) = a(A064216(n+1)). [Cf. the formula for A252463.]
Instead of the equation for a(2n+1) above, we may write a(A003961(n)) = a(n). - Peter Munn, May 21 2022
Other identities. For all n >= 1, the following holds:
For all w >= 0, a(p_{i} * p_{j} * ... * p_{k}) = a(p_{i+w} * p_{j+w} * ... * p_{k+w}).
For all n >= 2, A001222(a(n)) = A001222(n)-1. [a(n) has one less prime factor than n. Thus each semiprime (A001358) is mapped to some prime (A000040), etc.]
For all n >= 2, a(n) = A078898(A249817(n)).
For semiprimes n = p_i * p_j, j >= i, a(n) = A000040(1+A243055(n)) = p_{1+j-i}.
a(n) = floor(A348717(n)/2). - Antti Karttunen, Apr 30 2022
If n has prime factorization Product_{i=1..k} prime(x_i), then a(n) = Product_{i=2..k} prime(x_i-x_1+1). The opposite version is A358195, prime indices A358172, even bisection A241916. - Gus Wiseman, Dec 29 2022

A083221 Sieve of Eratosthenes arranged as an array and read by antidiagonals as A(1,1), A(1,2), A(2,1), A(1,3), A(2,2), A(3,1), ...

Original entry on oeis.org

2, 4, 3, 6, 9, 5, 8, 15, 25, 7, 10, 21, 35, 49, 11, 12, 27, 55, 77, 121, 13, 14, 33, 65, 91, 143, 169, 17, 16, 39, 85, 119, 187, 221, 289, 19, 18, 45, 95, 133, 209, 247, 323, 361, 23, 20, 51, 115, 161, 253, 299, 391, 437, 529, 29, 22, 57, 125, 203, 319, 377, 493, 551, 667
Offset: 2

Views

Author

Yasutoshi Kohmoto, Jun 05 2003

Keywords

Comments

This is permutation of natural numbers larger than 1.
From Antti Karttunen, Dec 19 2014: (Start)
If we assume here that a(1) = 1 (but which is not explicitly included because outside of the array), then A252460 gives an inverse permutation. See also A249741.
For navigating in this array:
A055396(n) gives the row number of row where n occurs, and A078898(n) gives its column number, both starting their indexing from 1.
A250469(n) gives the number immediately below n, and when n is an odd number >= 3, A250470(n) gives the number immediately above n. If n is a composite, A249744(n) gives the number immediately left of n.
First cube of each row, which is {the initial prime of the row}^3 and also the first number neither a prime or semiprime, occurs on row n at position A250474(n).
(End)
The n-th row contains the numbers whose least prime factor is the n-th prime: A020639(T(n,k)) = A000040(n). - Franklin T. Adams-Watters, Aug 07 2015

Examples

			The top left corner of the array:
   2,   4,   6,    8,   10,   12,   14,   16,   18,   20,   22,   24,   26
   3,   9,  15,   21,   27,   33,   39,   45,   51,   57,   63,   69,   75
   5,  25,  35,   55,   65,   85,   95,  115,  125,  145,  155,  175,  185
   7,  49,  77,   91,  119,  133,  161,  203,  217,  259,  287,  301,  329
  11, 121, 143,  187,  209,  253,  319,  341,  407,  451,  473,  517,  583
  13, 169, 221,  247,  299,  377,  403,  481,  533,  559,  611,  689,  767
  17, 289, 323,  391,  493,  527,  629,  697,  731,  799,  901, 1003, 1037
  19, 361, 437,  551,  589,  703,  779,  817,  893, 1007, 1121, 1159, 1273
  23, 529, 667,  713,  851,  943,  989, 1081, 1219, 1357, 1403, 1541, 1633
  29, 841, 899, 1073, 1189, 1247, 1363, 1537, 1711, 1769, 1943, 2059, 2117
  ...
		

Crossrefs

Transpose of A083140.
One more than A249741.
Inverse permutation: A252460.
Column 1: A000040, Column 2: A001248.
Row 1: A005843, Row 2: A016945, Row 3: A084967, Row 4: A084968, Row 5: A084969, Row 6: A084970.
Main diagonal: A083141.
First semiprime in each column occurs at A251717; A251718 & A251719 with additional criteria. A251724 gives the corresponding semiprimes for the latter. See also A251728.
Permutations based on mapping numbers between this array and A246278: A249817, A249818, A250244, A250245, A250247, A250249. See also: A249811, A249814, A249815.
Also used in the definition of the following arrays of permutations: A249821, A251721, A251722.

Programs

  • Mathematica
    lim = 11; a = Table[Take[Prime[n] Select[Range[lim^2], GCD[# Prime@ n, Product[Prime@ i, {i, 1, n - 1}]] == 1 &], lim], {n, lim}]; Flatten[Table[a[[i, n - i + 1]], {n, lim}, {i, n}]] (* Michael De Vlieger, Jan 04 2016, after Yasutoshi Kohmoto at A083140 *)

Extensions

More terms from Hugo Pfoertner, Jun 13 2003

A249745 Permutation of natural numbers: a(n) = (1 + A064989(A007310(n))) / 2.

Original entry on oeis.org

1, 2, 3, 4, 6, 7, 9, 10, 5, 12, 15, 8, 16, 19, 21, 22, 13, 24, 11, 27, 30, 17, 31, 34, 36, 18, 37, 40, 20, 42, 28, 26, 45, 49, 51, 52, 54, 55, 29, 33, 25, 14, 57, 64, 43, 66, 69, 39, 35, 70, 75, 44, 76, 48, 79, 82, 61, 84, 23, 87, 90, 47, 46, 91, 96, 97, 99, 58, 56, 60, 100, 62, 73, 72, 106, 112, 114, 115, 65, 117, 120, 38, 94, 121
Offset: 1

Views

Author

Antti Karttunen, Nov 23 2014

Keywords

Crossrefs

Inverse: A249746.
Row 2 of A251721.

Programs

  • Mathematica
    a249745[n_Integer] := Module[{f, p, a064989, a007310, a},
      f[x_] := Transpose@FactorInteger[x];
      p[x_] := Which[
        x == 1, 1,
        x == 2, 1,
        True, NextPrime[x, -1]];
      a064989[x_] := Times @@ Power[p /@ First[f[x]], Last[f[x]]];
      a007310[x_] := Select[Range[x], MemberQ[{1, 5}, Mod[#, 6]] &];
      a[x_] := (1 + a064989 /@ a007310[x])/2;
    a[n]]; a249745[252] (* Michael De Vlieger, Dec 18 2014, after Harvey P. Dale at A007310 *)
  • PARI
    A249745(n)=A064989(A007310(n))\2+1 \\ M. F. Hasler, Jan 19 2016
  • Scheme
    (define (A249745 n) (/ (+ 1 (A064989 (A007310 n))) 2))
    

Formula

a(n) = (A064989(A007310(n)) + 1) / 2.
a(n) = A048673(A249823(n)), as a composition of related permutations.
A007310(n) = A249735(a(n)) for all n >= 1. (This is the permutation which sorts the terms of A249735 into an ascending order, as they occur in A007310.)

A249821 Square array of permutations: A(row,col) = A246277(A083221(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, 4, 3, 2, 1, 5, 5, 3, 2, 1, 6, 4, 5, 3, 2, 1, 7, 7, 7, 5, 3, 2, 1, 8, 11, 11, 7, 5, 3, 2, 1, 9, 6, 13, 11, 7, 5, 3, 2, 1, 10, 13, 17, 13, 11, 7, 5, 3, 2, 1, 11, 17, 4, 17, 13, 11, 7, 5, 3, 2, 1, 12, 10, 19, 19, 17, 13, 11, 7, 5, 3, 2, 1, 13, 19, 23, 23, 19, 17, 13, 11, 7, 5, 3, 2, 1, 14, 9, 6, 29, 23, 19, 17, 13, 11, 7, 5, 3, 2, 1, 15, 8, 29, 31, 29, 23, 19, 17, 13, 11, 7, 5, 3, 2, 1
Offset: 1

Views

Author

Antti Karttunen, Nov 06 2014

Keywords

Comments

Permutation A249817 preserves the smallest prime factor of n, i.e., A055396(A249817(n)) = A055396(n), in other words, keeps all the terms that appear on any row of A246278 on the same row of A083221. Permutations in this table are induced by changes that A249817 does onto each row of the latter table, thus permutation on row r of this table can be used to sort row r of A246278 into ascending order. I.e., A246278(r, A(r,c)) = A083221(r,c) [the corresponding row in the Sieve of Eratosthenes, where each row appears in monotone order].
The multi-set of cycle-sizes of permutation A249817 is a disjoint union of cycle-sizes of all permutations in this array. For example, A249817 has a 7-cycle (33 39 63 57 99 81 45) which originates from the 7-cycle (6 7 11 10 17 14 8) of A064216, which occurs as the second row in this table.
On each row, 4 is the first composite number (and the first term less than previous, apart from row 1), and on row n it occurs in position A250474(n). This follows because A001222(A246277(n)) = A001222(n)-1 and because on each row of A083221 (see A083140) all terms between the square of prime (second term on each row) and the first cube (of the same prime, this cube mapping in this array to 4) are nonsquare semiprimes (A006881), this implies that the corresponding terms in this array must be primes.
Also, as the smaller prime factor of the terms on row n of A083221 is constant, A020639(n), and for all i < j: A246277(p_{i} * p_{j}) < A246277(p_i * p_{j+1}), the primes on any row appear in monotone order.

Examples

			The top left corner of the array:
1, 2, 3, 4, 5,  6,  7,  8,  9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, ...
1, 2, 3, 5, 4,  7, 11,  6, 13, 17, 10, 19,  9,  8, 23, 29, 14, 15, 31, ...
1, 2, 3, 5, 7, 11, 13, 17,  4, 19, 23,  6, 29, 31, 37, 41,  9, 43, 10, ...
1, 2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37,  4, 41, 43, 47, 53, 59, ...
...
		

Crossrefs

Inverse permutations can be found from table A249822.
Row k+1 is a left-to-right composition of the first k rows of A251721.
Row 1: A000027 (an identity permutation), Row 2: A064216, Row 3: A249823, Row 4: A249825.
The initial growing part of each row converges towards A008578.

Programs

Formula

A(row,col) = A246277(A083221(row,col)).
A001222(A(row,col)) = A001222(A083221(row,col)) - 1. [This follows directly from the properties of A246277.]

A250474 Number of times prime(n) occurs as the least prime factor among numbers 1 .. prime(n)^3: a(n) = A078898(A030078(n)).

Original entry on oeis.org

4, 5, 9, 14, 28, 36, 57, 67, 93, 139, 154, 210, 253, 272, 317, 396, 473, 504, 593, 658, 687, 792, 866, 979, 1141, 1229, 1270, 1356, 1397, 1496, 1849, 1947, 2111, 2159, 2457, 2514, 2695, 2880, 3007, 3204, 3398, 3473, 3828, 3904, 4047, 4121, 4583, 5061, 5228, 5309, 5474, 5743, 5832, 6269, 6543, 6816, 7107, 7197, 7488, 7686, 7784, 8295, 9029, 9248, 9354, 9568, 10351
Offset: 1

Views

Author

Antti Karttunen, Nov 23 2014

Keywords

Comments

Position of the first composite number (which is always 4) on row n of A249821. The fourth column of A249822.
Position of the first nonfixed term on row n of arrays of permutations A251721 and A251722.
According to the definition, this is the number of multiples of prime(n) below prime(n)^3 (and thus, the number of numbers below prime(n)^2) which do not have a smaller factor than prime(n). That is, the numbers remaining below prime(n)^2 after deleting all multiples of primes less than prime(n), as is done by applying the first n-1 steps of the sieve of Eratosthenes (when the first step is elimination of multiples of 2). This explains that the first differences are a(n+1)-a(n) = A050216(n)-1 for n>1, and a(n) = A054272(n)+2. - M. F. Hasler, Dec 31 2014

Examples

			prime(1) = 2 occurs as the least prime factor in range [1,8] for four times (all even numbers <= 8), thus a(1) = 4.
prime(2) = 3 occurs as the least prime factor in range [1,27] for five times (when n is: 3, 9, 15, 21, 27), thus a(2) = 5.
		

Crossrefs

One more than A250473. Two more than A054272.
Column 4 of A249822.
Cf. also A250477 (column 6), A250478 (column 8).

Programs

  • Mathematica
    f[n_] := Count[Range[Prime[n]^3], x_ /; Min[First /@ FactorInteger[x]] == Prime@ n]; Array[f, 16] (* Michael De Vlieger, Mar 30 2015 *)
  • PARI
    A250474(n) = 3 + primepi(prime(n)^2) - n; \\ Fast implementation.
    for(n=1, 5001, write("b250474.txt", n, " ", A250474(n)));
    \\ The following program reflects the given sum formula, but is far from the optimal solution:
    allocatemem(234567890);
    A002110(n) = prod(i=1, n, prime(i));
    A020639(n) = if(1==n,n,vecmin(factor(n)[,1]));
    A055396(n) = if(1==n,0,primepi(A020639(n)));
    A250474(n) = { my(p2 = prime(n)^2); sumdiv(A002110(n-1), d, moebius(d)*(p2\d)); };
    for(n=1, 23, print1(A250474(n),", "));
    
  • Scheme
    (define (A250474 n) (let loop ((k 2)) (if (not (prime? (A249821bi n k))) k (loop (+ k 1))))) ;; This is even slower. Code for A249821bi given in A249821.

Formula

a(n) = 3 + A000879(n) - n = A054272(n) + 2 = A250473(n) + 1.
a(n) = A078898(A030078(n)).
a(1) = 1, a(n) = Sum_{d|A002110(n-1)} moebius(d)*floor(prime(n)^2/d). [Follows when A030078(n), prime(n)^3 is substituted to the similar formula given for A078898(n). Here A002110(n) gives the product of the first n primes. Because the latter is always squarefree, one could use also Liouville's lambda (A008836) instead of Moebius mu (A008683)].
Other identities. For all n >= 1:
A249821(n, a(n)) = 4.

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

Views

Author

Antti Karttunen, Dec 07 2014

Keywords

Comments

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.

Examples

			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, ...
...
		

Crossrefs

Inverse permutations can be found from array A251721.
Row 1: A048673, Row 2: A249746, Row 3: A250476.

Programs

Formula

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

A249820 a(1) = 0 and for n > 1: a(n) = A249810(n) - A078898(n) = A078898(A003961(n)) - A078898(n).

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 1, 0, -1, 0, 2, 0, -1, 0, 6, 0, 4, 0, 1, 0, -4, 0, 11, 0, -4, 4, 3, 0, 3, 0, 25, -1, -7, 0, 20, 0, -7, -1, 12, 0, 7, 0, -2, 4, -8, 0, 44, 0, 0, -2, 0, 0, 36, 0, 22, -2, -13, 0, 23, 0, -12, 8, 90, 0, 0, 0, -5, -2, 4, 0, 77, 0, -16, 4, -3, 0, 4, 0, 55, 28, -19, 0, 41, 0, -19, -4, 15, 0, 43, 0, -2, -3, -20, 0, 155, 0, 12, 5, 24, 0
Offset: 1

Views

Author

Antti Karttunen, Dec 08 2014

Keywords

Comments

a(n) tells how many columns off A003961(n) is from the column where n is in square array A083221 (Cf. A083140, the sieve of Eratosthenes. The column index of n in that table is given by A078898(n)).

Examples

			For n = 8 = 2*2*2, A003961(8) = 27 (3*3*3), and while 8 is on row 1 and column 4 of A083221, 27 on the next row is in column 5, thus a(8) = 5 - 4 = 1.
For n = 10 = 2*5, A003961(10) = 21 (3*7), and while 10 is on row 1 and column 5 of A083221, 21 on the next row is in column 4, thus a(10) = 4 - 5 = -1.
		

Crossrefs

Programs

Formula

a(n) = A249810(n) - A078898(n) = A078898(A003961(n)) - A078898(n).
a(k) = 0 when k is a prime or square of prime, among some other numbers.

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

Views

Author

Antti Karttunen, Nov 23 2014

Keywords

Comments

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.

Crossrefs

One more than A054272.
One less than A250474.

Formula

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

A249810 a(1) = 0, a(n) = A078898(A003961(n)).

Original entry on oeis.org

0, 1, 1, 2, 1, 3, 1, 5, 2, 4, 1, 8, 1, 6, 3, 14, 1, 13, 1, 11, 4, 7, 1, 23, 2, 9, 9, 17, 1, 18, 1, 41, 5, 10, 3, 38, 1, 12, 6, 32, 1, 28, 1, 20, 12, 15, 1, 68, 2, 25, 7, 26, 1, 63, 4, 50, 8, 16, 1, 53, 1, 19, 19, 122, 5, 33, 1, 29, 10, 39, 1, 113, 1, 21, 17, 35, 3, 43, 1, 95, 42, 22, 1, 83, 6, 24, 11, 59, 1, 88, 4, 44, 13, 27, 7, 203
Offset: 1

Views

Author

Antti Karttunen, Dec 08 2014

Keywords

Crossrefs

Programs

Formula

a(1) = 0, a(n) = A078898(A003961(n)).
a(1) = 0, a(n) = A078898(n) + A249820(n).
Showing 1-10 of 12 results. Next