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

A355382 Number of divisors d of n such that bigomega(d) = omega(n).

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jul 02 2022

Keywords

Comments

The statistic omega = A001221 counts distinct prime factors (without multiplicity).
The statistic bigomega = A001222 counts prime factors with multiplicity.
If positive integers are regarded as arrows from the number of prime factors to the number of distinct prime factors, this sequence counts divisible composable pairs. Is there a nice choice of a composition operation making this into an associative category?

Examples

			The set of divisors of 180 satisfying the condition is {12, 18, 20, 30, 45}, so a(180) = 5.
		

Crossrefs

The version with multiplicity is A181591.
For partitions we have A355383, with multiplicity A339006.
The version for compositions is A355384.
Positions of first appearances are A355386.
A000005 counts divisors.
A001221 counts prime indices without multiplicity.
A001222 count prime indices with multiplicity.
A070175 gives representatives for bigomega and omega, triangle A303555.

Programs

  • Mathematica
    Table[Length[Select[Divisors[n],PrimeOmega[#]==PrimeNu[n]&]],{n,100}]

A355386 Position of first appearance of n in A355382, where A355382(m) = number of divisors d of m such that bigomega(d) = omega(m); or a(n) = -1 if n does not appear in A355382.

Original entry on oeis.org

1, 12, 36, 120, 180, 360, 840, 1260, 5400, 27000, 2520, 5040, 6300, 7560, 15120, 12600, 25200
Offset: 1

Views

Author

Gus Wiseman, Jul 02 2022

Keywords

Comments

The first position of -1 appears to be 18, pointed out by Amiram Eldar.
The terms are not always increasing.
The statistic omega = A001221 counts distinct prime factors (without multiplicity).
The statistic bigomega = A001222 counts prime factors with multiplicity.

Examples

			The terms together with their prime indices begin:
      1: {}
     12: {1,1,2}
     36: {1,1,2,2}
    120: {1,1,1,2,3}
    180: {1,1,2,2,3}
    360: {1,1,1,2,2,3}
    840: {1,1,1,2,3,4}
   1260: {1,1,2,2,3,4}
   5400: {1,1,1,2,2,2,3,3}
  27000: {1,1,1,2,2,2,3,3,3}
   2520: {1,1,1,2,2,3,4}
   5040: {1,1,1,1,2,2,3,4}
   6300: {1,1,2,2,3,3,4}
   7560: {1,1,1,2,2,2,3,4}
  15120: {1,1,1,1,2,2,2,3,4}
The terms together with their divisors satisfying the condition begin:
      1:   1
     12:   4,   6
     36:   4,   6,   9
    120:   8,  12,  20,  30
    180:  12,  18,  20,  30,  45
    360:   8,  12,  18,  20,  30,  45
    840:  24,  40,  56,  60,  84, 140, 210
   1260:  36,  60,  84,  90, 126, 140, 210, 315
   5400:   8,  12,  18,  20,  27,  30,  45,  50,  75
  27000:   8,  12,  18,  20,  27,  30,  45,  50,  75, 125
   2520:  24,  36,  40,  56,  60,  84,  90, 126, 140, 210, 315
   5040:  16,  24,  36,  40,  56,  60,  84,  90, 126, 140, 210, 315
   6300:  36,  60,  84,  90, 100, 126, 140, 150, 210, 225, 315, 350, 525
		

Crossrefs

These are the positions of first appearances in A355382, which is the version of A181591 without multiplicity.
A000005 counts divisors.
A001221 counts prime indices without multiplicity.
A001222 counts prime indices with multiplicity.
A070175 gives representatives for bigomega and omega, triangle A303555.
A355383 counts cmpsbl. pairs of partitions with containment, comps. A355384.

Programs

  • Mathematica
    tf=Table[Length[Select[Divisors[n],PrimeOmega[#]==PrimeNu[n]&]],{n,1000}];
    Table[Position[tf,n][[1,1]],{n,Select[Union[tf],SubsetQ[tf,Range[#]]&]}]

A355387 Number of ways to choose a distinct subsequence of an integer composition of n.

Original entry on oeis.org

1, 2, 5, 14, 37, 98, 259, 682, 1791, 4697, 12303, 32196, 84199, 220087, 575067, 1502176, 3923117, 10244069, 26746171, 69825070, 182276806, 475804961, 1241965456, 3241732629, 8461261457, 22084402087, 57640875725, 150442742575, 392652788250, 1024810764496
Offset: 0

Views

Author

Gus Wiseman, Jul 04 2022

Keywords

Comments

By "distinct" we mean equal subsequences are counted only once. For example, the pair (1,1)(1) is counted only once even though (1) is a subsequence of (1,1) in two ways. The version with multiplicity is A025192.

Examples

			The a(3) = 14 pairings of a composition with a chosen subsequence:
  (3)()     (3)(3)
  (21)()    (21)(1)   (21)(2)    (21)(21)
  (12)()    (12)(1)   (12)(2)    (12)(12)
  (111)()   (111)(1)  (111)(11)  (111)(111)
		

Crossrefs

For partitions we have A000712, composable A339006.
The homogeneous version is A011782, without containment A000302.
With multiplicity we have A025192, for partitions A070933.
The strict case is A032005.
The case of strict subsequences is A236002.
The composable case is A355384, homogeneous without containment A355388.
A075900 counts compositions of each part of a partition.
A304961 counts compositions of each part of a strict partition.
A307068 counts strict compositions of each part of a composition.
A336127 counts compositions of each part of a strict composition.

Programs

  • Mathematica
    Table[Sum[Length[Union[Subsets[y]]],{y,Join@@Permutations/@IntegerPartitions[n]}],{n,0,6}]
  • PARI
    lista(n)=my(f=sum(k=1,n,(x^k+x*O(x^n))/(1-x/(1-x)+x^k)));Vec((1-x)/((1-2*x)*(1-f))) \\ Christian Sievers, May 06 2025

Formula

G.f.: (1-x)/((1-2*x)*(1-f)) where f = Sum_{k>=1} x^k/(1-x/(1-x)+x^k) is the generating function for A331330. - Christian Sievers, May 06 2025

Extensions

a(16) and beyond from Christian Sievers, May 06 2025

A289621 Compound filter (omega & bigomega): a(1) = 0, for n > 1, a(n) = P(A001221(n), A001222(n)), where P(n,k) is sequence A000027 used as a pairing function.

Original entry on oeis.org

0, 1, 1, 2, 1, 5, 1, 4, 2, 5, 1, 8, 1, 5, 5, 7, 1, 8, 1, 8, 5, 5, 1, 12, 2, 5, 4, 8, 1, 13, 1, 11, 5, 5, 5, 12, 1, 5, 5, 12, 1, 13, 1, 8, 8, 5, 1, 17, 2, 8, 5, 8, 1, 12, 5, 12, 5, 5, 1, 18, 1, 5, 8, 16, 5, 13, 1, 8, 5, 13, 1, 17, 1, 5, 8, 8, 5, 13, 1, 17, 7, 5, 1, 18, 5, 5, 5, 12, 1, 18, 5, 8, 5, 5, 5, 23, 1, 8, 8
Offset: 1

Views

Author

Antti Karttunen, Jul 16 2017

Keywords

Crossrefs

Cf. A001221, A001222, A008966, A046660, A070012, A070013, A070014, A088529, A088530, A181591 (sequences with matching equivalence classes).

Programs

  • PARI
    A289621(n) = if(1==n,0,(1/2)*(2 + ((omega(n)+bigomega(n))^2) - omega(n) - 3*bigomega(n)));
    
  • Scheme
    (define (A289621 n) (if (= 1 n) 0 (* (/ 1 2) (+ (expt (+ (A001221 n) (A001222 n)) 2) (- (A001221 n)) (- (* 3 (A001222 n))) 2))))

Formula

a(1) = 0, for n > 1, a(n) = (1/2)*(2 + ((A001221(n)+A001222(n))^2) - A001221(n) - 3*A001222(n)).
Previous Showing 11-14 of 14 results.