2019hbcpc（第三届河北省大学生程序设计竞赛）

题解链接：http://icebound.cc/2019/05/25/hbcpc2019.html

热身赛 511

Description
给定n,m，请你计算[1,n!]范围内有多少数字能被m整除。由于答案很大，请输出答案对10e9+7取模的结果。

Input
第一行一个整数T，代表接下来有T组数据
接下来每组数据，有两个数字n,m
1<n,m<=1e6

Output
输出答案对10e9+7取模的结果。

Sample input
2
9 3
10 9

Sample output
120960
403200

题意为求 $(n!/m)modM$，可以利用费马小定理将原式化为：
$$(n!\ast m^{M-2})modM$$

再根据乘法的取模公式，进一步将式子化为：
$$(n!modM\ast m^{M-2}modM)modM$$

左边的 $n!modM$ 可以进一步拆分为
$$[nmodM\ast (n-1)modM\ast ......\ast 1modM]modM$$

右边的 $m^{M-2}modM$则可以利用快速幂解出

#include <iostream>
#include <algorithm>
using namespace std;
const int M = 1e9 + 7;
const int N = 1e6 + 5;
typedef long long ll;

ll res[N];

ll mod_pow(ll x, ll n, ll mod)
{
   
	ll res = 1; 
	while (n > 0) {
   
		if (n & 1) res = res * x % mod;
		x = x * x % mod;
		n >>=