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

A335738 Factorize each integer m >= 2 as the product of powers of nonunit squarefree numbers with distinct exponents that are powers of 2. The sequence lists m such that the factor with the largest exponent is a power of 2.

Original entry on oeis.org

2, 4, 8, 12, 16, 20, 24, 28, 32, 40, 44, 48, 52, 56, 60, 64, 68, 76, 80, 84, 88, 92, 96, 104, 112, 116, 120, 124, 128, 132, 136, 140, 144, 148, 152, 156, 160, 164, 168, 172, 176, 184, 188, 192, 204, 208, 212, 220, 224, 228, 232, 236, 240, 244, 248, 256, 260, 264, 268, 272
Offset: 1

Views

Author

Peter Munn, Jun 20 2020

Keywords

Comments

2 is the only term not divisible by 4. All powers of 2 are present. Every term divisible by an odd square is divisible by 16, the first such being 144.
The defined factorization is unique. Every positive number is a product of at most one squarefree number (A005117), at most one square of a squarefree number (A062503), at most one 4th power of a squarefree number (A113849), at most one 8th power of a squarefree number, and so on.
Iteratively map m using A000188, until 1 is reached, as A000188^k(m), for some k >= 1. m is in the sequence if and only if the preceding number, A000188^(k-1)(m), is 2. k can be shown to be A299090(m).
Closed under squaring, but not closed under multiplication: 12 = 3^1 * 2^2 and 432 = 3^1 * 3^2 * 2^4 are in the sequence, but 12 * 432 = 5184 = 3^4 * 2^6 = 2^2 * 6^4 is not.
The asymptotic density of this sequence is Sum_{k>=0} (d(2^(k+1)) - d(2^k))/2^(2^(k+1)-1) = 0.21363357193921052068..., where d(k) = 2^(k-1)/((2^k-1)*zeta(k)) is the asymptotic density of odd k-free numbers for k >= 2, and d(1) = 0. - Amiram Eldar, Feb 10 2024

Examples

			6 is a squarefree number, so its factorization for the definition (into powers of nonunit squarefree numbers with distinct exponents that are powers of 2) is the trivial "6^1". 6^1 is therefore the factor with the largest exponent, and is not a power of 2, so 6 is not in the sequence.
48 factorizes for the definition as 3^1 * 2^4. The factor with the largest exponent is 2^4, which is a power of 2, so 48 is in the sequence.
10^100 (a googol) factorizes in this way as 10^4 * 10^32 * 10^64. The factor with the largest exponent, 10^64, is a power of 10, not a power of 2, so 10^100 is not in the sequence.
		

Crossrefs

Complement within A020725 of A335740.
A000188, A299090 are used in a formula defining this sequence.
Powers of squarefree numbers: A005117(1), A144338(1), A062503(2), A113849(4).
Subsequences: A000079\{1}, A001749, A181818\{1}, A273798.
Numbers in the even bisection of A336322.
Row m of A352780 essentially gives the defined factorization of m.

Programs

  • Mathematica
    f[p_, e_] := p^Floor[e/2]; s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; Select[Range[2, 300], FixedPointList[s, #] [[-3]] == 2 &] (* Amiram Eldar, Nov 27 2020 *)
  • PARI
    is(n) = {my(e = valuation(n, 2), o = n >> e); if(e == 0, 0, if(o == 1, n > 1, floor(logint(e, 2)) > floor(logint(vecmax(factor(o)[,2]), 2))));} \\ Amiram Eldar, Feb 10 2024

Formula

{a(n)} = {m : m >= 2 and A000188^(k-1)(m) = 2, where k = A299090(m)}.
{a(n)} = {m : m >= 2 and A352780(m,e) = 2^(2^e), where e = A299090(m)-1}. - Peter Munn, Jun 24 2022

A340675 Exponential of Mangoldt function conjugated by Tek's flip: a(n) = A225546(A014963(A225546(n))).

Original entry on oeis.org

1, 2, 2, 4, 2, 2, 2, 1, 4, 2, 2, 1, 2, 2, 2, 16, 2, 1, 2, 1, 2, 2, 2, 1, 4, 2, 1, 1, 2, 2, 2, 1, 2, 2, 2, 4, 2, 2, 2, 1, 2, 2, 2, 1, 1, 2, 2, 1, 4, 1, 2, 1, 2, 1, 2, 1, 2, 2, 2, 1, 2, 2, 1, 1, 2, 2, 2, 1, 2, 2, 2, 1, 2, 2, 1, 1, 2, 2, 2, 1, 16, 2, 2, 1, 2, 2, 2, 1, 2, 1, 2, 1, 2, 2, 2, 1, 2, 1, 1, 4, 2, 2, 2, 1, 2
Offset: 1

Views

Author

Antti Karttunen and Peter Munn, Feb 01 2021

Keywords

Comments

Nonunit squarefree numbers take the value 2, other nonsquares take the value 1, and squares take the square of the value taken by their square root.

Crossrefs

Sequences used in a definition of this sequence: A014963, A048298, A225546, A267116, A297108, A340676.
Positions of 1's: {1} U A340681, 2's: A005117 \ {1}, of 4's: A062503 \ {1}, of 16's: A113849.
Positions of terms > 1: A340682, of terms > 2: A340674.
Sequences used to express relationship between terms of this sequence: A003961, A331590.

Programs

  • PARI
    A340675(n) = if(1==n,n,if(issquarefree(n), 2, if(!issquare(n), 1, A340675(sqrtint(n))^2)));

Formula

a(n) = 2^A048298(A267116(n)).
If A340673(n) = 1, then a(n) = 1, otherwise a(n) = 2^A297108(A340673(n)).
If A340676(n) = 0, then a(n) = 1, otherwise a(n) = 2^(2^(A340676(n)-1)).
If n = s^(2^k), s squarefree >= 2, k >= 0, then a(n) = 2^(2^k), otherwise a(n) = 1.
For n, k > 1, if a(n) = a(k) then a(A331590(n, k)) = a(n), otherwise a(A331590(n, k)) = 1.
a(n^2) = a(n)^2.
a(A003961(n)) = a(n).
a(A051144(n)) = 1.
a(n) = 1 if and only if A331591(n) <> 1, otherwise a(n) = 2^A051903(n).

A322554 Numbers whose product of prime indices is a power of a squarefree number (A072774).

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, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 36, 38, 40, 41, 42, 43, 44, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 62, 63, 64, 66, 67, 68, 72, 73, 76, 79, 80
Offset: 1

Views

Author

Gus Wiseman, Dec 15 2018

Keywords

Comments

The complement is {35, 37, 39, 45, 61, 65, ...}.
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. The multiset multisystem with MM-number n is formed by taking the multiset of prime indices of each part of the multiset of prime indices of n. For example, the prime indices of 78 are {1,2,6}, so the multiset multisystem with MM-number 78 is {{},{1},{1,2}}. This sequence lists all MM-numbers of regular multiset multisystems, where regularity means all vertex-degrees are equal.

Examples

			Most small numbers belong to this sequence. However, the sequence of multiset multisystems whose MM-numbers do not belong to this sequence begins:
  35: {{2},{1,1}}
  37: {{1,1,2}}
  39: {{1},{1,2}}
  45: {{1},{1},{2}}
  61: {{1,2,2}}
  65: {{2},{1,2}}
  69: {{1},{2,2}}
  70: {{},{2},{1,1}}
  71: {{1,1,3}}
  74: {{},{1,1,2}}
  75: {{1},{2},{2}}
  77: {{1,1},{3}}
  78: {{},{1},{1,2}}
  87: {{1},{1,3}}
  89: {{1,1,1,2}}
  90: {{},{1},{1},{2}}
  91: {{1,1},{1,2}}
  95: {{2},{1,1,1}}
  99: {{1},{1},{3}}
		

Crossrefs

Programs

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

A335324 Square part of 4th-power-free part of n.

Original entry on oeis.org

1, 1, 1, 4, 1, 1, 1, 4, 9, 1, 1, 4, 1, 1, 1, 1, 1, 9, 1, 4, 1, 1, 1, 4, 25, 1, 9, 4, 1, 1, 1, 1, 1, 1, 1, 36, 1, 1, 1, 4, 1, 1, 1, 4, 9, 1, 1, 1, 49, 25, 1, 4, 1, 9, 1, 4, 1, 1, 1, 4, 1, 1, 9, 4, 1, 1, 1, 4, 1, 1, 1, 36, 1, 1, 25, 4, 1, 1, 1, 1, 1, 1, 1, 4, 1
Offset: 1

Views

Author

Peter Munn, May 31 2020

Keywords

Comments

Equivalently, biquadratefree (4th-power-free) part of square part of n.
Multiplicative. The terms are squares of squarefree numbers (A062503).
Every positive integer n is the product of a unique subset S_n of the terms of A050376 (sometimes called Fermi-Dirac primes). a(n) is the product of the members of S_n that are squares of prime numbers (A001248).

Examples

			Removing the 4th powers from 192 = 2^6 * 3^1 gives 2^(6 - 4) * 3^1 = 2^2 * 3 = 12. So the 4th-power-free part of 192 is 12. The square part of 12 (largest square dividing 12) is 4. So a(192) = 4.
		

Crossrefs

A007913, A008833, A008835, A053165 are used in formulas defining the sequence.
Column 1 of A352780.
Range of values is A062503.
Positions of 1's: A252895.
Related to A038500 by A225546.
The formula section details how the sequence maps the terms of A003961, A331590.

Programs

  • Mathematica
    f[p_, e_] := p^(2*Floor[e/2] - 4*Floor[e/4]); a[n_] := Times @@ (f @@@ FactorInteger[n]); Array[a, 100] (* Amiram Eldar, Jun 01 2020 *)
  • PARI
    A053165(n)=my(f=factor(n)); f[, 2]=f[, 2]%4; factorback(f);
    a(n) = my(m=A053165(n)); m/core(m); \\ Michel Marcus, Jun 01 2020
    
  • Python
    from math import prod
    from sympy import factorint
    def A335324(n): return prod(p**(e&2) for p, e in factorint(n).items()) # Chai Wah Wu, Aug 07 2024

Formula

a(n) = A053165(A008833(n)) = A008833(A053165(n)).
a(n) = A053165(n) / A007913(n).
a(n) = A008833(n) / A008835(n).
n = A007913(n) * a(n) * A008835(n).
a(n) = A225546(A038500(A225546(n))).
a(n^2) = A007913(n)^2.
a(A003961(n)) = A003961(a(n)).
a(A331590(n, k)) = A331590(a(n), a(k)).
a(p^e) = p^(2*floor(e/2) - 4*floor(e/4)). - Amiram Eldar, Jun 01 2020
From Amiram Eldar, Sep 21 2023: (Start)
Dirichlet g.f.: zeta(s) * zeta(2*s-2) * zeta(4*s)/(zeta(2*s) * zeta(4*s-4)).
Sum_{k=1..n} a(k) ~ (4*zeta(3/2)*zeta(4))/(21*zeta(3)) * n^(3/2). (End)

A335740 Factorize each integer m >= 2 as the product of powers of nonunit squarefree numbers with distinct exponents that are powers of 2. The sequence lists m such that the factor with the largest exponent is not a power of 2.

Original entry on oeis.org

3, 5, 6, 7, 9, 10, 11, 13, 14, 15, 17, 18, 19, 21, 22, 23, 25, 26, 27, 29, 30, 31, 33, 34, 35, 36, 37, 38, 39, 41, 42, 43, 45, 46, 47, 49, 50, 51, 53, 54, 55, 57, 58, 59, 61, 62, 63, 65, 66, 67, 69, 70, 71, 72, 73, 74, 75, 77, 78, 79, 81, 82, 83, 85, 86, 87, 89, 90
Offset: 1

Views

Author

Peter Munn, Jun 20 2020

Keywords

Comments

Every missing number greater than 2 is a multiple of 4. Every power of 2 is missing. Every positive power of every squarefree number greater than 2 is present.
The defined factorization is unique. Every positive number is a product of at most one squarefree number (A005117), at most one square of a squarefree number (A062503), at most one 4th power of a squarefree number (A113849), at most one 8th power of a squarefree number, and so on.
Iteratively map m using A000188, until 1 is reached, as A000188^k(m), for some k >= 1. m is in the sequence if and only if the preceding number, A000188^(k-1)(m), is greater than 2. k can be shown to be A299090(m).
The asymptotic density of this sequence is 1 - Sum_{k>=0} (d(2^(k+1)) - d(2^k))/2^(2^(k+1)-1) = 0.78636642806078947931..., where d(k) = 2^(k-1)/((2^k-1)*zeta(k)) is the asymptotic density of odd k-free numbers for k >= 2, and d(1) = 0. - Amiram Eldar, Feb 10 2024

Examples

			6 is a squarefree number, so its factorization for the definition (into powers of nonunit squarefree numbers with distinct exponents that are powers of 2) is the trivial "6^1". 6^1 is therefore the factor with the largest exponent, and is not a power of 2, so 6 is in the sequence.
48 factorizes for the definition as 3^1 * 2^4. The factor with the largest exponent is 2^4, which is a power of 2, so 48 is not in the sequence.
10^100 (a googol) factorizes in this way as 10^4 * 10^32 * 10^64. The factor with the largest exponent, 10^64, is a power of 10, not a power of 2, so 10^100 is in the sequence.
		

Crossrefs

Complement within A020725 of A335738.
A000188, A299090 are used in a formula defining this sequence.
Powers of squarefree numbers: A005117(1), A144338(1), A062503(2), A113849(4).
Subsequences: A042968\{1,2}, A182853, A268390.
With {1}, numbers in the odd bisection of A336322.

Programs

  • Mathematica
    f[p_, e_] := p^Floor[e/2]; s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; Select[Range[2, 100], FixedPointList[s, #] [[-3]] > 2 &] (* Amiram Eldar, Nov 27 2020 *)
  • PARI
    is(n) = {my(e = valuation(n, 2), o = n >> e); if(e == 0, n > 1, if(o == 1, e < 1, floor(logint(e, 2)) <= floor(logint(vecmax(factor(o)[,2]), 2))));} \\ Amiram Eldar, Feb 10 2024

Formula

{a(n)} = {m : m >= 2 and A000188^(k-1)(m) > 2, where k = A299090(m)}.

A351011 Numbers k such that the k-th composition in standard order has even length and alternately equal and unequal parts, i.e., all run-lengths equal to 2.

Original entry on oeis.org

0, 3, 10, 36, 43, 58, 136, 147, 228, 235, 528, 547, 586, 676, 698, 904, 915, 2080, 2115, 2186, 2347, 2362, 2696, 2707, 2788, 2795, 3600, 3619, 3658, 3748, 3770, 8256, 8323, 8458, 8740, 8747, 8762, 9352, 9444, 9451, 10768, 10787, 10826, 11144, 11155, 14368
Offset: 1

Views

Author

Gus Wiseman, Feb 03 2022

Keywords

Comments

The k-th composition in standard order (graded reverse-lexicographic, A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again. This gives a bijective correspondence between nonnegative integers and integer compositions.

Examples

			The terms together with their binary expansions and standard compositions begin:
    0:           0  ()
    3:          11  (1,1)
   10:        1010  (2,2)
   36:      100100  (3,3)
   43:      101011  (2,2,1,1)
   58:      111010  (1,1,2,2)
  136:    10001000  (4,4)
  147:    10010011  (3,3,1,1)
  228:    11100100  (1,1,3,3)
  235:    11101011  (1,1,2,2,1,1)
  528:  1000010000  (5,5)
  547:  1000100011  (4,4,1,1)
  586:  1001001010  (3,3,2,2)
  676:  1010100100  (2,2,3,3)
  698:  1010111010  (2,2,1,1,2,2)
  904:  1110001000  (1,1,4,4)
  915:  1110010011  (1,1,3,3,1,1)
		

Crossrefs

The case of twins (binary weight 2) is A000120.
All terms are evil numbers A001969.
These compositions are counted by A003242 interspersed with 0's.
Partitions of this type are counted by A035457, any length A351005.
The Heinz numbers of these compositions are A062503.
Taking singles instead of twins gives A333489, complement A348612.
This is the anti-run case of A351010.
The strict case (distinct twins) is A351009, counted by A077957(n-2).
A011782 counts compositions.
A085207/A085208 represent concatenation of standard compositions.
A345167 ranks alternating compositions, counted by A025047.
A350355 ranks up/down compositions, counted by A025048.
A350356 ranks down/up compositions, counted by A025049.
A351014 counts distinct runs in standard compositions.
Selected statistics of standard compositions:
- Length is A000120.
- Sum is A070939.
- Heinz number is A333219.
- Number of distinct parts is A334028.
Selected classes of standard compositions:
- Partitions are A114994, strict A333256.
- Multisets are A225620, strict A333255.
- Strict compositions are A233564.
- Constant compositions are A272919.

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join@@ Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    Select[Range[0,1000],And@@(#==2&)/@Length/@Split[stc[#]]&]

A360902 Numbers with the same number of squarefree divisors and powerful divisors.

Original entry on oeis.org

1, 4, 9, 25, 36, 48, 49, 80, 100, 112, 121, 162, 169, 176, 196, 208, 225, 272, 289, 304, 361, 368, 405, 441, 464, 484, 496, 529, 567, 592, 656, 676, 688, 720, 752, 841, 848, 891, 900, 944, 961, 976, 1008, 1053, 1072, 1089, 1136, 1156, 1168, 1200, 1225, 1250, 1264
Offset: 1

Views

Author

Amiram Eldar, Feb 25 2023

Keywords

Comments

Numbers k such that A034444(k) = A005361(k).
Numbers whose squarefree kernel (A007947) and powerful part (A057521) have the same number of divisors (A000005).
If k and m are coprime terms, then k*m is also a term.
All the terms are exponentially 2^n-numbers (A138302).
The characteristic function of this sequence depends only on prime signature.
Numbers whose canonical prime factorization has exponents whose geometric mean is 2.
Equivalently, numbers of the form Product_{i=1..m} p_i^(2^k_i), where p_i are distinct primes, and Sum_{i=1..m} k_i = m (i.e., the exponents k_i have an arithmetic mean 1).
1 is the only squarefree (A005117) term.
Includes the squares of squarefree numbers (A062503), which are the powerful (A001694) terms of this sequence.
The squares of primes (A001248) are the only terms that are prime powers (A246655).
Numbers of the for m*p^(2^k), where m is squarefree, p is prime, gcd(m, p) = 1 and omega(m) = k - 1, are all terms. In particular, this sequence includes numbers of the form p^4*q, where p != q are primes (A178739), and numbers of the form p^8*q*r where p, q, and r are distinct primes (A179747).
The corresponding numbers of squarefree (or powerful) divisors are 1, 2, 2, 2, 4, 4, 2, 4, 4, 4, 2, 4, ... . The least term with 2^k squarefree divisors is A360903(k).

Examples

			4 is a term since it has 2 squarefree divisors (1 and 2) and 2 powerful divisors (1 and 4).
36 is a term since it has 4 squarefree divisors (1, 2, 3 and 6) and 4 powerful divisors (1, 4, 9 and 36).
		

Crossrefs

Programs

  • Mathematica
    q[n_] := Module[{e = FactorInteger[n][[;; , 2]]}, Times @@ e == 2^Length[e]]; q[1] = True; Select[Range[1300], q]
  • PARI
    is(k) = {my(e = factor(k)[,2]); prod(i = 1, #e, e[i]) == 2^#e; }

A370329 a(n) is the number of coreful divisors of the n-th powerful number that are also powerful numbers.

Original entry on oeis.org

1, 1, 2, 1, 3, 1, 2, 4, 1, 1, 5, 2, 3, 1, 2, 1, 2, 6, 3, 1, 1, 2, 4, 1, 4, 7, 4, 1, 3, 2, 1, 2, 3, 6, 1, 1, 2, 8, 1, 5, 3, 6, 2, 1, 5, 3, 4, 1, 8, 1, 1, 2, 4, 4, 9, 1, 2, 6, 1, 1, 9, 2, 2, 2, 1, 2, 1, 1, 4, 5, 1, 10, 1, 2, 1, 3, 8, 6, 3, 10, 1, 6, 2, 1, 7, 2, 3
Offset: 1

Views

Author

Amiram Eldar, Feb 15 2024

Keywords

Comments

A coreful divisor d of a number n is a divisor with the same set of distinct prime factors as n (see A307958).
The positive terms of A361430.

Crossrefs

Cf. A001694, A002117, A062503, A078434, A307958, A360908 (analogous with squares), A361430, A370328.

Programs

  • Mathematica
    f[p_, e_] := e - 1; s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; With[{max = 10^4}, s /@ Union@ Flatten@ Table[i^2*j^3, {j, 1, max^(1/3)}, {i, 1, Sqrt[max/j^3]}]]
  • PARI
    lista(kmax) = {my(e); for(k = 1, kmax, e = factor(k)[,2]; if(k == 1 || vecmin(e) > 1, print1(prod(i = 1, #e, e[i]-1), ", ")));}

Formula

a(n) = A361430(A001694(n)).
a(n) = 1 if and only if n is the square of a squarefree number (A062503).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = zeta(3/2) * zeta(3) * Product_{p prime} (1 + 2/p^2 + 2/p^(5/2) - 1/p^3 - 2/p^(7/2) - 2/p^4) = 6.91748056612108993003... . (The infinite product of primes is the value of f(1/2) in A361430).

A340673 If n is of the form s^(2^e), where s is a squarefree number, and e >= 0, then a(n) is the (1+e)-th prime, otherwise a(n) = 1.

Original entry on oeis.org

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

Views

Author

Antti Karttunen, Jan 31 2021

Keywords

Crossrefs

Positions of 1's: {1} U A340681, of 2's: A005117 \ {1}, of 3's: A062503 \ {1}, of 5's: A113849.
Positions of terms > 1: A340682, of terms > 2: A340674.

Programs

Formula

a(n) = A014963(A225546(n)).
a(n) = A225546(A340675(n)).
a(n) = A008578(1+A340676(n)).
If n is of the form s^(2^e), where s is a squarefree number, and e >= 0, then a(n) = A000040(1+e), otherwise a(n) = 1.

A340676 If n is of the form s^(2^e), where s is a squarefree number, and e >= 0, then a(n) = 1+e, otherwise a(n) = 0.

Original entry on oeis.org

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

Views

Author

Antti Karttunen, Feb 01 2021

Keywords

Crossrefs

Positions of zeros: {1} U A340681, of 1's: A005117 \ {1}, of 2's: A062503 \ {1}, of 3's: A113849.
Positions of nonzero terms: A340682, of terms > 1: A340674.

Programs

  • Mathematica
    a[1] = 0; a[n_] := If[Length[(u = Union[FactorInteger[n][[;; , 2]]])] == 1 && u[[1]] == 2^(e = IntegerExponent[u[[1]], 2]), e + 1, 0]; Array[a, 100] (* Amiram Eldar, Feb 10 2021 *)
  • PARI
    A001511(n) = 1+valuation(n,2);
    A209229(n) = (n && !bitand(n,n-1));
    A104117(n) = (A209229(n)*A001511(n));
    A267116(n) = if(n>1, fold(bitor, factor(n)[, 2]), 0);
    A340676(n) = if(1==n,0,A104117(A267116(n)));

Formula

a(n) = A297109(A225546(n)).
For n > 1, a(n) = A104117(A267116(n)). - Peter Munn, Feb 05 2021
Previous Showing 31-40 of 62 results. Next