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 31-40 of 77 results. Next

A366322 Heinz numbers of integer partitions containing at least one odd part. Numbers divisible by at least one prime of odd index.

Original entry on oeis.org

2, 4, 5, 6, 8, 10, 11, 12, 14, 15, 16, 17, 18, 20, 22, 23, 24, 25, 26, 28, 30, 31, 32, 33, 34, 35, 36, 38, 40, 41, 42, 44, 45, 46, 47, 48, 50, 51, 52, 54, 55, 56, 58, 59, 60, 62, 64, 65, 66, 67, 68, 69, 70, 72, 73, 74, 75, 76, 77, 78, 80, 82, 83, 84, 85, 86
Offset: 1

Views

Author

Gus Wiseman, Oct 14 2023

Keywords

Comments

The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). This gives a bijective correspondence between positive integers and integer partitions.

Examples

			The terms together with their prime indices begin:
    2: {1}
    4: {1,1}
    5: {3}
    6: {1,2}
    8: {1,1,1}
   10: {1,3}
   11: {5}
   12: {1,1,2}
   14: {1,4}
   15: {2,3}
   16: {1,1,1,1}
   17: {7}
   18: {1,2,2}
   20: {1,1,3}
   22: {1,5}
   23: {9}
   24: {1,1,1,2}
		

Crossrefs

The complement is A066207, counted by A035363.
For all odd parts we have A066208, counted by A000009.
Partitions of this type are counted by A086543.
For even instead of odd we have A324929, counted by A047967.
A031368 lists primes of odd index.
A112798 list prime indices, sum A056239.
A257991 counts odd prime indices, distinct A324966.

Programs

  • Mathematica
    Select[Range[100],Or@@OddQ/@PrimePi/@First/@FactorInteger[#]&]

Formula

A257991(a(n)) > 0.

A372591 Numbers whose binary weight (A000120) plus bigomega (A001222) is even.

Original entry on oeis.org

2, 6, 7, 8, 9, 10, 11, 13, 15, 19, 24, 28, 31, 32, 33, 34, 36, 37, 39, 40, 41, 42, 44, 46, 47, 50, 51, 52, 54, 57, 58, 59, 60, 61, 65, 67, 70, 73, 76, 77, 79, 85, 86, 90, 95, 96, 97, 98, 103, 106, 107, 109, 110, 111, 112, 117, 119, 123, 124, 126, 127, 128, 129
Offset: 1

Views

Author

Gus Wiseman, May 14 2024

Keywords

Comments

The odd version is A372590.

Examples

			The terms (center), their binary indices (left), and their weakly decreasing prime indices (right) begin:
          {2}   2  (1)
        {2,3}   6  (2,1)
      {1,2,3}   7  (4)
          {4}   8  (1,1,1)
        {1,4}   9  (2,2)
        {2,4}  10  (3,1)
      {1,2,4}  11  (5)
      {1,3,4}  13  (6)
    {1,2,3,4}  15  (3,2)
      {1,2,5}  19  (8)
        {4,5}  24  (2,1,1,1)
      {3,4,5}  28  (4,1,1)
  {1,2,3,4,5}  31  (11)
          {6}  32  (1,1,1,1,1)
        {1,6}  33  (5,2)
        {2,6}  34  (7,1)
        {3,6}  36  (2,2,1,1)
      {1,3,6}  37  (12)
    {1,2,3,6}  39  (6,2)
        {4,6}  40  (3,1,1,1)
      {1,4,6}  41  (13)
      {2,4,6}  42  (4,2,1)
		

Crossrefs

For sum (A372428, zeros A372427) we have A372587, complement A372586.
For minimum (A372437) we have A372440, complement A372439.
Positions of even terms in A372441, zeros A071814.
For maximum (A372442, zeros A372436) we have A372589, complement A372588.
The complement is A372590.
For just binary indices:
- length: A001969, complement A000069
- sum: A158704, complement A158705
- minimum: A036554, complement A003159
- maximum: A053754, complement A053738
For just prime indices:
- length: A026424 A028260 (count A027187), complement (count A027193)
- sum: A300061 (count A058696), complement A300063 (count A058695)
- minimum: A340933 (count A026805), complement A340932 (count A026804)
- maximum: A244990 (count A027187), complement A244991 (count A027193)
A019565 gives Heinz number of binary indices, adjoint A048675.
A029837 gives greatest binary index, least A001511.
A031215 lists even-indexed primes, odd A031368.
A048793 lists binary indices, length A000120, reverse A272020, sum A029931.
A070939 gives length of binary expansion.
A112798 lists prime indices, length A001222, reverse A296150, sum A056239.

Programs

  • Mathematica
    Select[Range[100],EvenQ[DigitCount[#,2,1]+PrimeOmega[#]]&]

A379307 Positive integers whose prime indices include no squarefree numbers.

Original entry on oeis.org

1, 7, 19, 23, 37, 49, 53, 61, 71, 89, 97, 103, 107, 131, 133, 151, 161, 173, 193, 197, 223, 227, 229, 239, 251, 259, 263, 281, 307, 311, 337, 343, 359, 361, 371, 379, 383, 409, 419, 427, 433, 437, 457, 463, 479, 497, 503, 521, 523, 529, 541, 569, 593, 613, 623
Offset: 1

Views

Author

Gus Wiseman, Dec 27 2024

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:
    1: {}
    7: {4}
   19: {8}
   23: {9}
   37: {12}
   49: {4,4}
   53: {16}
   61: {18}
   71: {20}
   89: {24}
   97: {25}
  103: {27}
  107: {28}
  131: {32}
  133: {4,8}
  151: {36}
  161: {4,9}
  173: {40}
		

Crossrefs

Partitions of this type are counted by A114374, strict A256012.
Positions of zero in A379306.
For a unique squarefree part we have A379316, counted by A379308 (strict A379309).
A000040 lists the primes, differences A001223.
A005117 lists the squarefree numbers, differences A076259.
A008966 is the characteristic function for the squarefree numbers.
A013929 lists the nonsquarefree numbers, differences A078147.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798, counted by A001222.
A061398 counts squarefree numbers between primes, zeros A068360.
A377038 gives k-th differences of squarefree numbers.
Other counts of prime indices:
- A330944 nonprime, see A000586, A000607, A076610, A330945.
- A379310 nonsquarefree, see A302478.
- A379311 old prime, see A204389, A320629, A379312-A379315.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[1000],Length[Select[prix[#],SquareFreeQ]]==0&]

A379310 Number of nonsquarefree prime indices of n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Dec 27 2024

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 prime indices of 39 are {2,6}, so a(39) = 0.
The prime indices of 70 are {1,3,4}, so a(70) = 1.
The prime indices of 98 are {1,4,4}, so a(98) = 2.
The prime indices of 294 are {1,2,4,4}, a(294) = 2.
The prime indices of 1911 are {2,4,4,6}, so a(1911) = 2.
The prime indices of 2548 are {1,1,4,4,6}, so a(2548) = 2.
		

Crossrefs

Positions of first appearances are A000420.
Positions of zero are A302478, counted by A073576 (strict A087188).
No squarefree parts: A379307, counted by A114374 (strict A256012).
One squarefree part: A379316, counted by A379308 (strict A379309).
A000040 lists the primes, differences A001223.
A005117 lists the squarefree numbers, differences A076259.
A008966 is the characteristic function for the squarefree numbers.
A013929 lists the nonsquarefree numbers, differences A078147.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798, counted by A001222.
A061398 counts squarefree numbers between primes, zeros A068360.
A377038 gives k-th differences of squarefree numbers.
Other counts of prime indices:
- A330944 nonprime, see A000586, A000607, A076610, A330945.
- A379311 old prime, see A204389, A320629, A379312-A379315.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Length[Select[prix[n],Not@*SquareFreeQ]],{n,100}]

Formula

Totally additive with a(prime(k)) = A107078(k) = 1 - A008966(k).

A019507 Droll numbers: numbers > 1 whose sum of even prime factors equals the sum of odd prime factors.

Original entry on oeis.org

72, 240, 672, 800, 2240, 4224, 5184, 6272, 9984, 14080, 17280, 33280, 39424, 48384, 52224, 57600, 93184, 116736, 161280, 174080, 192000, 247808, 304128, 373248, 389120, 451584, 487424, 537600, 565248, 585728, 640000, 718848, 1013760, 1089536, 1244160, 1384448
Offset: 1

Views

Author

Mario Velucchi (mathchess(AT)velucchi.it)

Keywords

Examples

			6272 = 2*2*2*2*2*2*2*7*7 is droll since 2+2+2+2+2+2+2 = 14 = 7+7.
		

Crossrefs

For count instead of sum we have A072978.
Partitions of this type are counted by A239261, without zero terms A249914.
For prime indices instead of factors we have A366748, zeros of A366749.
The LHS is A366839 with alternating zeros, for indices A366531, triangle A113686.
The RHS is A366840, for indices A366528, triangle A113685.
A000009 counts partitions into odd parts, ranks A066208.
A035363 counts partitions into even parts, ranks A066207.
A112798 lists prime indices, length A001222, sum A056239.
A257991 counts odd prime indices, even A257992.
A300061 lists numbers with even sum of prime indices, odd A300063.

Programs

  • Maple
    f:= proc(k, m) # numbers whose sum of prime factors >= m is k; m is prime
       local S,p,j;
       option remember;
       if k = 0 then return [1]
       elif m > k then return []
       fi;
       S:= NULL:
       p:= nextprime(m);
       for j from k by -m to 0 do
         S:= S, op(map(`*`,  procname(j,p) , m^((k-j)/m)))
       od;
       [S]
    end proc:
    g:= proc(N) local m,R;
      R:= NULL;
      for m from 1 while 2^m < N do
       R:= R, op(map(`*`,select(`<=`,f(2*m,3), N/2^m),2^m));
      od;
      sort([R])
    end proc:
    g(10^8); # Robert Israel, Feb 20 2025
  • Mathematica
    Select[Range[2, 2*10^6, 2], First[#] == Total[Rest[#]] & [Times @@@ FactorInteger[#]] &] (* Paolo Xausa, Feb 19 2025 *)
  • PARI
    isok(n) = {if (n % 2, return (0)); f = factor(n); return (2*f[1,2] == sum(i=2, #f~, f[i,1]*f[i,2]));} \\ Michel Marcus, Jun 21 2013

Formula

These are even numbers k such that A366839(k/2) = A366840(k). - Gus Wiseman, Oct 25 2023 (corrected Feb 19 2025)

Extensions

Name edited by Paolo Xausa, Feb 19 2025

A340933 Numbers whose least prime index is even. Heinz numbers of integer partitions whose last part is even.

Original entry on oeis.org

3, 7, 9, 13, 15, 19, 21, 27, 29, 33, 37, 39, 43, 45, 49, 51, 53, 57, 61, 63, 69, 71, 75, 77, 79, 81, 87, 89, 91, 93, 99, 101, 105, 107, 111, 113, 117, 119, 123, 129, 131, 133, 135, 139, 141, 147, 151, 153, 159, 161, 163, 165, 169, 171, 173, 177, 181, 183
Offset: 1

Views

Author

Gus Wiseman, Feb 12 2021

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. 1 has no prime indices so is not counted.

Examples

			The sequence of terms together with their prime indices begins:
      3: {2}         51: {2,7}         99: {2,2,5}
      7: {4}         53: {16}         101: {26}
      9: {2,2}       57: {2,8}        105: {2,3,4}
     13: {6}         61: {18}         107: {28}
     15: {2,3}       63: {2,2,4}      111: {2,12}
     19: {8}         69: {2,9}        113: {30}
     21: {2,4}       71: {20}         117: {2,2,6}
     27: {2,2,2}     75: {2,3,3}      119: {4,7}
     29: {10}        77: {4,5}        123: {2,13}
     33: {2,5}       79: {22}         129: {2,14}
     37: {12}        81: {2,2,2,2}    131: {32}
     39: {2,6}       87: {2,10}       133: {4,8}
     43: {14}        89: {24}         135: {2,2,2,3}
     45: {2,2,3}     91: {4,6}        139: {34}
     49: {4,4}       93: {2,11}       141: {2,15}
		

Crossrefs

These partitions are counted by A026805.
Looking at length or at maximum gives A028260/A244990, counted by A027187.
If all prime indices are even we get A066207, counted by A035363.
The complement is {1} \/ A340932, counted by A026804.
A001222 counts prime factors.
A005843 lists even numbers.
A031215 lists even-indexed primes.
A055396 selects least prime index.
A056239 adds up prime indices.
A058695 counts partitions of even numbers, ranked by A300061.
A061395 selects greatest prime index.
A112798 lists the prime indices of each positive integer.

Programs

  • Mathematica
    Select[Range[2,100],EvenQ[PrimePi[FactorInteger[#][[1,1]]]]&]

Formula

A055396(a(n)) belongs to A005843.
Closed under multiplication.

A379306 Number of squarefree prime indices of n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Dec 25 2024

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 prime indices of 39 are {2,6}, so a(39) = 2.
The prime indices of 70 are {1,3,4}, so a(70) = 2.
The prime indices of 98 are {1,4,4}, so a(98) = 1.
The prime indices of 294 are {1,2,4,4}, a(294) = 2.
The prime indices of 1911 are {2,4,4,6}, so a(1911) = 2.
The prime indices of 2548 are {1,1,4,4,6}, so a(2548) = 3.
		

Crossrefs

Positions of first appearances are A000079.
Positions of zero are A379307, counted by A114374 (strict A256012).
Positions of one are A379316, counted by A379308 (strict A379309).
A000040 lists the primes, differences A001223.
A005117 lists the squarefree numbers, differences A076259.
A008966 is the characteristic function for the squarefree numbers.
A013929 lists the nonsquarefree numbers, differences A078147.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798, counted by A001222.
A061398 counts squarefree numbers between primes, zeros A068360.
A377038 gives k-th differences of squarefree numbers.
Other counts of prime indices:
- A087436 postpositive, see A038550.
- A330944 nonprime, see A000586, A000607, A076610, A330945.
- A379310 nonsquarefree, see A302478.
- A379311 old prime, see A204389, A320629, A379312-A379315.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Length[Select[prix[n],SquareFreeQ]],{n,100}]

Formula

Totally additive with a(prime(k)) = A008966(k).

A379317 Positive integers with a unique even prime index.

Original entry on oeis.org

3, 6, 7, 12, 13, 14, 15, 19, 24, 26, 28, 29, 30, 33, 35, 37, 38, 43, 48, 51, 52, 53, 56, 58, 60, 61, 65, 66, 69, 70, 71, 74, 75, 76, 77, 79, 86, 89, 93, 95, 96, 101, 102, 104, 106, 107, 112, 113, 116, 119, 120, 122, 123, 130, 131, 132, 138, 139, 140, 141, 142
Offset: 1

Views

Author

Gus Wiseman, Dec 29 2024

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}
   6: {1,2}
   7: {4}
  12: {1,1,2}
  13: {6}
  14: {1,4}
  15: {2,3}
  19: {8}
  24: {1,1,1,2}
  26: {1,6}
  28: {1,1,4}
  29: {10}
  30: {1,2,3}
  33: {2,5}
  35: {3,4}
  37: {12}
  38: {1,8}
  43: {14}
  48: {1,1,1,1,2}
		

Crossrefs

Partitions of this type are counted by A038348 (strict A096911).
For all even parts we have A066207, counted by A035363 (strict A000700).
For no even parts we have A066208, counted by A000009 (strict A035457).
Positions of 1 in A257992.
A000040 lists the primes, differences A001223.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798, counted by A001222.
Other counts of prime indices:
- A330944 nonprime, see A000586, A000607, A076610, A330945.
- A379311 old prime, see A204389, A320629, A379312-A379315.

Programs

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

A332822 One part of a 3-way classification of the positive integers. Numbers n for which A048675(n) == 2 (mod 3).

Original entry on oeis.org

3, 4, 7, 10, 13, 18, 19, 22, 24, 25, 29, 32, 34, 37, 42, 43, 45, 46, 53, 55, 56, 60, 61, 62, 71, 78, 79, 80, 81, 82, 85, 89, 94, 98, 99, 101, 104, 105, 107, 108, 113, 114, 115, 118, 121, 131, 132, 134, 139, 140, 144, 146, 150, 151, 152, 153, 155, 163, 166, 173, 174, 176, 181, 182, 187, 189, 192, 193, 194, 195, 199, 200, 204
Offset: 1

Views

Author

Antti Karttunen and Peter Munn, Feb 25 2020

Keywords

Comments

The positive integers are partitioned between A332820, A332821 and this sequence.
For each prime p, the terms include exactly one of p and p^2. The primes alternate between this sequence and A332821. This sequence has the primes with even indexes, those in A031215.
The terms are the even numbers in A332820 halved. The terms are also the numbers m such that 5m is in A332820, and so on for alternate primes: 11, 17, 23 etc. Likewise, the terms are the numbers m such that 3m is in A332821, and so on for alternate primes: 7, 13, 19, 29 etc.
If we take each odd term of this sequence and replace each prime in its factorization by the next smaller prime, we get the same set of numbers as we get from halving the even terms of this sequence, and A332821 consists exactly of those numbers. The numbers that are one third of the terms that are multiples of 3 are in A332820, which consists exactly of those numbers. The numbers that are one fifth of the terms that are multiples of 5 constitute A332821, and for larger primes, an alternating pattern applies as described in the previous paragraph.
The product of any 2 terms of this sequence is in A332821, the product of any 3 terms is in A332820, and the product of a term of A332820 and a term of this sequence is in this sequence. So if a number k is present, k^2 is in A332821, k^3 is in A332820, and k^4 is in this sequence.
If k is an even number, exactly one of {k/2, k, 2k} is in the sequence (cf. A191257 / A067368 / A213258); and generally if k is a multiple of a prime p, exactly one of {k/p, k, k*p} is in the sequence.

Crossrefs

Positions of terms valued -1 in A332823; equivalently, numbers in row 3k-1 of A277905 for some k >= 1.
Subsequences: intersection of A026478 and A066207, A031215 (prime terms), A033430\{0}, A117642\{0}, A169604, A244727\{0}, A244729\{0}, A338910 (semiprime terms).

Programs

  • Mathematica
    Select[Range@ 204, Mod[Total@ #, 3] == 2 &@ Map[#[[-1]]*2^(PrimePi@ #[[1]] - 1) &, FactorInteger[#]] &] (* Michael De Vlieger, Mar 15 2020 *)
  • PARI
    isA332822(n) =  { my(f = factor(n)); (2==((sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2)%3)); };

Formula

{a(n) : n >= 1} = {2 * A332821(k) : k >= 1} U {A003961(A332821(k)) : k >= 1}.
{a(n) : n >= 1} = {A332821(k)^2 : k >= 1} U {A331590(2, A332821(k)) : k >= 1}.

A340786 Number of factorizations of 4n into an even number of even factors > 1.

Original entry on oeis.org

1, 1, 1, 3, 1, 2, 1, 3, 2, 2, 1, 4, 1, 2, 2, 6, 1, 3, 1, 4, 2, 2, 1, 6, 2, 2, 2, 4, 1, 4, 1, 7, 2, 2, 2, 7, 1, 2, 2, 6, 1, 4, 1, 4, 3, 2, 1, 10, 2, 3, 2, 4, 1, 4, 2, 6, 2, 2, 1, 8, 1, 2, 3, 12, 2, 4, 1, 4, 2, 4, 1, 10, 1, 2, 3, 4, 2, 4, 1, 10, 3, 2, 1, 8, 2, 2
Offset: 1

Views

Author

Gus Wiseman, Jan 31 2021

Keywords

Examples

			The a(n) factorizations for n = 6, 12, 24, 36, 60, 80, 500:
  4*6   6*8      2*48      2*72      4*60      4*80          40*50
  2*12  2*24     4*24      4*36      6*40      8*40          4*500
        4*12     6*16      6*24      8*30      10*32         8*250
        2*2*2*6  8*12      8*18      10*24     16*20         10*200
                 2*2*4*6   12*12     12*20     2*160         20*100
                 2*2*2*12  2*2*6*6   2*120     2*2*2*40      2*1000
                           2*2*2*18  2*2*2*30  2*2*4*20      2*2*10*50
                                     2*2*6*10  2*2*8*10      2*2*2*250
                                               2*4*4*10      2*10*10*10
                                               2*2*2*2*2*10
		

Crossrefs

Note: A-numbers of Heinz-number sequences are in parentheses below.
Positions of ones are 1 and A000040, or A008578.
A version for partitions is A027187 (A028260).
Allowing odd length gives A108501 (bisection of A340785).
Allowing odd factors gives A339846.
An odd version is A340102.
- Factorizations -
A001055 counts factorizations, with strict case A045778.
A316439 counts factorizations by product and length.
A340101 counts factorizations into odd factors.
A340653 counts balanced factorizations.
A340831/A340832 count factorizations with odd maximum/minimum.
- Even -
A027187 counts partitions of even maximum (A244990).
A058696 counts partitions of even numbers (A300061).
A067661 counts strict partitions of even length (A030229).
A236913 counts partitions of even length and sum (A340784).
A340601 counts partitions of even rank (A340602).

Programs

  • Maple
    g:= proc(n, m, p)
     option remember;
     local F,r,x,i;
     # number of factorizations of n into even factors > m with number of factors == p (mod 2)
     if n = 1 then if p = 0 then return 1 else return 0 fi fi;
     if m > n  or n::odd then return 0 fi;
     F:= sort(convert(select(t -> t > m and t::even, numtheory:-divisors(n)),list));
     r:= 0;
     for x in F do
       for i from 1 while n mod x^i = 0 do
         r:= r + procname(n/x^i, x, (p+i) mod 2)
     od od;
     r
    end proc:
    f:= n -> g(4*n, 1, 0):
    map(f, [$1..100]); # Robert Israel, Mar 16 2023
  • 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[4n],EvenQ[Length[#]]&&Select[#,OddQ]=={}&]],{n,100}]
  • PARI
    A340786aux(n, m=n, p=0) = if(1==n, (0==p), my(s=0); fordiv(n, d, if((d>1)&&(d<=m)&&!(d%2), s += A340786aux(n/d, d, 1-p))); (s));
    A340786(n) = A340786aux(4*n); \\ Antti Karttunen, Apr 14 2022
Previous Showing 31-40 of 77 results. Next