|
|
A360418
|
|
Numbers k such that, in a listing of all congruence classes of positive integers, the k-th congruence class contains k. Here the class r mod m (with r in {1,...,m}) precedes the class a' mod b' (with r' in {1,...,m'}) iff m < m' or r > r'.
|
|
4
|
|
|
1, 2, 3, 5, 13, 17, 20, 25, 41, 48, 53, 61, 85, 95, 102, 113, 145, 158, 167, 181, 221, 237, 248, 265, 313, 332, 345, 365, 421, 443, 458, 481, 545, 570, 587, 613, 685, 713, 732, 761, 841, 872, 893, 925
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
The sequence appears to be the interleaving of the four sequences A080856, A102083, A360416, A360417. This has been verified for values of k up to one million as of February 06 2023.
|
|
LINKS
|
|
|
EXAMPLE
|
The 1st congruence class in the list (with m=1 and r=1) is {1,2,3,...} which contains 1, so 1 is in the sequence. The 2nd congruence class (with m=2 and r=2) is {2,4,6,...} which contains 2, so 2 is in the sequence. The 3rd congruence class (with m=2 and r=1) is {1,3,5,...} which contains 3, so 3 is in the sequence. The 4th congruence class (with m=3 and r=3) is {3,6,9,...} which does not contain 4, so 4 is not in the sequence.
|
|
MATHEMATICA
|
mval[n_] := Floor[Sqrt[2 n] + 1/2]; (* A002024 *)
rval[n_] := (2 - 2 n + Round[Sqrt[2 n]] + Round[Sqrt[2 n]]^2)/2; (* A004736 *)
test[n_] := Mod[n - rval[n], mval[n]] == 0;
Select[Range[10000], test[#] &]
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|