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

A338899 Concatenated sequence of prime indices of squarefree semiprimes (A006881).

Original entry on oeis.org

1, 2, 1, 3, 1, 4, 2, 3, 2, 4, 1, 5, 1, 6, 2, 5, 1, 7, 3, 4, 1, 8, 2, 6, 1, 9, 2, 7, 3, 5, 2, 8, 1, 10, 1, 11, 3, 6, 2, 9, 1, 12, 4, 5, 1, 13, 3, 7, 1, 14, 2, 10, 4, 6, 2, 11, 1, 15, 3, 8, 1, 16, 2, 12, 3, 9, 1, 17, 4, 7, 1, 18, 2, 13, 2, 14, 4, 8, 1, 19, 2, 15
Offset: 1

Views

Author

Gus Wiseman, Nov 16 2020

Keywords

Comments

This is a triangle with two columns and strictly increasing rows, namely {A270650(n), A270652(n)}.
A squarefree semiprime is a product of any two distinct 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

			The sequence of terms together with their prime indices begins:
      6: {1,2}     57: {2,8}     106: {1,16}    155: {3,11}
     10: {1,3}     58: {1,10}    111: {2,12}    158: {1,22}
     14: {1,4}     62: {1,11}    115: {3,9}     159: {2,16}
     15: {2,3}     65: {3,6}     118: {1,17}    161: {4,9}
     21: {2,4}     69: {2,9}     119: {4,7}     166: {1,23}
     22: {1,5}     74: {1,12}    122: {1,18}    177: {2,17}
     26: {1,6}     77: {4,5}     123: {2,13}    178: {1,24}
     33: {2,5}     82: {1,13}    129: {2,14}    183: {2,18}
     34: {1,7}     85: {3,7}     133: {4,8}     185: {3,12}
     35: {3,4}     86: {1,14}    134: {1,19}    187: {5,7}
     38: {1,8}     87: {2,10}    141: {2,15}    194: {1,25}
     39: {2,6}     91: {4,6}     142: {1,20}    201: {2,19}
     46: {1,9}     93: {2,11}    143: {5,6}     202: {1,26}
     51: {2,7}     94: {1,15}    145: {3,10}    203: {4,10}
     55: {3,5}     95: {3,8}     146: {1,21}    205: {3,13}
		

Crossrefs

A270650 is the first column.
A270652 is the second column.
A320656 counts multiset partitions using these rows, or factorizations into squarefree semiprimes.
A338898 is the version including squares, with columns A338912 and A338913.
A338900 gives row differences.
A338901 gives the row numbers for first appearances.
A001221 and A001222 count distinct/all prime indices.
A001358 lists semiprimes.
A004526 counts 2-part partitions, with strict case shifted right once.
A005117 lists squarefree numbers.
A006881 lists squarefree semiprimes.
A046315 and A100484 list odd and even semiprimes.
A046388 lists odd squarefree semiprimes.
A166237 gives first differences of squarefree semiprimes.

Programs

  • Mathematica
    Join@@Cases[Select[Range[100],SquareFreeQ[#]&&PrimeOmega[#]==2&],k_:>PrimePi/@First/@FactorInteger[k]]

A338913 Greater prime index of the n-th semiprime.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Nov 20 2020

Keywords

Comments

A semiprime is a product of any two prime numbers. A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
After the first three terms, there appear to be no adjacent equal terms.

Examples

			The semiprimes are:
  2*2, 2*3, 3*3, 2*5, 2*7, 3*5, 3*7, 2*11, 5*5, 2*13, ...
so the greater prime factors are:
  2, 3, 3, 5, 7, 5, 7, 11, 5, 13, ...
with indices:
  1, 2, 2, 3, 4, 3, 4, 5, 3, 6, ...
		

Crossrefs

A115392 lists positions of first appearances of each positive integer.
A270652 is the squarefree case, with lesser part A270650.
A338898 has this as second column.
A338912 is the corresponding lesser prime index.
A001221 counts distinct prime indices.
A001222 counts prime indices.
A001358 lists semiprimes, with odd/even terms A046315/A100484.
A006881 lists squarefree semiprimes, with odd/even terms A046388/A100484.
A087794/A176504/A176506 are product/sum/difference of semiprime indices.
A338910/A338911 list products of pairs of odd/even-indexed primes.

Programs

  • Mathematica
    Table[Max[PrimePi/@First/@FactorInteger[n]],{n,Select[Range[100],PrimeOmega[#]==2&]}]

Formula

a(n) = A000720(A084127(n)).

A338898 Concatenated sequence of prime indices of semiprimes (A001358).

Original entry on oeis.org

1, 1, 1, 2, 2, 2, 1, 3, 1, 4, 2, 3, 2, 4, 1, 5, 3, 3, 1, 6, 2, 5, 1, 7, 3, 4, 1, 8, 2, 6, 1, 9, 4, 4, 2, 7, 3, 5, 2, 8, 1, 10, 1, 11, 3, 6, 2, 9, 1, 12, 4, 5, 1, 13, 3, 7, 1, 14, 2, 10, 4, 6, 2, 11, 1, 15, 3, 8, 1, 16, 2, 12, 3, 9, 1, 17, 4, 7, 5, 5, 1, 18, 2
Offset: 1

Views

Author

Gus Wiseman, Nov 15 2020

Keywords

Comments

This is a triangle with two columns and weakly increasing rows, namely {A338912(n), A338913(n)}.
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

			The sequence of semiprimes together with their prime indices begins:
      4: {1,1}     46: {1,9}      91: {4,6}     141: {2,15}
      6: {1,2}     49: {4,4}      93: {2,11}    142: {1,20}
      9: {2,2}     51: {2,7}      94: {1,15}    143: {5,6}
     10: {1,3}     55: {3,5}      95: {3,8}     145: {3,10}
     14: {1,4}     57: {2,8}     106: {1,16}    146: {1,21}
     15: {2,3}     58: {1,10}    111: {2,12}    155: {3,11}
     21: {2,4}     62: {1,11}    115: {3,9}     158: {1,22}
     22: {1,5}     65: {3,6}     118: {1,17}    159: {2,16}
     25: {3,3}     69: {2,9}     119: {4,7}     161: {4,9}
     26: {1,6}     74: {1,12}    121: {5,5}     166: {1,23}
     33: {2,5}     77: {4,5}     122: {1,18}    169: {6,6}
     34: {1,7}     82: {1,13}    123: {2,13}    177: {2,17}
     35: {3,4}     85: {3,7}     129: {2,14}    178: {1,24}
     38: {1,8}     86: {1,14}    133: {4,8}     183: {2,18}
     39: {2,6}     87: {2,10}    134: {1,19}    185: {3,12}
		

Crossrefs

A112798 restricted to rows of length 2 gives this triangle.
A115392 is the row number for the first appearance of each positive integer.
A176506 gives row differences.
A338899 is the squarefree version.
A338912 is column 1.
A338913 is column 2.
A001221 counts a number's distinct prime indices.
A001222 counts a number's prime indices.
A001358 lists semiprimes.
A004526 counts 2-part partitions.
A006881 lists squarefree semiprimes.
A037143 lists primes and semiprimes.
A046315 and A100484 list odd and even semiprimes.
A046388 and A100484 list odd and even squarefree semiprimes.
A065516 gives first differences of semiprimes.
A084126 and A084127 give the prime factors of semiprimes.
A270650 and A270652 give the prime indices of squarefree semiprimes.
A320655 counts factorizations into semiprimes.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Join@@primeMS/@Select[Range[100],PrimeOmega[#]==2&]

A338912 Lesser prime index of the n-th semiprime.

Original entry on oeis.org

1, 1, 2, 1, 1, 2, 2, 1, 3, 1, 2, 1, 3, 1, 2, 1, 4, 2, 3, 2, 1, 1, 3, 2, 1, 4, 1, 3, 1, 2, 4, 2, 1, 3, 1, 2, 3, 1, 4, 5, 1, 2, 2, 4, 1, 2, 1, 5, 3, 1, 3, 1, 2, 4, 1, 6, 2, 1, 2, 3, 5, 1, 2, 1, 4, 3, 1, 5, 2, 1, 3, 4, 1, 2, 6, 1, 3, 2, 6, 2, 5, 1, 4, 1, 3, 2, 1
Offset: 1

Views

Author

Gus Wiseman, Nov 20 2020

Keywords

Comments

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

Examples

			The semiprimes are:
  2*2, 2*3, 3*3, 2*5, 2*7, 3*5, 3*7, 2*11, 5*5, 2*13, ...
so the lesser prime factors are:
  2, 2, 3, 2, 2, 3, 3, 2, 5, 2, ...
with indices:
  1, 1, 2, 1, 1, 2, 2, 1, 3, 1, ...
		

Crossrefs

A084126 is the lesser prime factor (not index).
A084127 is the greater factor, with index A338913.
A115392 lists positions of ones.
A128301 lists positions of first appearances of each positive integer.
A270650 is the squarefree case, with greater part A270652.
A338898 has this as first column.
A001221 counts distinct prime indices.
A001222 counts prime indices.
A001358 lists semiprimes, with odds A046315 and evens A100484.
A006881 lists squarefree semiprimes, with odds A046388 and evens A100484.
A087794/A176504/A176506 are product/sum/difference of semiprime indices.
A338910/A338911 list products of pairs of odd/even-indexed primes.

Programs

  • Mathematica
    Table[Min[PrimePi/@First/@FactorInteger[n]],{n,Select[Range[100],PrimeOmega[#]==2&]}]

Formula

a(n) = A000720(A084126(n)).

A046337 Odd numbers with an even number of prime factors (counted with multiplicity).

Original entry on oeis.org

1, 9, 15, 21, 25, 33, 35, 39, 49, 51, 55, 57, 65, 69, 77, 81, 85, 87, 91, 93, 95, 111, 115, 119, 121, 123, 129, 133, 135, 141, 143, 145, 155, 159, 161, 169, 177, 183, 185, 187, 189, 201, 203, 205, 209, 213, 215, 217, 219, 221, 225, 235, 237, 247, 249, 253, 259
Offset: 1

Views

Author

Patrick De Geest, Jun 15 1998

Keywords

Crossrefs

Intersection of A005408 and A028260.
Setwise difference A005408 \ A067019.
Setwise difference A028260 \ A063745.
Union of A359161 and A359163.
Union of A327862 and A360110.
Subsequence of A345452, of A356312 and of A359371.
Positions of positive terms in A166698, positions of even terms in A327858 and A356299.
Subsequences: A002557, A046315 (odd semiprimes), A056913, A359596, A359607, A359608 (without its term 2).
Cf. A000035, A008836, A046338, A046470, A353557 (characteristic function), A358777.
Cf. also A036349, A297845.

Programs

  • Mathematica
    Select[Range[1,301,2],EvenQ[PrimeOmega[#]]&] (* Harvey P. Dale, Jul 25 2011 *)
  • PARI
    lista(nn) = {forstep(n=1, nn, 2, if (bigomega(n) % 2 == 0, print1(n, ", ")));} \\ Michel Marcus, Jul 04 2015

Formula

{k | A000035(k) > 0 and A008836(k) > 0}. - Antti Karttunen, Jan 13 2023

A176504 a(n) = m + k where prime(m)*prime(k) = semiprime(n).

Original entry on oeis.org

2, 3, 4, 4, 5, 5, 6, 6, 6, 7, 7, 8, 7, 9, 8, 10, 8, 9, 8, 10, 11, 12, 9, 11, 13, 9, 14, 10, 15, 12, 10, 13, 16, 11, 17, 14, 12, 18, 11, 10, 19, 15, 16, 12, 20, 17, 21, 11, 13, 22, 14, 23, 18, 13, 24, 12, 19, 25, 20, 15, 12, 26, 21, 27, 14, 16, 28, 13, 22, 29, 17, 15, 30, 23, 13, 31
Offset: 1

Views

Author

Juri-Stepan Gerasimov, Apr 19 2010

Keywords

Examples

			From _Gus Wiseman_, Dec 04 2020: (Start)
A semiprime (A001358) is a product of any two prime numbers. The sequence of all semiprimes together with their prime indices and weights begins:
   4: 1 + 1 = 2
   6: 1 + 2 = 3
   9: 2 + 2 = 4
  10: 1 + 3 = 4
  14: 1 + 4 = 5
  15: 2 + 3 = 5
  21: 2 + 4 = 6
  22: 1 + 5 = 6
  25: 3 + 3 = 6
  26: 1 + 6 = 7
(End)
		

Crossrefs

A056239 is the version for not just semiprimes.
A087794 gives the product of the same two indices.
A176506 gives the difference of the same two indices.
A338904 puts the n-th semiprime in row a(n).
A001358 lists semiprimes.
A006881 lists squarefree semiprimes.
A338898/A338912/A338913 give the prime indices of semiprimes.
A338899/A270650/A270652 give the prime indices of squarefree semiprimes, with product/sum/difference A339361/A339362/A338900.

Programs

  • Maple
    From R. J. Mathar, Apr 20 2010: (Start)
    isA001358 := proc(n) numtheory[bigomega](n) = 2 ; end proc:
    A001358 := proc(n) option remember ; if n = 1 then return 4 ; else for a from procname(n-1)+1 do if isA001358(a) then return a; end if; end do; end if; end proc:
    A084126 := proc(n) min(op(numtheory[factorset](A001358(n)))) ; end proc:
    A084127 := proc(n) max(op(numtheory[factorset](A001358(n)))) ; end proc:
    A176504 := proc(n) numtheory[pi](A084126(n)) + numtheory[pi](A084127(n)) ; end proc: seq(A176504(n),n=1..80) ; (End)
  • Mathematica
    Table[If[SquareFreeQ[n],Total[PrimePi/@First/@FactorInteger[n]],2*PrimePi[Sqrt[n]]],{n,Select[Range[100],PrimeOmega[#]==2&]}] (* Gus Wiseman, Dec 04 2020 *)

Formula

a(n) = A056239(A001358(n)) = A338912(n) + A338913(n). - Gus Wiseman, Dec 04 2020
sqrt(n/(log n log log n)) << a(n) << n/log log n. - Charles R Greathouse IV, Apr 17 2024

Extensions

Entries checked by R. J. Mathar, Apr 20 2010

A087794 Products of prime-indices of factors of semiprimes.

Original entry on oeis.org

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

Views

Author

Reinhard Zumkeller, Oct 09 2003

Keywords

Comments

A semiprime (A001358) 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. - Gus Wiseman, Dec 04 2020

Examples

			A001358(20)=57=3*19=A000040(2)*A000040(8), therefore a(20)=2*8=16.
From _Gus Wiseman_, Dec 04 2020: (Start)
The sequence of all semiprimes together with the products of their prime indices begins:
   4: 1 * 1 = 1
   6: 1 * 2 = 2
   9: 2 * 2 = 4
  10: 1 * 3 = 3
  14: 1 * 4 = 4
  15: 2 * 3 = 6
  21: 2 * 4 = 8
  22: 1 * 5 = 5
  25: 3 * 3 = 9
  26: 1 * 6 = 6
(End)
		

Crossrefs

A003963 is the version for not just semiprimes.
A176504 gives the sum of the same two indices.
A176506 gives the difference of the same two indices.
A339361 is the squarefree case.
A001358 lists semiprimes.
A006881 lists squarefree semiprimes.
A289182/A115392 list the positions of odd/even terms of A001358.
A338898/A338912/A338913 give the prime indices of semiprimes.
A338899/A270650/A270652 give the prime indices of squarefree semiprimes.
A338904 groups semiprimes by weight.

Programs

  • Mathematica
    Table[If[SquareFreeQ[n],Times@@PrimePi/@First/@FactorInteger[n],PrimePi[Sqrt[n]]^2],{n,Select[Range[100],PrimeOmega[#]==2&]}] (* Gus Wiseman, Dec 04 2020 *)

Formula

a(n) = A003963(A001358(n)) = A338912(n) * A338913(n). - Gus Wiseman, Dec 04 2020

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}]

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.

A108352 a(n) = primal code characteristic of n, which is the least positive integer, if any, such that (n o)^k = 1, otherwise equal to 0. Here "o" denotes the primal composition operator, as illustrated in A106177 and A108371 and (n o)^k = n o ... o n, with k occurrences of n.

Original entry on oeis.org

1, 0, 2, 2, 2, 0, 2, 2, 0, 0, 2, 0, 2, 0, 2, 2, 2, 0, 2, 3, 2, 0, 2, 3, 2, 0, 2, 3, 2, 0, 2, 2, 2, 0, 2, 0, 2, 0, 2, 0, 2, 0, 2, 3, 0, 0, 2, 3, 2, 0, 2, 3, 2, 0, 2, 3, 2, 0, 2, 0, 2, 0, 0, 2, 2, 0, 2, 3, 2, 0, 2, 0, 2, 0, 3, 3, 2, 0, 2, 3, 2, 0, 2, 0, 2, 0, 2, 3, 2, 0, 2, 3, 2, 0, 2, 3, 2, 0, 0, 2, 2, 0, 2, 3, 2
Offset: 1

Views

Author

Jon Awbrey, May 31 2005, revised Jun 01 2005

Keywords

Examples

			a(1) = 1 because (1 o)^1 = ({ } o)^1 = 1.
a(2) = 0 because (2 o)^k = (1:1 o)^k = 2, for all positive k.
a(3) = 2 because (3 o)^2 = (2:1 o)^2 = 1.
a(4) = 2 because (4 o)^2 = (1:2 o)^2 = 1.
a(5) = 2 because (5 o)^2 = (3:1 o)^2 = 1.
a(6) = 0 because (6 o)^k = (1:1 2:1 o)^k = 6, for all positive k.
a(7) = 2 because (7 o)^2 = (4:1 o)^1 = 1.
a(8) = 2 because (8 o)^2 = (1:3 o)^1 = 1.
a(9) = 0 because (9 o)^k = (2:2 o)^k = 9, for all positive k.
a(10) = 0 because (10 o)^k = (1:1 3:1 o)^k = 10, for all positive k.
Detail of calculation for compositional powers of 12:
(12 o)^2 = (1:2 2:1) o (1:2 2:1) = (1:1 2:2) = 18
(12 o)^3 = (1:1 2:2) o (1:2 2:1) = (1:2 2:1) = 12
Detail of calculation for compositional powers of 20:
(20 o)^2 = (1:2 3:1) o (1:2 3:1) = (3:2) = 25
(20 o)^3 = (3:2) o (1:2 3:1) = 1.
From _Antti Karttunen_, Nov 20 2019: (Start)
For n=718, because 718 = prime(1)^1 * prime(72)^1, its partial function primal code is (1:1 72:1), which, when composed with itself stays same (that is, A106177(718,718) = 718), thus, as 1 is never reached, a(718) = 0, like is true for all even nonsquare semiprimes.
For n=1804, as 1804 = prime(1)^2 * prime(5)^1 * prime(13)^1, its primal code is (1:2 5:1 13:1), which, when composed with itself yields 203401 = prime(5)^2 * prime(13)^2, i.e., primal code (5:2 13:2), which when composed with (1:2 5:1 13:1) yields 1, which happened on the second iteration, thus a(1804) = 2+1 = 3.
(End)
		

Crossrefs

Programs

  • PARI
    A106177sq(n,k) = { my(f = factor(k)); prod(i=1,#f~,f[i, 1]^valuation(n, prime(f[i, 2]))); }; \\ As in A106177.
    A108352(n) = { my(orgn=n,xs=Set([]), k=1); while(n>1, if(vecsearch(xs,n), return(0)); xs = setunion([n],xs); n = A106177sq(n,orgn); k++); (k); }; \\ Antti Karttunen, Nov 20 2019

Formula

a(A065091(n)) = 2 for all n, a(A001747(n)) = 0 for all n, except n=2, and a(A046315(n)) = 2 for n > 1. - Antti Karttunen, Nov 20 2019

Extensions

Links and cross-references added, Aug 19 2005
Term a(63) corrected and five more terms added (up to a(105)) by Antti Karttunen, Nov 20 2019
Previous Showing 11-20 of 121 results. Next