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

A356635 Numbers k that can be written as the sum of 7 divisors of k (not necessarily distinct).

Original entry on oeis.org

7, 8, 9, 10, 12, 14, 15, 16, 18, 20, 21, 22, 24, 27, 28, 30, 32, 33, 35, 36, 39, 40, 42, 44, 45, 48, 49, 50, 52, 54, 55, 56, 60, 63, 64, 66, 68, 70, 72, 75, 77, 78, 80, 81, 84, 88, 90, 91, 96, 98, 99, 100, 102, 104, 105, 108, 110, 112, 114, 117, 119, 120, 126, 128, 130
Offset: 1

Views

Author

Wesley Ivan Hurt, Aug 18 2022

Keywords

Comments

If k is in the sequence then so is k*m for positive m. - David A. Corneth, Aug 19 2022
Numbers that are divisible by at least one of 7, 8, 9, 10, 12, 15, 22, 33, 39, 52, 55, 68, 102, 114, 138. For proof, see link. - Robert Israel, Sep 02 2022
The asymptotic density of this sequence is 17819629/37182145 = 0.479252... . - Amiram Eldar, Aug 08 2023

Examples

			10 is in the sequence since 10 = 2+2+2+1+1+1+1, where each summand divides 10.
		

Crossrefs

Numbers k that can be written as the sum of j divisors of k (not necessarily distinct) for j=1..10: A000027 (j=1), A299174 (j=2), A355200 (j=3), A354591 (j=4), A355641 (j=5), A356609 (j=6), this sequence (j=7), A356657 (j=8), A356659 (j=9), A356660 (j=10).

Programs

  • Maple
    filter:= n -> ormap(t -> n mod t = 0, [7, 8, 9, 10, 12, 15, 22, 33, 39, 52, 55, 68, 102, 114, 138]):
    select(filter, [$1..200]); # Robert Israel, Sep 02 2022
  • Mathematica
    q[n_, k_] := AnyTrue[Tuples[Divisors[n], k], Total[#] == n &]; Select[Range[130], q[#, 7] &] (* Amiram Eldar, Aug 19 2022 *)
  • PARI
    isok(k) = my(d=divisors(k)); forpart(p=k, if (setintersect(d, Set(p)) == Set(p), return(1)), , [7,7]); \\ Michel Marcus, Aug 19 2022

A356657 Numbers k that can be written as the sum of 8 divisors of k (not necessarily distinct).

Original entry on oeis.org

8, 10, 12, 14, 16, 18, 20, 22, 24, 26, 28, 30, 32, 36, 40, 42, 44, 48, 50, 52, 54, 56, 60, 64, 66, 68, 70, 72, 76, 78, 80, 84, 88, 90, 96, 98, 100, 102, 104, 108, 110, 112, 114, 120, 126, 128, 130, 132, 136, 138, 140, 144, 150, 152, 154, 156, 160, 162, 168, 170, 174, 176
Offset: 1

Views

Author

Wesley Ivan Hurt, Aug 20 2022

Keywords

Comments

Terms are even. Proof by contradiction. Suppose m = a(n) is odd. Then each divisor is odd. Adding 8 odd numbers gives an even number. A contradiction. - David A. Corneth, Sep 02 2022

Examples

			14 is in the sequence since 14 = 2+2+2+2+2+2+1+1, where each summand divides 14.
		

Crossrefs

Numbers k that can be written as the sum of j divisors of k (not necessarily distinct) for j=1..10: A000027 (j=1), A299174 (j=2), A355200 (j=3), A354591 (j=4), A355641 (j=5), A356609 (j=6), A356635 (j=7), this sequence (j=8), A356659 (j=9), A356660 (j=10).

Programs

  • PARI
    isok(k) = my(d=divisors(k)); forpart(p=k, if (setintersect(d, Set(p)) == Set(p), return(1)), , [8,8]); \\ Michel Marcus, Aug 21 2022
    
  • PARI
    is(n) = if(n % 2 == 1, return(0)); my(d = divisors(n)); forvec(x = vector(8, i, [1, #d-1]), s=sum(i=1, #x, d [x[i]]); if(n == s, print(vector(#x, j, d[x[j]]));return(1)), 1); 0 \\ David A. Corneth, Aug 21 2022

A356659 Numbers k that can be written as the sum of 9 divisors of k (not necessarily distinct).

Original entry on oeis.org

9, 10, 12, 14, 15, 16, 18, 20, 21, 22, 24, 25, 26, 27, 28, 30, 32, 33, 35, 36, 39, 40, 42, 44, 45, 48, 50, 51, 52, 54, 55, 56, 57, 60, 63, 64, 65, 66, 68, 70, 72, 75, 76, 77, 78, 80, 81, 84, 85, 88, 90, 92, 96, 98, 99, 100, 102, 104, 105, 108, 110, 112, 114, 117, 120, 125
Offset: 1

Views

Author

Wesley Ivan Hurt, Aug 20 2022

Keywords

Comments

If k is in the sequence then so is k*m. - David A. Corneth, Oct 08 2022

Examples

			14 is in the sequence since 14 = 2+2+2+2+2+1+1+1+1, where each summand divides 14.
		

Crossrefs

Numbers k that can be written as the sum of j divisors of k (not necessarily distinct) for j=1..10: A000027 (j=1), A299174 (j=2), A355200 (j=3), A354591 (j=4), A355641 (j=5), A356609 (j=6), A356635 (j=7), A356657 (j=8), this sequence (j=9), A356660 (j=10).

Programs

  • PARI
    upto(n) = { my(v = vector(n,i,-1), t = 0); for(i = 1, n, if(v[i] == -1, print1(i", "); v[i] = is(i, 9); if(v[i] == 1, for(j = 2, n \ i, v[i*j] = 1; ) ) ); ); select(x->x >= 1, v, 1); }
    is(n, {qd = 10}) = { my(d = divisors(n)); d = d[^#d]; forvec(x = vector(qd-1, i, [1, #d]), s = sum(i = 1, qd-1, d[x[i]]); if(n - s >= d[x[qd - 1]], if(n % (n - s) == 0, return(1); ) ) , 1 ); 0 } \\ David A. Corneth, Oct 08 2022

Formula

a(n + t) = a(n) + s for some finite t and s. - David A. Corneth, Oct 08 2022

A356660 Numbers k that can be written as the sum of 10 divisors of k (not necessarily distinct).

Original entry on oeis.org

10, 12, 14, 16, 18, 20, 22, 24, 26, 28, 30, 32, 34, 36, 40, 42, 44, 48, 50, 52, 54, 56, 60, 64, 66, 68, 70, 72, 76, 78, 80, 84, 88, 90, 92, 96, 98, 100, 102, 104, 108, 110, 112, 114, 116, 120, 126, 128, 130, 132, 136, 138, 140, 144, 150, 152, 154, 156, 160, 162
Offset: 1

Views

Author

Wesley Ivan Hurt, Aug 20 2022

Keywords

Comments

From David A. Corneth, Oct 08 2022: (Start)
All terms are even. Proof: suppose a term is odd. Then all divisors are odd. Adding 10 odd numbers gives an even number. A contradiction.
If k is a term then so is k*m for m >= 1. Proof: Multiply each divisor in this sum of 10 divisors that give k with m. Then each term is a divisor of k*m and their sum is k*m. (End)

Examples

			14 is in the sequence since 14 = 2+2+2+2+1+1+1+1+1+1, where each summand divides 14.
		

Crossrefs

Numbers k that can be written as the sum of j divisors of k (not necessarily distinct) for j=1..10: A000027 (j=1), A299174 (j=2), A355200 (j=3), A354591 (j=4), A355641 (j=5), A356609 (j=6), A356635 (j=7), A356657 (j=8), A356659 (j=9), this sequence (j=10).

Programs

  • PARI
    upto(n) = { my(v = vector(n,i,-1), t = 0); forstep(i = 2, n, 2, if(v[i] == -1, v[i] = is(i); if(v[i] == 1, for(j = 2, n \ i, v[i*j] = 1; ) ) ); ); select(x->x >= 1, v, 1); }
    is(n, {qd = 10}) = { my(d = divisors(n), res = 0); d = d[^#d]; forvec(x = vector(qd-1, i, [1, #d]), s = sum(i = 1, qd-1, d[x[i]]); if(n - s >= d[x[qd - 1]], if(n % (n - s) == 0,  return(1); ) ) , 1 ); 0 } \\ David A. Corneth, Oct 08 2022
  • Python
    from sympy import divisors
    def t_sum_of_n_div(n, target):
        out, p = [], divisors(n)[::-1][1:]
        def dfs(t, divs,  index_s, kk):
            if len(out)!=0 or kk>target:return
            if kk == target and t == 0:
                out.append(divs)
                return
            for i in range(index_s, len(p)):
                if t >= p[i]:
                    temp_divs = divs.copy()
                    temp_divs.append(p[i])
                    dfs(t-p[i], temp_divs, i, kk+1)
        dfs(n, [], 0, 0)
        return out
    terms = [i for i in range(2, 200) if len(t_sum_of_n_div(i,10))!=0]
    print(terms) # Gleb Ivanov, Sep 02 2022
    

A358838 Minimum number of jumps needed to go from slab 0 to slab n in Jane Street's infinite sidewalk.

Original entry on oeis.org

0, 1, 2, 5, 3, 6, 9, 4, 7, 10, 10, 5, 8, 8, 11, 11, 11, 6, 14, 9, 9, 12, 12, 12, 15, 12, 7, 18, 15, 10, 10, 10, 13, 13, 13, 13, 16, 16, 13, 16, 8, 19, 19, 16, 11, 11, 11, 11, 19, 14, 14, 14, 14, 14, 22, 17, 17, 17, 14, 17, 17, 9, 20, 20, 20, 17, 17, 12, 12, 12
Offset: 0

Views

Author

Frederic Ruget, Dec 02 2022

Keywords

Comments

Slabs on the sidewalk are numbered n = 0, 1, 2,... and each has a label L(n) = 1, 1, 2, 2, 3, 3,...
At a given slab, a jump can be taken forward or backward by L(n) places (but not back before slab 0).
.
For every n >= 0,
let L(n) = 1 + floor(n/2) -- the label on slab n,
let forward(n) = n + L(n) -- jumping forward,
if n > 0, let backward(n) = n - L(n) -- jumping backward,
let lambda(n) = floor((2/3)*n),
let mu(n) = 1 + 2*n,
let nu(n) = 2 + 2*n.
Observe that given n >= 0, there are exactly two ways of landing onto slab n with a direct jump backwards:
backward-jumping from mu(n) to n, and
backward-jumping from nu(n) to n.
If n is a multiple of 3, there is no other ways of jumping onto slab n. But if n is not a multiple of 3, there is one additional way:
forward-jumping from lambda(n) to n.
(Note that L is A008619, forward(n) == A006999(n+1), lambda is A004523, mu is A005408, nu is A299174.)
.
Every slab n > 0 is reachable from slab 0, since there always exists some slab s < n which reaches n by one or more jumps:
if n != 0 (mod 3), then s = lambda(n) = floor((2/3)*n) takes one forward jump to n,
if n == 0 (mod 3) but n != 0 (mod 9), then s = lambda o lambda o mu(n) = floor((8/9)*n) takes two forward jumps and one backward jump to n,
if n == 0 (mod 9), then s = lambda o lambda o nu(n) = floor((8/9)*n + 6/9) takes two forward jumps and one backward jump to n.
This demonstrates that the sequence never stops.
This also gives the following bounds:
a(n) <= 1 + (4/3)*n,
a(n) <= 6*log(n)/log(9/8).
.
The sequence is a surjective mapping N -> N, since given any n >= 0:
a(forward^n(0)) == n.

Examples

			For n=0, a(0) = 0 since it takes zero jump to go from slab 0 to slab 0.
For n=3, a(3) = 5 jumps is the minimum needed to go from slab 0 to slab 3:
.
        1st   2nd      3rd           4th
        jump  jump     jump          jump
        ->-   ->-   ---->----   ------->-------
       /   \ /   \ /         \ /               \
n     0     1     2     3     4     5     6     7     8  ...
L(n)  1     1     2     2     3     3     4     4     5  ...
                         \                     /
                          ----------<----------
                                 5th jump (backwards)
		

Crossrefs

Always jumping forwards yields A006999.
In the COMMENTS section, L is A008619, forward(n) == A006999(n+1), lambda is A004523, mu is A005408, nu is A299174.
For related sequences, see A360744-A360746 and A360593-A360595.

Programs

  • Python
    def a(n: int) -> int:
        import itertools
        if n < 0: raise Exception("n must be a nonnegative integer")
        if n == 0: return 0
        if n == 1: return 1
        visited = {0, 1}  # the slabs we have visited so far
        rings = [{0}, {1}]  # the slabs ordered by depth (min length of path from 0)
        for depth in itertools.count(2):
            new_ring = set()
            for slab in rings[depth - 1]:
                label = (slab >> 1) + 1
                for next_slab in {slab - label, slab + label}:
                    if not next_slab in visited:
                        if next_slab == n: return depth
                        visited.add(next_slab)
                        new_ring.add(next_slab)
            rings.append(new_ring)

A177712 Even numbers that have a nontrivial odd divisor.

Original entry on oeis.org

6, 10, 12, 14, 18, 20, 22, 24, 26, 28, 30, 34, 36, 38, 40, 42, 44, 46, 48, 50, 52, 54, 56, 58, 60, 62, 66, 68, 70, 72, 74, 76, 78, 80, 82, 84, 86, 88, 90, 92, 94, 96, 98, 100, 102, 104, 106, 108, 110, 112, 114, 116, 118, 120, 122, 124, 126, 130, 132, 134, 136, 138, 140
Offset: 1

Views

Author

Keywords

Comments

Numbers which can be expressed as a sum of a set of positive consecutive even numbers: sum_{i=m..m+k} A005843(i), m>=1, k>=1.
Differs from A054741, which contains 105 for example.
These are the numbers that are not free of odd prime factors, but are not odd. Compare with A051144, nonsquarefree nonsquares. The self-inverse function defined by A225546 maps the members of either set 1:1 onto the other set. - Peter Munn, Jul 31 2020 with edit Feb 14 2022

Examples

			6=2+4. 10=4+6. 12=2+4+6. 14=6+8. 18=4+6+8. 20=2+4+6+8. 22=10+12. 24=6+8+10.
		

Crossrefs

Intersection of A057716 and A299174.
Related to A051144 via A225546.

Programs

  • Mathematica
    z=200;lst={};Do[c=a;Do[c+=b;If[c<=2*z,AppendTo[lst,c]],{b,a-2,1,-2}], {a,2,z,2}];Union@lst
  • PARI
    isA177712(n) = (!(n%2)&&(0<#select(x -> x%2,factor(n)[,1]))); \\ Antti Karttunen, Jul 31 2020
    
  • PARI
    isA177712(n) = (!(n%2)&&bitand(n,n-1)); \\ Antti Karttunen, Jul 31 2020
    
  • Python
    def A177712(n): return n+(m:=n.bit_length())+(n>=(1<Chai Wah Wu, Jun 30 2024

Formula

a(n) = 2 * A057716(n).

Extensions

Definition moved into a comment by R. J. Mathar, Aug 15 2010
New name from Peter Munn, Jul 31 2020

A119432 Numbers k such that 2*phi(k) <= k.

Original entry on oeis.org

2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 26, 28, 30, 32, 34, 36, 38, 40, 42, 44, 46, 48, 50, 52, 54, 56, 58, 60, 62, 64, 66, 68, 70, 72, 74, 76, 78, 80, 82, 84, 86, 88, 90, 92, 94, 96, 98, 100, 102, 104, 105, 106, 108, 110, 112, 114, 116, 118, 120, 122, 124, 126, 128, 130
Offset: 1

Views

Author

Keywords

Comments

Equivalently, numbers k such that totient(k) <= cototient(k).
Using the primes up to 23 it is possible to show that this sequence has (lower) density greater than 0.51. - Charles R Greathouse IV, Oct 26 2015
The asymptotic density of this sequence is in the interval (0.51120, 0.51176) (Kobayashi, 2016, improving the bounds 0.5105 and 0.5241 that were given by Wall, 1972). - Amiram Eldar, Oct 15 2020

Crossrefs

Disjoint union of A119434 and A299174. - Amiram Eldar, Oct 15 2020

Programs

Formula

Elements of A054741 together with all 2^n for n>0.

A334748 Let p be the smallest odd prime not dividing the squarefree part of n. Multiply n by p and divide by the product of all smaller odd primes.

Original entry on oeis.org

3, 6, 5, 12, 15, 10, 21, 24, 27, 30, 33, 20, 39, 42, 7, 48, 51, 54, 57, 60, 35, 66, 69, 40, 75, 78, 45, 84, 87, 14, 93, 96, 55, 102, 105, 108, 111, 114, 65, 120, 123, 70, 129, 132, 135, 138, 141, 80, 147, 150, 85, 156, 159, 90, 165, 168, 95, 174, 177, 28, 183, 186, 189
Offset: 1

Views

Author

Peter Munn, May 09 2020

Keywords

Comments

A permutation of A028983.
A007417 (which has asymptotic density 3/4) lists index n such that a(n) = 3n. The sequence maps the terms of A007417 1:1 onto A145204\{0}, defining a bijection between them.
Similarly, bijections are defined from the odd numbers (A005408) to the nonsquare odd numbers (A088828), from the positive even numbers (A299174) to A088829, from A003159 to the nonsquares in A003159, and from A325424 to the nonsquares in A036668. The latter two bijections are between sets where membership depends on whether a number's squarefree part divides by 2 and/or 3.

Examples

			84 = 21*4 has squarefree part 21 (and square part 4). The smallest odd prime absent from 21 = 3*7 is 5 and the product of all smaller odd primes is 3. So a(84) = 84*5/3 = 140.
		

Crossrefs

Permutation of A028983.
Row 3, and therefore column 3, of A331590. Cf. A334747 (row 2).
A007913, A034386, A225546, A284723 are used in formulas defining the sequence.
The formula section details how the sequence maps the terms of A003961, A019565, A070826; and how f(a(n)) relates to f(n) for f = A008833, A048675, A267116; making use of A003986.
Subsequences: A016051, A145204\{0}, A329575.
Bijections are defined that relate to A003159, A005408, A007417, A036668, A088828, A088829, A299174, A325424.

Programs

  • PARI
    a(n) = {my(c=core(n), m=n); forprime(p=3, , if(c % p, m*=p; break, m/=p)); m;} \\ Michel Marcus, May 22 2020

Formula

a(n) = n * p / (A034386(p-1)/2), where p = A284723(A007913(n)).
a(n) = A334747(A334747(n)).
a(n) = A331590(3, n) = A225546(4 * A225546(n)).
a(2*n) = 2 * a(n).
a(A019565(n)) = A019565(n+2).
a(k * m^2) = a(k) * m^2.
a(A003961(n)) = A003961(A334747(n)).
a(A070826(n)) = prime(n+1).
A048675(a(n)) = A048675(n) + 2.
A008833(a(n)) = A008833(n).
A267116(a(n)) = A267116(n) OR 1, where OR denotes the bitwise operation A003986.
a(A007417(n)) = A145204(n+1) = 3 * A007417(n).

A356066 Numbers with a prime index that is not a prime-power. Complement of A355743.

Original entry on oeis.org

2, 4, 6, 8, 10, 12, 13, 14, 16, 18, 20, 22, 24, 26, 28, 29, 30, 32, 34, 36, 37, 38, 39, 40, 42, 43, 44, 46, 47, 48, 50, 52, 54, 56, 58, 60, 61, 62, 64, 65, 66, 68, 70, 71, 72, 73, 74, 76, 78, 79, 80, 82, 84, 86, 87, 88, 89, 90, 91, 92, 94, 96, 98, 100, 101
Offset: 1

Views

Author

Gus Wiseman, Jul 31 2022

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:
    2: {1}
    4: {1,1}
    6: {1,2}
    8: {1,1,1}
   10: {1,3}
   12: {1,1,2}
   13: {6}
   14: {1,4}
   16: {1,1,1,1}
   18: {1,2,2}
   20: {1,1,3}
   22: {1,5}
   24: {1,1,1,2}
		

Crossrefs

The complement is A355743, counted by A023894.
The squarefree complement is A356065, counted by A054685.
Allowing prime index 1 gives A356064, complement A302492.
A000688 counts factorizations into prime-powers, strict A050361.
A001222 counts prime-power divisors.
A034699 gives the maximal prime-power divisor.
A246655 lists the prime-powers (A000961 includes 1), towers A164336.
A355742 chooses a prime-power divisor of each prime index.

Programs

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

Formula

Union of A299174 and A356064.

A043703 Numbers whose base-14 representation has an even number of runs.

Original entry on oeis.org

14, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 76, 77, 78
Offset: 1

Views

Author

Keywords

Crossrefs

Programs

  • Mathematica
    Select[Range[80],EvenQ[Length[Split[IntegerDigits[#,14]]]]&] (* Harvey P. Dale, Aug 30 2020 *)
  • PARI
    is(n) = qstreaks(n)%2 == 0
    qstreaks(n, {b = 14}) = {my(d = digits(n, b)); sum(i = 2, #d, d[i]!=d[i-1])+1} \\ David A. Corneth, Aug 31 2020
Previous Showing 11-20 of 29 results. Next