【洛谷】1593 因子和组合数学

最新推荐文章于 2025-11-16 17:14:54 发布

原创最新推荐文章于 2025-11-16 17:14:54 发布 · 389 阅读

1 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#数学

水题同时被 3 个专栏收录

110 篇文章

订阅专栏

洛谷

22 篇文章

订阅专栏

组合数学

5 篇文章

订阅专栏

本文介绍了一种利用正整数的素因子分解来高效计算其因子和的方法，并提供了一个具体的C++实现示例。通过应用数学公式和快速幂运算，可以解决在特定模数下计算因子和的问题。

题目传送门

好像这题并不是组合数学啊……只是一个套公式的题目罢了。

一个正整数 $n$ 有素因子分解 $=p_1^{a_1}*p_2^{a_2}*…*p_m^{a_m}$ ，则有以下两个定理：

1. n 的 因 子 和 = p a 1 + 1 1 - 1 p 1 - 1 * p a 2 + 1 2 - 1 p 2 - 1 * \dots * p a m + 1 m - 1 p m - 1 = \prod i = 1 m p a i + 1 i - 1 p i - 1

$1.n的因子和=\frac{p_1^{a_1+1}-1}{p_1-1}*\frac{p_2^{a_2+1}-1}{p_2-1}*…*\frac{p_m^{a_m+1}-1}{p_m-1}=\prod_{i=1}^m\frac{p_i^{a_i+1}-1}{p_i-1}$

2. n 的 因 子 个 数 = (a 1 + 1) * (a 2 + 1) * \dots * (a m + 1) = \prod i = 1 m (a i + 1)

$2.n的因子个数=(a_1+1)*(a_2+1)*…*(a_m+1)=\prod_{i=1}^m(a_i+1)$

然后这题就变成了直接套第一条定理的水题，直接对于每个素因子快速幂求贡献即可。

然而这题给出的模数9901是一个质数，可能被给出的包含，需要特判一下。

p.s.这题在poj上也有，题目要为1845，只是因为poj最近在维护所以不能交就去洛谷了。

附上AC代码：

#include <cstdio>
using namespace std;

typedef long long ll;
const int mod=9901;
ll a,b,c,ans;

inline ll ksm(ll x,ll y){
    ll sum=1;
    while (y){
        if (y&1) sum=sum*x%mod;
        x=x*x%mod,y>>=1;
    }
    return sum;
}

int main(void){
    scanf("%d%d",&a,&b),ans=1;
    for (int i=2; i*i<=a; ++i)
        if (a%i==0){
            c=0;while (a%i==0) a/=i,++c;
            if (i==9901) continue;
            ans=ans*(ksm(i,c*b+1)-1)%mod*ksm(i-1,mod-2)%mod;
        }
    if (a>1&&a!=9901) ans=ans*(ksm(a,b+1)-1)%mod*ksm(a-1,mod-2)%mod;
    printf("%lld",ans);
    return 0;
}