rt,可替换为:
$$ \left(\sum^r_{j=l}\left(\sum^j_{j=1}(\gcd(i, j))^k\right)\right)\bmod(10^9+7) $$
另外,这里的 %\%% 应该是指 mod\text{mod}mod 吧。