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

A046022 Primes together with 1 and 4.

Original entry on oeis.org

1, 2, 3, 4, 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
Offset: 1

Views

Author

Keywords

Comments

Also the numbers which are incrementally largest values of A002034. - validated by Franklin T. Adams-Watters, Jul 13 2012
Solutions to A000005(x) + A000010(x) - x - 1 = 0. - Labos Elemer, Aug 23 2001
Also numbers m such that m, phi(m) and tau(m) form an integer triangle, where phi=A000010 is the totient and tau=A000005 the number of divisors (see also A084820). - Reinhard Zumkeller, Jun 04 2003
Terms > 1 are n such that n does not divide (n-1)!. - Benoit Cloitre, Nov 12 2003
Terms > 1 are the sum of their prime factors; 4 (= 2+2) is the only such composite number. - Stuart Orford (sjorford(AT)yahoo.co.uk), Aug 04 2005
From Jonathan Vos Post, Aug 23 2010, Robert G. Wilson v, Aug 25 2010, proof by D. S. McNeil, Aug 29 2010: (Start)
Also the numbers n which divide A001414(n), or equivalently divide A075254(n). Proof:
Theorem: for a multiset of m >= 2 integers a_i, each a_i >= 2, Product_{i=1..m} a_i >= Sum_{i=1..m} a_i, with equality only at (a_1,a_2) = (2,2).
Lemma: For integers x,y >= 2, if x > 2 or y > 2, x*y > x + y. This follows from distributing (x-1)*(y-1) > 1.
[Proof of the theorem by induction on m:
first consider m=2. We have equality at (2,2) and for any product(a_i) > 4 there is some a_i > 2, so the lemma gives a_1*a_2 > a_1+a_2.
Then the induction m->m+1: Product_{i=1..m+1} a_i = a_(m+1)*Product_{i=1..m} a_i >= a_(m+1) * Sum_{i=1..m} a_i.
Since a_(m+1) >= 2 and the sum >= 4, the lemma applies, and we find a_(m+1) * Sum+{i=1..m} a_i > a_(m+1) + Sum_{i=1..m} a_i = Sum_{i=1..m+1} a_i and thus Product_{i=1..m+1} a_i > Sum_{i=1..m+1} a_i, QED.]
For composite n > 4, applying the theorem to the multiset of prime factors with multiplicity yields n > sopfr(n), so there are no composite numbers greater than 4 such that they divide sopfr(n).
(End)
Numbers k such that the k-th Fibonacci number is relatively prime to all smaller Fibonacci numbers. - Charles R Greathouse IV, Jul 13 2012
Numbers k such that (-1)^k*floor(d(k)*(-1)^k/2) = 1, where d(k) is the number of divisors of k. - Wesley Ivan Hurt, Oct 11 2013
Also, union of odd primes (A065091) and the divisors of 4. Also, union of A008578 and 4. - Omar E. Pol, Nov 04 2013
Numbers k such that sigma(k!) is divisible by sigma((k-1)!). - Altug Alkan, Jul 18 2016

Crossrefs

Programs

  • Haskell
    a046022 n = a046022_list !! (n-1)
    a046022_list = [1..4] ++ drop 2 a000040_list
    -- Reinhard Zumkeller, Apr 06 2014
    
  • Maple
    A046022:=n-> `if`((-1)^n*floor(numtheory[tau](n)*(-1)^n/2) = 1, n, NULL); seq(A046022(j), j=1..260); # Wesley Ivan Hurt, Oct 11 2013
  • Mathematica
    max = 0; a = {}; Do[m = FactorInteger[n]; w = Sum[m[[k]][[1]]*m[[k]][[2]], {k, 1, Length[m]}]; If[w > max, AppendTo[a, w]; max = w], {n, 1, 1000}]; a (* Artur Jasinski, Apr 06 2008 *)
  • PARI
    a(n)=if(n<6,n,prime(n-2)) \\ Charles R Greathouse IV, Apr 28 2015
    
  • Python
    from sympy import prime
    def A046022(n): return prime(n-2) if n>4 else n # Chai Wah Wu, Oct 17 2024

Formula

A141295(a(n)) = a(n). - Reinhard Zumkeller, Jun 23 2008
A018194(a(n)) = 1. - Reinhard Zumkeller, Mar 09 2012
A240471(a(n)) = 1. - Reinhard Zumkeller, Apr 06 2014

Extensions

Better description from Frank Ellermann, Jun 15 2001

A084821 Odd numbers m such that sigma(m) > m + phi(m), where sigma=A000203 is the divisor sum and phi=A000010 the totient.

Original entry on oeis.org

15, 45, 63, 75, 105, 135, 165, 189, 195, 225, 231, 255, 273, 285, 297, 315, 345, 357, 375, 399, 405, 429, 435, 441, 465, 483, 495, 525, 555, 567, 585, 609, 615, 645, 651, 675, 693, 705, 735, 765, 777, 795, 819, 825, 855, 861, 885, 891, 903, 915, 945, 975
Offset: 1

Views

Author

Reinhard Zumkeller, Jun 04 2003

Keywords

Comments

This sequence is closed under multiplication by odd values, and in particular with multiplication by itself. - Charles R Greathouse IV, Feb 19 2013

Crossrefs

Programs

  • Mathematica
    Select[Range[1, 1000, 2], DivisorSigma[1, #] > EulerPhi[#] + # &] (* Amiram Eldar, Sep 12 2019 *)
  • PARI
    is(n)=n%2 && sigma(n)>eulerphi(n)+n \\ Charles R Greathouse IV, Feb 19 2013

A268281 Numbers n such that n-tau(n), phi(n) and n form a Heronian triangle, where tau=A000005 is the number of divisors and phi=A000010 the totient.

Original entry on oeis.org

5, 34, 53, 90, 120, 440, 780, 1954, 120994, 140453, 28813276834
Offset: 1

Views

Author

Frank M Jackson, Jan 29 2016

Keywords

Comments

For all n, n > tau(n) and n > phi(n) and if n is prime then n-tau(n) = n-2 and phi(n) = n-1. So n = 5 gives the triangle {3, 4, 5} which is a primitive Pythagorean triangle and this is the only one. Other Pythagorean triangles are {30, 16, 34} and {756, 192, 780}, the remainder are only Heronian.
It is not known if this sequence is infinite. Prime numbers in the sequence are 5, 53 and 140453 and generate triangles {3, 4, 5}, {51, 52, 53} and {140451, 140452, 140453}.
If n = 2p where p is prime then n-tau(n) = n-4 and phi(n) = n/2-1. So n = 34 gives the triangle {16, 30, 34}. Similar numbers in this sequence are a(8), a(9) and a(11). See A272365 for generating Heronian triangles with sides n, n-4, n/2-1.
a(12) > 2*10^12. - Giovanni Resta, Apr 14 2016
Next prime value of a(n) after 140453 is > 2*10^5719. See A003500 for generating Heronian triangles with consecutive sides. - Frank M Jackson, Apr 19 2016
A003500(n)+1 is a member of this sequence iff it is prime. Also A272365(n) is a member of this sequence iff A272365(n)/2 is prime. - Frank M Jackson, Apr 29 2016

Examples

			a(2) = 34 because the triangle so formed has sides 30, 16, 34. It is Heronian with integer area 240 and is also Pythagorean. It is the second Heronian triangle.
The triangle corresponding to a(11) has sides n = 28813276834, n-tau(n) = 28813276830, phi(n) = 14406638416, and area 200960614753814018640.
		

Crossrefs

Programs

  • Mathematica
    triples[n_] := ({a, b, c}={n-DivisorSigma[0, n], EulerPhi[n], n}; s=(a+b+c)/2; If[a+b>c&&IntegerQ[Sqrt[s(s-a)(s-b)(s-c)]], {a, b, c}, {}]); lst={}; Do[If[triples[n]!={}, AppendTo[lst, Last[triples[n]]]], {n, 1, 200000}]; lst

Extensions

a(11) from Giovanni Resta, Apr 14 2016

A342144 Numbers m with integer solution to x^x == (x+1)^(x+1) (mod m) with x > 0.

Original entry on oeis.org

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

Views

Author

Owen C. Keith, Mar 01 2021

Keywords

Comments

Some values of m have multiple solutions.
For example, for m = 49, 25^25 == 26^26 (mod 49) and 37^37 == 38^38 (mod 49).
All terms are odd.
First differs from A334420 at a(70) which is 167 for this sequence and 165 for A334420.
First differs from A056911 at a(21) which is 49 for this sequence and 51 for A056911.

Examples

			3 is a term since 1^1 == 2^2 (mod 3).
5 is a term since 11^11 == 12^12 (mod 5).
		

Crossrefs

Programs

  • Mathematica
    seqQ[n_] := AnyTrue[Range[LCM[n, CarmichaelLambda[n]]+1], PowerMod[#, #, n] == PowerMod[# + 1, # + 1, n] &]; Select[Range[145], seqQ]
Showing 1-4 of 4 results.