OFFSET
0,2
COMMENTS
a(0) = 1 by convention. We observe that a(n) == 0 (mod 6) for n >= 2.
EXAMPLE
+----+----+---------------------------------+-----------------------------+
| n | k | divisors of k | D(k) mod M |
+----+----+---------------------------------+-----------------------------+
| 0 | 1 | {1} | x x x x x x x |
+----+----+---------------------------------+-----------------------------+
| 1 | 2 | {1,2} | {0,1} mod 2 |
+----+----+---------------------------------+-----------------------------+
| 2 | 6 | {1,2,3,6} | {0,1} mod 2, {0,1,2} mod 3 |
+----+----+---------------------------------+-----------------------------+
| 3 | 12 | {1,2,3,4,6,12} | {0,1} mod 2, {0,1,2} mod 3, |
| | | | {0,1,2,3} mod 4 |
+----+----+---------------------------------+-----------------------------+
| 4 | 84 | {1,2,3,4,6,7,12,14,21,28,42,84} | {0,1} mod 2, {0,1,2} mod 3, |
| | | {0,1,2,3} mod 4, |
| | | | {0,1,2,3,4,5,6} mod 7 |
+----+----+---------------------------------+-----------------------------+
| 5 | 60 | {1,2,3,4,5,6,10,12,15,20,30,60} | {0,1} mod 2, {0,1,2} mod 3, |
| | | | {0,1,2,3} mod 4, |
| | | | {0,1,2,3,4} mod 5, |
| | | | {0,1,2,3,4,5} mod 6 |
+----+----+---------------------------------+-----------------------------+
MAPLE
with(numtheory):
nn:=10^7:
for n from 0 to nn do:
ii:=0:
for k from 1 to nn while(ii=0) do:
x:=divisors(k):n0:=nops(x):ind:=0:
for M from 2 to k do:
lst:={}:
for i from 1 to n0 do:
lst:=lst union {irem(x[i], M)}:
od:
lst1:={}:jj:=0:n1:=nops(lst):
for j from 0 to n1-1 do:
if lst[j+1]=j then jj:=jj+1: elst fi:
od:
if jj=M
then ind:=ind+1:else fi:
od:
if ind=n then ii:=1:printf(`%d %d \n`, n, k):
else fi:
od:
od:
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Dec 28 2024
EXTENSIONS
More terms from Jinyuan Wang, Jan 11 2025
STATUS
approved
