Crash的数字表格

Description

今天的数学课上，Crash小朋友学习了最小公倍数(Least Common Multiple)。对于两个正整数a和b，LCM(a, b)表示能同时被a和b整除的最小正整数。例如，LCM(6, 8) = 24。回到家后，Crash还在想着课上学的东西，为了研究最小公倍数，他画了一张NM的表格。每个格子里写了一个数字，其中第i行第j列的那个格子里写着数为LCM(i, j)。一个4 5的表格如下： 1 2 3 4 5 2 2 6 4 10 3 6 3 12 15 4 4 12 4 20 看着这个表格，Crash想到了很多可以思考的问题。不过他最想解决的问题却是一个十分简单的问题：这个表格中所有数的和是多少。当N和M很大时，Crash就束手无策了，因此他找到了聪明的你用程序帮他解决这个问题。由于最终结果可能会很大，Crash只想知道表格里所有数的和mod 20101009的值。

Input

输入的第一行包含两个正整数，分别表示 $N$ 和 $M$。

Output

输出一个正整数，表示表格中所有数的和$mod \ 20101009$的值。

Sample Input

$4 \ 5$

Sample Output

$122$

【数据规模和约定】

100%的数据满足$N, M ≤ 10^7$。

由题可得：我们应该求的是 $\sum_{i =1} ^ n\ \sum_{j=1}^m\ lcm(i,j)$ （不妨设 $n<=m$）

先可以由原式化简：

$ans = \sum_{i=1}^n\ \sum_{j=1}^n\ \frac{ij}{gcd(i,j)}$

$ans = \sum_{d=1}^n\ \sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\ \sum_{j=1}^{\lfloor\frac{m}{d}\rfloor}\ [gcd(i,j)=1]\ ijd$

经过最开始的基本化简以后，既然这道题和 $gcd$ 有许许多多不可描述的关系，我们可以考虑莫比乌斯反演。
既然要用莫比乌斯反演，我们就应该来构造相应的 $f(x)$ 和 $F(x)$

我们设 $f(x, y, k)$ 表示 $\sum_{i=1}^x\ \sum_{j=1}^y\ [gcd[i,j]=k]\ ij$

再设 $F(x,y,t)$ 表示 $\sum_{i=1}^x\ \sum_{j=1}^y\ [t\mid gcd(i,j)]ij$

$\therefore \ ans=\sum_{d=1}^n\ d\ f(\lfloor\frac{n}{d}\rfloor,\lfloor\frac{m}{d}\rfloor,1)$

现在，我们来仔细观察一下 $F(x,y,t)$ 和 $f(x,y,k)$ 两个函数的关系

$F(x,y,t)=\sum_{t\mid d}\ f(x,y,d)$

由莫比乌斯反演后可得:

$f(x,y,d)=\sum_{d\mid t}\ \mu(\frac{t}{d})\ F(x,y,t)$

$\because\ d=1$

$\therefore\ f(x,y,1)=\sum_{t=1}^x\ \mu(t)\ F(x,y,t)$

现在，我们回过头了思考一下我们上面一系列操作的意义：
我们由题目要求推得我们需要求 $f(x,y,1)$
但是我们发现直接求并不好求，所以我们反演以后转化为去求 $F(x,y,t)$ 再进一步求得 $ans$
既然我们转化为 $F(x,y,t)$ ，那么这个函数应该要比较方面我们求值才可以达到我们的要求
所以我们来考虑一下 $F(x,y,t)$ 这个函数
仔细思考后可以发现：

$F(x,y,t)=\sum_{d=1}^x\ d^2\ (\sum_{a=1}^{\lfloor\frac{x}{d}\rfloor}\ a\ \sum_{b=1}^{\lfloor\frac{y}{d}\rfloor}\ b)$

设 $sum(x,y)=\sum_{a=1}^x\ \sum_{b=1}^y\ ab$
显然可以用高斯求和 $O(1)$ 求得

$\therefore \ F(x,y,t)=\sum_{d=1}^x\ d^2\ sum(\lfloor\frac{x}{d}\rfloor,\lfloor\frac{y}{d}\rfloor)$

那么整个分析过程就差不多完成了，综上所述：

$ans=\sum_{d=1}^n\ d\ f(\lfloor\frac{n}{d}\rfloor,\lfloor\frac{m}{d}\rfloor,1)$
我们可以分块处理 $f$ ，此处可以分块复杂度$O(\sqrt{n})$

$f(x,y,1)=\sum_{t=1}^x\ \mu(t)\ F(x,y,t)=\sum_{t=1}^x\ \mu(t)\ t^2\ sum(\lfloor\frac{x}{t}\rfloor,\lfloor\frac{y}{t}\rfloor)$
我们可以 $O(1)$ 预处理 $\mu(t)\ t^2$，计算 $sum$

所以总复杂度 $O(n\sqrt{n})$