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-8 of 8 results.

A111907 Numbers k such that the same number of primes, among primes <= the largest prime dividing k, divide k as do not.

Original entry on oeis.org

1, 3, 9, 14, 21, 27, 28, 35, 56, 63, 78, 81, 98, 112, 130, 147, 156, 175, 182, 189, 195, 196, 224, 234, 243, 245, 260, 273, 286, 312, 364, 392, 429, 441, 448, 455, 468, 520, 567, 570, 572, 585, 624, 650, 686, 702, 715, 728, 729, 784, 798, 819, 875, 896, 936
Offset: 1

Views

Author

Leroy Quet, Aug 19 2005

Keywords

Comments

Also numbers whose greatest prime index (A061395) is twice their number of distinct prime factors (A001221). - Gus Wiseman, Mar 19 2023

Examples

			28 is included because 7 is the largest prime dividing 28. And of the primes <= 7 (2,3,5,7), 2 and 7 (2 primes) divide 28 and 3 and 5 (also 2 primes) do not divide 28.
From _Gus Wiseman_, Mar 19 2023: (Start)
The terms together with their prime indices begin:
    1: {}
    3: {2}
    9: {2,2}
   14: {1,4}
   21: {2,4}
   27: {2,2,2}
   28: {1,1,4}
   35: {3,4}
   56: {1,1,1,4}
   63: {2,2,4}
   78: {1,2,6}
   81: {2,2,2,2}
   98: {1,4,4}
  112: {1,1,1,1,4}
  130: {1,3,6}
  147: {2,4,4}
  156: {1,1,2,6}
For example, 156 is included because it has prime indices {1,1,2,6}, with distinct parts {1,2,6} and distinct non-parts {3,4,5}, both of length 3. Alternatively, 156 has greatest prime index 6 and omega 3, and 6 = 2*3.
(End)
		

Crossrefs

For length instead of maximum we have A067801.
These partitions are counted by A239959.
A001222 (bigomega) counts prime factors, distinct A001221 (omega).
A061395 gives greatest prime index.
A112798 lists prime indices, sum A056239.
Comparing twice the number of distinct parts to greatest part:
less: A360254, ranks A111906
equal: A239959, ranks A111907
greater: A237365, ranks A111905
less or equal: A237363, ranks A361204
greater or equal: A361394, ranks A361395

Programs

  • Mathematica
    Select[Range[100],2*PrimeNu[#]==PrimePi[FactorInteger[#][[-1,1]]]&] (* Gus Wiseman, Mar 19 2023 *)
  • PARI
    {m=950;v=vector(m);for(n=1,m,f=factor(n)[,1]~;c=0;pc=0;forprime(p=2,vecmax(f), j=1;s=length(f);while(j<=s&&p!=f[j],j++);if(j<=s,c++);pc++);v[n]=sign(pc-2*c)); for(n=1,m,if(v[n]==0,print1(n,",")))} \\ Klaus Brockhaus, Aug 21 2005
    
  • Python
    from itertools import count, islice
    from sympy import sieve, factorint
    def a_gen():
        yield 1
        for k in count(3):
            f = [sieve.search(i)[0] for i in factorint(k)]
            if 2*len(f) == f[-1]:
                yield k
    A111907_list = list(islice(a_gen(), 100)) # John Tyler Rascoe, Jun 20 2024

Extensions

More terms from Klaus Brockhaus, Aug 21 2005

A361856 Positive integers whose prime indices satisfy (maximum) = 2*(median).

Original entry on oeis.org

12, 24, 42, 48, 60, 63, 72, 96, 126, 130, 140, 144, 189, 192, 195, 252, 266, 288, 308, 325, 330, 360, 378, 384, 399, 420, 432, 495, 546, 567, 572, 576, 588, 600, 630, 638, 650, 665, 756, 768, 819, 864, 882, 884, 931, 945, 957, 962, 975, 1122, 1134, 1152, 1190
Offset: 1

Views

Author

Gus Wiseman, Apr 02 2023

Keywords

Comments

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.
The median of a multiset is either the middle part (for odd length), or the average of the two middle parts (for even length).
These are Heinz numbers of partitions satisfying (maximum) = 2*(median).

Examples

			The terms together with their prime indices begin:
    12: {1,1,2}
    24: {1,1,1,2}
    42: {1,2,4}
    48: {1,1,1,1,2}
    60: {1,1,2,3}
    63: {2,2,4}
    72: {1,1,1,2,2}
    96: {1,1,1,1,1,2}
   126: {1,2,2,4}
   130: {1,3,6}
   140: {1,1,3,4}
   144: {1,1,1,1,2,2}
The prime indices of 126 are {1,2,2,4}, with maximum 4 and median 2, so 126 is in the sequence.
The prime indices of 308 are {1,1,4,5}, with maximum 5 and median 5/2, so 308 is in the sequence.
		

Crossrefs

The LHS (greatest prime index) is A061395.
The RHS (twice median) is A360005, distinct A360457.
These partitions are counted by A361849.
For mean instead of median we have A361855, counted by A361853.
For minimum instead of median we have A361908, counted by A118096.
For length instead of median we have A361909, counted by A237753.
A000975 counts subsets with integer median.
A001222 (bigomega) counts prime factors, distinct A001221 (omega).
A112798 lists prime indices, sum A056239.
A325347 counts partitions with integer median, complement A307683.
A359893 and A359901 count partitions by median.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{}, Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],Max@@prix[#]==2*Median[prix[#]]&]

Formula

A061395(a(n)) = 2*A360005(a(n)).

A361855 Numbers > 1 whose prime indices satisfy (maximum) * (length) = 2*(sum).

Original entry on oeis.org

28, 40, 78, 84, 171, 190, 198, 220, 240, 252, 280, 351, 364, 390, 406, 435, 714, 748, 756, 765, 777, 784, 814, 840, 850, 925, 988, 1118, 1197, 1254, 1330, 1352, 1419, 1425, 1440, 1505, 1564, 1600, 1638, 1716, 1755, 1794, 1802, 1820, 1950, 2067, 2204, 2254
Offset: 1

Views

Author

Gus Wiseman, Mar 29 2023

Keywords

Comments

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.
Also positive integers whose prime indices satisfy (maximum) = 2*(mean).
Also Heinz numbers of partitions of the same size as their complement (see example).

Examples

			The terms together with their prime indices begin:
   28: {1,1,4}
   40: {1,1,1,3}
   78: {1,2,6}
   84: {1,1,2,4}
  171: {2,2,8}
  190: {1,3,8}
  198: {1,2,2,5}
  220: {1,1,3,5}
  240: {1,1,1,1,2,3}
  252: {1,1,2,2,4}
  280: {1,1,1,3,4}
The prime indices of 84 are {1,1,2,4}, with maximum 4, length 4, and sum 8, and 4*4 = 2*8, so 84 is in the sequence.
The prime indices of 120 are {1,1,1,2,3}, with maximum 3, length 5, and sum 8, and 3*5 != 2*8, so 120 is not in the sequence.
The prime indices of 252 are {1,1,2,2,4}, with maximum 4, length 5, and sum 10, and 4*5 = 2*10, so 252 is in the sequence.
The partition (5,2,2,1) with Heinz number 198 has diagram:
  o o o o o
  o o . . .
  o o . . .
  o . . . .
Since the partition and its complement (shown in dots) both have size 10, 198 is in the sequence.
		

Crossrefs

These partitions are counted by A361853, strict A361854.
For median instead of mean we have A361856, counted by A361849.
For minimum instead of mean we have A361908, counted by A118096.
For length instead of mean we have A361909, counted by A237753.
A001222 (bigomega) counts prime factors, distinct A001221 (omega).
A061395 gives greatest prime index.
A112798 lists prime indices, sum A056239.
A326567/A326568 gives mean of prime indices.

Programs

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

Formula

A061395(a(n)) * A001222(a(n)) = 2*A056239(a(n)).

A361395 Positive integers k such that 2*omega(k) >= bigomega(k).

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 9, 10, 11, 12, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 28, 29, 30, 31, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 73, 74
Offset: 1

Views

Author

Gus Wiseman, Mar 16 2023

Keywords

Comments

Differs from A068938 in having 1 and 4 and lacking 80.
Includes all squarefree numbers.

Examples

			The prime indices of 80 are {1,1,1,1,3}, with 5 parts and 2 distinct parts, and 2*2 < 5, so 80 is not in the sequence.
		

Crossrefs

Complement of A360558.
Positions of nonnegative terms in A361205.
These partitions are counted by A361394.
A001222 (bigomega) counts prime factors, distinct A001221 (omega).
A112798 lists prime indices, sum A056239.
A360005 gives median of prime indices (times 2), distinct A360457.
Comparing twice the number of distinct parts to the number of parts:
less: A360254, ranks A360558
equal: A239959, ranks A067801
greater: A237365, ranks A361393
less or equal: A237363, ranks A361204
greater or equal: A361394, ranks A361395

Programs

  • Mathematica
    Select[Range[100],2*PrimeNu[#]>=PrimeOmega[#]&]

Formula

A001222(a(n)) <= 2*A001221(a(n)).

A361909 Positive integers > 1 whose prime indices satisfy: (maximum) = 2*(length).

Original entry on oeis.org

3, 14, 21, 35, 49, 52, 78, 117, 130, 152, 182, 195, 228, 273, 286, 325, 338, 342, 380, 429, 455, 464, 507, 513, 532, 570, 637, 696, 715, 798, 836, 845, 855, 950, 988, 1001, 1044, 1160, 1183, 1184, 1197, 1254, 1292, 1330, 1425, 1444, 1482, 1566, 1573, 1624
Offset: 1

Views

Author

Gus Wiseman, Apr 05 2023

Keywords

Comments

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 terms together with their prime indices begin:
     3: {2}
    14: {1,4}
    21: {2,4}
    35: {3,4}
    49: {4,4}
    52: {1,1,6}
    78: {1,2,6}
   117: {2,2,6}
   130: {1,3,6}
   152: {1,1,1,8}
   182: {1,4,6}
   195: {2,3,6}
   228: {1,1,2,8}
   273: {2,4,6}
   286: {1,5,6}
   325: {3,3,6}
   338: {1,6,6}
   342: {1,2,2,8}
		

Crossrefs

The LHS is A061395 (greatest prime index), least A055396.
Without multiplying by 2 in the RHS, we have A106529.
For omega instead of bigomega we have A111907, counted by A239959.
Partitions of this type are counted by A237753.
The RHS is A255201 (twice bigomega).
For mean instead of length we have A361855, counted by A361853.
For median instead of length we have A361856, counted by A361849.
For minimum instead of length we have A361908, counted by A118096.
A001221 (omega) counts distinct prime factors.
A001222 (bigomega) counts prime factors with multiplicity.
A112798 lists prime indices, sum A056239.
A316413 ranks partitions with integer mean, counted by A067538.
A326567/A326568 gives mean of prime indices.

Programs

  • Mathematica
    Select[Range[2,100],PrimePi[FactorInteger[#][[-1,1]]]==2*PrimeOmega[#]&]

A351414 Number of divisors of n that are either prime or have at least 1 square divisor > 1 and at least two distinct prime factors.

Original entry on oeis.org

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

Views

Author

Wesley Ivan Hurt, Feb 10 2022

Keywords

Examples

			a(24) = 4; 24 has divisors 2,3 (primes) and 12,24 (which both have at least 1 square divisor > 1 and at least two distinct prime factors).
a(36) = 5; 36 has divisors 2,3 (primes) and 12,18,36 (which all have at least 1 square divisor > 1 and at least two distinct prime factors).
		

Crossrefs

Cf. A001221 (omega), A008683 (mu), A048105, A361205.

Programs

  • Mathematica
    a[n_] := Module[{e = FactorInteger[n][[;;, 2]], d, nu, omega}, d = Times @@ (e+1); nu = Length[e]; omega = Total[e]; d - 2^nu - omega + 2*nu]; a[1] = 0; Array[a, 100] (* Amiram Eldar, Oct 06 2023 *)
  • PARI
    a(n) = {my(f = factor(n), d = numdiv(f), nu = omega(f), om = bigomega(f)); d - 2^nu - om + 2*nu;} \\ Amiram Eldar, Oct 06 2023

Formula

a(n) = Sum_{d|n} [[omega(d) = 1] = mu(d)^2], where [ ] is the Iverson bracket.
a(n) = A048105(n) + A361205(n). - Amiram Eldar, Oct 06 2023

A363222 Numbers whose multiset of prime indices satisfies (maximum) - (minimum) = (length).

Original entry on oeis.org

10, 21, 28, 42, 55, 70, 88, 91, 98, 99, 132, 165, 187, 198, 208, 220, 231, 247, 308, 312, 325, 330, 351, 363, 391, 455, 462, 468, 484, 520, 544, 550, 551, 585, 702, 713, 715, 726, 728, 770, 780, 816, 819, 833, 845, 975, 1073, 1078, 1092, 1144, 1170, 1210, 1216
Offset: 1

Views

Author

Gus Wiseman, May 29 2023

Keywords

Comments

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 terms together with their prime indices begin:
    10: {1,3}
    21: {2,4}
    28: {1,1,4}
    42: {1,2,4}
    55: {3,5}
    70: {1,3,4}
    88: {1,1,1,5}
    91: {4,6}
    98: {1,4,4}
    99: {2,2,5}
   132: {1,1,2,5}
   165: {2,3,5}
   187: {5,7}
   198: {1,2,2,5}
		

Crossrefs

The RHS is A001222.
Partitions of this type are counted by A237832.
The LHS (maximum minus minimum) is A243055.
A001221 (omega) counts distinct prime factors.
A112798 lists prime indices, sum A056239.
A360005 gives median of prime indices, distinct A360457.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],Max@@prix[#]-Min@@prix[#]==Length[prix[#]]&]

Formula

A061395(a(n)) - A055396(a(n)) = A001222(a(n)).

A385576 Numbers whose prime indices have the same number of distinct elements as maximal anti-runs.

Original entry on oeis.org

1, 2, 3, 5, 7, 11, 12, 13, 17, 18, 19, 20, 23, 28, 29, 31, 37, 41, 43, 44, 45, 47, 50, 52, 53, 59, 61, 63, 67, 68, 71, 73, 75, 76, 79, 83, 89, 92, 97, 98, 99, 101, 103, 107, 109, 113, 116, 117, 120, 124, 127, 131, 137, 139, 147, 148, 149, 151, 153, 157, 163
Offset: 1

Views

Author

Gus Wiseman, Jul 04 2025

Keywords

Comments

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.
These are also numbers with the same number of adjacent equal prime indices as adjacent unequal prime indices.

Examples

			The prime indices of 2640 are {1,1,1,1,2,3,5}, with 4 distinct parts {1,2,3,5} and 4 maximal anti-runs ((1),(1),(1),(2,3,5)), so 2640 is in the sequence.
The terms together with their prime indices begin:
   1: {}
   2: {1}
   3: {2}
   5: {3}
   7: {4}
  11: {5}
  12: {1,1,2}
  13: {6}
  17: {7}
  18: {1,2,2}
  19: {8}
  20: {1,1,3}
  23: {9}
  28: {1,1,4}
  29: {10}
  31: {11}
  37: {12}
  41: {13}
  43: {14}
  44: {1,1,5}
  45: {2,2,3}
  47: {15}
		

Crossrefs

The LHS is the rank statistic A001221, triangle counted by A116608.
The RHS is the rank statistic A375136, triangle counted by A133121.
These partitions are counted by A385574.
A034296 counts flat or gapless partitions, ranks A066311 or A073491.
A047993 counts partitions with max part = length, ranks A106529.
A356235 counts partitions with a neighborless singleton, ranks A356237.
A384877 gives lengths of maximal anti-runs of binary indices, firsts A384878.
A384893 counts subsets by maximal anti-runs, for partitions A268193, strict A384905.
A385572 counts subsets with the same number of runs as anti-runs, ranks A385575.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],#==1||PrimeNu[#]==Length[Split[prix[#],UnsameQ]]&]

Formula

A001221(a(n)) = A375136(a(n)).
Showing 1-8 of 8 results.