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

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.

Links

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.
Cf. A028260, A045931, A100367, A106529, A195017, A239241, A241638, A325698, A325700, A366533.

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

A366840 Sum of odd prime factors of n, counted with multiplicity.

Original entry on oeis.org

0, 0, 3, 0, 5, 3, 7, 0, 6, 5, 11, 3, 13, 7, 8, 0, 17, 6, 19, 5, 10, 11, 23, 3, 10, 13, 9, 7, 29, 8, 31, 0, 14, 17, 12, 6, 37, 19, 16, 5, 41, 10, 43, 11, 11, 23, 47, 3, 14, 10, 20, 13, 53, 9, 16, 7, 22, 29, 59, 8, 61, 31, 13, 0, 18, 14, 67, 17, 26, 12, 71, 6
Offset: 1

Views

Author

Gus Wiseman, Oct 27 2023

Keywords

Comments

Contains all positive integers except 1, 2, 4.

Examples

			The prime factors of 60 are {2,2,2,3,5}, of which the odd factors are {3,5}, so a(60) = 8.

Crossrefs

The compound version is A001414, triangle A331416.
For count instead of sum we have A087436, even version A007814.
Odd-indexed terms are A100005.
Positions of odd terms are A335657, even A036349.
For prime indices we have A366528, triangle A113685 (without zeros A365067)
The even version is A366839 = 2*A001511.
The partition triangle for this statistic is A366851, even version A116598.
A019507 lists numbers with (even factor sum) = (odd factor sum).
A066207 lists numbers with all even prime indices, counted by A035363.
A112798 lists prime indices, reverse A296150, length A001222, sum A056239.
A162641 counts even prime exponents, odd A162642.
A239261 counts partitions with (sum of odd parts) = (sum of even parts).
A257992 counts even prime indices, odd A257991.
Cf. A000009, A066208, A113686, A174713, A258117, A325698, A366531, A366850.

Programs

  • Mathematica
    Table[Total[Times@@@DeleteCases[If[n==1,{}, FactorInteger[n]],{2,_}]],{n,100}]
  • PARI
    a(n) = my(f=factor(n), j=if(n%2, 1, 2)); sum(i=j, #f~, f[i,1]*f[i,2]); \\ Michel Marcus, Oct 30 2023

Formula

a(n) = A100006(n) - A366839(n).
a(2n) = a(n).
a(2n-1) = A001414(2n-1) = A100005(n).
Completely additive with a(2^e) = 0 and a(p^e) = e*p for an odd prime p. - Amiram Eldar, Nov 03 2023

A366841 Least positive integer whose odd prime factors sum to n, starting with n = 5.

Original entry on oeis.org

5, 9, 7, 15, 27, 21, 11, 35, 13, 33, 105, 39, 17, 65, 19, 51, 195, 57, 23, 95, 171, 69, 285, 115, 29, 161, 31, 87, 483, 93, 261, 155, 37, 217, 465, 111, 41, 185, 43, 123, 555, 129, 47, 215, 387, 141, 645, 235, 53, 329, 705, 159, 987, 265, 59, 371, 61, 177
Offset: 5

Views

Author

Gus Wiseman, Oct 27 2023

Keywords

Comments

All terms are odd.
It seems that all composite terms not divisible by 3 form a supersequence of A292081. - Ivan N. Ianakiev, Oct 30 2023

Examples

			The terms together with their prime factors (which sum to n) begin:
    5 = 5
    9 = 3*3
    7 = 7
   15 = 3*5
   27 = 3*3*3
   21 = 3*7
   11 = 11
   35 = 5*7
   13 = 13
   33 = 3*11
  105 = 3*5*7

Links

Crossrefs

This is the odd case of A056240.
Positions of first appearances in A366840 (sum of odd prime factors).
The partition triangle for this statistic is A366851, even A116598.
A001414 adds up prime factors, triangle A331416.
A019507 lists numbers with (even factor sum) = (odd factor sum).
A027746 lists prime factors, length A001222.
A087436 counts odd prime factors, even A007814.
A366528 adds up odd prime indices, triangle A113685 (without zeros A365067).
Cf. A031368, A036349, A056239, A066208, A100005, A112798, A239261, A257991, A335657, A366839.

Programs

  • Mathematica
    nn=1000;
w=Table[Total[Times@@@DeleteCases[If[n==1,{},FactorInteger[n]],{2,_}]],{n,nn}];
spnm[y_]:=Max@@Select[Union[y],Function[i,Union[Select[y,#<=i&]]==Range[i]]];
Table[Position[w,k][[1,1]],{k,5,spnm[Join[{1,2,3,4},Take[w,nn]/.(0->1)]]}]
  • PARI
    f(n) = my(f=factor(n), j=if (n%2, 1, 2)); sum(i=j, #f~, f[i,1]*f[i,2]); \\ A366840
a(n) = my(k=1); while (f(k) != n, k++); k; \\ Michel Marcus, Nov 02 2023
Showing 1-3 of 3 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.