The answer is a=1680.
Since ab≡a(−a)=−a2(moda+b), we have a+b∣ab if and only if a+b∣a2. This shows a+b can be any divisor of a2 greater than a. Also, any such divisor corresponds to a unique positive integer b. Therefore, M(a) is the number of divisors of a2 greater than a. Since the positive divisors of a2 less than a are in one-to-one correspondence to the divisors of a2 greater than a, it remains to find a which has the maximum number of positive divisors.
For any a=p1c1p2c2⋯pscs where p1,p2,…,ps are distinct prime divisors of a, recall that d(a2)=(2c1+1)(2c2+1)⋯(2cs+1). WLOG we may assume c1≥c2≥⋯≥cs. Since 2×3×5×7×11=2310>2013, we have s≤4.
* If s=4, we must have c3=c4=1 since 2232527=6300>2013.
- If c2≥2, we must have c1=c2=2 since 23325⋅7=2520>2013. When a=p12p22p3p4, we have d(a2)=52⋅32=225.
- If c2=1, we have c1≤4 since 253⋅5⋅7=3360>2013. When a=p14p2p3p4, we have d(a2)=9⋅33=243.
• If s=3, we must have c3≤2 since 233353=27000>2013.
– If c3=2, we must have c2=2 since 233352=5400>2013. Then c1≤3 since 243252=3600>2013. When a=p13p22p32, we have d(a2)=7⋅52=175.
– If c3=1, we must have c2≤3 since 24345=6480>2013.
* If c2=3, we have c1=3 since 24335=2160>2013. When a=p13p23p3, we have d(a2)=72⋅3=147.
* If c2=2, we have c1≤5 since 26325=2880>2013. When a=p15p22p3, we have d(a2)=11⋅5⋅3=165.
* If c2=1, we have c1≤7 since 283⋅5=3840>2013. When a=p17p2p3, we have d(a2)=15⋅32=135.
• If s=2, we must have c2≤4 since 2535=7776>2013.
– If c2=4, we have c1=4 since 2534=2592>2013. When a=p14p24, we have d(a2)=92=81.
– If c2=3, we have c1≤6 since 2733=3456>2013. When a=p16p23, we have d(a2)=13⋅7=91.
– If c2=2, we have c1≤7 since 2832=2304>2013. When a=p17p22, we have d(a2)=15⋅5=75.
– If c2=1, we have c1≤9 since 2103=3072>2013. When a=p19p2, we have d(a2)=19⋅3=57.
• If s≤1, we have c1≤10 since 211=2048>2013. When a=p110, we have d(a2)=21.
Therefore, a2 has at most 243 positive divisors. This holds when a=p14p2p3p4 for some distinct primes p1,p2,p3,p4. As 342⋅3⋅5=2430>2013, we must have p1=2. Also, since 243⋅5⋅11=2640>2013, the only possibility is a=243⋅5⋅7=1680.