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 41-50 of 355 results. Next

A001047 a(n) = 3^n - 2^n.

Original entry on oeis.org

0, 1, 5, 19, 65, 211, 665, 2059, 6305, 19171, 58025, 175099, 527345, 1586131, 4766585, 14316139, 42981185, 129009091, 387158345, 1161737179, 3485735825, 10458256051, 31376865305, 94134790219, 282412759265, 847255055011, 2541798719465, 7625463267259, 22876524019505
Offset: 0

Views

Author

Keywords

Comments

a(n+1) is the sum of the elements in the n-th row of triangle pertaining to A036561. - Amarnath Murthy, Jan 02 2002
Number of 2 X n binary arrays with a path of adjacent 1's and no path of adjacent 0's from top row to bottom row. - R. H. Hardin, Mar 21 2002
With offset 1, partial sums of A027649. - Paul Barry, Jun 24 2003
Number of distinct lines through the origin in the n-dimensional lattice of side length 2. A049691 has the values for the 2-dimensional lattice of side length n. - Joshua Zucker, Nov 19 2003
a(n+1)/(n+1)=(3*3^n-2*2^n)/(n+1) is the second binomial transform of the harmonic sequence 1/(n+1). - Paul Barry, Apr 19 2005
a(n+1) is the sum of n-th row of A036561. - Reinhard Zumkeller, May 14 2006
The sequence gives the sum of the lengths of the segments in Cantor's dust generating sequence up to the i-th step. Measurement unit = length of the segment of i-th step. - Giorgio Balzarotti, Nov 18 2006
Let T be a binary relation on the power set P(A) of a set A having n = |A| elements such that for every element x, y of P(A), xTy if x is a proper subset of y. Then a(n) = |T|. - Ross La Haye, Dec 22 2006
From Alexander Adamchuk, Jan 04 2007: (Start)
a(n) is prime for n in A057468.
p divides a(p) - 1 for prime p.
Quotients (3^p - 2^p - 1)/p, where p = prime(n), are listed in A127071.
Numbers k such that k divides 3^k - 2^k - 1 are listed in A127072.
Pseudoprimes in A127072(n) include all powers of primes {2,3,7} and some composite numbers that are listed in A127073, which includes all Carmichael numbers A002997.
Numbers n such that n^2 divides 3^n - 2^n - 1 are listed in A127074.
5 divides a(2n).
5^2 divides a(2*5n).
5^3 divides a(2*5^2n).
5^4 divides a(2*5^3n).
7^2 divides a(6*7n).
13 divides a(4n).
13^2 divides a(4*13n).
19 divides a(3n).
19^2 divides a(3*19n).
23^2 divides a(11n).
23^3 divides a(11*23n).
23^4 divides a(11*23^2n).
29 divides a(7n).
p divides a((p-1)n) for prime p>3.
p divides a((p-1)/2) for prime p in A097934. Also primes p such that 6 is a square mod p, except {2,3}, A038876(n).
p^(k+1) divides a(p^k*(p-1)/2*n) for prime p in A097934.
p^(k+1) divides a(p^k*(p-1)*n) for prime p>3.
Note the exception that for p = 23, p^(k+2) divides a(p^k*(p-1)/2*n).
There are no more such exceptions for primes p up to 600000. (End)
a(n) divides a(q*(n+1)-1), for all q integer. Leonardo Sarasua, Apr 15 2024
Final digits of terms follow sequence 1,5,9,5. - Enoch Haga, Nov 26 2007
This is also the second column sequence of the Sheffer triangle A143494 (2-restricted Stirling2 numbers). See the e.g.f. given below. - Wolfdieter Lang, Oct 08 2011
Partial sums give A000392. - Jon Perry, Apr 05 2014
For n >= 1, this is also row 2 of A281890: when consecutive positive integers are written as a product of primes in nondecreasing order, "3" occurs in n-th position a(n) times out of every 6^n. - Peter Munn, May 17 2017
a(n) is the number of ternary sequences of length n which include the digit 2. For example, a(2)=5 since the sequences are 02,20,12,21,22. - Enrique Navarrete, Apr 05 2021
a(n-1) is the number of ways we can form disjoint unions of two nonempty subsets of [n] such that the union contains n. For example, for n = 3, a(2) = 5 since the disjoint unions are {1}U{3}, {1}U{2,3}, {2}U{3}, {2}U{1,3}, and {1,2}U{3}. Cf. A000392 if we drop the requirement that the union contains n. - Enrique Navarrete, Aug 24 2021
Configures as a composite Koch Snowflake Fractal (see illustration in links) based on the five-fold division of the Cantor Square/Cantor Dust Fractal of (9^n-4^n)/5 see my illustration in (A016153). - John Elias, Oct 13 2021
Number of pairs (A,B) where B is a subset of {1,2,...,n} and A is a proper subset of B. - Jianing Song, Jun 18 2022
From Manfred Boergens, Mar 29 2023: (Start)
With regard to the comments by Ross La Haye and Jianing Song: Omitting "proper" gives A000244.
Number of pairs (A,B) where B is a nonempty subset of {1,2,...,n} and A is a nonempty subset of B. For nonempty proper subsets see a(n+1) in A028243. (End)
a(n) is the number of n-digit numbers whose smallest decimal digit is 7. - Stefano Spezia, Nov 15 2023
a(n-1) is the number of all possible player-reduced binary games observed by each player in an nx2 game assuming the individual strategies of k < n - 1 players are fixed and the remaining n - k - 1 player will play as one, either maintaining their status quo strategies or jointly adopting an alternative strategy. - Ambrosio Valencia-Romero, Apr 11 2024

References

  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 86-87.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

a(n) = row sums of A091913, row 2 of A047969, column 1 of A090888 and column 1 of A038719.
Cf. partitions: A241766, A241759.
A diagonal of A262307.

Programs

  • Haskell
    a001047 n = a001047_list !! n
    a001047_list = map fst $ iterate (\(u, v) -> (3 * u + v, 2 * v)) (0, 1)
    -- Reinhard Zumkeller, Jun 09 2013
  • Magma
    [3^n - 2^n: n in [0..30]]; // Vincenzo Librandi, Jul 17 2011
    
  • Maple
    seq(3^n - 2^n, n=0..40); # Giorgio Balzarotti, Nov 18 2006
    A001047:=1/(3*z-1)/(2*z-1); # Simon Plouffe in his 1992 dissertation, dropping the initial zero
  • Mathematica
    Table[ 3^n - 2^n, {n, 0, 25} ]
    LinearRecurrence[{5, -6}, {0, 1}, 25] (* Harvey P. Dale, Aug 18 2011 *)
    Numerator@NestList[(3#+1)/2&,1/2,100] (* Zak Seidov, Oct 03 2011 *)
  • PARI
    {a(n) = 3^n - 2^n};
    
  • Python
    [3**n - 2**n for n in range(25)] # Ross La Haye, Aug 19 2005; corrected by David Radcliffe, Jun 26 2016
    
  • Sage
    [lucas_number1(n, 5, 6) for n in range(26)]  # Zerinvary Lajos, Apr 22 2009
    

Formula

G.f.: x/((1-2*x)*(1-3*x)).
a(n) = 5*a(n-1) - 6*a(n-2).
a(n) = 3*a(n-1) + 2^(n-1). - Jon Perry, Aug 23 2002
Starting 0, 0, 1, 5, 19, ... this is 3^n/3 - 2^n/2 + 0^n/6, the binomial transform of A086218. - Paul Barry, Aug 18 2003
a(n) = A083323(n)-1 = A056182(n)/2 = (A002783(n)-1)/2 = (A003063(n+2)-A003063(n+1))/2. - Ralf Stephan, Jan 12 2004
Binomial transform of A000225. - Ross La Haye, Feb 07 2005
a(n) = Sum_{k=0..n-1} binomial(n, k)*2^k. - Ross La Haye, Aug 20 2005
a(n) = 2^(2n) - A083324(n). - Ross La Haye, Sep 10 2005
a(n) = A112626(n, 1). - Ross La Haye, Jan 11 2006
E.g.f.: exp(3*x) - exp(2*x). - Mohammad K. Azarian, Jan 14 2009
a(n) = A217764(n,1). - Ross La Haye, Mar 27 2013
a(n) = 2*a(n-1) + 3^(n-1). - Toby Gottfried, Mar 28 2013
a(n) = A000244(n) - A000079(n). - Omar E. Pol, Mar 28 2013
a(n) = Sum_{k=0..2} Stirling1(2,k)*(k+1)^n = c_2^{(-n)}, poly-Cauchy numbers. - Takao Komatsu, Mar 28 2013
a(n) = A227048(n,A098294(n)). - Reinhard Zumkeller, Jun 30 2013
a(n+1) = Sum_{k=0..n} 2^k*3^(n-k). - J. M. Bergot, Mar 27 2018
Sum_{n>=1} 1/a(n) = A329064. - Amiram Eldar, Nov 20 2020
a(n) = (1/2)*Sum_{k=0..n} binomial(n, k)*(2^(n-k) + 2^k - 2).
a(n) = A001117(n) + 2*A000918(n) + 1. - Ambrosio Valencia-Romero, Mar 08 2022
a(n) = A000225(n) + A028243(n+1). - Ambrosio Valencia-Romero, Mar 09 2022
From Peter Bala, Jun 27 2025: (Start)
exp(Sum_{n >=1} a(2*n)/a(n)*x^n/n) = Sum_{n >= 0} a(n+1)*x^n.
exp(Sum_{n >=1} a(3*n)/a(n)*x^n/n) = 1 + 19*x + 247*x^2 + ... is the g.f. of A019443.
exp(Sum_{n >=1} a(4*n)/a(n)*x^n/n) = 1 + 65*x + 2743*x^2 + ... is the g.f. of A383754.
The following are all examples of telescoping series:
Sum_{n >= 1} 6^n/(a(n)*a(n+1)) = 2, since 6^n/(a(n)*a(n+1)) = b(n) - b(n+1), where b(n) = 2^n/a(n);
Sum_{n >= 1} 18^n/(a(n)*a(n+1)*a(n+2)) = 22/75, since 18^n/(a(n)*a(n+1)*a(n+2)) = c(n) - c(n+1), where c(n) = (5*6^n - 2*4^n)/(15*a(n)*a(n+1));
Sum_{n >= 1} 54^n/(a(n)*a(n+1)*a(n+2)*a(n+3)) = 634/48735 since 54^n/(a(n)*a(n+1)*a(n+2)*a(n+3)) = d(n) - d(n+1), where d(n) = (57*18^n - 38*12^n + 8*8^n)/(513*a(n)*a(n+1)*a(n+2)).
Sum_{n >= 1} 6^n/(a(n)*a(n+2)) = 14/25; Sum_{n >= 1} (-6)^n/(a(n)*a(n+2)) = -6/25.
Sum_{n >= 1} 6^n/(a(n)*a(n+3)) = 306/1805.
Sum_{n >= 1} 6^n/(a(n)*a(n+4)) = 4282/80275; Sum_{n >= 1} (-6)^n/(a(n)*a(n+4)) = -1698/80275. (End)

Extensions

Edited by Charles R Greathouse IV, Mar 24 2010

A003277 Cyclic numbers: k such that k and phi(k) are relatively prime; also k such that there is just one group of order k, i.e., A000001(k) = 1.

Original entry on oeis.org

1, 2, 3, 5, 7, 11, 13, 15, 17, 19, 23, 29, 31, 33, 35, 37, 41, 43, 47, 51, 53, 59, 61, 65, 67, 69, 71, 73, 77, 79, 83, 85, 87, 89, 91, 95, 97, 101, 103, 107, 109, 113, 115, 119, 123, 127, 131, 133, 137, 139, 141, 143, 145, 149, 151, 157, 159, 161, 163, 167, 173
Offset: 1

Views

Author

Keywords

Comments

Except for a(2)=2, all the terms in the sequence are odd. This is because of the existence of a non-cyclic dihedral group of order 2n for each n>1. - Ahmed Fares (ahmedfares(AT)my-deja.com), May 09 2001
Also gcd(n, A051953(n)) = 1. - Labos Elemer
n such that x^n == 1 (mod n) has no solution 2 <= x <= n. - Benoit Cloitre, May 10 2002
There is only one group (the cyclic group of order n) whose order is n. - Gerard P. Michon, Jan 08 2008 [This is a 1947 result of Tibor Szele. - Charles R Greathouse IV, Nov 23 2011]
Any divisor of a Carmichael number (A002997) must be odd and cyclic. Conversely, G. P. Michon conjectured (c. 1980) that any odd cyclic number has at least one Carmichael multiple (if the conjecture is true, each of them has infinitely many such multiples). In 2007, Michon & Crump produced explicit Carmichael multiples of all odd cyclic numbers below 10000 (see link, cf. A253595). - Gerard P. Michon, Jan 08 2008
Numbers n such that phi(n)^phi(n) == 1 (mod n). - Michel Lagneau, Nov 18 2012
Contains A000040, and all members of A006094 except 6. - Robert Israel, Jul 08 2015
Number m such that n^n == r (mod m) is solvable for any r. - David W. Wilson, Oct 01 2015
Numbers m such that A074792(m) = m + 1. - Thomas Ordowski, Jul 16 2017
Squarefree terms of A056867 (see McCarthy link p. 592 and similar comment with "cubefree" in A051532). - Bernard Schott, Mar 24 2022

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 840.
  • J. S. Rose, A Course on Group Theory, Camb. Univ. Press, 1978, see p. 7.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Subsequence of A051532. Intersection of A056867 and A005117.
Cf. A000010, A008966, A009195, A050384 (the same sequence but with the primes removed). Also A000001(a(n)) = 1.

Programs

  • Haskell
    import Data.List (elemIndices)
    a003277 n = a003277_list !! (n-1)
    a003277_list = map (+ 1) $ elemIndices 1 a009195_list
    -- Reinhard Zumkeller, Feb 27 2012
    
  • Magma
    [n: n in [1..200] | Gcd(n, EulerPhi(n)) eq 1]; // Vincenzo Librandi, Jul 09 2015
    
  • Maple
    select(t -> igcd(t, numtheory:-phi(t))=1, [$1..1000]); # Robert Israel, Jul 08 2015
  • Mathematica
    Select[Range[175], GCD[#, EulerPhi[#]] == 1 &] (* Jean-François Alcover, Apr 04 2011 *)
    Select[Range@175, FiniteGroupCount@# == 1 &] (* Robert G. Wilson v, Feb 16 2017 *)
    Select[Range[200],CoprimeQ[#,EulerPhi[#]]&] (* Harvey P. Dale, Apr 10 2022 *)
  • PARI
    isA003277(n) = gcd(n,eulerphi(n))==1 \\ Michael B. Porter, Feb 21 2010
    
  • Sage
    # Compare A050384.
    def isPrimeTo(n, m): return gcd(n, m) == 1
    def isCyclic(n): return isPrimeTo(n, euler_phi(n))
    [n for n in (1..173) if isCyclic(n)] # Peter Luschny, Nov 14 2018

Formula

n = p_1*p_2*...*p_k (for some k >= 0), where the p_i are distinct primes and no p_j-1 is divisible by any p_i.
A000001(a(n)) = 1.
Erdős proved that a(n) ~ e^gamma n log log log n, where e^gamma is A073004. - Charles R Greathouse IV, Nov 23 2011
A000005(a(n)) = 2^k. - Carlos Eduardo Olivieri, Jul 07 2015
A008966(a(n)) = 1. - Bernard Schott, Mar 24 2022

Extensions

More terms from Christian G. Bower

A087788 3-Carmichael numbers: Carmichael numbers equal to the product of 3 primes: k = p*q*r, where p < q < r are primes such that a^(k-1) == 1 (mod k) if a is prime to k.

Original entry on oeis.org

561, 1105, 1729, 2465, 2821, 6601, 8911, 10585, 15841, 29341, 46657, 52633, 115921, 162401, 252601, 294409, 314821, 334153, 399001, 410041, 488881, 512461, 530881, 1024651, 1152271, 1193221, 1461241, 1615681, 1857241, 1909001, 2508013
Offset: 1

Views

Author

Miklos Kristof, Oct 07 2003

Keywords

Comments

It is interesting that most of the numbers have the last digit 1. For example 530881, 3581761, 7207201, etc.
Granville & Pomerance conjecture that there are ~ c x^(1/3)/(log x)^3 terms of this sequence up to x. Heath-Brown proves that, for any e > 0, there are O(x^(7/20 + e)) terms of this sequence up to x. - Charles R Greathouse IV, Nov 19 2012

Examples

			a(6)=6601=7*23*41: 7-1|6601-1, 23-1|6601-1, 41-1|6601-1, i.e., 6|6600, 22|6600, 40|6600.
		

References

  • O. Ore, Number Theory and Its History, McGraw-Hill, 1948, Reprinted by Dover Publications, 1988, Chapter 14.

Crossrefs

Intersection of A002997 and A007304.
Cf. A162290.

Programs

  • PARI
    list(lim)=my(v=List());forprime(p=3,(lim)^(1/3), forprime(q=p+1, sqrt(lim\p),forprime(r=q+1,lim\(p*q),if((q*r-1)%(p-1)||(p*r-1)%(q-1)||(p*q-1)%(r-1),,listput(v,p*q*r)))));vecsort(Vec(v)) \\ Charles R Greathouse IV, Nov 19 2012

Formula

k is composite and squarefree and for p prime, p|k => p-1|k-1. A composite odd number k is a Carmichael number if and only if k is squarefree and p-1 divides k-1 for every prime p dividing k (Korselt, 1899) k = p*q*r, p-1|k-1, q-1|k-1, r-1|k-1.

Extensions

Minor edit to definition by N. J. A. Sloane, Sep 14 2009

A208728 Composite numbers n such that b^(n+1) == 1 (mod n) for every b coprime to n.

Original entry on oeis.org

15, 35, 255, 455, 1295, 2703, 4355, 6479, 9215, 10439, 11951, 16211, 23435, 27839, 44099, 47519, 47879, 62567, 63167, 65535, 93023, 94535, 104195, 120959, 131327, 133055, 141155, 142883, 157079, 170819, 196811, 207935, 260831, 283679, 430199, 560735, 576719
Offset: 1

Views

Author

Paolo P. Lava, Mar 01 2012

Keywords

Comments

GCD(b,n)=1 and b^(n+1) == 1 (mod n).
The sequence lists the squarefree composite numbers n such that every prime divisor p of n satisfies (p-1)|(n+1) (similar to Korselt's criterion).
The sequence can be considered as an extension of k-Knödel numbers to k negative, in this case equal to -1.
Numbers n > 3 such that b^(n+2) == b (mod n) for every integer b. Also, numbers n > 3 such that A002322(n) divides n+1. Are there infinitely many such numbers? It seems that such numbers n > 35 have at least three prime factors. - Thomas Ordowski, Jun 25 2017
Proof that 15 and 35 are the only numbers with only two prime factors (and so all others have >= 3): Since n is squarefree and composite, it has at least two prime factors, p and q. If these are the only two, n = p*q. Then the criterion is that (p-1)|(n+1) -> (p-1)|pq+1, and q-1|pq+1. Write pq+1 = j*(p-1) = k*(q-1). Rearranging, p*(j-q)=j+1 and q*(k-p)=k+1. Since j = (pq+1)/(p-1), for large n, j~=q, and k~=p. But we see that p divides j+1~=q, and q divides k+1~=p. For large n this is only possible if p and q are roughly equal, so j-q=k-p=1. Otherwise, we have j+1 >= 2*p and k+1 >= 2*q, and which puts upper bounds on p and q. Enumerating these gives (p,q)=(3,5) and (p,q)=(5,7) as the only solutions. - Alex Meiburg, Oct 03 2024

Examples

			6479 is part of the sequence because its prime factors are 11, 19 and 31: (6479+1)/(11-1)=648, (6479+1)/(19-1)=360 and (6479+1)/(31-1)=216.
		

Crossrefs

Programs

  • Maple
    with(numtheory); P:=proc(n) local d, ok, p;
    if issqrfree(n) then p:=factorset(n); ok:=1;
    for d from 1 to nops(p) do if frac((n+1)/(p[d]-1))>0 then ok:=0;
    break; fi; od; if ok=1 then n; fi; fi; end: seq(P(i),i=5..576719);
  • Mathematica
    Select[Range[2, 576719], SquareFreeQ[#] && ! PrimeQ[#] && Union[Mod[# + 1, Transpose[FactorInteger[#]][[1]] - 1]] == {0} &] (* T. D. Noe, Mar 05 2012 *)
  • PARI
    is(n)=if(isprime(n)||!issquarefree(n)||n<3, return(0)); my(f=factor(n)[, 1]); for(i=1, #f, if((n+1)%(f[i]-1), return(0))); 1 \\ Charles R Greathouse IV, Mar 05 2012

Extensions

Definition corrected by Thomas Ordowski, Jun 25 2017

A027760 Denominator of Sum_{p prime, p-1 divides n} 1/p.

Original entry on oeis.org

2, 6, 2, 30, 2, 42, 2, 30, 2, 66, 2, 2730, 2, 6, 2, 510, 2, 798, 2, 330, 2, 138, 2, 2730, 2, 6, 2, 870, 2, 14322, 2, 510, 2, 6, 2, 1919190, 2, 6, 2, 13530, 2, 1806, 2, 690, 2, 282, 2, 46410, 2, 66, 2, 1590, 2, 798, 2, 870, 2, 354, 2, 56786730, 2, 6, 2, 510, 2
Offset: 1

Views

Author

Keywords

Comments

The GCD of integers x^(n+1)-x, for all integers x. - Roger Cuculiere (cuculier(AT)imaginet.fr), Jan 19 2002
If each x in a ring satisfies x^(n+1)=x, the characteristic of the ring is a divisor of a(n) (Rosenblum 1977). - Daniel M. Rosenblum (DMRosenblum(AT)world.oberlin.edu), Sep 24 2008
The denominators of the Bernoulli numbers for n>0. B_n sequence begins 1, -1/2, 1/6, 0/2, -1/30, 0/2, 1/42, 0/2, ... This is an alternative version of A027642 suggested by the theorem of Clausen. To add a(0) = 1 has been proposed in A141056. - Peter Luschny, Apr 29 2009
For N > 1, a(n) is the greatest number k such that x*y^n ==y*x^n (mod k) for any integers x and y. Example: a(19) = 798 because x*y^19 ==y*x^19 (mod 798). - Michel Lagneau, Apr 21 2012
a(n) is the largest k such that b^(n+1) == b (mod k) for every integer b. - Mateusz Szymański, Feb 18 2016, corrected by Thomas Ordowski, Jul 01 2018
When n is even, a(n) is the product of the distinct primes dividing the denominator of zeta(1-n), where zeta(s) is the Riemann zeta function. - Griffin N. Macris, Jun 13 2016
If n+1 is prime, then A002322(a(n)) = n. Composite numbers n+1 such that A002322(a(n)) = n are in A317210. - Max Alekseyev and Thomas Ordowski, Jul 09 2018

Examples

			1/2, 5/6, 1/2, 31/30, 1/2, 41/42, 1/2, 31/30, 1/2, 61/66, 1/2, 3421/2730, 1/2, 5/6, 1/2, 557/510, ...
		

Crossrefs

Programs

  • Maple
    A027760 := proc(n) local s,p; s := 0 ; p := 2; while p <= n+1 do if n mod (p-1) = 0 then s := s+1/p; fi; p := nextprime(p) ; od: denom(s) ; end: # R. J. Mathar, Aug 12 2008
  • Mathematica
    clausen[n_] := Product[i, {i, Select[ Map[ # + 1 &, Divisors[n]], PrimeQ]}]
    Table[clausen[i], {i, 1, 20}] (* Peter Luschny, Apr 29 2009 *)
    f[n_] := Times @@ Select[Divisors@n + 1, PrimeQ]; Array[f, 56] (* Robert G. Wilson v, Apr 25 2012 *)
  • PARI
    a(n)=denominator(sumdiv(n,d,if(isprime(d+1),1/(d+1)))) \\ Charles R Greathouse IV, Jul 08 2011
    
  • PARI
    a(n)=my(pr=1);fordiv(n,d,if(isprime(d+1),pr*=d+1));pr \\ Charles R Greathouse IV, Jul 08 2011
    
  • Sage
    def A027760(n):
        return mul(filter(lambda s: is_prime(s), map(lambda i: i+1, divisors(n))))
    [A027760(n) for n in (1..56)]  # Peter Luschny, May 23 2013

Formula

a(2*k) = A091137(2*k)/A091137(2*k-1). - Paul Curtz, Aug 05 2008
a(n) = product_{p prime, p-1 divides n}. - Eric M. Schmidt, Aug 01 2013
a(2n-1) = 2. - Robert G. Wilson v, Jul 23 2018

Extensions

Formula submitted with A141417 added by R. J. Mathar, Nov 17 2010

A015919 Positive integers k such that 2^k == 2 (mod k).

Original entry on oeis.org

1, 2, 3, 5, 7, 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, 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, 341, 347, 349, 353, 359, 367
Offset: 1

Views

Author

Keywords

Comments

Includes 341 which is first pseudoprime to base 2 and distinguishes sequence from A008578.
First composite even term is a(14868) = 161038 = A006935(2). - Max Alekseyev, Feb 11 2015
If k is a term, then so is 2^k - 1. - Max Alekseyev, Sep 22 2016
Terms of the form 2^k - 2 correspond to k in A296104. - Max Alekseyev, Dec 04 2017
If 2^k - 1 is a term, then so is k. - Thomas Ordowski, Apr 27 2018

Crossrefs

Contains A002997 as a subsequence.
The odd terms form A176997.

Programs

  • Mathematica
    Prepend[ Select[ Range@370, PowerMod[2, #, #] == 2 &], {1, 2}] // Flatten (* Robert G. Wilson v, May 16 2018 *)
  • PARI
    is(n)=Mod(2,n)^n==2 \\ Charles R Greathouse IV, Mar 11 2014
    
  • Python
    def ok(n): return pow(2, n, n) == 2%n
    print([k for k in range(1, 400) if ok(k)]) # Michael S. Branicky, Jun 03 2022

Formula

Equals {1} U A000040 U A001567 U A006935 = A001567 U A006935 U A008578. - Ray Chandler, Dec 07 2003; corrected by Max Alekseyev, Feb 11 2015

A006931 Least Carmichael number with n prime factors, or 0 if no such number exists.

Original entry on oeis.org

561, 41041, 825265, 321197185, 5394826801, 232250619601, 9746347772161, 1436697831295441, 60977817398996785, 7156857700403137441, 1791562810662585767521, 87674969936234821377601, 6553130926752006031481761, 1590231231043178376951698401
Offset: 3

Views

Author

Keywords

Comments

Alford, Grantham, Hayman, & Shallue construct large Carmichael numbers, finding upper bounds for a(3)-a(19565220) and a(10333229505). - Charles R Greathouse IV, May 30 2013

References

  • J.-P. Delahaye, Merveilleux nombres premiers ("Amazing primes"), p. 269, Pour la Science, Paris 2000.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Mathematica
    (* Program not suitable to compute more than a few terms *)
    A2997 = Select[Range[1, 10^6, 2], CompositeQ[#] && Mod[#, CarmichaelLambda[#] ] == 1&];
    (First /@ Split[Sort[{PrimeOmega[#], #}& /@ A2997], #1[[1]] == #2[[1]]&])[[All, 2]] (* Jean-François Alcover, Sep 11 2018 *)
  • PARI
    Korselt(n,f=factor(n))=for(i=1,#f[,1],if(f[i,2]>1||(n-1)%(f[i,1]-1),return(0)));1
    a(n)=my(p=2,f);forprime(q=3,default(primelimit),forstep(k=p+2,q-2,2,f=factor(k);if(vecmax(f[,2])==1 && #f[,2]==n && Korselt(k,f), return(k)));p=q)
    \\ Charles R Greathouse IV, Apr 25 2012
    
  • PARI
    carmichael(A, B, k) = A=max(A, vecprod(primes(k+1))\2); (f(m, l, lo, k) = my(list=List()); my(hi=sqrtnint(B\m, k)); if(lo > hi, return(list)); if(k==1, lo=max(lo, ceil(A/m)); my(t=lift(1/Mod(m,l))); while(t < lo, t += l); forstep(p=t, hi, l, if(isprime(p), my(n=m*p); if((n-1)%(p-1) == 0, listput(list, n)))), forprime(p=lo, hi, if(gcd(m, p-1) == 1, list=concat(list, f(m*p, lcm(l, p-1), p+1, k-1))))); list); vecsort(Vec(f(1, 1, 3, k)));
    a(n) = if(n < 3, return()); my(x=vecprod(primes(n+1))\2,y=2*x); while(1, my(v=carmichael(x,y,n)); if(#v >= 1, return(v[1])); x=y+1; y=2*x); \\ Daniel Suteu, Feb 24 2023

Extensions

Corrected by Lekraj Beedassy, Dec 31 2002
More terms from Ralf Stephan, from the Pinch paper, Apr 16 2005
Edited by N. J. A. Sloane, May 16 2008 at the suggestion of R. J. Mathar.
Escape clause added by Jianing Song, Dec 12 2021

A063994 a(n) = Product_{primes p dividing n } gcd(p-1, n-1).

Original entry on oeis.org

1, 1, 2, 1, 4, 1, 6, 1, 2, 1, 10, 1, 12, 1, 4, 1, 16, 1, 18, 1, 4, 1, 22, 1, 4, 1, 2, 3, 28, 1, 30, 1, 4, 1, 4, 1, 36, 1, 4, 1, 40, 1, 42, 1, 8, 1, 46, 1, 6, 1, 4, 3, 52, 1, 4, 1, 4, 1, 58, 1, 60, 1, 4, 1, 16, 5, 66, 1, 4, 3, 70, 1, 72, 1, 4, 3, 4, 1, 78, 1, 2, 1, 82, 1, 16, 1, 4, 1, 88, 1, 36, 1
Offset: 1

Views

Author

N. J. A. Sloane, Sep 18 2001

Keywords

Comments

a(n) = number of bases b modulo n for which b^{n-1} == 1 (mod n).
a(A209211(n)) = 1. - Reinhard Zumkeller, Mar 02 2013
A049559(n) divides a(n) divides A000010(n). - Thomas Ordowski, Dec 14 2013
Note that a(n) = phi(n) iff n = 1 or n is prime or n is Carmichael number A002997. - Thomas Ordowski, Dec 17 2013

Crossrefs

Programs

  • Haskell
    a063994 n = product $ map (gcd (n - 1) . subtract 1) $ a027748_row n
    -- Reinhard Zumkeller, Mar 02 2013
  • Mathematica
    f[n_] := Times @@ GCD[n - 1, First /@ FactorInteger@ n - 1]; f[1] = 1; Array[f, 92] (* Robert G. Wilson v, Aug 08 2011 *)
  • PARI
    for (n=1, 1000, f=factor(n)~; a=1; for (i=1, length(f), a*=gcd(f[1, i] - 1, n - 1)); write("b063994.txt", n, " ", a) ) \\ Harry J. Smith, Sep 05 2009
    
  • PARI
    a(n)=my(f=factor(n)[,1]);prod(i=1,#f,gcd(f[i]-1,n-1)) \\ Charles R Greathouse IV, Dec 10 2013
    
  • Python
    def a(n):
        if n == 1: return 1
        return len([1 for witness in range(1,n) if pow(witness, n - 1, n) == 1])
    [a(n) for n in range(1, 100)]
    

Formula

a(p^m) = p-1 and a(2p^m) = 1 for prime p and integer m > 0. - Thomas Ordowski, Dec 15 2013
a(n) = Sum_{k=1..n}(floor((k^(n-1)-1)/n)-floor((k^(n-1)-2)/n)). - Anthony Browne, May 11 2016

Extensions

More terms from Robert G. Wilson v, Sep 21 2001

A324370 Product of all primes p not dividing n such that the sum of the base-p digits of n is at least p, or 1 if no such prime exists.

Original entry on oeis.org

1, 1, 2, 1, 6, 1, 6, 3, 10, 1, 6, 1, 210, 15, 2, 3, 30, 5, 210, 21, 110, 15, 30, 5, 546, 21, 14, 1, 30, 1, 462, 231, 1190, 105, 6, 1, 51870, 1365, 70, 21, 2310, 55, 2310, 105, 322, 105, 210, 35, 6630, 663, 286, 33, 330, 55, 798, 57, 290, 15, 30, 1, 930930, 15015, 1430, 2145, 1122, 85, 82110, 2415, 70, 3, 330, 55, 21111090, 285285
Offset: 1

Views

Author

Keywords

Comments

The product is finite, as the sum of the base-p digits of n is n if p > n.
a(198) = 2465 is the only term below 10^6 that is a Carmichael number (A002997).
It appears that a(n)=1 if and only if n is in A094960. - Robert Israel, Mar 30 2020
It turns out that a(n) equals the denominator of the first derivative of the Bernoulli polynomial B(n,x). So a(n)=1 if and only if n is in A094960, also impyling that n+1 is prime. A324370 is also involved in such formulas regarding higher derivatives. See Kellner 2023. - Bernd C. Kellner, Oct 12 2023

Examples

			For p = 2, 3, and 5, the sum of the base p digits of 7 is 1+1+1 = 3 >= 2, 2+1 = 3 >= 3, and 1+2 = 3 < 5, respectively, so a(7) = 2*3 = 6.
		

Crossrefs

Programs

  • Maple
    N:= 100: # for a(1)..a(N)
    V:= Vector(N,1):
    p:= 1:
    for iter from 1 do
       p:= nextprime(p);
       if p >= N then break fi;
       for n from p+1 to N do
         if n mod p <> 0 and convert(convert(n,base,p),`+`)>= p then
           V[n]:= V[n]*p
         fi
    od od:
    convert(V,list); # Robert Israel, Mar 30 2020
    # Alternatively, note that this formula is suggesting offset 0 and a(0) = 1:
    seq(denom(diff(bernoulli(n, x), x)), n = 1..51); # Peter Luschny, Oct 13 2023
  • Mathematica
    SD[n_, p_] := If[n < 1 || p < 2, 0, Plus @@ IntegerDigits[n, p]];
    DD2[n_] := Times @@ Select[Prime[Range[PrimePi[(n+1)/(2+Mod[n+1, 2])]]], !Divisible[n, #] && SD[n, #] >= # &];
    Table[DD2[n], {n, 1, 100}]
    (* From Bernd C. Kellner, Oct 12 2023 (Start) *)
    (* Denominator of first derivative of BP *)
    k = 1; Table[Denominator[Together[D[BernoulliB[n, x], {x, k}]]], {n, 1, 100}]
    (* End *)
  • Python
    from math import prod
    from sympy.ntheory import digits
    from sympy import primefactors, primerange
    def a(n):
        nonpf = set(primerange(1, n+1)) - set(primefactors(n))
        return prod(p for p in nonpf if sum(digits(n, p)[1:]) >= p)
    print([a(n) for n in range(1, 75)]) # Michael S. Branicky, Jul 03 2022

Formula

a(n) * A324369(n) = A195441(n-1) = denominator(Bernoulli_n(x) - Bernoulli_n).
a(n) * A324369(n) * A324371(n) = A144845(n-1) = denominator(Bernoulli_{n-1}(x)).
a(n+1) = A195441(n)/A324369(n+1) = A144845(n)/A007947(n+1) = A318256(n). Essentially the same as A318256. - Peter Luschny, Mar 05 2019
From Bernd C. Kellner, Oct 12 2023: (Start)
a(n) = denominator(Bernoulli_n(x)').
k-th derivative: let (n)_m be the falling factorial.
For n > k, a(n-k+1)/gcd(a(n-k+1), (n)_{k-1}) = denominator(Bernoulli_n(x)^(k)). Otherwise, the denominator equals 1. (End)

A121707 Numbers n > 1 such that n^3 divides Sum_{k=1..n-1} k^n = A121706(n).

Original entry on oeis.org

35, 55, 77, 95, 115, 119, 143, 155, 161, 187, 203, 209, 215, 221, 235, 247, 253, 275, 287, 295, 299, 319, 323, 329, 335, 355, 371, 377, 391, 395, 403, 407, 413, 415, 437, 455, 473, 475, 493, 497, 515, 517, 527, 533, 535, 539, 551, 559, 575, 581, 583, 589, 611
Offset: 1

Views

Author

Alexander Adamchuk, Aug 16 2006

Keywords

Comments

All terms belong to A038509 (Composite numbers with smallest prime factor >= 5). Many but not all terms belong to A060976 (Odd nonprimes, c, which divide Bernoulli(2*c)).
Many terms are semiprimes:
- the non-semiprimes are {275, 455, 475, 539, 575, 715, 775, 875, 935, ...}: see A321487;
- semiprime terms that are multiples of 5 have indices {7, 11, 19, 23, 31, 43, 47, 59, 67, 71, 79, 83, ...} = A002145 (Primes of form 4*k + 3, except 3, or k > 0; or Primes which are also Gaussian primes);
- semiprime terms that are multiples of 7 have indices {5, 11, 17, 23, 29, 41, 47, 53, 59, 71, 83, 89, ...} = A003627 (Primes of form 3*k - 1, except 2, or k > 1);
- semiprime terms that are multiples of 11 have indices {5, 7, 13, 17, 19, 23, 37, 43, 47, 53, 59, 67, 73, 79, 83, ...} = Primes of the form 4*k + 1 and 4*k - 1. [Edited by Michel Marcus, Jul 21 2018, M. F. Hasler, Nov 09 2018]
Conjecture: odd numbers n > 1 such that n divides Sum_{k=1..n-1} k^(n-1). - Thomas Ordowski and Robert Israel, Oct 09 2015. Professor Andrzej Schinzel (in a letter to me, dated Dec 29 2015) proved this conjecture. - Thomas Ordowski, Jul 20 2018
Note that n^2 divides Sum_{k=1..n-1} k^n for every odd number n > 1. - Thomas Ordowski, Oct 30 2015
Conjecture: these are "anti-Carmichael numbers" defined; n > 1 such that p - 1 does not divide n - 1 for every prime p dividing n. Equivalently, odd numbers n > 1 such that n is coprime to A027642(n-1). A number n > 1 is an "anti-Carmichael" if and only if gcd(n, b^n - b) = 1 for some integer b. - Thomas Ordowski, Jul 20 2018
It seems that these numbers are all composite terms of A317358. - Thomas Ordowski, Jul 30 2018
a(62) = 697 is the first term not in A267999: see A306097 for all these terms. - M. F. Hasler, Nov 09 2018
If the conjecture from Thomas Ordowski is true, then no term is a multiple of 2 or 3. - Jianing Song, Jan 27 2019
Conjecture: an odd number n > 1 is a term iff gcd(n, A027642(n-1)) = 1. - Thomas Ordowski, Jul 19 2019
Conjecture: Sequence consists of numbers n > 1 such that r = b^n + n - b will produce a prime for one or more integers b > 1. Only when n is in this sequence do one or more prime factors of n fail to divide r for all b. Also, n and b must be coprime for r to be prime. The above also applies to r = b^n - n - b, ignoring n=3, b=2. - Richard R. Forberg, Jul 18 2020
Odd numbers n > 1 such that Sum_{k(even)=2..n-1}2*k^(n-1) == 0 (mod n). - Davide Rotondo, Oct 28 2020
What is the asymptotic density of these numbers? The numbers A267999 have a slightly lower density. The difference between the densities is equal to the density of the numbers A306097. - Thomas Ordowski, Feb 15 2021
The asymptotic density of this sequence is in the interval (0.253, 0.265) (Ordowski, 2021). - Amiram Eldar, Feb 26 2021

Crossrefs

Cf. A000312, A002145, A002997, A027642, A031971, A038509, A060976, A121706, A267999 (probably a subsequence).
Cf. A306097 for terms of this sequence A121707 not in sequence A267999, A321487 for terms which are not semiprimes.
Cf. A191677 (n divides Sum_{k
Cf. A326478 for a conjectured connection with the Bernoulli numbers.

Programs

  • Maple
    filter:= n -> add(k &^ n mod n^3, k=1..n-1) mod n^3 = 0:
    select(filter, [$2..1000]); # Robert Israel, Oct 08 2015
  • Mathematica
    fQ[n_] := Mod[Sum[PowerMod[k, n, n^3], {k, n - 1}], n^3] == 0; Select[
    Range[2, 611], fQ] (* Robert G. Wilson v, Apr 04 2011 and slightly modified Aug 02 2018 *)
  • PARI
    is(n)=my(n3=n^3);sum(k=1,n-1,Mod(k,n3)^n)==0 \\ Charles R Greathouse IV, May 09 2013
    
  • PARI
    for(n=2, 1000, if(sum(k=1, n-1, k^n) % n^3 == 0, print1(n", "))) \\ Altug Alkan, Oct 15 2015
    
  • Sage
    # after Andrzej Schinzel
    def isA121707(n):
        if n == 1 or is_even(n): return False
        return n.divides(sum(k^(n-1) for k in (1..n-1)))
    [n for n in (1..611) if isA121707(n)] # Peter Luschny, Jul 18 2019

Extensions

Sequence corrected by Robert G. Wilson v, Apr 04 2011
Previous Showing 41-50 of 355 results. Next