|
|
A229303
|
|
Numbers n such that A031971(2*n) == n (mod 2*n).
|
|
19
|
|
|
1, 2, 4, 5, 7, 8, 11, 13, 14, 16, 17, 19, 22, 23, 25, 26, 28, 29, 31, 32, 34, 35, 37, 38, 41, 43, 44, 46, 47, 49, 52, 53, 56, 58, 59, 61, 62, 64, 65, 67, 68, 71, 73, 74, 76, 77, 79, 82, 83, 85, 86, 88, 89, 91, 92, 94, 95, 97, 98, 101, 103, 104, 106, 107, 109, 112, 113, 115, 116, 118, 119, 121, 122, 124, 125
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
The asymptotic density is in [0.583154, 0.58455].
The numbers k = 1, 2, 6, 42, 1806, 47058, 2214502422, 8490421583559688410706771261086 = A230311 are the only values of k such that the set {n: A031971(k*n) == n (mod k*n)} is nonempty. Its smallest element is n = 1, 1, 1, 1, 1, 5, 5, 39607528021345872635 = A231409. [Comment corrected and expanded by Jonathan Sondow, Dec 10 2013]
Up to (but excluding) the term 68 the exponents of even prime powers with squarefree neighbors. - Juri-Stepan Gerasimov, Apr 30 2016.
|
|
LINKS
|
|
|
MAPLE
|
a:= proc(n) option remember; local m;
for m from 1+`if`(n=1, 0, a(n-1)) do
if (t-> m=(add(k&^t mod t, k=1..t) mod t))(2*m)
then return m fi
od
end:
|
|
MATHEMATICA
|
g[n_] := Mod[Sum[PowerMod[i, n, n], {i, n}], n]; Select[Range[100], g[2*#] == # &]
|
|
PROG
|
(PARI) b(n)=sum(k=1, n, Mod(k, n)^n);
for(n=1, 200, if(b(2*n)==n, print1(n, ", ")));
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|