cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Previous Showing 11-18 of 18 results.

A215240 Sum of the numbers p such that phi(p) = n, where phi is Euler's totient function.

Original entry on oeis.org

3, 13, 0, 35, 0, 48, 0, 105, 0, 33, 0, 166, 0, 0, 0, 231, 0, 138, 0, 218, 0, 69, 0, 621, 0, 0, 0, 87, 0, 93, 0, 581, 0, 0, 0, 655, 0, 0, 0, 833, 0, 276, 0, 299, 0, 141, 0, 1514, 0, 0, 0, 159, 0, 243, 0, 377, 0, 177, 0, 1114, 0, 0, 0, 1315, 0, 201, 0, 0, 0, 213, 0
Offset: 1

Views

Author

T. D. Noe, Oct 12 2012

Keywords

Comments

These terms (greater than 0) are not unique. The first duplicate appears at a(256) = a(2236) = 6711.

Crossrefs

Cf. A002181 (smallest inverse), A006511 (largest inverse), A217842 (product of inverses).
Cf. A007617, A032447 (inverse of phi).

Programs

  • Mathematica
    Needs["CNT`"]; Table[Total[PhiInverse[n]], {n, 100}]
  • PARI
    a(n) = vecsum(invphi(n)); \\ Amiram Eldar, Nov 15 2024, using Max Alekseyev's invphi.gp

Formula

a(n) = 0 if and only if n is in A007617. - Amiram Eldar, Nov 15 2024

A276044 Least k such that phi(k) has exactly n divisors.

Original entry on oeis.org

1, 3, 5, 7, 17, 13, 85, 31, 37, 65, 1285, 61, 4369, 193, 185, 143, 65537, 181, 327685, 241, 577, 3281, 5570645, 403, 1297, 12289, 1057, 1037, 286331153, 779, 1431655765, 899, 9509, 197633, 5629, 1333, 137438953472, 786433, 42653, 1763, 2199023255552, 2993, 8796093022208, 15361, 3737, 12648641
Offset: 1

Views

Author

Altug Alkan, Aug 17 2016

Keywords

Comments

Least k such that A000005(A000010(k)) = n.
From Jon E. Schoenfield, Nov 13 2016: (Start)
For every n > 0, phi(2^n) = 2^(n-1) has exactly n divisors, so a(n) <= 2^n.
For every prime p, since phi(a(p)) has exactly p divisors, phi(a(p)) must be of the form q^(p-1), where q is a prime number. If q >= 3, we would have phi(a(p)) >= 3^(p-1), and since k > phi(k) for every k > 1, we would have a(p) >= 3^(p-1)+1, which would be contradicted by the upper bound a(p) <= 2^p (see above) unless 3^(p-1)+1 <= 2^p, which is true only for p = 2. Thus, for every prime p > 2, phi(a(p)) = 2^(p-1), so a(p) > 2^(p-1). In summary, we can state that, for every prime p > 2:
(1) a(p) is the least k such that phi(k) = 2^(p-1), and
(2) 2^(p-1) < a(p) <= 2^p.
After a(36)=1333, the next few known terms are a(38)=786433, a(39)=42653, a(40)=1763, and a(42)=2993; as shown above, known bounds on a(37) and a(41) are 2^36 < a(37) <= 2^37 and 2^40 < a(41) <= 2^41.
For prime p < 37, a(p) = A001317(p-1).
Observation: for prime p < 37, a(p) is the product of distinct Fermat primes 2^(2^j)+1 for j=0..4, i.e., 3, 5, 17, 257, and 65537 (see A019434), according to the locations of the 1-bits in p-1:
. p-1 in
p a(p) prime factorization of a(p) binary
== ========== =========================== ======
2 3 = 3 1
3 5 = 5 10
5 17 = 17 100
7 85 = 17 * 5 110
11 1285 = 257 * 5 1010
13 4369 = 257 * 17 1100
17 65537 = 65537 10000
19 327685 = 65537 * 5 10010
23 5570645 = 65537 * 17 * 5 10110
29 286331153 = 65537 * 257 * 17 11100
31 1431655765 = 65537 * 257 * 17 * 5 11110
.
This pattern does not continue to p=37, since 2^(2^5)+1 is not prime. (See also A038183 and the observation there from Arkadiusz Wesolowski.) (End)
As noted, for every prime p, phi(a(p))=2^(p-1), decompose a(p) = p_1^(e_1) *...* p_m^(e_m), then phi(a(p)) = p_1^(e_1-1)*(p_1 - 1) * ... * p_m^(e_m-1)*(p_m - 1). Thus a(p) is of the form 2^e * F_(a_1) *...* F_(a_l), where F_(a_i) = 2^(a_i) + 1 denote distinct Fermat primes. If e = 0, a_1 + ... + a_l = p - 1, while if e > 0, e + a_1 + ... + a_l = p. It can be deduced that a(p) = 2^p unless p-1 can be written as a_1 + ... a_l where 2^(a_i) + 1 are distinct Fermat primes. The only Fermat primes known have a_i in {1,2,4,8,16} and it is known that 2^a + 1 is composite for 16 < a < 2^33 (cf. A019434). It follows from the fact that 1 + 2 + 4 + 8 + 16 = 31 that a(p) = 2^p for primes p with 32 < p <= 2^33. - Pjotr Buys, Sep 18 2019

Examples

			a(5) = 17 because phi(17) = 16 has 5 positive divisors.
		

Crossrefs

Programs

  • Mathematica
    Table[k = 1; While[DivisorSigma[0, #] &@ EulerPhi@ k != n, k++]; k, {n, 28}] (* Michael De Vlieger, Aug 21 2016 *)
  • PARI
    a(n) = {my(k = 1); while(numdiv(eulerphi(k)) != n, k++); k; }

Formula

a(p) = 2^p for primes p with 32 < p <= 2^33. - Pjotr Buys, Sep 18 2019

Extensions

a(31)-a(36) from Michel Marcus and Jon E. Schoenfield, Nov 13 2016
a(37)-a(46) from Pjotr Buys, Sep 18 2019

A217842 Product of the numbers p such that phi(p) = n, where phi is Euler's totient function.

Original entry on oeis.org

2, 72, 1, 4800, 1, 15876, 1, 3456000, 1, 242, 1, 300500928, 1, 1, 1, 2130739200, 1, 1052676, 1, 119790000, 1, 1058, 1, 531598161669120000, 1, 1, 1, 1682, 1, 1922, 1, 20864198246400, 1, 1, 1, 1159208596538496, 1, 1, 1, 265804426800000000, 1, 17757796, 1
Offset: 1

Views

Author

T. D. Noe, Oct 12 2012

Keywords

Comments

It appears that all terms greater than 1 are distinct. This is true for all n <= 10^6.

Crossrefs

Cf. A002181 (smallest inverse), A006511 (largest inverse), A215240 (sum of inverses).
Cf. A032447 (inverse of phi).

Programs

  • Mathematica
    Needs["CNT`"]; Table[Times @@ PhiInverse[n], {n, 100}]
  • PARI
    a(n) = vecprod(invphi(n)); \\ Amiram Eldar, Nov 15 2024, using Max Alekseyev's invphi.gp

A297475 Numbers n such that phi(x) = n for more than one value of x, and the smallest such x divides the largest.

Original entry on oeis.org

1, 2, 8, 10, 22, 28, 30, 44, 46, 52, 54, 56, 58, 66, 70, 78, 82, 92, 102, 104, 106, 110, 116, 126, 128, 130, 136, 138, 140, 148, 150, 164, 166, 172, 178, 190, 196, 198, 204, 210, 212, 222, 226, 228, 238, 250, 260, 262, 268, 270, 282, 292, 294, 296, 306, 310, 316, 330, 332, 342, 344, 346, 356, 358, 366, 368, 372
Offset: 1

Views

Author

Torlach Rush, Dec 30 2017

Keywords

Comments

The larger endpoint is always twice the value of the smaller endpoint.
Conjecture 1: The number of solutions, excluding endpoints is always 0, or an odd number. (known to n = 2 * 10^5)
Conjecture 2: If both endpoints are divisible by 5, then the number of solutions (excluding terms of A007366) is of the form 4k + 1. (known to n = 2 * 10^5)
A007366 is contained in this sequence and the number of solutions, excluding endpoints is always 0.
Terms of this sequence are totients with a single odd totient inverse.

Examples

			2 is in the sequence because {phi^-1(2)} = {3,4,6}, and 2 = 6 / 3.
8 is in the sequence because {phi^-1(8)} = {15,...,30}, and 2 = 30 / 15.
10 is in the sequence because {phi^-1(10)} = {11,22}, and 2 = 22 / 11.
		

Crossrefs

Programs

  • Mathematica
    With[{nn = 67}, Take[#, nn] &@ Keys@ Select[KeySort@ PositionIndex@ Array[EulerPhi, nn^2], IntegerQ[#2/#1] & @@ {First@ #, Last@ #} &]] (* Michael De Vlieger, Dec 31 2017 *)
  • PARI
    isok(n) = my(vx = invphi(n)); (#vx > 1) && ((vecmax(vx) % vecmin(vx)) == 0); \\ Michel Marcus, Jul 18 2018

Formula

2 = max({phi^-1(n)}) / min({phi^-1(n)}).
0 = A006511(n) mod A002181(n).

A143423 Least even number k such that phi(k) = n, where n runs through the values (A002202) taken by phi.

Original entry on oeis.org

2, 4, 8, 14, 16, 22, 26, 32, 38, 44, 46, 52, 58, 62, 64, 74, 82, 86, 92, 94, 104, 106, 162, 116, 118, 122, 128, 134, 142, 146, 158, 164, 166, 172, 178, 188, 194, 202, 206, 212, 214, 218, 242, 226, 236, 244, 254, 256, 262, 268, 274, 278, 284, 292, 298, 302, 314
Offset: 1

Views

Author

T. D. Noe, Aug 14 2008

Keywords

Comments

Such an even number always exists.

Crossrefs

Cf. A002181 (least k such that phi(k)=n), A006511 (largest k such that phi(k)=n).

Programs

  • Maple
    f:= proc(m) local L;
          L:= numtheory:-invphi(m);
          if L = [] then NULL
          else min(select(type,L,even))
          fi
    end proc:
    map(f, [1,seq(2*k,k=1..1000)]); # Robert Israel, Oct 07 2015
  • Mathematica
    f[m_] := Module[{L}, L = invphi[m]; If[L == {}, Nothing, Min[Select[L, EvenQ]]]];
    f /@ Join [{1}, 2 Range[1000]] (* Jean-François Alcover, Aug 28 2020, using Maxim Rytin's invphi function *)

A143511 Least number k such that phi(k) = n, where n runs through the values in A143510.

Original entry on oeis.org

33817088, 67634176, 101451264, 135268352, 169085440, 202902528, 236719616, 270536704, 8589934592
Offset: 1

Views

Author

T. D. Noe, Aug 21 2008

Keywords

References

Crossrefs

Cf. A002181 (least k such that phi(k)=n).

A196079 Difference between the largest and smallest inverse of totient function.

Original entry on oeis.org

1, 3, 7, 11, 15, 11, 29, 43, 35, 41, 23, 55, 29, 31, 69, 89, 109, 55, 69, 47, 145, 53, 81, 87, 59, 137, 155, 67, 71, 197, 79, 207, 83, 165, 187, 141, 323, 149, 103, 159, 107, 269, 121, 235, 177, 319, 127, 255, 131, 253, 137, 139, 213, 445, 149, 151
Offset: 1

Views

Author

Franz Vrabec, Sep 27 2011

Keywords

Comments

No terms are zero if Carmichael's conjecture is true.
Even terms are rare: e.g., all inverses of 257*2^16 are even [Foster], so the difference between the largest and smallest inverse is even.

Examples

			Let n=3. The largest inverse of A002202(3)=4 is A006511(3)=12, the smallest inverse is A002181(3)=5, so a(3)=12-5=7.
		

Crossrefs

Programs

  • Mathematica
    max = 300; inversePhi[?OddQ] = {}; inversePhi[1] = {1, 2}; inversePhi[m] := Module[{p, nmax, n, nn}, p = Select[Divisors[m] + 1, PrimeQ]; nmax = m*Times @@ (p/(p - 1)); n = m; nn = Reap[While[n <= nmax, If[EulerPhi[n] == m, Sow[n]]; n++]] // Last; If[nn == {}, {}, First[nn]]]; Join[{2}, Reap[For[n = 2, n <= max, n = n + 2, nn = inversePhi[n] ; If[ nn != {} , Sow[Max[nn] - Min[nn]]]]] // Last // First] (* Jean-François Alcover, Nov 21 2013 *)

Formula

a(n) = A006511(n) - A002181(n).

Extensions

a(1) corrected by the editors, Nov 23 2013
a(1) in b-file corrected by Andrew Howroyd, Feb 22 2018

A281627 a(n) is the smallest number k such that sigma(phi(k)) = A062402(k) is the n-th Mersenne prime (A000668(n)), or 0 if no such k exists.

Original entry on oeis.org

3, 5, 17, 85, 4369, 65537, 327685, 1431655765, 2305843009213693952, 618970019642690137449562112, 162259276829213363391578010288128, 170141183460469231731687303715884105728
Offset: 1

Views

Author

Jaroslav Krizek, Feb 11 2017

Keywords

Comments

Conjecture 1: If A062402(n) = A000203(A000010(n)) = sigma(phi(n)) is a prime p for some n, then p is Mersenne prime (A000668); a(n) > 0 for all n.
Conjecture 2: a(n) = the smallest number k such that phi(k) has exactly A000043(n)-1 divisors; see A276044.
Conjecture 3: a(n) = the smallest number k such that phi(k) has exactly A000043(n)-1 prime factors (counted with multiplicity); see A275969.
a(n) <= A000668(n) for n = 1-18; conjecture: a(n) <= A000668(n) for all n.
Equals A002181 (index in A002202 of (intersection of A023194 and A002202)). - Michel Marcus, Feb 12 2017

Crossrefs

Cf. A053576 (includes the first 13 known terms of this sequence).

Programs

  • Magma
    A281627:=func; Set(Sort([A281627(n): n in [SumOfDivisors(EulerPhi(n)): n in[1..1000000] | IsPrime(SumOfDivisors(EulerPhi(n)))]]));
    
  • PARI
    terms() = {v = readvec("b023194.txt"); for(i=1, #v, if (istotient(v[i], &n), print1(n/2, ", ")););} \\ Michel Marcus, Feb 12 2017
    
  • PARI
    f(p) = {my(s = invsigma(p), t, minv = 0); for(i = 1 ,#s, t = invphi(s[i]); for(j = 1, #t, if(minv == 0, minv = t[j]); if(t[j] < minv, minv = t[j]))); minv;} \\ using Max Alekseyev's invphi.gp
    list(pmax) = forprime(p = 1, pmax, if(isprime(2^p-1), print1(f(2^p-1), ", "))); \\ Amiram Eldar, Dec 23 2024

Extensions

a(8) from Michel Marcus, Feb 12 2017
a(9)-a(12) from Amiram Eldar, Dec 23 2024
Previous Showing 11-18 of 18 results.