|
|
A336505
|
|
5-practical numbers: numbers m such that the polynomial x^m - 1 has a divisor of every degree <= m in the prime field F_5[x].
|
|
4
|
|
|
1, 2, 3, 4, 5, 6, 8, 10, 12, 15, 16, 18, 20, 24, 25, 30, 32, 35, 36, 39, 40, 42, 44, 45, 48, 50, 52, 54, 55, 56, 60, 62, 64, 65, 66, 70, 72, 75, 78, 80, 84, 88, 90, 93, 96, 100, 104, 105, 108, 110, 112, 117, 120, 124, 125, 126, 128, 130, 132, 135, 140, 144, 150
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
For a rational prime number p, a "p-practical number" is a number m such that the polynomial x^m - 1 has a divisor of every degree <= m in F_p[x], the prime field of order p.
A number m is 5-practical if and only if every number 1 <= k <= m can be written as Sum_{d|m} A007736(d) * n_d, where A007736(d) is the multiplicative order of 5 modulo the largest divisor of d not divisible by 5, and 0 <= n_d <= phi(d)/A007736(d).
The number of terms not exceeding 10^k for k = 1, 2, ... are 7, 46, 286, 2179, 16847, 141446, 1223577, ...
|
|
LINKS
|
|
|
MATHEMATICA
|
rep[v_, c_] := Flatten @ Table[ConstantArray[v[[i]], {c[[i]]}], {i, Length[c]}]; mo[n_, p_] := MultiplicativeOrder[p, n/p^IntegerExponent[n, p]]; ppQ[n_, p_] := Module[{d = Divisors[n]}, m = mo[#, p] & /@ d; ns = EulerPhi[d]/m; r = rep[m, ns]; Min @ Rest @ CoefficientList[Series[Product[1 + x^r[[i]], {i, Length[r]}], {x, 0, n}], x] > 0]; Select[Range[200], ppQ[#, 5] &]
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|