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

A366716 a(n) = phi(12^n+1), where phi is Euler's totient function (A000010).

Original entry on oeis.org

1, 12, 112, 1296, 20416, 229680, 2306304, 32916240, 400515072, 3863116800, 47825825600, 685853880624, 8732596764672, 97509650382144, 990242755633152, 11148606564480000, 184883057981234176, 2047145911595946000, 20281543142263603200, 294779525244632305920
Offset: 0

Views

Author

Sean A. Irvine, Oct 17 2023

Keywords

Crossrefs

Programs

  • Mathematica
    EulerPhi[12^Range[0,19] + 1] (* Paul F. Marrero Romero, Oct 27 2023 *)
  • PARI
    {a(n) = eulerphi(12^n+1)}

Formula

a(n) = A000010(A178248(n)). - Paul F. Marrero Romero, Oct 27 2023

A369291 Array read by antidiagonals: T(n,k) = phi(k^n-1)/n, where phi is Euler's totient function (A000010), n >= 1, k >= 2.

Original entry on oeis.org

1, 1, 1, 2, 2, 2, 2, 4, 4, 2, 4, 4, 12, 8, 6, 2, 12, 20, 32, 22, 6, 6, 8, 56, 48, 120, 48, 18, 4, 18, 36, 216, 280, 288, 156, 16, 6, 16, 144, 160, 1240, 720, 1512, 320, 48, 4, 30, 96, 432, 1120, 5040, 5580, 4096, 1008, 60, 10, 16, 216, 640, 5400, 6048, 31992, 14976, 15552, 2640, 176
Offset: 1

Views

Author

Andrew Howroyd, Jan 28 2024

Keywords

Comments

For k a prime power, T(n,k) is the number of primitive polynomials of degree n over GF(k). See A011260, A027385 for additional information.

Examples

			Array begins:
n\k|  2   3    4     5      6      7      8       9 ...
---+---------------------------------------------------
 1 |  1   1    2     2      4      2      6       4 ...
 2 |  1   2    4     4     12      8     18      16 ...
 3 |  2   4   12    20     56     36    144      96 ...
 4 |  2   8   32    48    216    160    432     640 ...
 5 |  6  22  120   280   1240   1120   5400    5280 ...
 6 |  6  48  288   720   5040   6048  23328   27648 ...
 7 | 18 156 1512  5580  31992  37856 254016  340704 ...
 8 | 16 320 4096 14976 139968 192000 829440 1966080 ...
  ...
		

Crossrefs

Rows n=1..3 and 5 are A000010(k-1), A319210, A319213, A319214.
Cf. A319183.

Programs

  • Mathematica
    A369291[n_, k_] := EulerPhi[k^n - 1]/n;
    Table[A369291[k, n-k+2], {n, 15}, {k, n}] (* Paolo Xausa, Jun 17 2024 *)
  • PARI
    T(n,k) = eulerphi(k^n-1)/n

A112954 Number of numbers m such that phi(m) = n*tau(m), with phi=A000010 and tau=A000005.

Original entry on oeis.org

7, 9, 10, 9, 7, 17, 4, 17, 14, 15, 7, 19, 2, 16, 20, 21, 0, 29, 0, 29, 9, 13, 7, 32, 7, 11, 23, 21, 7, 39, 0, 19, 17, 4, 11, 44, 2, 0, 11, 41, 7, 24, 2, 19, 30, 11, 0, 55, 4, 23, 7, 21, 7, 46, 9, 27, 4, 11, 0, 61, 0, 0, 27, 29, 9, 30, 2, 10, 19, 31, 0, 57, 2, 9, 27, 4, 4, 30, 2, 50, 29, 9
Offset: 1

Views

Author

Reinhard Zumkeller, Oct 07 2005

Keywords

Crossrefs

Extensions

More terms from Max Alekseyev, Mar 01 2010

A132191 Square array a(m,n) read by antidiagonals, defined by A000010(n)*a(m,n) = Sum_{k=1..n, gcd(k,n)=1} m^{ Sum_{d|n} A000010(d)/ (multiplicative order of k modulo d) }.

Original entry on oeis.org

1, 1, 2, 1, 4, 3, 1, 6, 9, 4, 1, 12, 18, 16, 5, 1, 12, 54, 40, 25, 6, 1, 40, 72, 160, 75, 36, 7, 1, 28, 405, 280, 375, 126, 49, 8, 1, 96, 390, 2176, 825, 756, 196, 64, 9, 1, 104, 1944, 2800, 8125, 2016, 1372, 288, 81, 10, 1, 280, 3411, 17920, 13175, 23976, 4312, 2304, 405
Offset: 1

Views

Author

N. J. A. Sloane, Dec 01 2007, based on email from Max Alekseyev, Nov 08 2007

Keywords

Comments

From Andrew Howroyd, Apr 22 2017: (Start)
Number of step shifted (decimated) sequences of length n using a maximum of m different symbols. See A056371 for an explanation of step shifts. -
Number of mappings with domain {0..n-1} and codomain {1..m} up to equivalence. Mappings A and B are equivalent if there is a d, prime to n, such that A(i) = B(i*d mod n) for i in {0..n-1}. (End)

Examples

			Array begins:
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, ...
2, 4, 6, 12, 12, 40, 28, 96, 104, 280, 216, 1248, 704, 2800, 4344, 8928, 8232, 44224, 29204, 136032, ...
3, 9, 18, 54, 72, 405, 390, 1944, 3411, 14985, 17802, 139968, 133104, 798525, 1804518, 5454378, 8072532, 64599849, 64573626, 437732424, ...
4, 16, 40, 160, 280, 2176, 2800, 17920, 44224, 263296, 419872, 4280320, 5594000, 44751616, 134391040, 539054080, 1073758360, 11453771776, 15271054960, 137575813120, ...
5, 25, 75, 375, 825, 8125, 13175, 103125, 327125, 2445625, 4884435, 61640625, 101732425, 1017323125, 3816215625, 19104609375, 47683838325, 635787765625, 1059638680675, 11924780390625, ...
		

Crossrefs

Row m=2 is A056371
Row m=3 is A056372
Row m=4 is A056373
Row m=5 is A056374
Row m=6 is A056375
Column n=2 is A000290
Column n=3 is A002411
Column n=4 is A019582

Programs

  • Mathematica
    a[m_, n_] := (1/EulerPhi[n])*Sum[If[GCD[k, n]==1, m^DivisorSum[n, EulerPhi[#] / MultiplicativeOrder[k, #]&], 0], {k, 1, n}]; Table[a[m-n+1, n], {m, 1, 15}, {n, m, 1, -1}] // Flatten (* Jean-François Alcover, Dec 01 2015 *)
  • PARI
    for(i=1,15,for(m=1,i,n=i-m+1; print1(sum(k=1, n, if(gcd(k,n)==1, m^sumdiv(n,d,eulerphi(d)/znorder(Mod(k,d))),0))/eulerphi(n)","))) \\ Herman Jamke (hermanjamke(AT)fastmail.fm), Apr 26 2008

Extensions

More terms from Herman Jamke (hermanjamke(AT)fastmail.fm), Apr 26 2008
Offset corrected by Andrew Howroyd, Apr 20 2017

A295500 a(n) = phi(3^n-1), where phi is Euler's totient function (A000010).

Original entry on oeis.org

1, 4, 12, 32, 110, 288, 1092, 2560, 9072, 26400, 84700, 165888, 797160, 2384928, 6019200, 15728640, 64533700, 141087744, 580765248, 1246080000, 4823425152, 14758128000, 46070066188, 85996339200, 385087175000, 1270928131200, 3474144608256, 8810420097024
Offset: 1

Views

Author

Seiichi Manyama, Nov 22 2017

Keywords

Crossrefs

phi(k^n-1): A053287 (k=2), this sequence (k=3), A295501 (k=4), A295502 (k=5), A366623 (k=6), A366635 (k=7), A366654 (k=8), A366663 (k=9), A295503 (k=10), A366685 (k=11), A366711 (k=12).

Programs

  • Mathematica
    EulerPhi[3^Range[30] - 1] (* Paolo Xausa, Jun 18 2024 *)
  • PARI
    {a(n) = eulerphi(3^n-1)}

Formula

a(n) = n*A027385(n).
a(n) = A000010(A024023(n)). - Michel Marcus, Jun 18 2024

A366690 a(n) = phi(11^n+1), where phi is Euler's totient function (A000010).

Original entry on oeis.org

1, 4, 60, 432, 7320, 53680, 803520, 6495720, 100874752, 764738496, 12756110400, 89493288192, 1568774615040, 11278053084480, 180228847518720, 1310982643872000, 22974417331646464, 168479281019744640, 2521788545778163200, 20190830281379049600
Offset: 0

Views

Author

Sean A. Irvine, Oct 16 2023

Keywords

Crossrefs

Programs

  • Mathematica
    EulerPhi[11^Range[0,19] + 1] (* Paul F. Marrero Romero, Nov 10 2023 *)
  • PARI
    {a(n) = eulerphi(11^n+1)}

Formula

a(n) = A000010(A034524(n)). - Paul F. Marrero Romero, Nov 10 2023

A366711 a(n) = phi(12^n-1), where phi is Euler's totient function (A000010).

Original entry on oeis.org

10, 120, 1560, 13440, 226200, 2021760, 32518360, 274391040, 4534807680, 51953616000, 646094232960, 4662793175040, 97266341877120, 1070382142166400, 13666309113600000, 109897747141754880, 2016918439151095000, 17518491733377024000, 290436363064202660760
Offset: 1

Views

Author

Sean A. Irvine, Oct 17 2023

Keywords

Crossrefs

phi(k^n-1): A053287 (k=2), A295500 (k=3), A295501 (k=4), A295502 (k=5), A366623 (k=6), A366635 (k=7), A366654 (k=8), A366663 (k=9), A295503 (k=10), A366685 (k=11), this sequence (k=12).

Programs

  • Mathematica
    EulerPhi[12^Range[30] - 1]
  • PARI
    {a(n) = eulerphi(12^n-1)}

A061026 Smallest number m such that phi(m) is divisible by n, where phi = Euler totient function A000010.

Original entry on oeis.org

1, 3, 7, 5, 11, 7, 29, 15, 19, 11, 23, 13, 53, 29, 31, 17, 103, 19, 191, 25, 43, 23, 47, 35, 101, 53, 81, 29, 59, 31, 311, 51, 67, 103, 71, 37, 149, 191, 79, 41, 83, 43, 173, 69, 181, 47, 283, 65, 197, 101, 103, 53, 107, 81, 121, 87, 229, 59, 709, 61, 367, 311, 127, 85
Offset: 1

Views

Author

Melvin J. Knight (knightmj(AT)juno.com), May 25 2001

Keywords

Comments

Conjecture: a(n) is odd for all n. Verified up to n <= 3*10^5. - Jianing Song, Feb 21 2021
The conjecture above is false because a(16842752) = 33817088; see A002181 and A143510. - Flávio V. Fernandes, Oct 08 2023

Examples

			a(48) = 65 because phi(65) = phi(5)*phi(13) = 4*12 = 48 and no smaller integer m has phi(m) divisible by 48.
		

Crossrefs

Cf. A233516, A233517 (records).
Cf. A005179 (analog for number of divisors), A070982 (analog for sum of divisors).

Programs

  • Mathematica
    a = ConstantArray[1, 64]; k = 1; While[Length[vac = Rest[Flatten[Position[a, 1]]]] > 0, k++; a[[Intersection[Divisors[EulerPhi[k]], vac]]] *= k]; a  (* Ivan Neretin, May 15 2015 *)
  • PARI
    a(n) = my(s=1); while(eulerphi(s)%n, s++); s;
    vector(100, n, a(n))
    
  • Python
    from sympy import totient as phi
    def a(n):
      k = 1
      while phi(k)%n != 0: k += 1
      return k
    print([a(n) for n in range(1, 65)]) # Michael S. Branicky, Feb 21 2021

Formula

Sequence is unbounded; a(n) <= n^2 since phi(n^2) is always divisible by n.
If n+1 is prime then a(n) = n+1.
a(n) = min{ k : phi(k) == 0 (mod n) }.
a(n) = a(2n) for odd n > 1. - Jianing Song, Feb 21 2021

A066674 Least number m such that phi(m) = A000010(m) is divisible by the n-th prime.

Original entry on oeis.org

3, 7, 11, 29, 23, 53, 103, 191, 47, 59, 311, 149, 83, 173, 283, 107, 709, 367, 269, 569, 293, 317, 167, 179, 389, 607, 619, 643, 1091, 227, 509, 263, 823, 557, 1193, 907, 1571, 653, 2339, 347, 359, 1087, 383, 773, 3547, 797, 2111, 2677, 5449, 2749, 467
Offset: 1

Views

Author

Labos Elemer, Dec 19 2001

Keywords

Comments

All terms seem to be primes of the form a(n) = k*prime(n)+1 for some k.
Is this a duplicate of A035095? - R. J. Mathar, Dec 13 2008
For the first 5*10^6 terms, a(n) = A035095(n). - Donovan Johnson, Oct 21 2011
Comments on the relationship between A035095, A066674, A125878, added by N. J. A. Sloane, Jan 07 2013: (Start)
Let a(n) = A066674(n), b(n) = A035095(n), c(n) = A125878(n).
It is immediate from the definitions that a(n) <= b(n) and a(n) <= c(n).
Bjorn Poonen (Jan 06 2013) makes the following observations:
1) A prime p divides phi(m) if and only if p^2 | m or p | q-1 for some prime q | m. Thus the smallest m for p is either p^2 or the smallest prime q = 1 (mod p). In other words, a(n) = min(b(n),p(n)^2).
2) In particular, the m in the definition of a(n) is at most p(n)^2, so phi(m)/p(n) < p(n), so p(n) is the largest prime dividing phi(m), and phi(m)/(2 p(n)) < p(n)/2 < p(n-1), so p(n-1) does not divide phi(m)/2.
Thus c(n) = a(n).
Further comments from Eric Bach, Jan 07 2013: (Start)
As others have pointed out, the possible equivalence of a(n) and b(n) is basically the question of how quickly the least prime q == 1 mod p grows, as a function of p. In particular, if q < p^2, the two sequences are the same.
Here are some remarks connected with this.
1. There are probabilistic arguments suggesting that q = O(p (log p)^2). See Heath-Brown (1978), Wagstaff (1979), Bach and Huelsbergen (1993). Using the sieve of Eratosthenes, I found no exceptions to q < p^2 below p = 1254767. So it seems likely that a(n) and b(n) are the same.
2. If ERH holds, then q = O(p log p)^2, see Heath-Brown (1990), (1992). Explicitly, on the same hypothesis, q < 2(p log p)^2, see Bach and Sorenson (1996).
3. By Linnik's theorem, q = O(p^c) for some c > 0. This is unconditional, but the best known value of c, equal to 5.18 -- see Xylouris (2011) -- is nowhere near 2. Heath-Brown (1992) mentions the conjecture (generalized to Linnik's theorem) that q <= p^2. If true, a(n) and b(n) are identical, since p^2 cannot be 1 mod p. (End)
Don Reble (Jan 07 2013) observes that A074884 and A117673 are related to these questions.
Summary: A066674 and A125878 are the same, and A035095 is probably also the same, but this is an open question.
(End)

References

  • E. Bach and J. Shallit, Algorithmic Number Theory, Vol. 1: Efficient Algorithms, MIT Press, Cambridge, MA, 1996.

Crossrefs

Programs

  • Mathematica
    f[n_] := Block[{m = p = Prime@ n}, While[ Mod[ EulerPhi@ m, p] != 0, m += 2]; m]; f[1] = 3; Array[f, 60] (* Robert G. Wilson v, Dec 27 2014 *)

Formula

a(n) = min{m : phi(m) = 0 mod prime(n) = 0}.

Extensions

a(2) corrected by R. J. Mathar, Dec 13 2008

A179186 Numbers k such that phi(k) = phi(k+4), with Euler's totient function phi = A000010.

Original entry on oeis.org

8, 14, 16, 20, 35, 52, 64, 91, 140, 148, 244, 292, 403, 455, 616, 628, 772, 801, 1011, 1024, 1108, 1144, 1252, 1270, 1295, 1456, 1588, 1684, 1820, 1828, 2030, 2164, 2452, 2623, 2644, 2660, 2692, 2932, 3028, 3216, 3321, 3508, 3988, 4264, 4340, 4372, 4612, 4804, 4852, 4948
Offset: 1

Views

Author

M. F. Hasler, Jan 05 2011

Keywords

Comments

Is there some k > 5 such that phi(k) = phi(k+3)?
None up to 500000. - Harvey P. Dale, Feb 16 2011
No further solutions to the phi(k) = phi(k+3) problem less than 10^12. On the other hand, this sequence has 267797240 terms under 10^12. - Jud McCranie, Feb 13 2012
No reason is known that would prevent other solutions of phi(k) = phi(k+3), see Graham, Holt, & Pomerance. - Jud McCranie, Jan 03 2013
If a(n) is even then a(n)/2 is in A001494 - see comment at A217139. - Jud McCranie, Dec 31 2012

References

  • S. W. Graham, J. J. Holt, and C. Pomerance, "On the solutions to phi(n)=phi(n+k)", Number Theory in Progress, Proc. Intern. Conf. in Honor of 60th Birthday of A. Schinzel, Poland, 1997. Walter de Gruyter, 1999, pp. 867-82.

Crossrefs

Programs

  • Magma
    [n: n in [1..5000] | EulerPhi(n) eq EulerPhi(n+4)]; // Vincenzo Librandi, Sep 08 2016
  • Mathematica
    Select[Range[5000],EulerPhi[#]==EulerPhi[#+4]&]  (* Harvey P. Dale, Feb 16 2011 *)
    SequencePosition[EulerPhi[Range[5000]],{x_,,,_,x_}][[;;,1]] (* Harvey P. Dale, Sep 12 2024 *)
  • PARI
    {op=vector(N=4); for( n=1,1e4,if( op[n%N+1]+0==op[n%N+1]=eulerphi(n),print1(n-N,",")))}
    
Previous Showing 41-50 of 4244 results. Next