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 13 results. Next

A272214 Square array read by antidiagonals upwards in which T(n,k) is the product of the n-th prime and the sum of the divisors of k, n >= 1, k >= 1.

Original entry on oeis.org

2, 3, 6, 5, 9, 8, 7, 15, 12, 14, 11, 21, 20, 21, 12, 13, 33, 28, 35, 18, 24, 17, 39, 44, 49, 30, 36, 16, 19, 51, 52, 77, 42, 60, 24, 30, 23, 57, 68, 91, 66, 84, 40, 45, 26, 29, 69, 76, 119, 78, 132, 56, 75, 39, 36, 31, 87, 92, 133, 102, 156, 88, 105, 65, 54, 24, 37, 93, 116, 161, 114, 204, 104, 165, 91, 90, 36, 56
Offset: 1

Views

Author

Omar E. Pol, Apr 28 2016

Keywords

Comments

From Omar E. Pol, Dec 21 2021: (Start)
Also triangle read by rows: T(n,j) = A000040(n-j+1)*A000203(j), 1 <= j <= n.
For a visualization of T(n,j) first consider a tower (a polycube) in which the terraces are the symmetric representation of sigma(j), for j = 1 to n, starting from the top, and the heights of the terraces are the first n prime numbers respectively starting from the base. Then T(n,j) can be represented with a set of A237271(j) right prisms of height A000040(n-j+1) since T(n,j) is also the total number of cubes that are exactly below the parts of the symmetric representation of sigma(j) in the tower.
The sum of the n-th row of triangle is A086718(n) equaling the volume of the tower whose largest side of the base is n and its total height is the n-th prime.
The tower is an member of the family of the stepped pyramids described in A245092 and of the towers described in A221529. That is an infinite family of symmetric polycubes whose volumes represent the convolution of A000203 with any other integer sequence. (End)

Examples

			The corner of the square array begins:
   2,  6,   8,  14,  12,  24,  16,  30,  26,  36, ...
   3,  9,  12,  21,  18,  36,  24,  45,  39,  54, ...
   5, 15,  20,  35,  30,  60,  40,  75,  65,  90, ...
   7, 21,  28,  49,  42,  84,  56, 105,  91, 126, ...
  11, 33,  44,  77,  66, 132,  88, 165, 143, 198, ...
  13, 39,  52,  91,  78, 156, 104, 195, 169, 234, ...
  17, 51,  68, 119, 102, 204, 136, 255, 221, 306, ...
  19, 57,  76, 133, 114, 228, 152, 285, 247, 342, ...
  23, 69,  92, 161, 138, 276, 184, 345, 299, 414, ...
  29, 87, 116, 203, 174, 348, 232, 435, 377, 522, ...
  ...
From _Omar E. Pol_, Dec 21 2021: (Start)
Written as a triangle the sequence begins:
   2;
   3,  6;
   5,  9,  8;
   7, 15, 12,  14;
  11, 21, 20,  21,  12;
  13, 33, 28,  35,  18,  24;
  17, 39, 44,  49,  30,  36, 16;
  19, 51, 52,  77,  42,  60, 24,  30;
  23, 57, 68,  91,  66,  84, 40,  45, 26;
  29, 69, 76, 119,  78, 132, 56,  75, 39, 36;
  31, 87, 92, 133, 102, 156, 88, 105, 65, 54, 24;
...
Row sums give A086718. (End)
		

Crossrefs

Rows 1-4 of the square array: A074400, A272027, A274535, A319527.
Columns 1-5 of the square array: A000040, A001748, A001749, A138636, A272470.
Main diagonal of the square array gives A272211.
Cf. A086718 (antidiagonal sums of the square array, row sums of the triangle).

Programs

  • Mathematica
    Table[Prime[#] DivisorSigma[1, k] &@(n - k + 1), {n, 12}, {k, n}] // Flatten (* Michael De Vlieger, Apr 28 2016 *)

Formula

T(n,k) = prime(n)*sigma(k) = A000040(n)*A000203(k), n >= 1, k >= 1.
T(n,k) = A272400(n+1,k).

A306330 Squarefree n with >= 3 factors that admit idempotent factorizations n = p*q.

Original entry on oeis.org

30, 42, 66, 78, 102, 105, 114, 130, 138, 165, 170, 174, 182, 186, 195, 210, 222, 246, 255, 258, 266, 273, 282, 285, 290, 318, 330, 345, 354, 366, 370, 390, 399, 402, 410, 426, 434, 435, 438, 455, 462, 465, 474, 498, 510, 518, 530, 534, 546, 555, 570, 582, 602
Offset: 1

Views

Author

Barry Fagin, Feb 07 2019

Keywords

Comments

An idempotent factorization of n is a way of writing n = p*q such that b^(k(p-1)(q-1)+1) is congruent to b mod n for any integer k >= 0 and any b in Z_n. For example, p = 19, q = 15 is an idempotent factorization of n = 285. All factorizations of semiprimes are idempotent, so this sequence is restricted to n with >= 3 factors. Idempotent factorizations have the property that p and q generate correctly functioning RSA keys, even if one or both are composite.
We show in the reference below that a bipartite factorization of a squarefree integer n = pq is idempotent if and only if lambda(pq) divides (p-1)(q-1).
(p and q are not required to be primes. - N. J. A. Sloane, Feb 08 2019)

Examples

			30 = 5 * 6, 42 = 7 * 6, 66 = 11 * 6, 78 = 13 * 6, 102 = 17 * 6, 105 = 7 * 15, 114 = 19 * 6, 130 = 13 * 10 are the idempotent factorizations for the first 8 terms in the sequence. 210 = 10 * 21 is the smallest n with a fully composite idempotent factorization, one in which both p and q are composite. The number n = p * 6 is idempotent for any prime p >= 5.
		

Crossrefs

Subsequence of A120944 (composite squarefree numbers).

Programs

  • PARI
    isok3(p, q, n) = frac((p-1)*(q-1)/lcm(znstar(n)[2])) == 0;
    isok(n) = {if (issquarefree(n) && omega(n) >= 3, my(d = divisors(n)); for (k=1, #d\2, if ((d[k] != 1) && isok3(d[k], n/d[k], n), return (1););););} \\ Michel Marcus, Feb 22 2019

Extensions

Edited by N. J. A. Sloane, Feb 08 2019

A272470 7 times the primes.

Original entry on oeis.org

14, 21, 35, 49, 77, 91, 119, 133, 161, 203, 217, 259, 287, 301, 329, 371, 413, 427, 469, 497, 511, 553, 581, 623, 679, 707, 721, 749, 763, 791, 889, 917, 959, 973, 1043, 1057, 1099, 1141, 1169, 1211, 1253, 1267, 1337, 1351, 1379, 1393, 1477, 1561, 1589, 1603, 1631, 1673, 1687, 1757, 1799, 1841, 1883, 1897
Offset: 1

Views

Author

Omar E. Pol, Apr 30 2016

Keywords

Crossrefs

Column 4 of A272214.
k times the primes (k=1..6): A000040, A100484, A001748, A001749, A001750, A138636.

Programs

  • Mathematica
    7 Prime@ Range@ 58 (* Michael De Vlieger, May 01 2016 *)
  • PARI
    a(n) = 7*prime(n); \\ Michel Marcus, May 01 2016
    
  • Python
    from sympy import prime
    for n in range(1,1000):print(7*prime(n),end=", ") # Soumil Mandal, May 08 2016

Formula

a(n) = 7*prime(n) = 7*A000040(n).

A286592 Compound filter (prime signature & deficiency/abundance): a(n) = P(A046523(n), A286449(n)), where P(n,k) is sequence A000027 used as a pairing function.

Original entry on oeis.org

1, 3, 5, 10, 8, 42, 17, 36, 40, 27, 30, 183, 47, 34, 51, 136, 57, 243, 80, 288, 72, 177, 122, 765, 194, 72, 308, 117, 192, 1020, 212, 528, 142, 259, 196, 1576, 255, 111, 196, 1059, 302, 1020, 327, 103, 202, 471, 380, 2823, 500, 832, 306, 132, 498, 765, 672, 1564, 747, 786, 668, 4620, 743, 282, 337, 2080, 502, 1020, 782, 165, 441, 696, 822, 6288, 905, 747, 1047, 202
Offset: 1

Views

Author

Antti Karttunen, May 21 2017

Keywords

Comments

The lowermost conspicuous horizontal line in the scatter plot (at about log 3) is caused by value 1020, which corresponds to the prime signature 30 (p*q*r) and deficiency -12 packed together with the pairing function (as A002260(1020) = 30 and A004736(1020) = 16, A286449(24) = 16 and A033879(24) = -12). This value occurs in this sequence (at least) in the positions given by A138636, from its third term 30 onward.

Crossrefs

Programs

Formula

a(n) = (1/2)*(2 + ((A046523(n)+A286449(n))^2) - A046523(n) - 3*A286449(n)).

A306508 Squarefree numbers that have fully composite idempotent factorizations.

Original entry on oeis.org

210, 462, 570, 1155, 1302, 1330, 1365, 1785, 2210, 2310, 2730, 3003, 3410, 3710, 3990, 4305, 4515, 4758, 4810, 5005, 5187, 5474, 5610, 5642, 6006, 6105, 6118, 6270, 6510, 6622, 6630, 7410, 7770, 8265, 8385, 8463, 8645, 9282, 9471, 9870, 10010, 10101, 10230, 10374, 10545, 10582
Offset: 1

Views

Author

Barry Fagin, Feb 20 2019

Keywords

Comments

Fully composite idempotent factorizations are bipartite factorizations n=p*q such that p and q are composite numbers with the property that for any b in Z_n, b^(k(p-1)(q-1)+1) is congruent to b mod n for any integer k >= 0. Idempotent factorizations have the property that p and q generate correctly functioning RSA keys, even if one or both are composite.
2730 has more than one fully composite idempotent factorization (10*273, 21*130). It is the smallest positive integer with that property. 7770 and 8463 are similar.

Examples

			210=10*21, 462=22*21, 570=10*57, 1155=21*55, 1302=6*217, 1330=10*133, 1365=15*91 and 1785=21*85 are the fully composite idempotent factorizations for the first eight terms.
		

Crossrefs

Subsequence of A120944 (composite squarefree numbers).
Subsequence of A306330 (composite squarefree numbers with idempotent factorizations).

Programs

  • PARI
    isokc(p, q, n) = (p != 1) && !isprime(p) && !isprime(q) && (frac((p-1)*(q-1)/lcm(znstar(n)[2])) == 0);
    isok(n) = {if (issquarefree(n) && omega(n) >= 3, my(d = divisors(n)); for (k=1, #d\2, if (isokc(d[k], n/d[k], n), return (1););););} \\ Michel Marcus, Feb 22 2019
  • Python
    for n in range(2,max_n):
        factor_list = numbthy.factor(n)
        numFactors = len(factor_list)
        if numFactors <= 3:
            continue
        if not bsflib.is_composite_and_square_free_with_list(n,factor_list):
            continue
        fciFactorizations = bsflib.fullyCompositeIdempotentFactorizations(n,factor_list)
        numFCIFs = len(fciFactorizations)
        if numFCIPs > 0:
            fcIdempotents += 1
        print(n)
    

A309131 Triangle read by rows: T(n, k) = (n - k)*prime(1 + k), with 0 <= k < n.

Original entry on oeis.org

2, 4, 3, 6, 6, 5, 8, 9, 10, 7, 10, 12, 15, 14, 11, 12, 15, 20, 21, 22, 13, 14, 18, 25, 28, 33, 26, 17, 16, 21, 30, 35, 44, 39, 34, 19, 18, 24, 35, 42, 55, 52, 51, 38, 23, 20, 27, 40, 49, 66, 65, 68, 57, 46, 29, 22, 30, 45, 56, 77, 78, 85, 76, 69, 58, 31
Offset: 1

Views

Author

Stefano Spezia, Jul 14 2019

Keywords

Comments

T(n, k) is the k-superdiagonal sum of an n X n Toeplitz matrix M(n) whose first row consists of successive prime numbers prime(1), ..., prime(n).
The h-th subdiagonal of the triangle T gives the primes multiplied by (h + 1).
The k-th column of the triangle T gives the multiples of prime(1 + k).
Also array A(n, k) = n*prime(1 + k) read by ascending antidiagonals, with 0 <= k < n. - Michel Marcus, Jul 15 2019

Examples

			The triangle T(n, k) begins:
---+-----------------------------------------------------
n\k|    0     1     2     3     4     5     6     7     8
---+-----------------------------------------------------
1  |    2
2  |    4     3
3  |    6     6     5
4  |    8     9    10     7
5  |   10    12    15    14    11
6  |   12    15    20    21    22    13
7  |   14    18    25    28    33    26    17
8  |   16    21    30    35    44    39    34    19
9  |   18    24    35    42    55    52    51    38    23
...
For n = 3 the matrix M(3) is
          2,         3,         5
    M_{2,1},         2,         3
    M_{3,1},   M_{3,2},         2
and therefore T(3, 0) = 2 + 2 + 2 = 6, T(3, 1) = 3 + 3 = 6, and T(3, 2) = 5.
		

Crossrefs

Cf. A000040: diagonal; A001747: 1st subdiagonal; A001748: 2nd subdiagonal; A001749: 3rd subdiagonal; A001750: 4th subdiagonal; A005843: 0th column; A008585: 1st column; A008587: 2nd column; A008589: 3rd column; A008593: 4th column; A008595: 5th column; A008599: 6th column; A008601: 7th column; A014148: row sums; A138636: 5th subdiagonal; A272470: 6th subdiagonal.

Programs

  • Magma
    [[(n-k)*NthPrime(1+k): k in [0..n-1]]: n in [1..11]]; // triangle output
    
  • Maple
    a:=(n, k)->(n-k)*ithprime(1+k): seq(seq(a(n, k), k=0..n-1), n=1..11);
  • Mathematica
    Flatten[Table[(n-k)*Prime[1+k],{n,1,11},{k,0,n-1}]]
  • PARI
    T(n, k) = (n - k)*prime(1 + k);
    tabl(nn) = for(n=1, nn, for(k=0, n-1, print1(T(n, k), ", ")); print); \\ triangle output
    
  • Sage
    [[(n-k)*Primes().unrank(k) for k in (0..n-1)] for n in (1..11)] # triangle output

Formula

T(n, k) = A025581(n, k)*A000040(1 + k).

A160684 Numbers k >= 1 such that A000045(k)/A000005(k) is an integer.

Original entry on oeis.org

1, 3, 4, 6, 12, 24, 30, 36, 42, 54, 60, 64, 66, 72, 78, 84, 96, 102, 108, 114, 120, 132, 138, 156, 168, 174, 180, 186, 192, 204, 216, 222, 228, 240, 246, 252, 258, 264, 276, 282, 288, 300, 312, 318, 348, 354, 360, 366, 372, 384, 396, 400, 402, 405, 408, 420, 426
Offset: 1

Views

Author

Ctibor O. Zizka, May 23 2009

Keywords

Crossrefs

Cf. A000045, A000005, A377032. Includes A138636 except for 18.

Programs

  • Maple
    A000045 := proc(n) combinat[fibonacci](n) ; end: A000005 := proc(n) numtheory[tau](n) ; end: for n from 1 to 580 do if A000045(n) mod A000005(n) = 0 then printf("%d,",n) ; fi; od: # R. J. Mathar, May 25 2009

Extensions

71 removed and sequence extended by R. J. Mathar, May 25 2009

A306812 Maximally idempotent integers with three or more factors.

Original entry on oeis.org

273, 455, 1729, 2109, 2255, 2387, 3367, 3515, 4433, 4697, 4921, 5673, 6643, 6935, 7667, 8103, 8723, 8729, 9139, 9455, 10235, 10787, 11543, 13237, 13505, 14497, 16211, 16385, 16523, 17507, 18031, 18907, 20033, 20801, 21437, 22649, 23579, 24583
Offset: 1

Views

Author

Barry Fagin, Mar 11 2019

Keywords

Comments

An integer n has an idempotent factorization n=pq if b^(k(p-1)(q-1)+1) is congruent to b mod n for any integer k >= 0 and any b in Z_n (see A306330). An integer is maximally idempotent if all its bipartite factorizations n=pq are idempotent.
There are 15506 maximally idempotent integers less than 2^30. 15189 have three factors, 315 have four, two have five. The smallest maximally idempotent integer with four factors is 63973=7*13*19*37, a Carmichael number. The two with five factors are 13*19*37*73*109 and 11*31*41*101*151. The smallest maximally idempotent integer with six factors is 11*31*41*61*101*151.

Examples

			273 is the smallest maximally idempotent integer.  Factorization is (3,7,13).  Bipartite factorizations are (3,91), (7,39), (13,21).  Lambda(273) = 12, (2*90),(6*38) and (12*20) are all divisible by 12, thus 273 is maximally idempotent.  The same is true for 455 = 5*7*13.  The next entry in the sequence, 1729=7*13*19, is a Carmichael number, but most Carmichael numbers are not maximally idempotent.
		

Crossrefs

Subsequence of A120944 (composite squarefree numbers). Subsequence of A306330 (squarefree numbers that admit idempotent factorizations). Members of the sequence with >= 4 factors for a subsequence of A306508 (squarefree integers with fully composite idempotent factorizations).

Programs

  • Python
    ## This uses a custom library of number theory functions and the numbthy library.
    ## Hopefully the names of the functions make the process clear.
    for n in range(2,max_n):
        factor_list = numbthy.factor(n)
        numFactors = len(factor_list)
        if numFactors <= 2: # skip primes and semiprimes
            continue
        if not bsflib.is_composite_and_square_free_with_list(n,factor_list):
            continue
        ipList = bsflib.idempotentPartitions(n, factor_list)
        if len(ipList) == 2**(numFactors-1)-1:
            print(n)

A348508 a(n) = A003959(n) - 2*n, where A003959 is multiplicative with a(p^e) = (p+1)^e.

Original entry on oeis.org

-1, -1, -2, 1, -4, 0, -6, 11, -2, -2, -10, 12, -12, -4, -6, 49, -16, 12, -18, 14, -10, -8, -22, 60, -14, -10, 10, 16, -28, 12, -30, 179, -18, -14, -22, 72, -36, -16, -22, 82, -40, 12, -42, 20, 6, -20, -46, 228, -34, 8, -30, 22, -52, 84, -38, 104, -34, -26, -58, 96, -60, -28, 2, 601, -46, 12, -66, 26, -42, 4, -70, 288
Offset: 1

Views

Author

Antti Karttunen, Oct 30 2021

Keywords

Crossrefs

Programs

  • Mathematica
    f[p_, e_] := (p + 1)^e; a[1] = -1; a[n_] := Times @@ f @@@ FactorInteger[n] - 2*n; Array[a, 100] (* Amiram Eldar, Oct 30 2021 *)
  • PARI
    A003959(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1]++); factorback(f); };
    A348508(n) = (A003959(n) - 2*n);

Formula

a(n) = A003959(n) - 2*n.
a(n) = A348507(n) - n.
a(n) = A348029(n) - A033879(n).
From Antti Karttunen, Dec 05 2021: (Start)
a(n) = A168036(n) + A348970(n).
For all n >= 1, a(A138636(n)) = 12.
(End)
a(p) = 1 - p if p prime. - Bernard Schott, Feb 17 2022

A219742 Bernoulli denominators with 8 divisors in increasing order (without repetitions).

Original entry on oeis.org

30, 42, 66, 138, 282, 354, 498, 642, 1002, 1074, 1362, 1434, 1578, 2082, 2154, 2298, 2658, 2802, 2874, 3018, 3378, 3522, 3882, 3954, 4314, 4962, 5034, 5178, 5322, 5898, 6114, 7122, 7338, 7554, 7698, 7842, 7914, 8202, 8634, 8922, 8994, 9138, 9714, 10722
Offset: 1

Views

Author

Arkadiusz Wesolowski, Nov 29 2012

Keywords

Comments

Let m, n >= 1 and let f(m) denote number of Bernoulli numbers less than or equal to 10^m having denominator divisible by a(n). For any n, f(m) = floor(10^m/(a(n)/6 - 1)). It appears that the fraction of even Bernoulli numbers with denominator 6 is not so close to 1/6.

Crossrefs

Programs

  • Mathematica
    6*Prime@Flatten@Position[Table[p = Prime[n]; Length@Select[Divisors[p - 1] + 1, PrimeQ], {n, 277}], 3]

Formula

a(n) = 6*A092307(n).
A002445 INTERSECT A138636.
Showing 1-10 of 13 results. Next