Wannafly 模拟赛5 E 思维+数论

题目链接

---

参考博客（感谢）

用自己的语言简单总结一下…加深理解。

题意：
Aria正面临算设课程的考试。
设
$$F(n)=\prod _{i=1}^n\prod _{j=1}^{i}j=\prod _{i=1}^{n}i!$$
对于给定的n,m（其中n为质数），求F($n^m$)中质因子n的出现次数，即求一个最大的非负整数e满足$n^e$整除F($n^m$)。
由于是算设课的考试，答案当然是对$10^9+7$取模的。

思路：
首先此题需要基于一个重要的数论知识点：
对于一个数$n$，所包含的某个因子$m$的次数是：

int calc(int n,int m){
    int ans = 0;
    while(n){
        ans += n/m;
        n = n/m;
    }
    return ans;
}

证明：
令 $k=n/m$
因为：
$n!=1\ast 2\ast 3\ast ...n$
试着将因子m提出，则一定有：(C为一个不含因子m的数)
$n!=(m\ast (2\ast m)\ast (3\ast m)...(k\ast m))\ast C$
$=k!\ast m^k\ast C$

所以$k$一定是最后所求个数的一部分。
另一部分则是$k!$中所包含因子$m$的个数。递归处理即可。
证毕。

---

记$calc(n,m)$为$n!$中所含因子$m$的个数。
因为$n$为质数，则只需要单独考虑$1-n^m$中每个数所含因子$n$个数的和。

$$ans=\sum _{i=1}^{n^m}calc(i,n)$$

由上面的证明的数论知识，我们易得：
对于$i\in [kn,kn+n-1]$ ，$calc(i,n)$相同
故我们考虑将$1-n^m$中的所有数每$n$个分为一组，则一组中所有数对答案的贡献相同。
对于$n^m$则需要特殊考虑，故答案可以化简为：

$$ans=calc(n^m,n)+\sum _{i=1}^{n^m-1}calc(i,n)$$
$$=\sum _{i=0}^{m-1}{n}^i+n\ast \sum _{i=1}^{n^{m-1}-1}calc(i\ast n,n)$$

因为：$calc(i\ast n,n)=i+calc(i,n)$
所以：
$$ans=\frac{n^m-1}{n-1}+n\ast (\sum _{i=1}^{n^{m-1}-1}i+\sum _{i=1}^{n^{m-1}-1}calc(i,n))$$

对于括号里面的第二项，可以视为一个子问题递归处理，即得：
$$ans=\frac{n^m-1}{n-1}+\frac{1}{2}\ast (n^m(n^{m-1}-1)+n^m(n^{m-2}-1)+...+n^m(n-1))$$
$$=\frac{n^m-1}{n-1}+\frac{1}{2}\ast n^m\ast \sum _{i=1}^{m-1}(n^i-1)$$
$$=\frac{n^m-1}{n-1}+\frac{n^m}{2}(\frac{n^m-1}{n-1}-m)$$

快速幂+逆元即可。

代码：

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
using namespace std;
typedef long long ll;
 
const ll mod = 1e9 + 7;
 
ll fast_pow(ll n,ll m){
    ll res = 1;
    while(m){
        if(m&1) res = res*n%mod;
        n = n*n%mod;
        m >>= 1;
    }
    return res;
}
 
int main(){
    ll n,m;
    scanf("%lld%lld",&n,&m);
 
    ll res = (fast_pow(n,m)-1)*fast_pow(n-1,mod-2)%mod;
    ll ans = (res + (fast_pow(n,m)*fast_pow(2,mod-2)%mod*(res-m)%mod))%mod;
    if(ans<0) ans += mod;
    printf("%lld\n",ans);
    return 0;
}