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 61-70 of 485 results. Next

A062069 a(n) = sigma(d(n)), where d(k) is the number of divisors function (A000005) and sigma(k) is the sum of divisors function (A000203).

Original entry on oeis.org

1, 3, 3, 4, 3, 7, 3, 7, 4, 7, 3, 12, 3, 7, 7, 6, 3, 12, 3, 12, 7, 7, 3, 15, 4, 7, 7, 12, 3, 15, 3, 12, 7, 7, 7, 13, 3, 7, 7, 15, 3, 15, 3, 12, 12, 7, 3, 18, 4, 12, 7, 12, 3, 15, 7, 15, 7, 7, 3, 28, 3, 7, 12, 8, 7, 15, 3, 12, 7, 15, 3, 28, 3, 7, 12, 12, 7, 15, 3, 18, 6, 7, 3, 28, 7, 7, 7, 15, 3
Offset: 1

Views

Author

Amarnath Murthy, Jun 13 2001

Keywords

Comments

a(1) = 1, a(p) = 3 for p = primes (A000040), a(pq) = 7 for pq = product of two distinct primes (A006881), a(pq...z) = 2^(k+1)-1 = A000225(k+1) for pq...z = product of k (k > 2) distinct primes p,q,...,z (A120944), a(p^k) = sigma(k+1) = A000203(k+1) for p^k = prime powers (A000961(n) for n > 1). Sequence {1,3,4,12} is finite sequence of numbers n such that sigma(tau(n)) = n. [Jaroslav Krizek, Jul 16 2009]
For semiprime n, a(n) is either 4 or 7. Also a(n) = d(n) + omega(n) + mu(n), the sum of three core sequences A000005, A001221 and A008683. When n is semiprime, a(n) is completely defined by the Mobius function as: a(n) = 4 + 3*mu(n). a(n) also has the fractal-like identities a(d(n)) = d(n) and a(n) = sigma(a(d(n))). - Wesley Ivan Hurt, Sep 02 2013
If n is a triprime (A014612), d(n) is 4, 6, or 8 and a(n) = sigma(d(n)) is 7, 12, or 15 respectively. Then a(n) = -d(n)^2/4 + 5*d(n) - 9. - Wesley Ivan Hurt, Sep 08 2013

Examples

			sigma(d(12)) = sigma(6) = 12.
		

Crossrefs

Programs

  • Maple
    A062069:= (n-> numtheory[sigma](numtheory[tau](n))):
    seq (A062069(n), n=1..40); # Jani Melik, Jan 25 2011
  • Mathematica
    Table[DivisorSigma[1, DivisorSigma[0, n]], {n, 1, 80}] (* Carl Najafi, Aug 16 2011 *)
  • PARI
    v=[]; for(n=1,150,v=concat(v, sigma(numdiv(n)))); v
    
  • PARI
    { for (n=1, 1000, write("b062069.txt", n, " ", sigma(numdiv(n))) ) } \\ Harry J. Smith, Jul 31 2009

Formula

a(n) = A000203(A000005(n)). - Wesley Ivan Hurt, Sep 09 2013

Extensions

More terms from Jason Earls, Jun 19 2001

A338904 Irregular triangle read by rows where row n lists all semiprimes whose prime indices sum to n.

Original entry on oeis.org

4, 6, 9, 10, 14, 15, 21, 22, 25, 26, 33, 35, 34, 39, 49, 55, 38, 51, 65, 77, 46, 57, 85, 91, 121, 58, 69, 95, 119, 143, 62, 87, 115, 133, 169, 187, 74, 93, 145, 161, 209, 221, 82, 111, 155, 203, 247, 253, 289, 86, 123, 185, 217, 299, 319, 323, 94, 129, 205
Offset: 2

Views

Author

Gus Wiseman, Nov 28 2020

Keywords

Comments

A semiprime is a product of any two prime numbers. A prime index of n is a number m such that the m-th prime number divides n. The multiset of prime indices of n is row n of A112798.

Examples

			Triangle begins:
   4
   6
   9  10
  14  15
  21  22  25
  26  33  35
  34  39  49  55
  38  51  65  77
  46  57  85  91 121
  58  69  95 119 143
  62  87 115 133 169 187
  74  93 145 161 209 221
  82 111 155 203 247 253 289
  86 123 185 217 299 319 323
  94 129 205 259 341 361 377 391
		

Crossrefs

A004526 gives row lengths.
A024697 gives row sums.
A087112 is a different triangle of semiprimes.
A098350 has antidiagonals with the same distinct terms as these rows.
A338905 is the squarefree case, with row sums A025129.
A338907/A338906 are the union of odd/even rows.
A339114/A339115 are the row minima/maxima.
A001358 lists semiprimes, with odd/even terms A046315/A100484.
A006881 lists squarefree semiprimes, with odd/even terms A046388/A100484.
A014342 is the self-convolution of primes.
A037143 lists primes and semiprimes.
A056239 gives sum of prime indices (Heinz weight).
A062198 gives partial sums of semiprimes.
A084126 and A084127 give the prime factors of semiprimes.
A289182/A115392 list the positions of odd/even terms in A001358.
A332765 gives the greatest squarefree semiprime of weight n.
A338898, A338912, and A338913 give the prime indices of semiprimes, with product A087794, sum A176504, and difference A176506.
A338899, A270650, and A270652 give the prime indices of squarefree semiprimes, with difference A338900.

Programs

  • Mathematica
    Table[Sort[Table[Prime[k]*Prime[n-k],{k,n/2}]],{n,2,10}]

A167171 Squarefree semiprimes together with primes.

Original entry on oeis.org

2, 3, 5, 6, 7, 10, 11, 13, 14, 15, 17, 19, 21, 22, 23, 26, 29, 31, 33, 34, 35, 37, 38, 39, 41, 43, 46, 47, 51, 53, 55, 57, 58, 59, 61, 62, 65, 67, 69, 71, 73, 74, 77, 79, 82, 83, 85, 86, 87, 89, 91, 93, 94, 95, 97, 101, 103, 106, 107, 109, 111, 113, 115, 118, 119, 122, 123
Offset: 1

Views

Author

Juri-Stepan Gerasimov, Oct 29 2009

Keywords

Comments

Numbers such that d(n)=2*omega(n), where d = A000005 is the number of divisors.
Numbers n such that half of number of divisors of n is equal to number of distinct primes dividing n.
Numbers p*q such that p is 1 or a prime and q is a prime greater than p.

Examples

			a(1)=2 (d(2)=2*omega(2)); a(2)=3 (d(3)=2*omega(3)).
		

Crossrefs

Programs

  • Maple
    omega := proc(n) if n = 1 then 0 ; else nops( numtheory[factorset](n)) ; end if; end proc: isA167171 := proc(n) numtheory[tau](n) = 2*omega(n) ; end proc: for n from 1 to 300 do if isA167171(n) then printf("%d,",n) ; end if ; end do: # R. J. Mathar, Oct 31 2009
  • Mathematica
    a = {}; Do[If[1 <= PrimeOmega[n] <= 2 && SquareFreeQ[n], AppendTo[a, n]], {n, 123}]; a (* L. Edson Jeffery, Jan 01 2015 *)
  • PARI
    for(n=1, 1e3, if(numdiv(n)==2*omega(n), print1(n, ", "))) \\ Felix Fröhlich, Aug 11 2014

Formula

Equals A037143 \ A000290 = A006881 union A000040. - V. Raman, Sep 13 2012
a(n) ~ n log n/log log n. - Charles R Greathouse IV, Apr 05 2017

Extensions

Corrected by R. J. Mathar, Oct 31 2009
New name from Charles R Greathouse IV, Apr 05 2017

A338907 Semiprimes whose prime indices sum to an odd number.

Original entry on oeis.org

6, 14, 15, 26, 33, 35, 38, 51, 58, 65, 69, 74, 77, 86, 93, 95, 106, 119, 122, 123, 141, 142, 143, 145, 158, 161, 177, 178, 185, 201, 202, 209, 214, 215, 217, 219, 221, 226, 249, 262, 265, 278, 287, 291, 299, 302, 305, 309, 319, 323, 326, 327, 329, 346, 355
Offset: 1

Views

Author

Gus Wiseman, Nov 28 2020

Keywords

Comments

All terms are squarefree (A005117).
A semiprime is a product of any two prime numbers. A prime index of n is a number m such that the m-th prime number divides n. The multiset of prime indices of n is row n of A112798.
The semiprimes in A300063; the semiprimes in A332820. - Peter Munn, Dec 25 2020

Examples

			The sequence of terms together with their prime indices begins:
      6: {1,2}      95: {3,8}     202: {1,26}
     14: {1,4}     106: {1,16}    209: {5,8}
     15: {2,3}     119: {4,7}     214: {1,28}
     26: {1,6}     122: {1,18}    215: {3,14}
     33: {2,5}     123: {2,13}    217: {4,11}
     35: {3,4}     141: {2,15}    219: {2,21}
     38: {1,8}     142: {1,20}    221: {6,7}
     51: {2,7}     143: {5,6}     226: {1,30}
     58: {1,10}    145: {3,10}    249: {2,23}
     65: {3,6}     158: {1,22}    262: {1,32}
     69: {2,9}     161: {4,9}     265: {3,16}
     74: {1,12}    177: {2,17}    278: {1,34}
     77: {4,5}     178: {1,24}    287: {4,13}
     86: {1,14}    185: {3,12}    291: {2,25}
     93: {2,11}    201: {2,19}    299: {6,9}
		

Crossrefs

A031368 looks at primes instead of semiprimes.
A098350 has this as union of odd-indexed antidiagonals.
A300063 looks at all numbers (not just semiprimes).
A338904 has this as union of odd-indexed rows.
A338906 is the even version.
A001358 lists semiprimes, with odd/even terms A046315/A100484.
A006881 lists squarefree semiprimes, with odd/even terms A046388/A100484.
A056239 gives the sum of prime indices (Heinz weight).
A084126 and A084127 give the prime factors of semiprimes.
A087112 groups semiprimes by greater factor.
A289182/A115392 list the positions of odd/even terms in A001358.
A338898, A338912, and A338913 give the prime indices of semiprimes, with product A087794, sum A176504, and difference A176506.
A338899, A270650, and A270652 give the prime indices of squarefree semiprimes, with difference A338900.
A338908 lists squarefree semiprimes of even weight.
A339114/A339115 give the least/greatest semiprime of weight n.
Subsequence of A332820.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],PrimeOmega[#]==2&&OddQ[Total[primeMS[#]]]&]
  • Python
    from math import isqrt
    from sympy import primepi, primerange
    def A338907(n):
        def bisection(f,kmin=0,kmax=1):
            while f(kmax) > kmax: kmax <<= 1
            kmin = kmax >> 1
            while kmax-kmin > 1:
                kmid = kmax+kmin>>1
                if f(kmid) <= kmid:
                    kmax = kmid
                else:
                    kmin = kmid
            return kmax
        def f(x): return n+x-sum((primepi(x//p)-a>>1) for a,p in enumerate(primerange(isqrt(x)+1)))
        return bisection(f,n,n) # Chai Wah Wu, Apr 03 2025

Formula

Complement of A338906 in A001358.

A339560 Number of integer partitions of n that can be partitioned into distinct pairs of distinct parts, i.e., into a set of edges.

Original entry on oeis.org

1, 0, 0, 1, 1, 2, 2, 4, 5, 8, 8, 13, 17, 22, 28, 39, 48, 62, 81, 101, 127, 167, 202, 253, 318, 395, 486, 608, 736, 906, 1113, 1353, 1637, 2011, 2409, 2922, 3510, 4227, 5060, 6089, 7242, 8661, 10306, 12251, 14503, 17236, 20345, 24045, 28334, 33374, 39223, 46076
Offset: 0

Views

Author

Gus Wiseman, Dec 10 2020

Keywords

Comments

Naturally, such a partition must have an even number of parts. Its multiplicities form a graphical partition (A000569, A320922), and vice versa.

Examples

			The a(3) = 1 through a(11) = 13 partitions (A = 10):
  (21)  (31)  (32)  (42)  (43)    (53)    (54)    (64)    (65)
              (41)  (51)  (52)    (62)    (63)    (73)    (74)
                          (61)    (71)    (72)    (82)    (83)
                          (3211)  (3221)  (81)    (91)    (92)
                                  (4211)  (3321)  (4321)  (A1)
                                          (4221)  (5221)  (4322)
                                          (4311)  (5311)  (4331)
                                          (5211)  (6211)  (4421)
                                                          (5321)
                                                          (5411)
                                                          (6221)
                                                          (6311)
                                                          (7211)
For example, the partition y = (4,3,3,2,1,1) can be partitioned into a set of edges in two ways:
  {{1,2},{1,3},{3,4}}
  {{1,3},{1,4},{2,3}},
so y is counted under a(14).
		

Crossrefs

A338916 allows equal pairs (x,x).
A339559 counts the complement in even-length partitions.
A339561 gives the Heinz numbers of these partitions.
A339619 counts factorizations of the same type.
A000070 counts non-multigraphical partitions of 2n, ranked by A339620.
A000569 counts graphical partitions, ranked by A320922.
A001358 lists semiprimes, with squarefree case A006881.
A002100 counts partitions into squarefree semiprimes.
A058696 counts partitions of even numbers, ranked by A300061.
A209816 counts multigraphical partitions, ranked by A320924.
A320655 counts factorizations into semiprimes.
A320656 counts factorizations into squarefree semiprimes.
A339617 counts non-graphical partitions of 2n, ranked by A339618.
A339655 counts non-loop-graphical partitions of 2n, ranked by A339657.
A339656 counts loop-graphical partitions, ranked by A339658.
A339659 counts graphical partitions of 2n into k parts.
The following count partitions of even length and give their Heinz numbers:
- A027187 has no additional conditions (A028260).
- A096373 cannot be partitioned into strict pairs (A320891).
- A338914 can be partitioned into strict pairs (A320911).
- A338915 cannot be partitioned into distinct pairs (A320892).
- A338916 can be partitioned into distinct pairs (A320912).
- A339559 cannot be partitioned into distinct strict pairs (A320894).

Programs

  • Mathematica
    strs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[strs[n/d],Min@@#>d&]],{d,Select[Rest[Divisors[n]],And[SquareFreeQ[#],PrimeOmega[#]==2]&]}]];
    Table[Length[Select[IntegerPartitions[n],strs[Times@@Prime/@#]!={}&]],{n,0,15}]

Formula

A027187(n) = a(n) + A339559(n).

Extensions

More terms from Jinyuan Wang, Feb 14 2025

A339561 Products of distinct squarefree semiprimes.

Original entry on oeis.org

1, 6, 10, 14, 15, 21, 22, 26, 33, 34, 35, 38, 39, 46, 51, 55, 57, 58, 60, 62, 65, 69, 74, 77, 82, 84, 85, 86, 87, 90, 91, 93, 94, 95, 106, 111, 115, 118, 119, 122, 123, 126, 129, 132, 133, 134, 140, 141, 142, 143, 145, 146, 150, 155, 156, 158, 159, 161, 166
Offset: 1

Views

Author

Gus Wiseman, Dec 13 2020

Keywords

Comments

First differs from A320911 in lacking 36.
A squarefree semiprime (A006881) is a product of any two distinct prime numbers.
The following are equivalent characteristics for any positive integer n:
(1) the prime factors of n can be partitioned into distinct strict pairs (a set of edges);
(2) n can be factored into distinct squarefree semiprimes;
(3) the prime signature of n is graphical.

Examples

			The sequence of terms together with their prime indices begins:
      1: {}        55: {3,5}         91: {4,6}
      6: {1,2}     57: {2,8}         93: {2,11}
     10: {1,3}     58: {1,10}        94: {1,15}
     14: {1,4}     60: {1,1,2,3}     95: {3,8}
     15: {2,3}     62: {1,11}       106: {1,16}
     21: {2,4}     65: {3,6}        111: {2,12}
     22: {1,5}     69: {2,9}        115: {3,9}
     26: {1,6}     74: {1,12}       118: {1,17}
     33: {2,5}     77: {4,5}        119: {4,7}
     34: {1,7}     82: {1,13}       122: {1,18}
     35: {3,4}     84: {1,1,2,4}    123: {2,13}
     38: {1,8}     85: {3,7}        126: {1,2,2,4}
     39: {2,6}     86: {1,14}       129: {2,14}
     46: {1,9}     87: {2,10}       132: {1,1,2,5}
     51: {2,7}     90: {1,2,2,3}    133: {4,8}
For example, the number 1260 can be factored into distinct squarefree semiprimes in two ways, (6*10*21) or (6*14*15), so 1260 is in the sequence. The number 69300 can be factored into distinct squarefree semiprimes in seven ways:
  (6*10*15*77)
  (6*10*21*55)
  (6*10*33*35)
  (6*14*15*55)
  (6*15*22*35)
  (10*14*15*33)
  (10*15*21*22),
so 69300 is in the sequence. A complete list of all strict factorizations of 24 is: (2*3*4), (2*12), (3*8), (4*6), (24), all of which contain at least one number that is not a squarefree semiprime, so 24 is not in the sequence.
		

Crossrefs

A309356 is a kind of universal embedding.
A320894 is the complement in A028260.
A320911 lists all (not just distinct) products of squarefree semiprimes.
A339560 counts the partitions with these Heinz numbers.
A339661 has nonzero terms at these positions.
A001358 lists semiprimes, with squarefree case A006881.
A005117 lists squarefree numbers.
A320656 counts factorizations into squarefree semiprimes.
The following count vertex-degree partitions and give their Heinz numbers:
- A058696 counts partitions of 2n (A300061).
- A000070 counts non-multigraphical partitions of 2n (A339620).
- A209816 counts multigraphical partitions (A320924).
- A320921 counts connected graphical partitions (A320923).
- A339655 counts non-loop-graphical partitions of 2n (A339657).
- A339656 counts loop-graphical partitions (A339658).
- A339617 counts non-graphical partitions of 2n (A339618).
- A000569 counts graphical partitions (A320922).
The following count partitions of even length and give their Heinz numbers:
- A027187 has no additional conditions (A028260).
- A096373 cannot be partitioned into strict pairs (A320891).
- A338914 can be partitioned into strict pairs (A320911).
- A338915 cannot be partitioned into distinct pairs (A320892).
- A338916 can be partitioned into distinct pairs (A320912).
- A339559 cannot be partitioned into distinct strict pairs (A320894).
- A339560 can be partitioned into distinct strict pairs (A339561 [this sequence]).

Programs

  • Mathematica
    sqs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[sqs[n/d],Min@@#>d&]],{d,Select[Divisors[n],SquareFreeQ[#]&&PrimeOmega[#]==2&]}]];
    Select[Range[100],sqs[#]!={}&]

Formula

A339741 Products of distinct primes or squarefree semiprimes.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 10, 11, 12, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 26, 28, 29, 30, 31, 33, 34, 35, 36, 37, 38, 39, 41, 42, 43, 44, 45, 46, 47, 50, 51, 52, 53, 55, 57, 58, 59, 60, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 73, 74, 75, 76, 77, 78, 79, 82, 83, 84
Offset: 1

Views

Author

Gus Wiseman, Dec 23 2020

Keywords

Comments

First differs from A212167 in lacking 1080, with prime indices {1,1,1,2,2,2,3}.
First differs from A335433 in lacking 72 (see example).
A squarefree semiprime (A006881) is a product of any two distinct prime numbers.
The following are equivalent characteristics for any positive integer n:
(1) the prime factors of n can be partitioned into distinct singletons and strict pairs, i.e., into a set of half-loops and edges;
(2) n can be factored into distinct primes or squarefree semiprimes;
(3) the prime signature of n is half-loop-graphical.

Examples

			The sequence of terms together with their prime indices begins:
       1: {}           20: {1,1,3}        39: {2,6}
       2: {1}          21: {2,4}          41: {13}
       3: {2}          22: {1,5}          42: {1,2,4}
       5: {3}          23: {9}            43: {14}
       6: {1,2}        26: {1,6}          44: {1,1,5}
       7: {4}          28: {1,1,4}        45: {2,2,3}
      10: {1,3}        29: {10}           46: {1,9}
      11: {5}          30: {1,2,3}        47: {15}
      12: {1,1,2}      31: {11}           50: {1,3,3}
      13: {6}          33: {2,5}          51: {2,7}
      14: {1,4}        34: {1,7}          52: {1,1,6}
      15: {2,3}        35: {3,4}          53: {16}
      17: {7}          36: {1,1,2,2}      55: {3,5}
      18: {1,2,2}      37: {12}           57: {2,8}
      19: {8}          38: {1,8}          58: {1,10}
For example, we have 36 = (2*3*6), so 36 is in the sequence. On the other hand, a complete list of all strict factorizations of 72 is: (2*3*12), (2*4*9), (2*36), (3*4*6), (3*24), (4*18), (6*12), (8*9), (72). Since none of these consists of only primes or squarefree semiprimes, 72 is not in the sequence. A complete list of all factorizations of 1080 into primes or squarefree semiprimes is:
  (2*2*2*3*3*3*5)
  (2*2*2*3*3*15)
  (2*2*3*3*3*10)
  (2*2*3*3*5*6)
  (2*2*3*6*15)
  (2*3*3*6*10)
  (2*3*5*6*6)
  (2*6*6*15)
  (3*6*6*10)
  (5*6*6*6)
Since none of these is strict, 1080 is not in the sequence.
		

Crossrefs

See link for additional cross-references.
Allowing only primes gives A013929.
Not allowing primes gives A339561.
Complement of A339740.
Positions of positive terms in A339742.
Allowing squares of primes gives the complement of A339840.
Unlabeled multiset partitions of this type are counted by A339888.
A001055 counts factorizations.
A001358 lists semiprimes, with squarefree case A006881.
A002100 counts partitions into squarefree semiprimes.
A339841 have exactly one factorization into primes or semiprimes.

Programs

  • Mathematica
    sqps[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[sqps[n/d],Min@@#>d&]],{d,Select[Divisors[n],PrimeQ[#]||SquareFreeQ[#]&&PrimeOmega[#]==2&]}]];
    Select[Range[100],sqps[#]!={}&]

A080257 Numbers having at least two distinct or a total of at least three prime factors.

Original entry on oeis.org

6, 8, 10, 12, 14, 15, 16, 18, 20, 21, 22, 24, 26, 27, 28, 30, 32, 33, 34, 35, 36, 38, 39, 40, 42, 44, 45, 46, 48, 50, 51, 52, 54, 55, 56, 57, 58, 60, 62, 63, 64, 65, 66, 68, 69, 70, 72, 74, 75, 76, 77, 78, 80, 81, 82, 84, 85, 86, 87, 88, 90, 91, 92, 93, 94, 95, 96, 98, 99, 100
Offset: 1

Views

Author

Reinhard Zumkeller, Feb 10 2003

Keywords

Comments

Complement of A000430; A080256(a(n)) > 3.
A084114(a(n)) > 0, see also A084110.
Also numbers greater than the square of their smallest prime-factor: a(n)>A020639(a(n))^2=A088377(a(n));
a(n)>A000430(k) for n<=13, a(n) < A000430(k) for n>13.
Numbers with at least 4 divisors. - Franklin T. Adams-Watters, Jul 28 2006
Union of A024619 and A033942; A211110(a(n)) > 2. - Reinhard Zumkeller, Apr 02 2012
Also numbers > 1 that are neither prime nor a square of a prime. Also numbers whose omega-sequence (A323023) has sum > 3. Numbers with omega-sequence summing to m are: A000040 (m = 1), A001248 (m = 3), A030078 (m = 4), A068993 (m = 5), A050997 (m = 6), A325264 (m = 7). - Gus Wiseman, Jul 03 2019
Numbers n such that sigma_2(n)*tau(n) = A001157(n)*A000005(n) >= 4*n^2. Note that sigma_2(n)*tau(n) >= sigma(n)^2 = A072861 for all n. - Joshua Zelinsky, Jan 23 2025

Examples

			8=2*2*2 and 10=2*5 are terms; 4=2*2 is not a term.
From _Gus Wiseman_, Jul 03 2019: (Start)
The sequence of terms together with their prime indices begins:
   6: {1,2}
   8: {1,1,1}
  10: {1,3}
  12: {1,1,2}
  14: {1,4}
  15: {2,3}
  16: {1,1,1,1}
  18: {1,2,2}
  20: {1,1,3}
  21: {2,4}
  22: {1,5}
  24: {1,1,1,2}
  26: {1,6}
  27: {2,2,2}
  28: {1,1,4}
  30: {1,2,3}
  32: {1,1,1,1,1}
(End)
		

Crossrefs

Programs

  • Haskell
    a080257 n = a080257_list !! (n-1)
    a080257_list = m a024619_list a033942_list where
       m xs'@(x:xs) ys'@(y:ys) | x < y  = x : m xs ys'
                               | x == y = x : m xs ys
                               | x > y  = y : m xs' ys
    -- Reinhard Zumkeller, Apr 02 2012
    
  • Mathematica
    Select[Range[100],PrimeNu[#]>1||PrimeOmega[#]>2&] (* Harvey P. Dale, Jul 23 2013 *)
  • PARI
    is(n)=omega(n)>1 || isprimepower(n)>2
    
  • PARI
    is(n)=my(k=isprimepower(n)); if(k, k>2, !isprime(n)) \\ Charles R Greathouse IV, Jan 23 2025

Formula

a(n) = n + O(n/log n). - Charles R Greathouse IV, Sep 14 2015

Extensions

Definition clarified by Harvey P. Dale, Jul 23 2013

A129521 Numbers of the form p*q, p and q prime with q=2*p-1.

Original entry on oeis.org

6, 15, 91, 703, 1891, 2701, 12403, 18721, 38503, 49141, 79003, 88831, 104653, 146611, 188191, 218791, 226801, 269011, 286903, 385003, 497503, 597871, 665281, 721801, 736291, 765703, 873181, 954271, 1056331, 1314631, 1373653, 1537381
Offset: 1

Views

Author

Reinhard Zumkeller, Apr 19 2007

Keywords

Comments

All terms are Fermat 4-pseudoprimes, i.e., satisfy 4^n == 4 (mod n). See A020136 and A122781.

Crossrefs

Subsequence of A006881, A129510, and A122781.
Intersection of A000384 and A001358, "hexagonal semiprimes". - Wesley Ivan Hurt, Jul 04 2013

Programs

  • Haskell
    a129521 n = p * (2 * p - 1) where p = a005382 n
    -- Reinhard Zumkeller, Nov 10 2013
  • Magma
    [2*n^2-n: n in [0..1000]|IsPrime(n) and IsPrime(2*n-1)]; // Vincenzo Librandi, Dec 27 2010
    
  • Mathematica
    p = Select[Prime[Range[155]], PrimeQ[2# - 1] &]; p (2p - 1) (* Robert G. Wilson v, Sep 11 2011 *)
  • PARI
    forprime(p=2,10000,q=2*p-1;if(isprime(q),print1(p*q,", ")))
    

Formula

a(n) = A005382(n)*A005383(n).

A320893 Numbers with an even number of prime factors (counted with multiplicity) that can be factored into squarefree semiprimes (A320911) but cannot be factored into distinct semiprimes (A320892).

Original entry on oeis.org

1296, 7776, 10000, 12960, 18144, 19440, 21600, 27216, 28512, 33696, 36000, 38416, 42336, 42768, 44064, 46656, 48600, 49248, 50544, 50625, 59616, 60000, 66096, 73872, 75168, 77760, 80352, 89424, 95256, 95904, 98784, 100000
Offset: 1

Views

Author

Gus Wiseman, Oct 23 2018

Keywords

Comments

A semiprime (A001358) is a product of any two not necessarily distinct primes.

Crossrefs

Programs

  • Mathematica
    sqfsemfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[sqfsemfacs[n/d],Min@@#>=d&]],{d,Select[Rest[Divisors[n]],And[SquareFreeQ[#],PrimeOmega[#]==2]&]}]];
    strsemfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[strsemfacs[n/d],Min@@#>d&]],{d,Select[Rest[Divisors[n]],PrimeOmega[#]==2&]}]];
    Select[Range[10000],And[EvenQ[PrimeOmega[#]],strsemfacs[#]=={},sqfsemfacs[#]!={}]&]
Previous Showing 61-70 of 485 results. Next