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 11-20 of 37 results. Next

A001414 Integer log of n: sum of primes dividing n (with repetition). Also called sopfr(n).

Original entry on oeis.org

0, 2, 3, 4, 5, 5, 7, 6, 6, 7, 11, 7, 13, 9, 8, 8, 17, 8, 19, 9, 10, 13, 23, 9, 10, 15, 9, 11, 29, 10, 31, 10, 14, 19, 12, 10, 37, 21, 16, 11, 41, 12, 43, 15, 11, 25, 47, 11, 14, 12, 20, 17, 53, 11, 16, 13, 22, 31, 59, 12, 61, 33, 13, 12, 18, 16, 67, 21, 26, 14, 71, 12, 73, 39, 13, 23, 18, 18
Offset: 1

Views

Author

Keywords

Comments

MacMahon calls this the potency of n.
Downgrades the operators in a prime decomposition. E.g., 40 factors as 2^3 * 5 and sopfr(40) = 2 * 3 + 5 = 11.
Consider all ways of writing n as a product of zero, one, or more factors; sequence gives smallest sum of terms. - Amarnath Murthy, Jul 07 2001
a(n) <= n for all n, and a(n) = n iff n is 4 or a prime.
Look at the graph of this sequence. At the lower edge of the logarithmic scatterplot there is a set of fuzzy but unmistakable diagonal fringes, sloping toward the southeast. Their spacing gradually increases, and their slopes gradually decrease; they are more distinct toward the lower edge of the range. Is any explanation known? - Allan C. Wechsler, Oct 11 2015
For n >= 2, the glb and lub are: 3 * log(n) / log(3) <= a(n) <= n, where the lub occurs when n = 3^k, k >= 1. (Jakimczuk 2012) - Daniel Forgues, Oct 12 2015
Except for the initial term, row sums of A027746. - M. F. Hasler, Feb 08 2016
Atanassov proves that a(n) <= A065387(n) - n. - Charles R Greathouse IV, Dec 06 2016
From Robert G. Wilson v, Aug 15 2022: (Start)
Differs from A337310 beginning with n at 64, 192, 256, 320, 448, 512, ..., .
The number of terms which equal k is A000607(k).
The first occurrence of k>1 is A056240(k).
The last occurrence of k>1 is A000792(k).
The Amarnath Murthy comment of Jul 07 2001 is a result of the fundamental theorem of arithmetic.
(End)

Examples

			a(24) = 2+2+2+3 = 9.
a(30) = 10: 30 can be written as 30, 15*2, 10*3, 6*5, 5*3*2. The corresponding sums are 30, 17, 13, 11, 10. Among these 10 is the least.
		

References

  • K. Atanassov, New integer functions, related to ψ and σ functions. IV., Bull. Number Theory Related Topics 12 (1988), pp. 31-35.
  • Amarnath Murthy, Generalization of Partition function and introducing Smarandache Factor Partition, Smarandache Notions Journal, Vol. 11, 1-2-3, Spring-2000.
  • Amarnath Murthy and Charles Ashbacher, Generalized Partitions and Some New Ideas on Number Theory and Smarandache Sequences, Hexis, Phoenix; USA 2005. See Section 1.4.
  • Joe Roberts, Lure of the Integers, Math. Assoc. America, 1992, p. 89.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

A000607(n) gives the number of values of k for which A001414(k) = n.
Cf. A036349 (indices of even terms), A356163 (their char. function), A335657 (indices of odd terms), A289142 (of multiples of 3), A373371 (their char. function).
For sum of squares of prime factors see A067666, for cubes see A224787.
Other completely additive sequences with primes p mapped to a function of p include: A001222 (with a(p)=1), A056239 (with a(p)=primepi(p)), A059975 (with a(p)=p-1), A064097 (with a(p)=1+a(p-1)), A113177 (with a(p)=Fib(p)), A276085 (with a(p)=p#/p), A341885 (with a(p)=p*(p+1)/2), A373149 (with a(p)=prevprime(p)), A373158 (with a(p)=p#).
For other completely additive sequences see the cross-references in A104244.

Programs

  • Haskell
    a001414 1 = 0
    a001414 n = sum $ a027746_row n
    -- Reinhard Zumkeller, Feb 27 2012, Nov 20 2011
    
  • Magma
    [n eq 1 select 0 else (&+[j[1]*j[2]: j in Factorization(n)]): n in [1..100]]; // G. C. Greubel, Jan 10 2019
  • Maple
    A001414 := proc(n) add( op(1,i)*op(2,i),i=ifactors(n)[2]) ; end proc:
    seq(A001414(n), n=1..100); # Peter Luschny, Jan 17 2011
  • Mathematica
    a[n_] := Plus @@ Times @@@ FactorInteger@ n; a[1] = 0; Array[a, 78] (* Ray Chandler, Nov 12 2005 *)
  • PARI
    a(n)=local(f); if(n<1,0,f=factor(n); sum(k=1,matsize(f)[1],f[k,1]*f[k,2]))
    
  • PARI
    A001414(n) = (n=factor(n))[,1]~*n[,2] \\ M. F. Hasler, Feb 07 2009
    
  • Python
    from sympy import factorint
    def A001414(n):
        return sum(p*e for p,e in factorint(n).items()) # Chai Wah Wu, Jan 08 2016
    
  • Sage
    [sum(factor(n)[j][0]*factor(n)[j][1] for j in range(0,len(factor(n)))) for n in range(1,79)] # Giuseppe Coppoletta, Jan 19 2015
    

Formula

If n = Product p_j^k_j then a(n) = Sum p_j * k_j.
Dirichlet g.f. f(s)*zeta(s), where f(s) = Sum_{p prime} p/(p^s-1) = Sum_{k>0} primezeta(k*s-1) is the Dirichlet g.f. for A120007. Totally additive with a(p^e) = p*e. - Franklin T. Adams-Watters, Jun 02 2006
For n > 1: a(n) = Sum_{k=1..A001222(n)} A027746(n,k). - Reinhard Zumkeller, Aug 27 2011
Sum_{n>=1} (-1)^a(n)/n^s = ((2^s + 1)/(2^s - 1))*zeta(2*s)/zeta(s), if Re(s)>1 and 0 if s=1 (Alladi and Erdős, 1977). - Amiram Eldar, Nov 02 2020
a(n) >= k*log(n), where k = 3/log(3). This bound is sharp. - Charles R Greathouse IV, Jul 28 2025

A000792 a(n) = max{(n - i)*a(i) : i < n}; a(0) = 1.

Original entry on oeis.org

1, 1, 2, 3, 4, 6, 9, 12, 18, 27, 36, 54, 81, 108, 162, 243, 324, 486, 729, 972, 1458, 2187, 2916, 4374, 6561, 8748, 13122, 19683, 26244, 39366, 59049, 78732, 118098, 177147, 236196, 354294, 531441, 708588, 1062882, 1594323, 2125764, 3188646, 4782969, 6377292
Offset: 0

Views

Author

Keywords

Comments

Numbers of the form 3^k, 2*3^k, 4*3^k with a(0) = 1 prepended.
If a set of positive numbers has sum n, this is the largest value of their product.
In other words, maximum of products of partitions of n: maximal value of Product k_i for any way of writing n = Sum k_i. To find the answer, take as many of the k_i's as possible to be 3 and then use one or two 2's (see formula lines below).
a(n) is also the maximal size of an Abelian subgroup of the symmetric group S_n. For example, when n = 6, one of the Abelian subgroups with maximal size is the subgroup generated by (123) and (456), which has order 9. [Bercov and Moser] - Ahmed Fares (ahmedfares(AT)my-deja.com), Apr 19 2001
Also the maximum number of maximal cliques possible in a graph with n vertices (cf. Capobianco and Molluzzo). - Felix Goldberg (felixg(AT)tx.technion.ac.il), Jul 15 2001 [Corrected by Jim Nastos and Tanya Khovanova, Mar 11 2009]
Every triple of alternate terms {3*k, 3*k+2, 3*k+4} in the sequence forms a geometric progression with first term 3^k and common ratio 2. - Lekraj Beedassy, Mar 28 2002
For n > 4, a(n) is the least multiple m of 3 not divisible by 8 for which omega(m) <= 2 and sopfr(m) = n. - Lekraj Beedassy, Apr 24 2003
Maximal number of divisors that are possible among numbers m such that A080256(m) = n. - Lekraj Beedassy, Oct 13 2003
Or, numbers of the form 2^p*3^q with p <= 2, q >= 0 and 2p + 3q = n. Largest number obtained using only the operations +,* and () on the parts 1 and 2 of any partition of n into these two summands where the former exceeds the latter. - Lekraj Beedassy, Jan 07 2005
a(n) is the largest number of complexity n in the sense of A005520 (A005245). - David W. Wilson, Oct 03 2005
a(n) corresponds also to the ultimate occurrence of n in A001414 and thus stands for the highest number m such that sopfr(m) = n, for n >= 2. - Lekraj Beedassy, Apr 29 2002
a(n) for n >= 1 is a paradigm shift sequence with procedural length p = 0, in the sense of A193455. - Jonathan T. Rowell, Jul 26 2011
a(n) = largest term of n-th row in A212721. - Reinhard Zumkeller, Jun 14 2012
For n >= 2, a(n) is the largest number whose prime divisors (with multiplicity) add to n, whereas the smallest such number (resp. smallest composite number) is A056240(n) (resp. A288814(n)). - David James Sycamore, Nov 23 2017
For n >= 3, a(n+1) = a(n)*(1 + 1/s), where s is the smallest prime factor of a(n). - David James Sycamore, Apr 10 2018

Examples

			a{8} = 18 because we have 18 = (8-5)*a(5) = 3*6 and one can verify that this is the maximum.
a(5) = 6: the 7 partitions of 5 are (5), (4, 1), (3, 2), (3, 1, 1), (2, 2, 1), (2, 1, 1, 1), (1, 1, 1, 1, 1) and the corresponding products are 5, 4, 6, 3, 4, 2 and 1; 6 is the largest.
G.f. = 1 + x + 2*x^2 + 3*x^3 + 4*x^4 + 6*x^5 + 9*x^6 + 12*x^7 + 18*x^8 + ...
		

References

  • B. R. Barwell, Cutting String and Arranging Counters, J. Rec. Math., 4 (1971), 164-168.
  • B. R. Barwell, Journal of Recreational Mathematics, "Maximum Product": Solution to Prob. 2004;25(4) 1993, Baywood, NY.
  • M. Capobianco and J. C. Molluzzo, Examples and Counterexamples in Graph Theory, p. 207. North-Holland: 1978.
  • S. L. Greitzer, International Mathematical Olympiads 1959-1977, Prob. 1976/4 pp. 18;182-3 NML vol. 27 MAA 1978
  • J. L. Gross and J. Yellen, eds., Handbook of Graph Theory, CRC Press, 2004; p. 396.
  • P. R. Halmos, Problems for Mathematicians Young and Old, Math. Assoc. Amer., 1991, pp. 30-31 and 188.
  • L. C. Larson, Problem-Solving Through Problems. Problem 1.1.4 pp. 7. Springer-Verlag 1983.
  • D. J. Newman, A Problem Seminar. Problem 15 pp. 5;15. Springer-Verlag 1982.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

See A007600 for a left inverse.
Cf. array A064364, rightmost (nonvanishing) numbers in row n >= 2.
See A056240 and A288814 for the minimal numbers whose prime factors sums up to n.
A000792, A178715, A193286, A193455, A193456, and A193457 are closely related as paradigm shift sequences for (p = 0, ..., 5 respectively).
Cf. A202337 (subsequence).

Programs

  • Haskell
    a000792 n = a000792_list !! n
    a000792_list = 1 : f [1] where
       f xs = y : f (y:xs) where y = maximum $ zipWith (*) [1..] xs
    -- Reinhard Zumkeller, Dec 17 2011
    
  • Magma
    I:=[1,1,2,3,4]; [n le 5 select I[n] else 3*Self(n-3): n in [1..45]]; // Vincenzo Librandi, Apr 14 2015
  • Maple
    A000792 := proc(n)
        m := floor(n/3) ;
        if n mod 3 = 0 then
            3^m ;
        elif n mod 3 = 1 then
            4*3^(m-1) ;
        else
            2*3^m ;
        end if;
        floor(%) ;
    end proc: # R. J. Mathar, May 26 2013
  • Mathematica
    a[1] = 1; a[n_] := 4* 3^(1/3 *(n - 1) - 1) /; (Mod[n, 3] == 1 && n > 1); a[n_] := 2*3^(1/3*(n - 2)) /; Mod[n, 3] == 2; a[n_] := 3^(n/3) /; Mod[n, 3] == 0; Table[a[n], {n, 0, 40}]
    CoefficientList[Series[(1 + x + 2x^2 + x^4)/(1 - 3x^3), {x, 0, 50}], x] (* Harvey P. Dale, May 01 2011 *)
    f[n_] := Max[ Times @@@ IntegerPartitions[n, All, Prime@ Range@ PrimePi@ n]]; f[1] = 1; Array[f, 43, 0] (* Robert G. Wilson v, Jul 31 2012 *)
    a[ n_] := If[ n < 2, Boole[ n > -1], 2^Mod[-n, 3] 3^(Quotient[ n - 1, 3] + Mod[n - 1, 3] - 1)]; (* Michael Somos, Jan 23 2014 *)
    Join[{1, 1}, LinearRecurrence[{0, 0, 3}, {2, 3, 4}, 50]] (* Jean-François Alcover, Jan 08 2019 *)
    Join[{1,1},NestList[#+Divisors[#][[-2]]&,2,41]] (* James C. McMahon, Aug 09 2024 *)
  • PARI
    {a(n) = floor( 3^(n - 4 - (n - 4) \ 3 * 2) * 2^( -n%3))}; /* Michael Somos, Jul 23 2002 */
    
  • PARI
    lista(nn) = {print1("1, 1, "); print1(a=2, ", "); for (n=1, nn, a += a/divisors(a)[2]; print1(a, ", "););} \\ Michel Marcus, Apr 14 2015
    
  • PARI
    A000792(n)=if(n>1,3^((n-2)\3)*(2+(n-2)%3),1) \\ M. F. Hasler, Jan 19 2019
    

Formula

G.f.: (1 + x + 2*x^2 + x^4)/(1 - 3*x^3). - Simon Plouffe in his 1992 dissertation.
a(3n) = 3^n; a(3*n+1) = 4*3^(n-1) for n > 0; a(3*n+2) = 2*3^n.
a(n) = 3*a(n-3) if n > 4. - Henry Bottomley, Nov 29 2001
a(n) = n if n <= 2, otherwise a(n-1) + Max{gcd(a(i), a(j)) | 0 < i < j < n}. - Reinhard Zumkeller, Feb 08 2002
A007600(a(n)) = n; Andrew Chi-Chih Yao attributes this observation to D. E. Muller. - Vincent Vatter, Apr 24 2006
a(n) = 3^(n - 2 - 2*floor((n - 1)/3))*2^(2 - (n - 1) mod 3) for n > 1. - Hieronymus Fischer, Nov 11 2007
From Kiyoshi Akima (k_akima(AT)hotmail.com), Aug 31 2009: (Start)
a(n) = 3^floor(n/3)/(1 - (n mod 3)/4), n > 1.
a(n) = 3^(floor((n - 2)/3))*(2 + ((n - 2) mod 3)), n > 1. (End)
a(n) = (2^b)*3^(C - (b + d))*(4^d), n > 1, where C = floor((n + 1)/3), b = max(0, ((n + 1) mod 3) - 1), d = max(0, 1 - ((n + 1) mod 3)). - Jonathan T. Rowell, Jul 26 2011
G.f.: 1 / (1 - x / (1 - x / (1 + x / (1 - x / (1 + x / (1 + x^2 / (1 + x))))))). - Michael Somos, May 12 2012
3*a(n) = 2*a(n+1) if n > 1 and n is not divisible by 3. - Michael Somos, Jan 23 2014
a(n) = a(n-1) + largest proper divisor of a(n-1), n > 2. - Ivan Neretin, Apr 13 2015
a(n) = max{a(i)*a(n-i) : 0 < i < n} for n >= 4. - Jianing Song, Feb 15 2020
a(n+1) = a(n) + A038754(floor( (2*(n-1) + 1)/3 )), for n > 1. - Thomas Scheuerle, Oct 27 2022

Extensions

More terms and better description from Therese Biedl (biedl(AT)uwaterloo.ca), Jan 19 2000

A288814 a(n) is the smallest composite number whose prime divisors (with multiplicity) sum to n.

Original entry on oeis.org

4, 6, 8, 10, 15, 14, 21, 28, 35, 22, 33, 26, 39, 52, 65, 34, 51, 38, 57, 76, 95, 46, 69, 92, 115, 184, 161, 58, 87, 62, 93, 124, 155, 248, 217, 74, 111, 148, 185, 82, 123, 86, 129, 172, 215, 94, 141, 188, 235, 376, 329, 106, 159, 212, 265, 424, 371, 118, 177, 122, 183, 244, 305, 488, 427, 134, 201, 268, 335, 142
Offset: 4

Views

Author

David James Sycamore, Jun 16 2017

Keywords

Comments

Agrees with A056240(n) if n is composite (but not if n is prime).
For n prime, let P_n = greatest prime < n such that A056240(n-P_n) = A288313(m) for some m; then a(n) = Min{q*a(n-q): q prime, n-1 > q >= P_n}.
In most cases q is the greatest prime < p, but there are exceptions; e.g., p=211 is the smallest prime for which q (=197) is the second prime removed from 211, not the first. 541 is the next prime with this property (q=521). The same applies to p=16183, for which q=16139, the second prime removed from p. These examples all arise with q being the lesser of a prime pair.
For p prime, a(p) = q*a(p-q) for some prime q < p as described above. Then a(p-q) = 2,4,8 or 3*r for some prime r.
The subsequence of terms (4, 6, 8, 10, 14, 21, 22, 26, 34, ...), where for all m > n, a(m) > a(n) is the same as sequence A088686, and the sequence of its indices (4, 5, 6, 7, 9, 10, 13, 19, ...) is the same as A088685. - David James Sycamore, Jun 30 2017
Records are in A088685. - Robert G. Wilson v, Feb 26 2018
Number of terms less than 10^k, k=1,2,3,...: 3, 32, 246, 2046, 17053, 147488, ..., . - Robert G. Wilson v, Feb 26 2018

Examples

			a(5) = 6 = 2*3 is the smallest composite number whose prime divisors add to 5.
a(7) = 10 = 2*5 is the smallest composite number whose prime divisors add to 7.
12 = 2 * 2 * 3 is not in the sequence, since the sum of its prime divisors is 7, a value already obtained by the lesser 10. - _David A. Corneth_, Jun 22 2017
		

Crossrefs

Programs

A064364 Positive integers sorted by A001414(n), the sum of their prime divisors, as the major key and n as the minor key.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 8, 9, 7, 10, 12, 15, 16, 18, 14, 20, 24, 27, 21, 25, 30, 32, 36, 11, 28, 40, 45, 48, 54, 35, 42, 50, 60, 64, 72, 81, 13, 22, 56, 63, 75, 80, 90, 96, 108, 33, 49, 70, 84, 100, 120, 128, 135, 144, 162, 26, 44, 105, 112, 125, 126, 150, 160, 180, 192, 216, 243
Offset: 1

Views

Author

Howard A. Landman, Sep 25 2001

Keywords

Comments

This is a permutation of the positive integers.
a(1) could be taken as 0 because 1 is not a member of A001414 and one could start with a(0)=1 (see the W. Lang link).
The row length sequence of this array is A000607(n), n>=2.
If the array is [1,0,2,3,4,5,6,6,...] with offset 0 then the row length sequence is A000607(n), n>=0.
From David James Sycamore, May 11 2018: (Start)
For n > 1, a(n) is the smallest number not yet seen such that sopfr(a(n)) is the least possible integer. The sequence lists in increasing order elements of the finite sets S(k) = {x: sopfr(x)=k}, k >= 0, where sopfr(x) = 0 iff x = 1. When a(n) = A056240(k) for some k >= 2, then sopfr(a(n)) = k and a(n) is the first of A000607(k) terms, all of which have sopfr = k. (A000607(k) is the number of partitions of k into prime parts.) Consequently the sequence follows a sawtooth profile, rising from a(n) = A056240(k) to A000792(k), the greatest number with sopfr = k, then starting over with A056240(k+1) for the next larger value of sopfr. (End) [Edited by M. F. Hasler, Jan 19 2019]

Examples

			The triangle reads:
1,
(0,) (see comment in link to "first 16 rows" by W. Lang)
2,
3,
4,
5,  6,
8,  9,
7,  10, 12,
15, 16, 18,
14, 20, 24, 27,
21, 25, 30, 32, 36,
11, 28, 40, 45, 48, 54,
35, 42, 50, 60, 64, 72, 81,
13, 22, 56, 63, 75, 80, 90, 96, 108,
...
		

Crossrefs

Cf. A001414.
Cf. A000607 (row lengths), A002098 (row sums), A056240 (least = first term in the n-th row), A000792 (greatest term in the n-th row).
Cf. A257815 (inverse).

Programs

  • Haskell
    import Data.List (partition, union)
    a064364 n k = a064364_tabf !! (n-1) !! (k-1)
    a064364_row n = a064364_tabf !! (n-1)
    a064364_tabf = [1] : tail (f 1 [] 1 (map a000792 [2..])) where
       f k pqs v (w:ws) = (map snd pqs') :
         f (k + 1) (union pqs'' (zip (map a001414 us) us )) w ws where
           us = [v + 1 .. w]
           (pqs', pqs'') = partition ((== k) . fst) pqs
    a064364_list = concat a064364_tabf
    -- Reinhard Zumkeller, Jun 11 2015
    
  • Mathematica
    terms = 1000; nmax0 = 100000 (* a rough estimate of max sopfr *);
    sopfr[n_] := sopfr[n] = Total[Times @@@ FactorInteger[n]];
    f[n1_, n2_] := Which[t1 = sopfr[n1]; t2 = sopfr[n2]; t1 < t2, True, t1 == t2, n1 <= n2, True, False];
    Clear[g];
    g[nmax_] := g[nmax] = Sort[Range[nmax], f][[1 ;; terms]];
    g[nmax = nmax0];
    g[nmax += nmax0];
    While[g[nmax] != g[nmax - nmax0], Print[nmax]; nmax += nmax0];
    A064364 = g[nmax] (* Jean-François Alcover, Mar 13 2019 *)
  • PARI
    lista(nn) = {nmax = A000792(nn); v = vector(nmax, k, A001414(k)); for (n=1, nn, vn = select(x->x==n, v, 1); for (k = 1, #vn, print1(vn[k], ", ")))} \\ Michel Marcus, May 01 2018
    
  • PARI
    A064364_vec(N, k=6, L=9)={vector(N, i, if(i<7, N=i, until(A001414(N+=1)==k, ); NA056240(k)-1))} \\ To compute terms up to a given value of k=sopfr(n) and/or for large N >> 1000, it is more efficient to use code similar to lista() above, with "for(k...)" replaced by "a=concat(a, vn)". - M. F. Hasler, Jan 19 2019

Formula

If a(n) = A056240(k) for some k then a(n+A000607(k)-1) = A000792(k). - David James Sycamore, May 11 2018

Extensions

More terms from Vladeta Jovovic, Sep 27 2005

A295185 a(n) is the smallest composite number whose prime divisors (with multiplicity) sum to prime(n); n >= 3.

Original entry on oeis.org

6, 10, 28, 22, 52, 34, 76, 184, 58, 248, 148, 82, 172, 376, 424, 118, 488, 268, 142, 584, 316, 664, 1335, 388, 202, 412, 214, 436, 3729, 508, 1048, 274, 2919, 298, 1208, 1256, 652, 1336, 1384, 358, 3801, 382, 772, 394, 6501, 7385, 892, 454, 916, 1864, 478, 5061, 2008, 2056, 2104, 538, 2168, 1108, 562, 5943, 9669
Offset: 3

Views

Author

David James Sycamore, Nov 16 2017

Keywords

Comments

Sequence is undefined for n=1,2 since no composites exist whose prime divisors sum to 2, 3. For n >= 3, a(n) = A288814(prime(n)) = prime(n-k)*B(prime(n) - prime(n-k)) where B=A056240, and k >= 1 is the "type" of prime(n), indicated as prime(n)~k(g1,g2,...,gk) where gi = prime(n-(i-1)) - prime(n-i); 1 <= i <= k. Thus: 5~1(2), 211~2(12,2), 4327~3(30,8,6) etc. The sequence relates to gaps between odd primes, and in particular to the sequence of k prime gaps below prime(n). The even-indexed terms of B are relevant, as are those of subsequences:
C=A288313, 2,4 plus terms B(n) where n-3 is prime (A298252),
D=A297150, terms B(n) where n-5 is prime and n-3 is composite (A297925) and
E=A298615, terms B(n) where both n-3 and n-5 are composite (A298366).
The above sequences of indices 2m form a partition of the even numbers and the corresponding terms B(2m) form a partition of the even-indexed terms of A056240. The union of D and E is the sequence A292081 = B-C.
Let g(n,t) = prime(n) - prime(n-t), t < n, and h(n,t) = g(n,t) - g(n,1), 1 < t < n. If g1=g(n,1) is a term in A298252 (g1-3 is prime), then B(g1) is a term in C, so k=1. If g1 belongs to A297925 or A298366 then B(g1) is a term in D or E and the value of k depends on subsequent gaps below prime(n), within a range dependent on g1.
Let range R1(g1) = u - g(n,1) where u is the index in B of the greatest term in C such that C(u) < B(g1). Let range R2(g1) = v-g(n,1) where v is the index in B of the greatest term in D such that D(v) <= B(g1). For all n, R2 < R1, and if g1 is a term in D then R2(g1)=0. Examples: R1(12)=2, R2(12)=0, R1(30)=26, R2(30)=6.
k >= 1 is the smallest integer such that B(g(n,k)) <= B(g(n,t)) for all t satisfying g1 <= g(n,t) <= g1 + R1(g1). For g1-3 prime, k=1. If g1-3 is composite, let z be least integer > 1 such that g(n,z)-3 is prime, and let w be least integer >= 1 such that g(n,w)-5 is prime. Then z "complies" if h(n,z) <= R1, and w "complies" if h(n,w) <= R2. If g1-5 is prime then R2=w=0 and only z is relevant.
B(g1) must belong to C,D or E. If in C (g1-3 is prime) then k=1. If in D (g1-5 is prime), k=z if z complies, otherwise k=1. If B(g1) is in E and z complies but not w then k=z, or if w complies but not z then k=w. If B(g1) is in E and z,w both comply then k=z if 3*(g(n,z)-3) < 5*(g(n,w)-5), otherwise k=w. If neither z nor w comply, then k=1.
Conjecture: For all n >= 3, a(n) >= A288189(n).

Examples

			5=prime(3), g(3,1)=5-3=2, a term in C; k=1, and a(3)=3*B(5-3)=3*2=6; 5~1(2).
17=prime(7), g(7,1)=17-13=4, a term in C; k=1, a(7)=13*B(17-13)=13*4=52; 17~1(4).
211=prime(47); g(47,1)=12, a term in D, R1=2, R2=0, k=z=2, a(47)=197*b(211-197)=197*33=6501; 211~2(12,2), and 211 is first prime of type k=2.
8923=prime(1109); g(1109,1)=30, a term in E. R1=26, R2=6, z=3 and w=2 both comply  but 3*(g(n,3)-3)=159 > 5*(g(n,2)-5)=155, so k=w=2. Therefore a(1109)=8887*b(8923-8887)=8887*b(36)=8887*155=1377485; 8923~2(30,6).
40343=prime(4232); g(4232,1)=54, a term in E. R1=58, R2=12,z=6 and w=3, both comply, 3*(g(n,z)-3)=309 and 5*(g(n,w)-5)=305 therefore k=w=3 and a(4232) = 40277*b(40343-40277)=40277*b(66)=40277*305=12284485; 40343~3(54,6,6).
81611=prime(7981); g(81611,1)=42, a term in D, R1=22, R2=0; z complies, k=z=6, a(7981)=81547*b(81611-81547)=81546*b(64)=81546*183=14923101; 81611~6(42,6,4,6,2,4) and is the first prime of type k=6.
If p is the greater of twin/cousin primes then p~1(2), p~1(4), respectively.
		

Crossrefs

Programs

  • Mathematica
    b[n_] := b[n] = Total[Times @@@ FactorInteger[n]];
    a[n_] := For[k = 2, True, k++, If[CompositeQ[k], If[b[k] == Prime[n], Return[k]]]];
    Table[a[n], {n, 3, 63}] (* Jean-François Alcover, Feb 23 2018 *)
  • PARI
    a(n) = { my(p=prime(n)); forcomposite(x=6, , my(f=factor(x)); if(f[, 1]~*f[, 2]==p, return(x))); } \\ Iain Fox, Dec 08 2017

Formula

a(n) = A288814(prime(n)) = prime(n-k)*A056240(prime(n) - prime(n-k)) for some k >= 1 and prime(n-k) = gpf(A288814(prime(n)).
a(n) >= A288189(n).

A297925 Even numbers k such that k - 5 is prime but k - 3 is not prime.

Original entry on oeis.org

12, 18, 24, 28, 36, 42, 48, 52, 58, 66, 72, 78, 84, 88, 94, 102, 108, 114, 118, 132, 136, 144, 156, 162, 168, 172, 178, 186, 198, 204, 216, 228, 234, 238, 246, 256, 262, 268, 276, 282, 288, 298, 312, 318, 322, 336, 342, 354, 358, 364, 372, 378, 384, 388, 394, 402, 406, 414, 426, 438, 444, 448, 454
Offset: 1

Views

Author

David James Sycamore, Jan 08 2018

Keywords

Comments

Even numbers that are the sum of 5 and another prime, but not the sum of 3 and another prime. For n >= 1, a(n) - 5 = A049591(n), a(n) - 3 = A107986(n+1).
Let r(n) = a(n) - 5, Then r(n) is the greatest prime < a(n), and therefore A056240(a(n)) = 5*r(n). Furthermore, since r(n) + 2 must be composite, A056240(a(n)) = 5*A049591(n).
The terms in this sequence, combined with those in A298366 and A298252 form a partition of A005843(n);n>=3 (nonnegative even numbers>=6). This is because any even integer n>=6 satisfies either (i) n-3 is prime, (ii) n-5 is prime but n-3 is composite, or (iii) both n-5 and n-3 are composite.

Examples

			12 is a term because 12 - 5 = 7 is prime, and 12 - 3 = 9 is composite. Also A049591(1)+5=7+5=12 and A107986(2)+3=9+3=12.
18 is a term because 18 - 5 = 13 is prime, and 18 - 3 = 15 is composite.
16 is not a term because 16 - 5 = 11 and 16 - 3 = 13 are both prime.
		

Crossrefs

Similar to A130038. Subsequence of A175222.

Programs

  • GAP
    Filtered([8..500], k-> IsPrime(k-5) and not IsPrime(k-3) and (k mod 2)=0); # G. C. Greubel, May 21 2019
  • Magma
    [n: n in [3..500] | IsPrime(n-5) and not IsPrime(n-3) and (n mod 2) eq 0]; // G. C. Greubel, May 21 2019
    
  • Maple
    N:=100
    for n from 8 to N by 2 do
    if isprime(n-5) and not isprime(n-3) then print (n);
    end if
    end do
  • Mathematica
    Select[Range[6, 500, 2], And[PrimeQ[# - 5], ! PrimeQ[# - 3]] &] (* Michael De Vlieger, Jan 10 2018 *)
    Select[Range[6, 500, 2], Boole[PrimeQ[# -{5, 3}]] == {1, 0} &] (* Harvey P. Dale, Jan 30 2024 *)
  • PARI
    isok(n) = !(n % 2) && isprime(n-5) && !isprime(n-3); \\ Michel Marcus, Jan 09 2018
    
  • Sage
    [n for n in (3..500) if is_prime(n-5) and not is_prime(n-3) and (mod(n, 2)==0)] # G. C. Greubel, May 21 2019
    

Formula

a(n) = A049591(n) + 5 = A107986(n+1) + 3 for all n >= 1.

A064502 Smallest m such that sum of distinct primes dividing m equals n, or 0 if no such number exists (as at n=1,4,6).

Original entry on oeis.org

1, 0, 2, 3, 0, 5, 0, 7, 15, 14, 21, 11, 35, 13, 33, 26, 39, 17, 65, 19, 51, 38, 57, 23, 95, 46, 69, 285, 115, 29, 161, 31, 87, 62, 93, 741, 155, 37, 217, 74, 111, 41, 185, 43, 123, 86, 129, 47, 215, 94, 141, 645, 235, 53, 329, 106, 159, 987, 265, 59, 371, 61, 177, 122
Offset: 0

Views

Author

Robert G. Wilson v, Oct 05 2001

Keywords

Comments

a(n) = first occurrence of n in A008472 (sum of prime factors of n without repetition).
Note that for all primes p, p = a(p); if n is composite then a(n) must be a composite and the only zeros are 1, 4 and 6.

Examples

			n = 217 = 3*17*197: sum = 3 + 17 + 197 = 217 = n.
		

Crossrefs

Programs

  • Mathematica
    t = Table[0, {100} ]; Do[a = Apply[Plus, Transpose[ FactorInteger[n]] [[1]]]; If[ a < 101 && t[[a]] == 0, t[[a]] = n], {n, 2, 10^5} ]; Append[t, 0]
  • PARI
    sopf(n)= { local(f,s=0); f=factor(n); for(i=1, matsize(f)[1], s+=f[i, 1]); return(s) } { for (n=1, 1000, if (n==1 || n==4 || n==6, m=0, if (isprime(n), m=n, m=1; while(sopf(m) != n, m++))); write("b064502.txt", n, " ", m) ) } \\ Harry J. Smith, Sep 16 2009

Formula

a(n) = Min{x : A008472[x]=n}.

Extensions

a(0)=1 prepended by Alois P. Heinz, May 24 2023

A298252 Even integers n such that n-3 is prime.

Original entry on oeis.org

6, 8, 10, 14, 16, 20, 22, 26, 32, 34, 40, 44, 46, 50, 56, 62, 64, 70, 74, 76, 82, 86, 92, 100, 104, 106, 110, 112, 116, 130, 134, 140, 142, 152, 154, 160, 166, 170, 176, 182, 184, 194, 196, 200, 202, 214, 226, 230, 232, 236, 242, 244, 254, 260, 266, 272, 274, 280
Offset: 1

Views

Author

David James Sycamore, Jan 15 2018

Keywords

Comments

Subsequence of A005843, same as A113935 with first term (5) excluded, since it is odd, not even. Index in A056240 of terms in A288313 (except for first two terms 2,4 of latter).
The terms in this sequence, combined with those in A297925 and A298366 form a partition of A005843(n); n>=3 (nonnegative numbers>=6). This is because any even integer n>=6 satisfies either(i) n-3 is prime, (ii) n-5 prime but n-3 composite, or (iii) n-5 and n-3 both composite.
a(n) is the smallest even number e > prime(n+1) such that e has a Goldbach partition containing prime(n+1). - Felix Fröhlich, Aug 18 2019

Examples

			a(1)=6 because 6-3=3; prime, and no smaller even number has this property; also a(1)=A113935(2)=6.  a(2)=8 because 8-3=5 is prime; also A113935(3)=8.
12 is not in the sequence because 12-3 = 9, composite.
		

Crossrefs

Programs

  • GAP
    Filtered([1..300],n->IsEvenInt(n) and IsPrime(n-3)); # Muniru A Asiru, Mar 23 2018
    
  • Magma
    [NthPrime(n+1) +3: n in [1..70]]; // G. C. Greubel, May 21 2019
    
  • Maple
    N:=200
      for n from 6 to N by 2 do
    if isprime(n-3) then print(n);
    end if
    end do
  • Mathematica
    Select[2 Range@125, PrimeQ[# - 3] &] (* Robert G. Wilson v, Jan 15 2018 *)
    Select[Prime[Range[100]]+3,EvenQ] (* Harvey P. Dale, Mar 07 2022 *)
  • PARI
    a(n) = prime(n + 1) + 3 \\ David A. Corneth, Mar 23 2018
    
  • Sage
    [nth_prime(n+1) +3 for n in (1..70)] # G. C. Greubel, May 21 2019

Formula

a(n) = A113935(n+1), n>=1.
A056240(a(n)) = A288313(n+2).
a(n) = prime(n + 1) + 3 = A113935(n + 1). - David A. Corneth, Mar 23 2018

A298366 Even numbers n such that n-5 and n-3 are both composite.

Original entry on oeis.org

30, 38, 54, 60, 68, 80, 90, 96, 98, 120, 122, 124, 126, 128, 138, 146, 148, 150, 158, 164, 174, 180, 188, 190, 192, 206, 208, 210, 212, 218, 220, 222, 224, 240, 248, 250, 252, 258, 264, 270, 278, 290, 292, 294, 300, 302, 304, 306, 308, 324, 326, 328, 330, 332, 338, 344, 346, 348, 360, 366, 368, 374, 380
Offset: 1

Views

Author

David James Sycamore, Jan 17 2018

Keywords

Comments

The sequence displays runs of consecutive even integers, whose frequency and length are related to gaps between successive primes local to these numbers. Where primes are rare (large gaps), the runs of consecutive even integers are longer (run length proportional to gap size). Let p < q be consecutive primes such that g = q-p >= 6. A string of r consecutive terms differing by 2 will start at p+7, and continue to q+1, where r = (g-4)/2. Thus at prime gap 8 a string of 2 consecutive terms differing by 2 will occur, at gap 10 there will be 3, and at gap 30 there will be 13; and so on. As the gap size increases by 2 so the run length of consecutive even terms increases by 1. The first occurrence of run length m occurs at the term corresponding to 7 + A000230(m/2).
The terms in this sequence, combined with those in A297925 and A298252 form a partition of A005843(n); n >= 3; (nonnegative even numbers >= 6). This is because any even integer n >= 6 satisfies either: (i). n-3 is prime, (ii). n-5 is prime and n-3 is composite, or (iii). both n-5 and n-3 are composite.
For any n >= 1, A056240(a(n)) = A298615(n).

Examples

			30 is included because 30-5 = 25, and 30-3 = 27; both composite, and 30 is the smallest even number with this property, hence a(1)=30. Also, A056240(a(1)) = A056240(30) = 161 = A298615(1). 24 is not included because although 24 - 3 = 21, composite; 24 - 5 = 19, prime. 210 is in this sequence, since 205 and 207 are both composite. 113 is the first prime to have a gap 14 ahead of it. Therefore we would expect a run of (14 - 4)/2 = 5 consecutive terms to start at 7 + A000230(7) = 113 + 7 = 120; thus: 120,122,124,126,128. Likewise the first occurrence of run length 7 occurs at gap m = 2*7 + 4 = 18, namely the term corresponding to 7 + A000230(9) = 523 + 7 = 530; thus: 530,532,534,536,538,540,542.
		

Crossrefs

Programs

  • Magma
    [2*n: n in [8..200] | not IsPrime(2*n-5) and not IsPrime(2*n-3)]; // Vincenzo Librandi, Nov 16 2018
    
  • Maple
    N:=300:
    for n from 8 to N by 2 do
    if not isprime(n-5) and not isprime(n-3) then print(n);
    end if
    end do
  • Mathematica
    Rest[2 Select[Range[250], !PrimeQ[2 # - 5] && !PrimeQ[2 # - 3] &]] (* Vincenzo Librandi, Nov 16 2018 *)
    Select[Range[2,400,2],AllTrue[#-{3,5},CompositeQ]&] (* Harvey P. Dale, Jul 01 2025 *)
  • PARI
    select( is_A298366(n)=!(isprime(n-5)||isprime(n-3)||bitand(n,1)||n<9), [5..200]*2) \\ Last 2 conditions aren't needed if n > 4 and even. - M. F. Hasler, Nov 19 2018 and Apr 07 2020 after edit by Michel Marcus, Apr 04 2020

Formula

a(n) = A061673(n) + 4 = A269345(n) + 5. - M. F. Hasler, Nov 19 2018

A288189 a(n) is the smallest composite number whose sum of prime divisors (with multiplicity) is divisible by prime(n).

Original entry on oeis.org

4, 8, 6, 10, 28, 22, 52, 34, 76, 184, 58, 213, 148, 82, 172, 309, 424, 118, 393, 268, 142, 584, 316, 664, 573, 388, 202, 412, 214, 436, 753, 508, 813, 274, 1465, 298, 933, 974, 652, 1336, 1384, 358, 1137, 382, 772, 394, 1257, 1329, 892, 454, 916, 1864, 478, 1497, 1538, 1569
Offset: 1

Views

Author

David James Sycamore, Jul 01 2017

Keywords

Comments

In most cases a(n) = A288814(prime(n)) but there are exceptions, e.g., a(37)=213, whereas A288814(37)=248. Other exceptions include a(53), a(67), a(127), a(137), etc. These examples occur when there is a number r such that A001414(r*p) is less than A288814(p).
The strictly increasing subsequence of terms (10, 22, 34, 58, 82, 118, 142, 202, 214, 274, 298, ...) where for all m>n, a(m)>a(n) gives the semiprimes with prime sum of prime factors, A108605. The sequence of the indices of this subsequence (5, 7, 13, 19, 31, 43, 61, 73, 103, 109, 139, 151, ...) gives the greater of twin primes, A006512.

Examples

			a(5)=6 because 6 = 2*3 is the smallest number whose sum of prime divisors (2+3 = 5) is divisible by 5.
a(37) = 213 = A288814(74) = A288814(2*37).
		

Crossrefs

Programs

Previous Showing 11-20 of 37 results. Next