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

A327405 Quotient of n over the maximum divisor of n that is 1 or whose prime indices have a common divisor > 1.

Original entry on oeis.org

1, 2, 1, 4, 1, 2, 1, 8, 1, 2, 1, 4, 1, 2, 3, 16, 1, 2, 1, 4, 1, 2, 1, 8, 1, 2, 1, 4, 1, 6, 1, 32, 3, 2, 5, 4, 1, 2, 1, 8, 1, 2, 1, 4, 5, 2, 1, 16, 1, 2, 3, 4, 1, 2, 5, 8, 1, 2, 1, 12, 1, 2, 1, 64, 1, 6, 1, 4, 3, 10, 1, 8, 1, 2, 3, 4, 7, 2, 1, 16, 1, 2, 1, 4, 5
Offset: 1

Views

Author

Gus Wiseman, Sep 21 2019

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. Numbers whose prime indices have a common divisor > 1 are listed in A318978.

Examples

			The divisors of 90 that are 1 or whose prime indices have a common divisor > 1 are {1, 3, 5, 9}, so a(90) = 90/9 = 10.
		

Crossrefs

See link for additional cross-references.

Programs

  • Mathematica
    Table[n/Max[Select[Divisors[n],GCD@@PrimePi/@First/@FactorInteger[#]!=1&]],{n,100}]
  • PARI
    A327405(n) = (n / vecmax(select(d -> (1==d)||(gcd(apply(primepi,factor(d)[, 1]~))>1), divisors(n)))); \\ Antti Karttunen, Dec 06 2021

Formula

a(n) = n/A327656(n).

A327657 Number of divisors of n that are 1 or whose prime indices have a common divisor > 1.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Sep 21 2019

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. Numbers whose prime indices have a common divisor > 1 are listed in A318978.

Examples

			The divisors of 90 that are 1 or whose prime indices have a common divisor > 1 are {1, 3, 5, 9}, so a(90) = 4.
		

Crossrefs

See link for additional cross-references.

Programs

  • Mathematica
    Table[Length[Select[Divisors[n],GCD@@PrimePi/@First/@FactorInteger[#]!=1&]],{n,100}]
  • PARI
    A327657(n) = sumdiv(n, d, (1==d)||(gcd(apply(x->primepi(x), factor(d)[, 1]))>1)); \\ Antti Karttunen, Dec 05 2021

Formula

a(n) = A000005(n) - A318979(n). - Antti Karttunen, Dec 05 2021

Extensions

Data section extended up to 105 terms by Antti Karttunen, Dec 05 2021

A327658 Number of factorizations of n that are empty or whose factors have a common divisor > 1.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 3, 2, 1, 1, 2, 1, 1, 1, 5, 1, 2, 1, 2, 1, 1, 1, 4, 2, 1, 3, 2, 1, 1, 1, 7, 1, 1, 1, 4, 1, 1, 1, 4, 1, 1, 1, 2, 2, 1, 1, 7, 2, 2, 1, 2, 1, 4, 1, 4, 1, 1, 1, 3, 1, 1, 2, 11, 1, 1, 1, 2, 1, 1, 1, 7, 1, 1, 2, 2, 1, 1, 1, 7, 5, 1, 1, 3, 1, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Sep 21 2019

Keywords

Comments

First differs from A319786 at a(900) = 11, A319786(900) = 12.
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. Numbers whose prime indices have a common divisor > 1 are listed in A318978.

Examples

			The a(120) = 7 factorizations:
  (120)
  (2*60)
  (4*30)
  (6*20)
  (10*12)
  (2*2*30)
  (2*6*10)
		

Crossrefs

See link for additional cross-references.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[facs[n],#=={}||GCD@@#!=1&]],{n,100}]

A327685 Nonprime numbers whose prime indices have a common divisor > 1.

Original entry on oeis.org

9, 21, 25, 27, 39, 49, 57, 63, 65, 81, 87, 91, 111, 115, 117, 121, 125, 129, 133, 147, 159, 169, 171, 183, 185, 189, 203, 213, 235, 237, 243, 247, 259, 261, 267, 273, 289, 299, 301, 303, 305, 319, 321, 325, 333, 339, 343, 351, 361, 365, 371, 377, 387, 393, 399
Offset: 1

Views

Author

Gus Wiseman, Sep 22 2019

Keywords

Comments

First differs from A322336 in lacking 2535 = prime(2)*prime(3)*prime(6)^2.
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. Heinz numbers of integer partitions with a common divisor > 1 are A318978, and the enumeration of these partitions by sum is A108572.

Examples

			The sequence of terms together with their prime indices begins:
    9: {2,2}
   21: {2,4}
   25: {3,3}
   27: {2,2,2}
   39: {2,6}
   49: {4,4}
   57: {2,8}
   63: {2,2,4}
   65: {3,6}
   81: {2,2,2,2}
   87: {2,10}
   91: {4,6}
  111: {2,12}
  115: {3,9}
  117: {2,2,6}
  121: {5,5}
  125: {3,3,3}
  129: {2,14}
  133: {4,8}
  147: {2,4,4}
		

Crossrefs

Programs

  • Mathematica
    Select[Range[100],#>1&&!PrimeQ[#]&&GCD@@PrimePi/@First/@FactorInteger[#]>1&]

Formula

Complement of A000040 in A318978.

A338552 Non-powers of primes whose prime indices have a common divisor > 1.

Original entry on oeis.org

21, 39, 57, 63, 65, 87, 91, 111, 115, 117, 129, 133, 147, 159, 171, 183, 185, 189, 203, 213, 235, 237, 247, 259, 261, 267, 273, 299, 301, 303, 305, 319, 321, 325, 333, 339, 351, 365, 371, 377, 387, 393, 399, 417, 427, 441, 445, 453, 477, 481, 489, 497, 507
Offset: 1

Views

Author

Gus Wiseman, Nov 03 2020

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 Heinz numbers of non-constant, non-relatively prime partitions. The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), giving a bijective correspondence between positive integers and integer partitions.

Examples

			The sequence of terms together with their prime indices begins:
     21: {2,4}      183: {2,18}       305: {3,18}
     39: {2,6}      185: {3,12}       319: {5,10}
     57: {2,8}      189: {2,2,2,4}    321: {2,28}
     63: {2,2,4}    203: {4,10}       325: {3,3,6}
     65: {3,6}      213: {2,20}       333: {2,2,12}
     87: {2,10}     235: {3,15}       339: {2,30}
     91: {4,6}      237: {2,22}       351: {2,2,2,6}
    111: {2,12}     247: {6,8}        365: {3,21}
    115: {3,9}      259: {4,12}       371: {4,16}
    117: {2,2,6}    261: {2,2,10}     377: {6,10}
    129: {2,14}     267: {2,24}       387: {2,2,14}
    133: {4,8}      273: {2,4,6}      393: {2,32}
    147: {2,4,4}    299: {6,9}        399: {2,4,8}
    159: {2,16}     301: {4,14}       417: {2,34}
    171: {2,2,8}    303: {2,26}       427: {4,18}
		

Crossrefs

A318978 allows prime powers, counted by A018783, with complement A289509.
A327685 allows nonprime prime powers.
A338330 is the coprime instead of relatively prime version.
A338554 counts the partitions with these Heinz numbers.
A338555 is the complement.
A000740 counts relatively prime compositions.
A000961 lists powers of primes, with complement A024619.
A051424 counts pairwise coprime or singleton partitions.
A108572 counts nontrivial periodic partitions, with Heinz numbers A001597.
A291166 ranks relatively prime compositions, with complement A291165.
A302696 gives the Heinz numbers of pairwise coprime partitions.
A327516 counts pairwise coprime partitions, with Heinz numbers A302696.

Programs

  • Mathematica
    Select[Range[100],!(#==1||PrimePowerQ[#]||GCD@@PrimePi/@First/@FactorInteger[#]==1)&]

Formula

Equals A024619 /\ A318978.
Complement of A000961 \/ A289509.

A338555 Numbers that are either a power of a prime or have relatively prime prime indices.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Nov 03 2020

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 Heinz numbers of partitions either constant or relatively prime (A338553). The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), giving a bijective correspondence between positive integers and integer partitions.

Crossrefs

A327534 uses primes instead of prime powers.
A338331 is the pairwise coprime version, with complement A338330.
A338552 is the complement.
A338553 counts the partitions with these Heinz numbers.
A000837 counts relatively prime partitions, with Heinz numbers A289509.
A000961 lists powers of primes.
A018783 counts partitions whose prime indices are not relatively prime, with Heinz numbers A318978.
A051424 counts pairwise coprime or singleton partitions.
A291166 ranks relatively prime compositions, with complement A291165.
A327516 counts pairwise coprime partitions, with Heinz numbers A302696.

Programs

  • Mathematica
    Select[Range[100],#==1||PrimePowerQ[#]||GCD@@PrimePi/@First/@FactorInteger[#]==1&]

Formula

Equals A000961 \/ A289509.
Complement of A024619 /\ A318978.

A320439 Number of factorizations of n into factors > 1 where each factor's prime indices are relatively prime. Number of factorizations of n using elements of A289509.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jan 08 2019

Keywords

Comments

Also the number of multiset partitions of the multiset of prime indices of n using multisets each of which is relatively prime.
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.
Two or more numbers are relatively prime if they have no common divisor > 1. A single number is not considered to be relatively prime unless it is equal to 1.

Examples

			The a(72) = 6 factorizations are (2*2*18), (2*6*6), (2*36), (4*18), (6*12), (72).
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    facsrp[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#,d]&)/@Select[facsrp[n/d],Min@@#>=d&],{d,Select[Rest[Divisors[n]],GCD@@primeMS[#]==1&]}]];
    Table[Length[facsrp[n]],{n,100}]
  • PARI
    A320439(n, m=n) = if(1==n, 1, my(s=0); fordiv(n, d, if((d<=m)&&(1==gcd(apply(x->primepi(x), factor(d)[, 1]))), s += A320439(n/d, d))); (s)); \\ Antti Karttunen, Dec 06 2021

A327515 Number of steps to reach a fixed point starting with n and repeatedly taking the quotient by the maximum divisor that is 1, 2, or a nonprime number whose prime indices are pairwise coprime (A327512, A327514).

Original entry on oeis.org

0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 1, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 1, 1, 1, 1, 0, 1, 0, 1, 0, 1, 0, 1, 1, 1, 0, 1, 0, 1, 1, 1, 0, 1, 1, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 1, 1, 0, 1, 0, 1, 1, 1, 1, 1, 0, 1, 0, 1, 0, 1, 1, 1, 0, 1, 0, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 0, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Sep 19 2019

Keywords

Comments

Positions of zeros are A289509.
First term > 1 is a(225) = 2.
First zero not in A318978 is a(17719) = 0.
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. Numbers that are 1, 2, or a nonprime number whose prime indices are pairwise coprime are listed in A302696.

Examples

			We have 50625 -> 3375 -> 225 ->  15 -> 1, so a(50625) = 4.
		

Crossrefs

See link for additional cross-references.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Length[FixedPointList[#/Max[Select[Divisors[#],#==1||CoprimeQ@@primeMS[#]&]]&,n]]-2,{n,100}]
  • PARI
    isA302696(n) = if(isprimepower(n), !(n%2), if(!issquarefree(n>>valuation(n, 2)), 0, my(pis=apply(primepi, factor(n)[, 1])); (lcm(pis)==factorback(pis))));
    A327512(n) = vecmax(select(isA302696, divisors(n)));
    A327515(n) = for(k=0,oo,my(nextn=n/A327512(n)); if(nextn==n,return(k)); n = nextn); \\ Antti Karttunen, Jan 28 2025

Formula

a(15^n) = n.

Extensions

Data section extended to a(105) and secondary offset added by Antti Karttunen, Jan 28 2025

A335241 Numbers whose prime indices are not pairwise coprime, where a singleton is not coprime unless it is {1}.

Original entry on oeis.org

1, 3, 5, 7, 9, 11, 13, 17, 18, 19, 21, 23, 25, 27, 29, 31, 36, 37, 39, 41, 42, 43, 45, 47, 49, 50, 53, 54, 57, 59, 61, 63, 65, 67, 71, 72, 73, 75, 78, 79, 81, 83, 84, 87, 89, 90, 91, 97, 98, 99, 100, 101, 103, 105, 107, 108, 109, 111, 113, 114, 115, 117, 121
Offset: 1

Views

Author

Gus Wiseman, May 30 2020

Keywords

Comments

We use the Mathematica definition for CoprimeQ, so a singleton is not considered coprime unless it is (1).
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 sequence of terms together with their prime indices begins:
    1: {}          31: {11}          61: {18}
    3: {2}         36: {1,1,2,2}     63: {2,2,4}
    5: {3}         37: {12}          65: {3,6}
    7: {4}         39: {2,6}         67: {19}
    9: {2,2}       41: {13}          71: {20}
   11: {5}         42: {1,2,4}       72: {1,1,1,2,2}
   13: {6}         43: {14}          73: {21}
   17: {7}         45: {2,2,3}       75: {2,3,3}
   18: {1,2,2}     47: {15}          78: {1,2,6}
   19: {8}         49: {4,4}         79: {22}
   21: {2,4}       50: {1,3,3}       81: {2,2,2,2}
   23: {9}         53: {16}          83: {23}
   25: {3,3}       54: {1,2,2,2}     84: {1,1,2,4}
   27: {2,2,2}     57: {2,8}         87: {2,10}
   29: {10}        59: {17}          89: {24}
		

Crossrefs

The complement is A302696.
The version for relatively prime instead of coprime is A318978.
The version for standard compositions is A335239.
These are the Heinz numbers of the partitions counted by A335240.
Singleton or pairwise coprime partitions are counted by A051424.
Singleton or pairwise coprime sets are ranked by A087087.
Primes and numbers with pairwise coprime prime indices are A302569.
Numbers whose binary indices are pairwise coprime are A326675.
Coprime standard composition numbers are A333227.

Programs

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

A338554 Number of non-constant integer partitions of n whose parts have a common divisor > 1.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 1, 0, 2, 1, 5, 0, 9, 0, 13, 6, 18, 0, 33, 0, 40, 14, 54, 0, 87, 5, 99, 27, 133, 0, 211, 0, 226, 55, 295, 18, 443, 0, 488, 100, 637, 0, 912, 0, 1000, 198, 1253, 0, 1775, 13, 1988, 296, 2434, 0, 3358, 59, 3728, 489, 4563, 0, 6241, 0, 6840, 814
Offset: 0

Views

Author

Gus Wiseman, Nov 07 2020

Keywords

Examples

			The a(6) = 2 through a(15) = 6 partitions (empty columns indicated by dots, A = 10, B = 11, C = 12):
  (42)  .  (62)   (63)  (64)    .  (84)     .  (86)      (96)
           (422)        (82)       (93)        (A4)      (A5)
                        (442)      (A2)        (C2)      (C3)
                        (622)      (633)       (644)     (663)
                        (4222)     (642)       (662)     (933)
                                   (822)       (842)     (6333)
                                   (4422)      (A22)
                                   (6222)      (4442)
                                   (42222)     (6422)
                                               (8222)
                                               (44222)
                                               (62222)
                                               (422222)
		

Crossrefs

A046022 lists positions of zeros.
A082023(n) - A059841(n) is the 2-part version, n > 2.
A303280(n) - 1 is the strict case (n > 1).
A338552 lists the Heinz numbers of these partitions.
A338553 counts the complement, with Heinz numbers A338555.
A000005 counts constant partitions, with Heinz numbers A000961.
A000837 counts relatively prime partitions, with Heinz numbers A289509.
A018783 counts non-relatively prime partitions (ordered: A178472), with Heinz numbers A318978.
A282750 counts relatively prime partitions by sum and length.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],!SameQ@@#&&GCD@@#>1&]],{n,0,30}]

Formula

For n > 0, a(n) = A018783(n) - A000005(n) + 1.
Previous Showing 11-20 of 23 results. Next