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 21-25 of 25 results.

A032196 Number of necklaces with 11 black beads and n-11 white beads.

Original entry on oeis.org

1, 1, 6, 26, 91, 273, 728, 1768, 3978, 8398, 16796, 32066, 58786, 104006, 178296, 297160, 482885, 766935, 1193010, 1820910, 2731365, 4032015, 5864750, 8414640, 11920740, 16689036, 23107896, 31666376, 42975796
Offset: 11

Views

Author

Keywords

Comments

The g.f. is Z(C_11,x)/x^11, the 11-variate cycle index polynomial for the cyclic group C_11, with substitution x[i]->1/(1-x^i), i=1..11. By Polya enumeration, a(n+11) is the number of cyclically inequivalent 11-necklaces whose 11 beads are labeled with nonnegative integers such that the sum of labels is n, for n=0,1,2,... See A102190 for Z(C_11,x). See the comment in A032191 on the equivalence of this problem with the one given in the `Name' line. - Wolfdieter Lang, Feb 15 2005

Crossrefs

Programs

  • Mathematica
    k = 11; Table[Apply[Plus, Map[EulerPhi[ # ]Binomial[n/#, k/# ] &, Divisors[GCD[n, k]]]]/n, {n, k, 30}] (* Robert A. Russell, Sep 27 2004 *)
    DeleteCases[CoefficientList[Series[(x^11) (1 - 9 x + 41 x^2 - 109 x^3 + 191 x^4 - 229 x^5 + 191 x^6 - 109 x^7 + 41 x^8 - 9 x^9 + x^10)/((1 - x)^10 (1 - x^11)), {x, 0, 39}], x], 0] (* Michael De Vlieger, Oct 10 2016 *)

Formula

"CIK[ 11 ]" (necklace, indistinct, unlabeled, 11 parts) transform of 1, 1, 1, 1...
G.f.: (x^11) * (1 - 9*x + 41*x^2 - 109*x^3 + 191*x^4 - 229*x^5 + 191*x^6 - 109*x^7 + 41*x^8 - 9*x^9 + x^10) / ((1-x)^10 * (1-x^11)).
a(n) = ceiling(binomial(n, 11)/n) (conjecture Wolfdieter Lang).
From Herbert Kociemba, Oct 11 2016: (Start)
This conjecture indeed is true.
Sketch of proof:
There are binomial(n,11) ways to place the 11 black beads in the necklace with n beads. If n is not divisible by 11 there are no necklaces with a rotational symmetry. So exactly n necklaces are equivalent up to rotation and there are binomial(n,11)/n = ceiling(binomial(n,11)/n) equivalence classes.
If n is divisible by 11 the only way to get a necklace with rotational symmetry is to space out the 11 black beads evenly. There are n/11 ways to do this and all ways are equivalent up to rotation. So there are binomial(n,11) - n/11 unsymmetric necklaces which give binomial(n,11)/n - 1/11 equivalence classes. If we add the single symmetric equivalence class we get Binomial(n,11)/n - 1/11 + 1 which also is ceiling(binomial(n,11)/n). (End)
G.f.: (10/(1 - x^11) + 1/(1 - x)^11)*x^11/11. - Herbert Kociemba, Oct 16 2016

A032197 Number of necklaces with 12 black beads and n-12 white beads.

Original entry on oeis.org

1, 1, 7, 31, 116, 364, 1038, 2652, 6310, 14000, 29414, 58786, 112720, 208012, 371516, 643856, 1086601, 1789515, 2883289, 4552275, 7056280, 10752060, 16128424, 23841480, 34769374, 50067108, 71250060, 100276894, 139672312
Offset: 12

Views

Author

Keywords

Comments

The g.f. is Z(C_12,x)/x^12, the 12-variate cycle index polynomial for the cyclic group C_12, with substitution x[i]->1/(1-x^i), i=1,...,12. Therefore by Polya enumeration a(n+12) is the number of cyclically inequivalent 12-necklaces whose 12 beads are labeled with nonnegative integers such that the sum of labels is n, for n=0,1,2,... See A102190 for Z(C_12,x). See the comment in A032191 on the equivalence of this problem with the one given in the `Name' line. - Wolfdieter Lang, Feb 15 2005

Crossrefs

Programs

  • Mathematica
    k = 12; Table[Apply[Plus, Map[EulerPhi[ # ]Binomial[n/#, k/# ] &, Divisors[GCD[n, k]]]]/n, {n, k, 30}] (* Robert A. Russell, Sep 27 2004 *)

Formula

"CIK[ 12 ]" (necklace, indistinct, unlabeled, 12 parts) transform of 1, 1, 1, 1...
G.f.: (x^12)*(1-3*x+7*x^2+9*x^3+18*x^4+38*x^5+72*x^6+92*x^7+168*x^8+160*x^9+238*x^10+230*x^11+296*x^12+234*x^13+330*x^14+248*x^15+284*x^16+238*x^17+230*x^18+166*x^19+172*x^20+78*x^21+80*x^22+38*x^23+21*x^24+7*x^25+3*x^26+x^27) /((1+x)*(1-x)*(1-x^2)*(1-x^3)*(1-x)^5*(1+x+x^2)*(1-x^4)^2*(1-x^6)*(1-x^12)). - Wolfdieter Lang, Feb 15 2005 (see comment)
G.f.: 1/12 x^12 ((1 - x)^-12 + (1 - x^2)^-6 + 2 (1 - x^3)^-4 + 2 (1 - x^4)^-3 + 2 (1 - x^6)^-2 + 4 (1 - x^12)^-1). - Herbert Kociemba, Oct 22 2016

A328857 Numbers with a record number of divisors that have the same value of the Euler totient function (A000010).

Original entry on oeis.org

1, 2, 12, 120, 240, 2520, 10920, 21840, 32760, 65520, 600600, 900900, 1801800, 3603600, 4455360, 8910720, 17821440, 46060560, 69090840, 92121120, 126977760, 138181680, 245044800, 414545040, 490089600, 507911040, 1015822080, 1266665400, 1523733120, 2533330800
Offset: 1

Views

Author

Amiram Eldar, Oct 28 2019

Keywords

Comments

The corresponding record values are 1, 2, 3, 4, 5, 7, 8, 9, 10, 11, 13, 14, 16, 17, 18, 20, 21, 23, 24, 25, 26, 27, 29, 30, 31, 32, 33, 34, 35, 36, ...

Examples

			2 has 2 divisors with the same value of phi: phi(1) = phi(2) = 1.
12 has 3 divisors with the same value of phi: phi(3) = phi(4) = phi(6) = 2.
		

Crossrefs

Programs

  • Mathematica
    f[n_] := Max[Tally[EulerPhi /@ Divisors[n]][[;; , 2]]]; fm = 0; s = {}; Do[f1 = f[n]; If[f1 > fm, fm = f1; AppendTo[s, n]], {n, 1, 10^6}]; s

A373530 Numbers k such that k, k+1 and k+2 all have at least three divisors with the same value of the Euler totient function (A000010).

Original entry on oeis.org

14052608, 83025998, 87703714, 93978520, 117345150, 163338174, 213589088, 218539880, 294321950, 369698434, 401177798, 463425920, 470217824, 497434040, 529524918, 539318438, 554556078, 559474838, 581302358, 584754848, 608842934, 612448640, 617445814, 625591966
Offset: 1

Views

Author

Amiram Eldar, Jun 08 2024

Keywords

Comments

Numbers k such that k, k+1 and k+2 are all in A359565.
There must be 3 or more divisors of k that have the same Euler totient value, and ditto for k+1 and k+2, but those values may differ as among k, k+1, and k+2. - Harvey P. Dale, Sep 01 2024

Crossrefs

Subsequence of A359565 and A373529.

Programs

  • Mathematica
    q[n_] := Max[Tally[EulerPhi[Divisors[n]]][[;; , 2]]] > 2; seq[kmax_] := Module[{s = {}, q1 = 0, q2 = 0, q3}, Do[q3 = q[k]; If[q1 && q2 && q3, AppendTo[s, k-2]]; q1=q2; q2=q3, {k, 3, kmax}]; s]; seq[10^8]
    SequencePosition[Table[If[Max[Tally[EulerPhi[Divisors[n]]][[;;,2]]]>2,1,0],{n,88*10^6}],{1,1,1}] [[;;,1]] (* The program generates the first 3 terms of the sequence. *) (* Harvey P. Dale, Sep 01 2024 *)
  • PARI
    is(k) = vecmax(matreduce(apply(x->eulerphi(x), divisors(k)))[,2]) > 2;
    lista(kmax) = {my(q1 = 0, q2 = 0, q3); for(k = 3, kmax, q3 = is(k); if(q1 && q2 && q3, print1(k-2, ", ")); q1 = q2; q2 = q3);}

A363320 a(n) is the product of the frequencies of the distinct values obtained when the Euler totient function is applied to the divisors of n.

Original entry on oeis.org

1, 2, 1, 2, 1, 4, 1, 2, 1, 4, 1, 6, 1, 4, 1, 2, 1, 8, 1, 4, 1, 4, 1, 12, 1, 4, 1, 4, 1, 16, 1, 2, 1, 4, 1, 12, 1, 4, 1, 6, 1, 16, 1, 4, 1, 4, 1, 24, 1, 8, 1, 4, 1, 16, 1, 4, 1, 4, 1, 54, 1, 4, 2, 2, 1, 16, 1, 4, 1, 16, 1, 24, 1, 4, 1, 4, 1, 16, 1, 12, 1, 4, 1, 36, 1, 4, 1, 4, 1, 64
Offset: 1

Views

Author

Juri-Stepan Gerasimov, May 27 2023

Keywords

Comments

The product of the multiplicities of distinct values of the n-th row of A102190.

Examples

			The divisors of 12 are {1, 2, 3, 4, 6, 12} and their phi values are {1, 1, 2, 2, 2, 4} whose sum is also 12. The set of distinct values are {1, 2, 4} which occur with multiplicities {2, 3, 1} respectively. Therefore, a(12) = 2*3*1 = 6.
		

Crossrefs

Programs

  • Mathematica
    a[n_] := Times @@ Tally[EulerPhi[Divisors[n]]][[;; , 2]]; Array[a, 100] (* Amiram Eldar, May 27 2023 *)
  • PARI
    a(n)=my(f=vector(n)); fordiv(n,d,f[eulerphi(d)]++); vecprod([t | t<-f, t>0]) \\ Andrew Howroyd, May 27 2023
Previous Showing 21-25 of 25 results.