看了那么多提示取模的帖子还是错，具体死因 原先((i * (solve(b, i) - solve(a - 1, i)) + mod) % mod + mod) % mod中的i * (solve(b, i) - solve(a - 1, i))有可能会爆炸，所以要改成现加模数在模模数，最后乘i

改后的代码((i * (((solve(b, i) - solve(a - 1, i)) + mod) % mod)) + mod) % mod

还有最好每一个地方能取模的就取模