在 C++ 中如何用时间复杂度低的算法求 $\sum_{i=1}^{n}i^k$ 对 $10^9+7$ 取模的值?

其中 $1 \le n \le 10^9,1 \le k \le 10^6$。