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.

Previous Showing 31-40 of 56 results. Next

A268675 Self-inverse permutation of natural numbers: a(1) = 1; a(n) = A000079(A193231(A007814(n))) * A250469(a(A268674(n))).

Original entry on oeis.org

1, 2, 3, 8, 5, 6, 7, 4, 21, 10, 11, 24, 13, 14, 15, 32, 17, 42, 19, 40, 9, 22, 23, 12, 55, 26, 27, 56, 29, 30, 31, 16, 69, 34, 35, 168, 37, 38, 39, 20, 41, 18, 43, 88, 93, 46, 47, 96, 91, 110, 123, 104, 53, 54, 25, 28, 117, 58, 59, 120, 61, 62, 63, 64, 65, 138, 67, 136, 33, 70, 71, 84, 73, 74, 75, 152, 77, 78, 79, 160
Offset: 1

Views

Author

Antti Karttunen, Feb 11 2016

Keywords

Crossrefs

Formula

a(1) = 1, and for n > 1, a(n) = A000079(A193231(A007814(n))) * A250469(a(A268674(n))).
Other identities. For all n >= 1:
A000035(a(n)) = A000035(n). [This permutation preserves the parity of n.]
A020639(a(n)) = A020639(n). [More generally, it preserves the smallest prime dividing n.]
A055396(a(n)) = A055396(n).

A269359 Self-inverse permutation of natural numbers: a(1)=1, a(A269360(n)) = A250469(1+a(n)), a(A250469(1+n)) = A269360(a(n)).

Original entry on oeis.org

1, 3, 2, 9, 6, 5, 26, 11, 4, 27, 8, 65, 66, 25, 16, 15, 120, 71, 36, 169, 76, 33, 74, 41, 14, 7, 10, 81, 86, 185, 206, 215, 22, 195, 50, 19, 330, 515, 196, 75, 24, 337, 186, 49, 46, 45, 348, 247, 44, 35, 358, 213, 116, 353, 290, 143, 106, 507, 536, 295, 1266, 1345, 226, 99, 12, 13, 512, 2321, 220, 123, 18, 1285, 306, 23, 40, 21
Offset: 1

Views

Author

Antti Karttunen, Mar 13 2016

Keywords

Crossrefs

Similar or related permutations: A244319, A269863, A269864, A269865, A269866, A269867.

Formula

a(1) = 1, after which for even n, a(n) = A250469(1+a(A268674(n-1))), for odd n, a(n) = A269360(a(A268674(n)-1)).
The declarative form can be expressed in terms of A250469 only:
a(1)=1, a(1+A250469(n)) = A250469(1+a(n)), a(A250469(1+n)) = 1+A250469(a(n)).

A346473 a(n) = A250469(n) - sigma(n).

Original entry on oeis.org

0, 0, 1, 2, 1, 3, 3, 6, 12, 9, 1, 5, 3, 15, 11, 14, 1, 12, 3, 15, 23, 27, 5, 9, 18, 33, 25, 25, 1, 15, 5, 30, 37, 45, 29, 14, 3, 51, 39, 27, 1, 27, 3, 45, 37, 63, 5, 17, 64, 54, 53, 55, 5, 39, 19, 45, 65, 81, 1, 9, 5, 87, 51, 62, 35, 51, 3, 75, 79, 63, 1, 18, 5, 105, 61, 85, 47, 63, 3, 51, 84, 117, 5, 25, 25, 123, 95
Offset: 1

Views

Author

Antti Karttunen, Jul 28 2021

Keywords

Comments

The first negative term is a(120) = -3.

Crossrefs

Programs

  • Mathematica
    Block[{g}, g[n_] := If[n == 1, 0, PrimePi@FactorInteger[n][[1, 1]]]; Function[s, MapIndexed[Lookup[s, g[First@ #2] + 1][[#1]] - Boole[First@ #2 == 1] - DivisorSigma[1, First@ #2] &, #] &@ Map[Position[Lookup[s, g@ #], #][[1, 1]] &, Range@ 120]]@ PositionIndex@ Array[g, 10^4]] (* Michael De Vlieger, Oct 18 2021 *)
  • PARI
    A346473(n) = A250469(n)-sigma(n);

Formula

a(n) = A250469(n) - A000203(n).
a(n) = A286385(n) - A280692(n).
a(A001359(n)) = 1 for all n >= 1.

A347377 Möbius transform of A280692, A003961(n) - A250469(n).

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 6, 0, -6, 0, 12, 0, -6, 0, 30, 0, 24, 0, 12, 0, -24, 0, 48, 0, -24, 60, 24, 0, 24, 0, 114, -20, -42, 0, 84, 0, -42, -10, 60, 0, 48, 0, 12, 60, -48, 0, 168, 0, 6, -30, 24, 0, 132, 0, 108, -30, -78, 0, 96, 0, -72, 120, 390, 0, 44, 0, 12, -30, 36, 0, 288, 0, -96, 60, 24, 0, 58, 0, 228, 360, -114
Offset: 1

Views

Author

Antti Karttunen, Sep 01 2021

Keywords

Crossrefs

Programs

Formula

a(n) = Sum_{d|n} A008683(n/d) * A280692(d).
a(n) = A003972(n) - A347376(n).

A347378 a(n) = A250469(n) - n.

Original entry on oeis.org

0, 1, 2, 5, 2, 9, 4, 13, 16, 17, 2, 21, 4, 25, 20, 29, 2, 33, 4, 37, 34, 41, 6, 45, 24, 49, 38, 53, 2, 57, 6, 61, 52, 65, 42, 69, 4, 73, 56, 77, 2, 81, 4, 85, 70, 89, 6, 93, 72, 97, 74, 101, 6, 105, 36, 109, 88, 113, 2, 117, 6, 121, 92, 125, 54, 129, 4, 133, 106, 137, 2, 141, 6, 145, 110, 149, 66, 153, 4, 157, 124
Offset: 1

Views

Author

Antti Karttunen, Sep 01 2021

Keywords

Crossrefs

Programs

Formula

a(n) = A250469(n) - n.
a(n) = n - A346476(n).
a(n) = A336853(n) - A280692(n).
For all n >= 1, a(A000040(n)) = A001223(n).

A349631 Dirichlet convolution of A003961 with A346479, which is Dirichlet inverse of A250469.

Original entry on oeis.org

1, 0, 0, 0, 0, 0, 0, 6, 0, -6, 0, 12, 0, -6, 0, 18, 0, 24, 0, 24, 0, -24, 0, 0, 0, -24, 60, 36, 0, 48, 0, 42, -20, -42, 0, -12, 0, -42, -10, 12, 0, 72, 0, 60, 60, -48, 0, -24, 0, 42, -30, 72, 0, -84, 0, 12, -30, -78, 0, -120, 0, -72, 120, 126, 0, 180, 0, 96, -30, 132, 0, -48, 0, -96, 60, 108, 0, 174, 0, -84, 120
Offset: 1

Views

Author

Antti Karttunen, Nov 27 2021

Keywords

Comments

Note that for n = 2..36, a(n) = -A349632(n).
Dirichlet convolution of this sequence with A347376 is A003972.

Crossrefs

Cf. A003961, A250469, A346479, A349632 (Dirichlet inverse).
Cf. also A003972, A347376, A349381.
Cf. also arrays A083221, A246278, A249821, A249822 and permutations A250245, A250246.

Programs

  • PARI
    up_to = 20000;
    A020639(n) = if(1==n,n,vecmin(factor(n)[, 1]));
    ordinal_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), pt); for(i=1, length(invec), if(mapisdefined(om,invec[i]), pt = mapget(om, invec[i]), pt = 0); outvec[i] = (1+pt); mapput(om,invec[i],(1+pt))); outvec; };
    v078898 = ordinal_transform(vector(up_to,n,A020639(n)));
    A078898(n) = v078898[n];
    A250469(n) = if(1==n,n,my(spn = nextprime(1+A020639(n)), c = A078898(n), k = 0); while(c, k++; if((1==k)||(A020639(k)>=spn),c -= 1)); (k*spn));
    DirInverseCorrect(v) = { my(u=vector(#v)); u[1] = (1/v[1]); for(n=2, #v, u[n] = (-u[1]*sumdiv(n, d, if(dA250469(n)));
    A346479(n) = v346479[n];
    A003961(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ From A003961
    A349631(n) = sumdiv(n,d,A003961(d)*A346479(n/d));

Formula

a(n) = Sum_{d|n} A003961(d) * A346479(n/d).

A349632 Dirichlet convolution of A250469 with A346234, which is Dirichlet inverse of A003961.

Original entry on oeis.org

1, 0, 0, 0, 0, 0, 0, -6, 0, 6, 0, -12, 0, 6, 0, -18, 0, -24, 0, -24, 0, 24, 0, 0, 0, 24, -60, -36, 0, -48, 0, -42, 20, 42, 0, 12, 0, 42, 10, -12, 0, -72, 0, -60, -60, 48, 0, 24, 0, -42, 30, -72, 0, 84, 0, -12, 30, 78, 0, 120, 0, 72, -120, -90, 0, -180, 0, -96, 30, -132, 0, 48, 0, 96, -60, -108, 0, -174, 0, 12, -120
Offset: 1

Views

Author

Antti Karttunen, Nov 27 2021

Keywords

Comments

Note that for n = 2..36, a(n) = -A349631(n).
Dirichlet convolution of this sequence with A003972 is A347376.

Crossrefs

Cf. A003961, A250469, A346234, A349631 (Dirichlet inverse).
Cf. also A003972, A347376, A349382.
Cf. also arrays A083221, A246278, A249821, A249822 and permutations A250245, A250246.

Programs

  • PARI
    up_to = 20000;
    A020639(n) = if(1==n,n,vecmin(factor(n)[, 1]));
    ordinal_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), pt); for(i=1, length(invec), if(mapisdefined(om,invec[i]), pt = mapget(om, invec[i]), pt = 0); outvec[i] = (1+pt); mapput(om,invec[i],(1+pt))); outvec; };
    v078898 = ordinal_transform(vector(up_to,n,A020639(n)));
    A078898(n) = v078898[n];
    A250469(n) = if(1==n,n,my(spn = nextprime(1+A020639(n)), c = A078898(n), k = 0); while(c, k++; if((1==k)||(A020639(k)>=spn),c -= 1)); (k*spn));
    A346234(n) = (moebius(n)*A003961(n));
    A349632(n) = sumdiv(n,d,A250469(n/d)*A346234(d));

Formula

a(n) = Sum_{d|n} A250469(d) * A346234(n/d).

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

Original entry on oeis.org

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

Views

Author

Keywords

Comments

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)

Examples

			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.
		

References

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

Crossrefs

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

Programs

  • Haskell
    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
    
  • Maple
    a:= n-> mul(nextprime(i[1])^i[2], i=ifactors(n)[2]):
    seq(a(n), n=1..80);  # Alois P. Heinz, Sep 13 2017
  • Mathematica
    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 *)
  • PARI
    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]))
    
  • PARI
    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
    
  • Perl
    use ntheory ":all";  sub a003961 { vecprod(map { next_prime($) } factor(shift)); }  # _Dana Jacobsen, Mar 06 2016
    
  • Python
    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

Formula

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

A055396 Smallest prime dividing n is a(n)-th prime (a(1)=0).

Original entry on oeis.org

0, 1, 2, 1, 3, 1, 4, 1, 2, 1, 5, 1, 6, 1, 2, 1, 7, 1, 8, 1, 2, 1, 9, 1, 3, 1, 2, 1, 10, 1, 11, 1, 2, 1, 3, 1, 12, 1, 2, 1, 13, 1, 14, 1, 2, 1, 15, 1, 4, 1, 2, 1, 16, 1, 3, 1, 2, 1, 17, 1, 18, 1, 2, 1, 3, 1, 19, 1, 2, 1, 20, 1, 21, 1, 2, 1, 4, 1, 22, 1, 2, 1, 23, 1, 3, 1, 2, 1, 24, 1, 4, 1, 2, 1, 3, 1
Offset: 1

Views

Author

Henry Bottomley, May 15 2000

Keywords

Comments

Grundy numbers of the game in which you decrease n by a number prime to n, and the game ends when 1 is reached. - Eric M. Schmidt, Jul 21 2013
a(n) = the smallest part of the partition having Heinz number n. We define the Heinz number of a partition p = [p_1, p_2, ..., p_r] as Product(p_j-th prime, j=1...r) (concept used by Alois P. Heinz in A215366 as an "encoding" of a partition). For example, for the partition [1, 1, 2, 4, 10] we get 2*2*3*7*29 = 2436. Example: a(21) = 2; indeed, the partition having Heinz number 21 = 3*7 is [2,4]. - Emeric Deutsch, Jun 04 2015
a(n) is the number of numbers whose largest proper divisor is n, i.e., for n>1, number of occurrences of n in A032742. - Stanislav Sykora, Nov 04 2016
For n > 1, a(n) gives the number of row where n occurs in arrays A083221 and A246278. - Antti Karttunen, Mar 07 2017

Examples

			a(15) = 2 because 15=3*5, 3<5 and 3 is the 2nd prime.
		

References

  • John H. Conway, On Numbers and Games, 2nd Edition, p. 129.

Crossrefs

Programs

  • Haskell
    a055396 = a049084 . a020639  -- Reinhard Zumkeller, Apr 05 2012
    
  • Maple
    with(numtheory):
    a:= n-> `if`(n=1, 0, pi(min(factorset(n)[]))):
    seq(a(n), n=1..100);  # Alois P. Heinz, Aug 03 2013
  • Mathematica
    a[1] = 0; a[n_] := PrimePi[ FactorInteger[n][[1, 1]] ]; Table[a[n], {n, 1, 96}](* Jean-François Alcover, Jun 11 2012 *)
  • PARI
    a(n)=if(n==1,0,primepi(factor(n)[1,1])) \\ Charles R Greathouse IV, Apr 23 2015
    
  • Python
    from sympy import primepi, isprime, primefactors
    def a049084(n): return primepi(n)*(1*isprime(n))
    def a(n): return 0 if n==1 else a049084(min(primefactors(n))) # Indranil Ghosh, May 05 2017

Formula

From Reinhard Zumkeller, May 22 2003: (Start)
a(n) = A049084(A020639(n)).
A000040(a(n)) = A020639(n); a(n) <= A061395(n).
(End)
From Antti Karttunen, Mar 07 2017: (Start)
A243055(n) = A061395(n) - a(n).
a(A276086(n)) = A257993(n).
(End)

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
Previous Showing 31-40 of 56 results. Next