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-10 of 12 results. Next

A005867 a(0) = 1; for n > 0, a(n) = (prime(n)-1)*a(n-1).

Original entry on oeis.org

1, 1, 2, 8, 48, 480, 5760, 92160, 1658880, 36495360, 1021870080, 30656102400, 1103619686400, 44144787456000, 1854081073152000, 85287729364992000, 4434961926979584000, 257227791764815872000, 15433667505888952320000
Offset: 0

Views

Author

Keywords

Comments

Local minima of Euler's phi function. - Walter Nissen
Number of potential primes in a modulus primorial(n+1) sieve. - Robert G. Wilson v, Nov 20 2000
Let p=prime(n) and let p# be the primorial (A002110), then it can be shown that any p# consecutive numbers have exactly a(n-1) numbers whose lowest prime factor is p. For a proof, see the "Proofs Regarding Primorial Patterns" link. For example, if we let p=7 and consider the interval [101,310] containing 210 numbers, we find the 8 numbers 119, 133, 161, 203, 217, 259, 287, 301. - Dennis Martin (dennis.martin(AT)dptechnology.com), Jul 16 2006
From Gary W. Adamson, Apr 21 2009: (Start)
Equals (-1)^n * (1, 1, 1, 2, 8, 48, ...) dot (-1, 2, -3, 5, -7, 11, ...).
a(6) = 480 = (1, 1, 1, 2, 8, 48) dot (-1, 2, -3, 5, -7, 11) = (-1, 2, -3, 10, -56, 528). (End)
It can be proved that there are at least T prime numbers less than N, where the recursive function T is: T = N - N*Sum_{i=0..T(sqrt(N))} A005867(i)/A002110(i). This can show for example that at least 0.16*N numbers are primes less than N for 29^2 > N > 23^2. - Ben Paul Thurston, Aug 23 2010
First column of A096294. - Eric Desbiaux, Jun 20 2013
Conjecture: The g.f. for the prime(n+1)-rough numbers (A000027, A005408, A007310, A007775, A008364, A008365, A008366, A166061, A166063) is x*P(x)/(1-x-x^a(n)+x^(a(n)+1)), where P(x) is an order a(n) polynomial with symmetric coefficients (i.e., c(0)=c(n), c(1)=c(n-1), ...). - Benedict W. J. Irwin, Mar 18 2016
a(n)/A002110(n+1) (primorial(n+1)) is the ratio of natural numbers whose smallest prime factor is prime(n+1); i.e., prime(n+1) coprime to A002110(n). So the ratio of even numbers to natural numbers = 1/2; odd multiples of 3 = 1/6; multiples of 5 coprime to 6 (A084967) = 2/30 = 1/15; multiples of 7 coprime to 30 (A084968) = 8/210 = 4/105; etc. - Bob Selcoe, Aug 11 2016
The 2-adic valuation of a(n) is A057773(n), being sum of the 2-adic valuations of the product terms here. - Kevin Ryde, Jan 03 2023
For n > 1, a(n) is the number of prime(n+1)-rough numbers in [1, primorial(prime(n))]. - Alexandre Herrera, Aug 29 2023

Examples

			a(3): the mod 30 prime remainder set sieve representation yields the remainder set: {1, 7, 11, 13, 17, 19, 23, 29}, 8 elements.
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A057773 (2-adic valuation).
Column 1 of A281890.

Programs

  • Haskell
    a005867 n = a005867_list !! n
    a005867_list = scanl (*) 1 a006093_list
    -- Reinhard Zumkeller, May 01 2013
  • Maple
    A005867 := proc(n)
        mul(ithprime(j)-1,j=1..n) ;
    end proc: # Zerinvary Lajos, Aug 24 2008, R. J. Mathar, May 03 2017
  • Mathematica
    Table[ Product[ EulerPhi[ Prime[ j ] ], {j, 1, n} ], {n, 1, 20} ]
    RecurrenceTable[{a[0]==1,a[n]==(Prime[n]-1)a[n-1]},a,{n,20}] (* Harvey P. Dale, Dec 09 2013 *)
    EulerPhi@ FoldList[Times, 1, Prime@ Range@ 18] (* Michael De Vlieger, Mar 18 2016 *)
  • PARI
    for(n=0, 22, print1(prod(k=1,n, prime(k)-1), ", "))
    

Formula

a(n) = phi(product of first n primes) = A000010(A002110(n)).
a(n) = Product_{k=1..n} (prime(k)-1) = Product_{k=1..n} A006093(n).
Sum_{n>=0} a(n)/A002110(n+1) = 1. - Bob Selcoe, Jan 09 2015
a(n) = A002110(n)-((1/A000040(n+1) - A038110(n+1)/A038111(n+1))*A002110(n+1)). - Jamie Morken, Mar 27 2019
a(n) = |Sum_{k=0..n} A070918(n,k)|. - Alois P. Heinz, Aug 18 2019
a(n) = A058251(n)/A060753(n+1). - Jamie Morken, Apr 25 2022
a(n) = A002110(n) - A016035(A002110(n)) - 1 for n >= 1. - David James Sycamore, Sep 07 2024
Sum_{n>=0} 1/a(n) = A345974. - Amiram Eldar, Jun 26 2025

Extensions

Offset changed to 0, Name changed, and Comments and Examples sections edited by T. D. Noe, Apr 04 2010

A007775 Numbers not divisible by 2, 3 or 5.

Original entry on oeis.org

1, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 49, 53, 59, 61, 67, 71, 73, 77, 79, 83, 89, 91, 97, 101, 103, 107, 109, 113, 119, 121, 127, 131, 133, 137, 139, 143, 149, 151, 157, 161, 163, 167, 169, 173, 179, 181, 187, 191, 193, 197, 199, 203, 209
Offset: 1

Views

Author

Keywords

Comments

Also numbers n such that the sum of the 4th powers of the first n positive integers is divisible by n, or A000538(n) = n*(n+1)(2*n+1)(3*n^2+3*n-1)/30 is divisible by n. - Alexander Adamchuk, Jan 04 2007
Also the 7-rough numbers: positive integers that have no prime factors less than 7. - Michael B. Porter, Oct 09 2009
a(n) mod 3 has period 8, repeating [1,1,2,1,2,1,2,2] = (n mod 2) + floor(((n-1) mod 8)/7) - floor(((n-2) mod 8)/7) + 1. floor(a(n)/3) is the set of numbers k such that k is congruent to {0,2,3,4,5,6,7,9} mod 10 = floor((5*n-2)/4)-floor((n mod 8)/6). - Gary Detlefs, Jan 08 2012
Numbers k such that C(k+3,3)==1 (mod k) and C(k+5,5)==1 (mod k). - Gary Detlefs, Sep 15 2013
a(n) mod 30 has period 8 repeating [1, 7, 11, 13, 17, 19, 23, 29]. The mean of these 8 numbers is 120/8 = 15. (a(n)-15) mod 30 has period 8 repeating [-14, -8, -4, -2, 2, 4, 8, 14]. One half of the absolute value produces the symmetric sequence [7, 4, 2, 1, 1, 2, 4, 7] = A061501(((n-1) + 16) mod 8). - Gary Detlefs, Sep 24 2013
a(n) are exactly those positive integers m such that the sequence b(n) = n*(n + m)*(n + 2*m)*(n + 3*m)(n + 4*m)/120 is integral. Cf. A007310. - Peter Bala, Nov 13 2015
The asymptotic density of this sequence is 4/15. - Amiram Eldar, Sep 30 2020
If a(n) + a(n+1) = 0 (mod 30), then a(n-j) + a(n+j+1) = a(n) + a(n+1) for each j in [1, n-1]. - Alexandre Herrera, Jun 27 2023

Crossrefs

Cf. A000538, A054403, A145011 (first differences).
For k-rough numbers with other values of k, see A000027, A005408, A007310, A007775, A008364, A008365, A008366, A166061, A166063.
Complement is A080671.
For digital root of Fibonacci numbers indexed by this sequence, see A227896.

Programs

  • Haskell
    a007775 n = a007775_list !! (n-1)
    a007775_list = 1 : filter ((> 5) . a020639) [1..]
    -- Reinhard Zumkeller, Jan 06 2013
    
  • Magma
    I:=[1, 7, 11, 13, 17, 19, 23, 29, 31]; [n le 9 select I[n] else Self(n-1) +Self(n-8) - Self(n-9): n in [1..80]]; // G. C. Greubel, Oct 22 2018
    
  • Maple
    for i from 1 to 500 do if gcd(i,30) = 1 then print(i); fi; od;
    for k from 1 to 300 do if ((k^2 mod 48=1) or (k^2 mod 48=25)) and ((k^2 mod 120=1) or (k^2 mod 120=49)) then print(k) fi od. # Gary Detlefs, Dec 30 2011
  • Mathematica
    Select[ Range[ 300 ], GCD[ #1, 30 ]==1& ]
    Select[Range[250], Mod[#, 2]>0&&Mod[#, 3]>0&&Mod[#, 5]>0&] (* Vincenzo Librandi, Feb 08 2014 *)
    a[ n_] := Quotient[ n, 8, 1] 30 + {1, 7, 11, 13, 17, 19, 23, 29}[[Mod[n, 8, 1]]]; (* Michael Somos, Jun 02 2014 *)
    LinearRecurrence[{1, 0, 0, 0, 0, 0, 0, 1, -1}, {1, 7, 11, 13, 17, 19, 23, 29, 31}, 100] (* Mikk Heidemaa, Dec 07 2017 *)
    Cases[Range@1000, x_ /; NoneTrue[Array[Prime, 3], Divisible[x, #] &]] (* Mikk Heidemaa, Dec 07 2017 *)
    CoefficientList[ Series[(x^8 + 6x^7 + 4x^6 + 2x^5 + 4x^4 + 2x^3 + 4x^2 + 6x + 1)/((x - 1)^2 (x^7 + x^6 + x^5 + x^4 + x^3 + x^2 + x + 1)), {x, 0, 55}], x] (* Robert G. Wilson v, Dec 07 2017 *)
  • PARI
    isA007775(n) = gcd(n,30)==1 \\ Michael B. Porter, Oct 09 2009
    
  • PARI
    {a(n) = n\8 * 30 + [ -1, 1, 7, 11, 13, 17, 19, 23][n%8 + 1]} /* Michael Somos, Feb 05 2011 */
    
  • PARI
    {a(n) = n\8 * 6 + 9 + 3 * (n+1)\2 * 2 - max(5, (n-2)%8) * 2} /* Michael Somos, Jun 02 2014 */
    
  • PARI
    Vec(x*(1+6*x+4*x^2+2*x^3+4*x^4+2*x^5+4*x^6+6*x^7+x^8)/((1+x)*(x^2+1)*(x^4+1)*( x-1)^2) + O(x^100)) \\ Altug Alkan, Nov 16 2015
    
  • Python
    def A007775(n): return ((m:=n-1)<<2|1)-(m>>2&-2)+(2,0,-2,0)[m-1>>1&3] # Chai Wah Wu, Feb 02 2025
  • Sage
    a = lambda n: ((((n-1)<< 2)-((n-1)>>2))|1) + ((((n-1)<<1)-((n-1)>> 1)) & 2)
    print([a(n) for n in (1..56)]) # after Andrew Lelechenko, Peter Luschny, Jul 08 2017
    

Formula

A141256(a(n)) = n+1. - Reinhard Zumkeller, Jun 17 2008
From R. J. Mathar, Feb 27 2009: (Start)
a(n+8) = a(n) + 30.
a(n) = a(n-1) + a(n-8) - a(n-9).
G.f.: x*(1 + 6*x + 4*x^2 + 2*x^3 + 4*x^4 + 2*x^5 + 4*x^6 + 6*x^7 + x^8)/((1 + x)*(x^2 + 1)*(x^4 + 1)*(x-1)^2). (End)
a(n) = 4*n - 3 - 2*floor((n-1)/8) + (1 + (-1)^floor((n-2)/2))*(-1)^floor((n-2)/4), n >= 1. - Timothy Hopper, Mar 14 2010
a(1 - n) = -a(n). - Michael Somos, Feb 05 2011
Numbers k such that ((k^2 mod 48=1) or (k^2 mod 48=25)) and ((k^2 mod 120=1) or (k^2 mod 120=49)). - Gary Detlefs, Dec 30 2011
Numbers k such that k^2 mod 30 is 1 or 19. - Gary Detlefs, Dec 31 2011
a(n) = 3*(floor((5*n-2)/4) - floor((n mod 8)/6)) + (n mod 2) + floor(((n-1) mod 8)/7) - floor(((n-2) mod 8)/7) + 1. - Gary Detlefs, Jan 08 2012
a(n) = 4*n - 3 + 2*(floor((n+6)/8) - floor((n+4)/8) - floor((n+2)/8) + floor(n/8) - floor((n-1)/8)), n >= 1. From the o.g.f. given above by R. J. Mathar (with the denominator written as (1-x^8)*(1-x)), and a two-step reduction of the floor functions. Compare with Hopper's and Detlefs's formulas above. - Wolfdieter Lang, Jan 26 2012
a(n) = (6*f(n) - 3 + (-1)^f(n))/2, where f(n)= n + floor(n/4)+ floor(((n+4) mod 8)/6). - Gary Detlefs, Sep 15 2013
a(n) = 30*floor((n-1)/8) + 15 + 2*f((n-1) mod 8 + 16)*(-1)^floor(((n+3) mod 8)/4), where f(n) = (n*(n+1)/2+1) mod 10. - Gary Detlefs, Sep 24 2013
a(n) = 3*n + 6*floor(n/8) + (n mod 2) - 2*floor(((n-2) mod 8)/6) - 2*floor(((n-2) mod 8)/7) + 1. - Gary Detlefs, Jun 01 2014
a(n+1) = ((n << 2 - n >> 2) || 1) + ((n << 1 - n >> 1) && 2), where << and >> are bitwise left and right shifts, || and && are bitwise "or" and "and". - Andrew Lelechenko, Jul 08 2017
a(n) = 2*n + 2*floor(1/2 + (7*n)/8) + 2*(91 mod (2 - ((3*n)/4 + n^2/4) mod 2)) - 3 (n > 0). - Mikk Heidemaa, Dec 06 2017
Sum_{n>=1} (-1)^(n+1)/a(n) = sqrt(23 + sqrt(5) - sqrt(6*(5 + sqrt(5))))*Pi/15. - Amiram Eldar, Dec 13 2021

A008364 11-rough numbers: not divisible by 2, 3, 5 or 7.

Original entry on oeis.org

1, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 113, 121, 127, 131, 137, 139, 143, 149, 151, 157, 163, 167, 169, 173, 179, 181, 187, 191, 193, 197, 199, 209, 211, 221, 223, 227, 229, 233, 239, 241, 247
Offset: 1

Views

Author

Keywords

Comments

The first A005867(4) = 48 terms give the reduced residue system for the 4th primorial number 210 = A002110(4).
This sequence is closed under multiplication: any product of terms is also a term. - Labos Elemer, Feb 26 2003
Conjecture: these are numbers n such that (Sum_{k=1..n} k^4) mod n = 0 and (Sum_{k=1..n} k^6) mod n = 0. - Gary Detlefs, Dec 20 2011
From Peter Bala, May 03 2018: (Start)
The above conjecture is true. Let m be even and let the m-th Bernoulli number be written in reduced form as Bernoulli(m) = N(m)/D(m). Apply Ireland and Rosen, Proposition 15.2.2, to show the congruence D(m)*( Sum_{k = 1..n} k^m )/n = N(m) (mod n) holds for all n >= 1. It follows easily from this congruence that ( Sum_{k = 1..n} k^m )/n is integral iff n is coprime to D(m). Now Bernoulli(4) = -1/(2*3*5) and Bernoulli(6) = 1/(2*3*7) so the numbers n such that both (Sum_{k=1..n} k^4) mod n = 0 and (Sum_{k=1..n} k^6) mod n = 0 are exactly those numbers coprime to the primes 2, 3, 5 and 7, that is, the 11-rough numbers. (End)
Conjecture: these are numbers n such that (n^6 mod 210 = 1) or (n^6 mod 210 = 169). - Gary Detlefs, Dec 30 2011
The second Detlefs conjecture above is true and extremely easy to verify with some basic properties of congruences: take the terms of this sequence up to 209 and compute their sixth powers modulo 210: there should only be 1's and 169's there. Then take the complement of this sequence up to 210, where you will see no instances of 1 or 169. - Alonso del Arte, Jan 12 2014
It is well-known that the product of 7 consecutive integers is divisible by 7!. Conjecture: This sequence is exactly the set of positive values of r such that ( Product_{k = 0..6} n + k*r )/7! is an integer for all n. - Peter Bala, Nov 14 2015
From Ruediger Jehn, Nov 05 2020: (Start)
This conjecture is true. The first part of the proof deals with numbers not in A008364, i.e., numbers which are divisible by p (p either 2, 3, 5, 7). Let r = p*s and n = 1, then (Product_{k = 0..6} n + k*r) is not divisible by p, because none of the factors 1 + k*p*s are divisible by p. Hence dividing the product by 7! does not return an integer.
The second part deals with numbers in A008364. If r and q are coprime, then for any i < q there exists k < q with (k*r mod q) = i. From this, it also follows that for any n there exists k < q with ((n + k*r) mod q) = 0. But this means that Product_{k = 0..6} n + k*r is divisible by all numbers from 2 to 7 because there is always a factor that is divisible. We still have to show that the product is also divisible by 2 times 3 times 4 times 6. If the k_1 with ((n + k_1*r) mod 4) = 0 is even, then (n mod 2) = ((n + 2*r) mod 2) = ((n + 4*r) mod 2) = ((n + 6*r) mod 2) = 0. If this k_1 is odd, then ((n + r) mod 2) = ((n + 3*r) mod 2) = ((n + 5*r) mod 2) = 0. In both cases there are at least 2 other factors divisible by 2. If the k_2 with ((n + k_2*r) mod 6) = 0 is smaller than 4, then ((n + (k_2 + 3)*r) mod 3) = 0. Otherwise, ((n + (k_2 - 3)*r) mod 3) = 0. In both cases there is at least 1 other factor divisible by 3. And therefore Product_{k = 0..6} n + k*r is divisible by 7! for any n.
(End)

References

  • Diatomic sequence of 4th prime: A. de Polignac (1849), J. Dechamps (1907).
  • Dickson L. E., History of the Theory of Numbers, Vol. 1, p. 439, Chelsea, 1952.
  • K. Ireland and M. Rosen, A Classical Introduction to Modern Number Theory, Springer-Verlag, 1980.

Crossrefs

First differences give A049296. Cf. A002110, A048597.
For k-rough numbers with other values of k, see A000027, A005408, A007310, A007775, A008364, A008365, A008366, A166061, A166063. - Michael B. Porter, Oct 10 2009
Cf. A005867, A092695, A210679, A080672 (complement).

Programs

  • Haskell
    a008364 n = a008364_list !! (n-1)
    a008364_list = 1 : filter ((> 7) . a020639) [1..]
    -- Reinhard Zumkeller, Mar 26 2012
  • Maple
    for i from 1 to 500 do if gcd(i,210) = 1 then print(i); fi; od;
    t1:=[]; for i from 1 to 1000 do if gcd(i,210) = 1 then t1:=[op(t1),i]; fi; od: t1;
    S:= (j,n)-> sum(k^j,k=1..n): for n from 1 to 247 do if (S(4,n) mod n = 0) and (S(6,n) mod n = 0) then print(n) fi od; # Gary Detlefs, Dec 20 2011
  • Mathematica
    Select[ Range[ 300 ], GCD[ #1, 210 ] == 1 & ]
    Select[Range[250], Mod[#, 2]>0 && Mod[#, 3]>0 && Mod[#, 5]>0 && Mod[#, 7]>0 &] (* Vincenzo Librandi, Nov 16 2015 *)
    Cases[Range@1000, x_ /; NoneTrue[Array[Prime, 4], Divisible[x, #] &]] (* Mikk Heidemaa, Dec 07 2017 *)
    Select[Range[250],Union[Divisible[#,{2,3,5,7}]]=={False}&] (* Harvey P. Dale, Sep 24 2021 *)
  • PARI
    isA008364(n) = gcd(n,210)==1 \\ Michael B. Porter, Oct 10 2009
    

Formula

Starting with a(49) = 211, a(n) = a(n-48) + 210. - Zak Seidov, Apr 11 2011
a(n) = a(n-1) + a(n-48) - a(n-49). - Charles R Greathouse IV, Dec 21 2011
A020639(a(n)) > 7. - Reinhard Zumkeller, Mar 26 2012
G.f.: x*(x^48 + 10*x^47 + 2*x^46 + 4*x^45 + 2*x^44 + 4*x^43 + 6*x^42 + 2*x^41 + 6*x^40 + 4*x^39 + 2*x^38 + 4*x^37 + 6*x^36 + 6*x^35 + 2*x^34 + 6*x^33 + 4*x^32 + 2*x^31 + 6*x^30 + 4*x^29 + 6*x^28 + 8*x^27 + 4*x^26 + 2*x^25 + 4*x^24 + 2*x^23 + 4*x^22 + 8*x^21 + 6*x^20 + 4*x^19 + 6*x^18 + 2*x^17 + 4*x^16 + 6*x^15 + 2*x^14 + 6*x^13 + 6*x^12 + 4*x^11 + 2*x^10 + 4*x^9 + 6*x^8 + 2*x^7 + 6*x^6 + 4*x^5 + 2*x^4 + 4*x^3 + 2*x^2 + 10*x + 1) / (x^49 - x^48 - x + 1). - Colin Barker, Sep 27 2013
a(n) = 35*n/8 + O(1). - Charles R Greathouse IV, Sep 14 2015
A007775 INTERSECT A206547. - R. J. Mathar, Apr 10 2024

Extensions

New name from Charles R Greathouse IV, Dec 21 2011 based on comment from Michael B. Porter, Oct 10 2009

A297845 Encoded multiplication table for polynomials in one indeterminate with nonnegative integer coefficients. Symmetric square array T(n, k) read by antidiagonals, n > 0 and k > 0. See comment for details.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 4, 5, 4, 1, 1, 5, 9, 9, 5, 1, 1, 6, 7, 16, 7, 6, 1, 1, 7, 15, 25, 25, 15, 7, 1, 1, 8, 11, 36, 11, 36, 11, 8, 1, 1, 9, 27, 49, 35, 35, 49, 27, 9, 1, 1, 10, 25, 64, 13, 90, 13, 64, 25, 10, 1, 1, 11, 21, 81, 125, 77, 77, 125, 81
Offset: 1

Views

Author

Rémy Sigrist, Jan 10 2018

Keywords

Comments

For any number n > 0, let f(n) be the polynomial in a single indeterminate x where the coefficient of x^e is the prime(1+e)-adic valuation of n (where prime(k) denotes the k-th prime); f establishes a bijection between the positive numbers and the polynomials in a single indeterminate x with nonnegative integer coefficients; let g be the inverse of f; T(n, k) = g(f(n) * f(k)).
This table has many similarities with A248601.
For any n > 0 and m > 0, f(n * m) = f(n) + f(m).
Also, f(1) = 0 and f(2) = 1.
The function f can be naturally extended to the set of positive rational numbers: if r = u/v (not necessarily in reduced form), then f(r) = f(u) - f(v); as such, f is a homomorphism from the multiplicative group of positive rational numbers to the additive group of polynomials of a single indeterminate x with integer coefficients.
See A297473 for the main diagonal of T.
As a binary operation, T(.,.) is related to A306697(.,.) and A329329(.,.). When their operands are terms of A050376 (sometimes called Fermi-Dirac primes) the three operations give the same result. However the rest of the multiplication table for T(.,.) can be derived from these results because T(.,.) distributes over integer multiplication (A003991), whereas for A306697 and A329329, the equivalent derivation uses distribution over A059896(.,.) and A059897(.,.) respectively. - Peter Munn, Mar 25 2020
From Peter Munn, Jun 16 2021: (Start)
The operation defined by this sequence can be extended to be the multiplicative operator of a ring over the positive rationals that is isomorphic to the polynomial ring Z[x]. The extended function f (described in the author's original comments) is the isomorphism we use, and it has the same relationship with the extended operation that exists between their unextended equivalents.
Denoting this extension of T(.,.) as t_Q(.,.), we get t_Q(n, 1/k) = t_Q(1/n, k) = 1/T(n, k) and t_Q(1/n, 1/k) = T(n, k) for positive integers n and k. The result for other rationals is derived from the distributive property: t_Q(q, r*s) = t_Q(q, r) * t_Q(q, s); t_Q(q*r, s) = t_Q(q, s) * t_Q(r, s). This may look unusual because standard multiplication of rational numbers takes on the role of the ring's additive group.
There are many OEIS sequences that can be shown to be a list of the integers in an ideal of this ring. See the cross-references.
There are some completely additive sequences that similarly define by extension completely additive functions on the positive rationals that can be shown to be homomorphisms from this ring onto the integer ring Z, and these functions relate to some of the ideals. For example, the extended function of A048675, denoted A048675_Q, maps i/j to A048675(i) - A048675(j) for positive integers i and j. For any positive integer k, the set {r rational > 0 : k divides A048675_Q(r)} forms an ideal of the ring; for k=2 and k=3 the integers in this ideal are listed in A003159 and A332820 respectively.
(End)

Examples

			Array T(n, k) begins:
  n\k|  1   2   3    4    5    6    7     8    9    10
  ---+------------------------------------------------
    1|  1   1   1    1    1    1    1     1    1     1  -> A000012
    2|  1   2   3    4    5    6    7     8    9    10  -> A000027
    3|  1   3   5    9    7   15   11    27   25    21  -> A003961
    4|  1   4   9   16   25   36   49    64   81   100  -> A000290
    5|  1   5   7   25   11   35   13   125   49    55  -> A357852
    6|  1   6  15   36   35   90   77   216  225   210  -> A191002
    7|  1   7  11   49   13   77   17   343  121    91
    8|  1   8  27   64  125  216  343   512  729  1000  -> A000578
    9|  1   9  25   81   49  225  121   729  625   441
   10|  1  10  21  100   55  210   91  1000  441   550
From _Peter Munn_, Jun 24 2021: (Start)
The encoding, n, of polynomials, f(n), that is used for the table is further described in A206284. Examples of encoded polynomials:
   n      f(n)        n           f(n)
   1         0       16              4
   2         1       17            x^6
   3         x       21        x^3 + x
   4         2       25           2x^2
   5       x^2       27             3x
   6     x + 1       35      x^3 + x^2
   7       x^3       36         2x + 2
   8         3       49           2x^3
   9        2x       55      x^4 + x^2
  10   x^2 + 1       64              6
  11       x^4       77      x^4 + x^3
  12     x + 2       81             4x
  13       x^5       90   x^2 + 2x + 1
  15   x^2 + x       91      x^5 + x^3
(End)
		

Crossrefs

Row n: n=1: A000012, n=2: A000027, n=3: A003961, n=4: A000290, n=5: A357852, n=6: A191002, n=8: A000578.
Main diagonal: A297473.
Functions f satisfying f(T(n,k)) = f(n) * f(k): A001222, A048675 (and similarly, other rows of A104244), A195017.
Powers of k: k=3: A000040, k=4: A001146, k=5: A031368, k=6: A007188 (see also A066117), k=7: A031377, k=8: A023365, k=9: main diagonal of A329050.
Integers in the ideal of the related ring (see Jun 2021 comment) generated by S: S={3}: A005408, S={4}: A000290\{0}, S={4,3}: A003159, S={5}: A007310, S={5,4}: A339690, S={6}: A325698, S={6,4}: A028260, S={7}: A007775, S={8}: A000578\{0}, S={8,3}: A191257, S={8,6}: A332820, S={9}: A016754, S={10,4}: A340784, S={11}: A008364, S={12,8}: A145784, S={13}: A008365, S={15,4}: A345452, S={15,9}: A046337, S={16}: A000583\{0}, S={17}: A008366.
Equivalent sequence for polynomial composition: A326376.
Related binary operations: A003991, A306697/A059896, A329329/A059897.

Programs

  • PARI
    T(n,k) = my (f=factor(n), p=apply(primepi, f[, 1]~), g=factor(k), q=apply(primepi, g[, 1]~)); prod (i=1, #p, prod(j=1, #q, prime(p[i]+q[j]-1)^(f[i, 2]*g[j, 2])))

Formula

T is completely multiplicative in both parameters:
- for any n > 0
- and k > 0 with prime factorization Prod_{i > 0} prime(i)^e_i:
- T(prime(n), k) = T(k, prime(n)) = Prod_{i > 0} prime(n + i - 1)^e_i.
For any m > 0, n > 0 and k > 0:
- T(n, k) = T(k, n) (T is commutative),
- T(m, T(n, k)) = T(T(m, n), k) (T is associative),
- T(n, 1) = 1 (1 is an absorbing element for T),
- T(n, 2) = n (2 is an identity element for T),
- T(n, 2^i) = n^i for any i >= 0,
- T(n, 4) = n^2 (A000290),
- T(n, 8) = n^3 (A000578),
- T(n, 3) = A003961(n),
- T(n, 3^i) = A003961(n)^i for any i >= 0,
- T(n, 6) = A191002(n),
- A001221(T(n, k)) <= A001221(n) * A001221(k),
- A001222(T(n, k)) = A001222(n) * A001222(k),
- A055396(T(n, k)) = A055396(n) + A055396(k) - 1 when n > 1 and k > 1,
- A061395(T(n, k)) = A061395(n) + A061395(k) - 1 when n > 1 and k > 1,
- T(A000040(n), A000040(k)) = A000040(n + k - 1),
- T(A000040(n)^i, A000040(k)^j) = A000040(n + k - 1)^(i * j) for any i >= 0 and j >= 0.
From Peter Munn, Mar 13 2020 and Apr 20 2021: (Start)
T(A329050(i_1, j_1), A329050(i_2, j_2)) = A329050(i_1+i_2, j_1+j_2).
T(n, m*k) = T(n, m) * T(n, k); T(n*m, k) = T(n, k) * T(m, k) (T distributes over multiplication).
A104244(m, T(n, k)) = A104244(m, n) * A104244(m, k).
For example, for m = 2, the above formula is equivalent to A048675(T(n, k)) = A048675(n) * A048675(k).
A195017(T(n, k)) = A195017(n) * A195017(k).
A248663(T(n, k)) = A048720(A248663(n), A248663(k)).
T(n, k) = A306697(n, k) if and only if T(n, k) = A329329(n, k).
A007913(T(n, k)) = A007913(T(A007913(n), A007913(k))) = A007913(A329329(n, k)).
(End)

Extensions

New name from Peter Munn, Jul 17 2021

A084970 Numbers whose smallest prime factor is 13.

Original entry on oeis.org

13, 169, 221, 247, 299, 377, 403, 481, 533, 559, 611, 689, 767, 793, 871, 923, 949, 1027, 1079, 1157, 1261, 1313, 1339, 1391, 1417, 1469, 1651, 1703, 1781, 1807, 1937, 1963, 2041, 2119, 2171, 2197, 2249, 2327, 2353, 2483, 2509, 2561, 2587, 2743, 2873
Offset: 1

Views

Author

Robert G. Wilson v, Jun 15 2003

Keywords

Examples

			a(2) = 13*13, a(3) = 13*17.
		

Crossrefs

Sixth row of A083140.
Cf. A084967 (5), A084968 (7), A084969 (11), A332799 (17), A332798 (19), A332797 (23), A008365 (13-rough numbers).

Programs

Formula

a(n) = a(n-480) + 30030 = a(n-1) + a(n-480) - a(n-481). - Charles R Greathouse IV, Nov 19 2014
Lim_{n->infinity} a(n)/n = A038111(6)/A038110(6) = 1001/16 = 62.5625. - Vladimir Shevelev, Jan 20 2015
a(n) = 13*A008365(n).

Extensions

More terms from David Wasserman, Oct 19 2004

A008366 Smallest prime factor is >= 17.

Original entry on oeis.org

1, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 113, 127, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233, 239, 241, 251, 257, 263, 269, 271, 277, 281, 283
Offset: 1

Views

Author

Keywords

Comments

Also the 17-rough numbers: positive integers that have no prime factors less than 17. - Michael B. Porter, Oct 10 2009
a(n) - (1001/192) n is periodic with period 5760. - Robert Israel, Mar 18 2016
From Peter Bala, May 12 2018: (Start)
The product of two 17-rough numbers is a 17-rough number and the prime factors of a 17-rough number are 17-rough numbers.
Let k equal either 13, 14, 15 or 16. Then the product of k numbers n*(n + a)*(n + 2*a)*...*(n + (k-1)*a) in arithmetical progression is divisible by k! for all integer n if and only if a is a 17-rough number.
The sequence terms satisfy the congruence x^60 = 1 (mod 30030), where 30030 = 2*3*5*7*11*13. (End)
The asymptotic density of this sequence is 192/1001. - Amiram Eldar, Sep 30 2020

Crossrefs

For k-rough numbers with other values of k, see A000027 A005408 A007310 A007775 A008364 A008365 A008366 A166061 A166063.
Cf. A005867.

Programs

  • Maple
    for i from 1 to 500 do if gcd(i,30030) = 1 then print(i); fi; od;
  • Mathematica
    Select[ Range[ 300 ], GCD[ #1, 30030 ]==1& ]
    Join[{1},Select[Range[300],FactorInteger[#][[1,1]]>=17&]] (* Harvey P. Dale, Mar 28 2020 *)
  • PARI
    isA008366(n) = gcd(n,30030)==1 \\ Michael B. Porter, Oct 10 2009

Formula

Numbers n > 1 such that ((Sum_{k=1..n} k^10) mod n = 0) and ((Sum_{k=1..n} k^12) mod n = 0) (conjecture). - Gary Detlefs, Dec 27 2011
a(n) = a(n-1) + a(n-5760) - a(n-5761). - Vaclav Kotesovec, Mar 18 2016
G.f: x*P(x)/(1 - x - x^5760 + x^5761) where P(x) is a polynomial of degree 5760. - Benedict W. J. Irwin, Mar 23 2016
a(n) = (1001/192)*n + O(1), where the O(1) term is bounded by +/- 19. - Charles R Greathouse IV, Oct 13 2022
A008365 SETMINUS A084970 . - R. J. Mathar, Nov 05 2024

A166061 19-rough numbers: positive integers that have no prime factors less than 19.

Original entry on oeis.org

1, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 113, 127, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233, 239, 241, 251, 257, 263, 269, 271
Offset: 1

Views

Author

Michael B. Porter, Oct 05 2009

Keywords

Comments

Or, positive integers relatively prime to 510510 = 2*3*5*7*11*13*17.

Examples

			437 = 19 * 23 is in the sequence since the two prime factors, 19 and 23, are not less than 19.
		

Crossrefs

Programs

Formula

a(n) = k*n + O(1) where k = 17017/3072 = 5.539388.... In particular, k*n - 31 < a(n) < k*n + 25. - Charles R Greathouse IV, Sep 24 2018

A166063 23-rough numbers: positive integers that have no prime factors less than 23.

Original entry on oeis.org

1, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 113, 127, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233, 239, 241, 251, 257, 263, 269, 271, 277, 281, 283, 293, 307, 311, 313, 317, 331, 337, 347, 349, 353, 359, 367, 373, 379, 383, 389, 397, 401, 409, 419, 421, 431, 433, 439, 443, 449, 457, 461, 463, 467, 479, 487, 491, 499
Offset: 1

Views

Author

Michael B. Porter, Oct 05 2009

Keywords

Comments

Or, positive integers relatively prime to 9699690 = 2*3*5*7*11*13*17*19.
First composite term is 529 = 23^2.

Examples

			667 = 23 * 29 is in the sequence since the two prime factors, 23 and 29, are not less than 23.
		

Crossrefs

Cf. A332797 (subsequence).

Programs

  • Maple
    A166063 := proc(n)
        option remember;
        local a;
        if n =1 then
            1;
        else
            for a from procname(n-1)+1 do
                numtheory[factorset](a) ;
                if min(op(%)) >= 23 then
                    return a;
                end if;
            end do:
        end if;
    end proc:
    seq(A166063(n),n=1..80) ; # R. J. Mathar, Nov 05 2024
  • Mathematica
    Select[Range[500],FactorInteger[#][[1,1]]>22&] (* Harvey P. Dale, Nov 22 2010 *)
  • PARI
    isA166063(n) = gcd(n,9699690)==1 \\ Michael B. Porter, Oct 10 2009

Formula

a(n) = k*n + O(1) where k = 323323/55296 = 5.8471.... In particular, k*n - 51 < a(n) < k*n + 45. - Charles R Greathouse IV, Sep 21 2018
A166061 SETMINUS A332798 - R. J. Mathar, Nov 05 2024

Extensions

Additional terms provided provided by Harvey P. Dale, Nov 22 2010

A322274 Smallest multiplication factors f, prime or 1, for all b (mod 9240), coprime to 9240 (= 4*11#), so that b*f is a square mod 8, mod 3, mod 5, mod 7, and mod 11.

Original entry on oeis.org

1, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 113, 127, 131, 113, 19, 29, 79, 157, 67, 167, 1, 173, 179, 181, 71, 193, 197, 31, 211, 389, 103, 83, 181, 233, 239, 241, 463, 59, 257, 263, 269, 271, 277, 281, 283, 1, 173, 131, 283, 311, 97, 53, 443, 331, 193, 107, 61, 257, 239, 1, 103, 277
Offset: 1

Views

Author

Hans Ruegg, Dec 01 2018

Keywords

Comments

See sequence A322269 for further explanations. This sequence is related to A322269(5).
The sequence is periodic, repeating itself after phi(9240) terms. Its largest term is 1873, which is A322269(5). In order to satisfy the conditions, both f and b must be coprime to 9240. Otherwise, the product would be zero mod a prime <= 11.
The b(n) corresponding to each a(n) is A008365(n).
The first 28 entries are trivial: f=b, and then the product b*f naturally is a square modulo everything.

Examples

			The 30th number coprime to 9240 is 139. a(30) is 19, because 19 is the smallest prime by which we can multiply 139, so that the product (139*19 = 2641) is a square mod 8, and modulo all primes up to 11.
		

Crossrefs

Programs

  • PARI
    QresCode(n, nPrimes) = {
      code = bitand(n,7)>>1;
      for (j=2, nPrimes,
        x = Mod(n,prime(j));
        if (issquare(x), code += (1<A322271, sequence(3) returns A322272, ... sequence(6) returns A322275.

A334391 Numbers whose only palindromic divisor is 1.

Original entry on oeis.org

1, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 103, 107, 109, 113, 127, 137, 139, 149, 157, 163, 167, 169, 173, 179, 193, 197, 199, 211, 221, 223, 227, 229, 233, 239, 241, 247, 251, 257, 263, 269, 271, 277, 281, 283, 289, 293, 299, 307
Offset: 1

Views

Author

Bernard Schott, Apr 26 2020

Keywords

Comments

Equivalent: Numbers such that the LCM of their palindromic divisors (A087999) is 1, or,
Numbers such that the number of palindromic divisors (A087990) is 1.
All terms are odd.
The 1st family consists of non-palindromic primes that form the subsequence A334321.
The 2nd family consists of {p^k, p prime, k >= 2} such that p^j for 1 <= j <= k is not a palindrome {169 = 13^2, 289 = 17^2, 361 = 19^2, ..., 2197 = 13^3, ...} (see examples).
The 3rd family consists of products p_1^q_1 * ... * p_k^q_k with k >= 2, all of whose divisors are nonpalindromic {221 = 13 * 27, 247 = 13 * 19, 299 = 13 * 23, 377 = 13 * 29, 391 = 17 * 23, 403 = 13 * 31, 481 = 13 * 37, ...}.
Also, equivalent: numbers all of whose divisors > 1 are nonpalindromic (A029742). - Bernard Schott, Jul 14 2022

Examples

			49 = 7^2, the divisor 7 is a palindrome so 49 is not a term.
169 = 13^2, divisors of 169 are {1, 13, 169} and 169 is a term.
391 = 17*23, divisors of 391 are {1,17,23,391} and 391 is a term.
307^2 = 94249 that is palindrome, so 94249 is not a term.
		

Crossrefs

A334321 is a subsequence.

Programs

  • Maple
    notpali:= proc(n) local L;
      L:= convert(n,base,10);
      L <> ListTools:-Reverse(L)
    end proc:
    filter:= proc(n) option remember; andmap(notpali,numtheory:-divisors(n) minus {1}) end proc:
    select(filter, [seq(i,i=1..400,2)]); # Robert Israel, Apr 28 2020
  • Mathematica
    Select[Range[300], !AnyTrue[Rest @ Divisors[#], PalindromeQ] &] (* Amiram Eldar, Apr 26 2020 *)
  • PARI
    ispal(n) = my(d=digits(n)); d == Vecrev(d);
    isok(n) = fordiv(n, d, if (d>1 && ispal(d), return(0))); return(1); \\ Michel Marcus, Apr 26 2020
    
  • Python
    from sympy.ntheory import divisors, is_palindromic
    def ok(n): return not any(is_palindromic(d) for d in divisors(n)[1:])
    print(list(filter(ok, range(1, 308, 2)))) # Michael S. Branicky, May 08 2021

Formula

A087990(a(n)) = 1.
A087999(a(n)) = 1.
Showing 1-10 of 12 results. Next