AcWing - 97 - 约数之和(分治、因数和)

2021-03-13 15:27

阅读:602

标签:return   air   algorithm   algo   names   print   clu   sync   memset   

题目链接
??我们首先要知道怎么来求A的约数之和。首先,把A分解质因数,可得:\(A = q1^{k1}\times q2^{k2} ... \times \ qn^{kn}\)
然后我们用乘法的分配律可得A的因数之和为\(F(A) = (q1^0 + q1^1 + ... + q1^{k1}) \times (q2^0 + q2^1 + ... + q2^{k2})\ ... \times \ (qn^0 + qn^1 + ... + qn^{kn})\)
??那么对于\(A^B\)来说,\(A^B = q1^{k1\times B}\times q2^{k2\times B} ... \times \ qn^{kn\times B}\)\(F(A^B) = (q1^0 + q1^B + ... + q1^{k1\times B}) \times (q2^0 + q2^B + ... + q2^{k2\times B})\ ... \times \ (qn^0 + qn^B + ... + qn^{kn\times B})\)
所以说,我们只要求出来所有因数的所有幂次之和这道题就能很容易的求出来了\(q^0 + q^B + ... + q^{k\times B}\),但是我们发现\(k\times B\)看起来并不好求,用暴力的方法显然会超时,关于这个式子的求法可以看这里

//https://www.cnblogs.com/shuitiangong/
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define endl '\n'
#define rtl rt b ? a : b)
#define minn(a, b) (a  P;
typedef pair P2;
const double pi = acos(-1.0);
const double eps = 1e-7;
const ll MOD =  9901;
const int INF = 0x3f3f3f3f;
const int _NAN = -0x3f3f3f3f;
const double EULC = 0.5772156649015328;
const int NIL = -1;
template void read(T &x){
    x = 0;char ch = getchar();ll f = 1;
    while(!isdigit(ch)){if(ch == '-')f*=-1;ch=getchar();}
    while(isdigit(ch)){x = x*10+ch-48;ch=getchar();}x*=f;
}
const int maxn = 1e4+10;
P fac[maxn];
ll a, b; int kase;
void solve(int a) { //分解因数
    int t = a; kase = 0;
    for (int i = 2; i*i1) fac[kase++] = P(a, b); //如果a本身就是素数,那么a的因数就只有1和自己
}
ll solve2(ll x, int y) { //快速幂
    ll ans = 1; x %= MOD;
    while(y) {
        if (y&1) ans = ans*x%MOD;
        x = x*x%MOD;
        y >>= 1;
    }
    return ans;
}
ll solve3(ll a, ll b) { //计算每个因数所有幂次的累加和
    if (!b) return 1;
    ll res = 1;
    if (b&1) res = res*(1+solve2(a, b/2+1))%MOD*solve3(a, b/2)%MOD;
    else res = res*(((1+solve2(a, b/2))*solve3(a, b/2-1)%MOD + solve2(a, b))%MOD)%MOD;
    return res;
}
int main(void) {
    while(cin >> a >> b) {
        if (!a) {
            cout 

AcWing - 97 - 约数之和(分治、因数和)

标签:return   air   algorithm   algo   names   print   clu   sync   memset   

原文地址:https://www.cnblogs.com/shuitiangong/p/12538482.html


评论


亲,登录后才可以留言!