关闭

poj 1845 Sumdiv(同余模公式)

169人阅读 评论(0) 收藏 举报
分类:

http://poj.org/problem?id=1845

2^3 = 8. 
The natural divisors of 8 are: 1,2,4,8. Their sum is 15. 
15 modulo 9901 is 15 (that should be output). 

题意:求A^B,由于A和B很大,所以这里用到同余模公式

思路:根据唯一分解定理A = p1^a1*p2^a2.....*pn^an,

则A^B=(p1^a1*p2^a2.....*pn^an)^B = p1^(a1*B)*p2^(a2*B).....*pn^(an*B)

所有正约数和S = (1+p1+p1^2+...p1^(a1*B)) + (1+p2+p2^2+...p2^(a2*B)) + ... (1+pn+pn^2+...pn^(an*B))

可见ss = (1+p+p^2+...p^n)是等比数列的和,递归二分求该等比数列的和

当n为奇数时,一共有偶数项

ss = (1+p^(n/2+1)) + p* (1+p^(n/2+1)) + ... p^(n/2) *  (1+p^(n/2+1));

    =  (1+p+p^2+...p^(n/2)) * (1+p^(n/2+1))

当n为偶数时,一共有奇数项

ss = (1+p^(n/2)+1) * p * (1+p^(n/2)+1) * ... p^(n/2-1) * (1+p^(n/2)+1) + p^(n/2);

    = (1+p+p^2+...p^(n/2-1)) * (1+p^(n/2+1)) + (p^(n/2))

同余模公式:

(a+b)%mod = (a%mod+b%mod)%mod; 

(a*b)%mod = (a%mod*b%mod)%mod; 

根据同余模公式即可求的结果

#include <iostream>
#include <queue>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cstdlib>
#include <limits>
#include <stack>
#include <vector>
#include <map>

using namespace std;

#define N 1002000
#define INF 0xfffffff
#define PI acos (-1.0)
#define EPS 1e-8
#define met(a, b) memset (a, b, sizeof (a))

typedef long long LL;

LL a;
int isprim[N] = {1, 1}, prim[N], cnt = 0, t, b;

struct node
{
    LL p, a;
} stu[N];

void Init ()
{///素数打表
    for (int i=2; i<1001000; i++)
    {
        if (!isprim[i])
        {
            prim[cnt++] = i;

            for (int j=i+i; j<1001000; j+=i)
                isprim[j] = 1;
        }
    }
}

void Fenjie ()
{///唯一分解定理
    t = 0;

    for (int i=0; i<cnt; i++)
    {
        int k = 0;
        if (a%prim[i]==0)
        {
            while (a%prim[i]==0)
            {
                k++;
                a /= prim[i];
            }

            stu[t].a = k;
            stu[t++].p = prim[i]%9901;
        }
        if (a==1) break;
    }

    if (a!=1)
        stu[t].p = a%9901, stu[t++].a = 1;
}

LL Quick_Pow (int m, int n)
{///快速幂
    LL temp = 1;
    while (n)
    {
        if (n&1)
            temp = temp * m % 9901;
        n >>= 1;
        m = m * m % 9901;
    }
    return temp;
}


LL sum (int p, LL n)
{///等比数列递归二分求和
    if (n==0) return 1;

    if (n%2) return (sum (p, n/2) * (1+Quick_Pow(p, n/2+1)))%9901;
    else return (sum (p, n/2-1) * (1+Quick_Pow(p, n/2+1)) + Quick_Pow(p, n/2))%9901;
}

int main ()
{
    Init();

    while (scanf ("%I64d %d", &a, &b) != EOF)
    {
        met (stu, 0);

        if (a==1 || !a)
        {
            puts ("1");
            continue;
        }
        
        Fenjie();
        LL ans = 1;

        for (int i=0; i<t; i++)
            ans = (ans * sum(stu[i].p, stu[i].a*b)%9901) % 9901;

        printf ("%I64d\n", ans%9901);
    }
    return 0;
}

1
0
查看评论

POJ 1845:Sumdiv 快速幂+逆元

Sumdiv Time Limit: 1000MS   Memory Limit: 30000K Total Submissions: 16466   Accepted: 4101 Description Consi...
  • u010885899
  • u010885899
  • 2015-09-17 10:41
  • 1332

POJ1845-Sumdiv

转载请注明出处:優YoU  http://user.qzone.qq.com/289065406/blog/1309237394   大致题意: 求A^B的所有约数(即因子)之和,并对其取模 9901再输出。   解题思路: 要求有较强 数学思维 的题 ...
  • lyy289065406
  • lyy289065406
  • 2011-07-31 16:29
  • 5723

POJ 1845 Sumdiv (快速幂+质因数+约数和公式+同余模)

Sumdiv Time Limit: 1000MS   Memory Limit: 30000K Total Submissions: 16109   Accepted: 3992 Description Consider two natural number...
  • Grit_ICPC
  • Grit_ICPC
  • 2015-08-17 20:51
  • 568

poj 1845 Sumdiv 数论--等比数列和(逆元或者递归)

逆元求分数取模代码: #include #include #include #include using namespace std; const int mod=9901; int pow_mod(int a,int b,int n) { a=a%n; int s=1; while...
  • a601025382s
  • a601025382s
  • 2013-10-02 12:10
  • 1478

poj 1845 Sumdiv

点击打开链接poj 1845 思路:数学+二分 分析: 1 题目要求的是A^B的所有因子的和对9901取模 2 先看几个数学定理 1:整数的唯一分解定理(如果A本身就是素数的话,那么本身就是分解式) 任意正整数都有且只有一种方式写出其素因子的乘积表达式。 A = (p1^k1)...
  • cgl1079743846
  • cgl1079743846
  • 2012-11-11 21:38
  • 560

【POJ 1845】Sumdiv

首先可以把A质因数分解,然后用等比数列求和公式求出每个质因数可能的贡献。 现在问题来了,在模9901的情况下如何做除法。 一个显然的想法就是求出逆元,因为9901是一个质数,所以直接乘上逆元即可。 我一开始也是这么想的,但是疯狂WA啊。 后来网上找到了另一种方法:二分求等比数列和。大概的意思...
  • Ripped
  • Ripped
  • 2017-03-14 10:44
  • 99

poj 1845 Sumdiv --- 因数分解

http://blog.csdn.net/lyy289065406/article/details/6648539 讲的非常详细,业界良心博主 #include #include #include #include #include using namespace std; #define in...
  • u010489389
  • u010489389
  • 2014-02-23 19:50
  • 451

POJ 1845 Sumdiv 笔记

自然数A和B,求A ^ B的所有自然除数的总和。 
  • woniupengpeng
  • woniupengpeng
  • 2017-06-21 19:34
  • 53

POJ 1845 Sumdiv (快速幂)

题目链接 : http://poj.org/problem?id=1845 题意 : 求∑(x^n的所有因子) mod 9901; 思路 : 如果 x = p1^a1 * p2^a2 * p3^a3...pn^an, (p1, p2, ...都是x的质因子), 那么所求解S = (1 + ...
  • u010619381
  • u010619381
  • 2014-01-19 20:30
  • 653

POJ 1845 Sumdiv

Sol:求A^B所有约数和% MOD的结果。 根据唯一分解定理将A进行因式分解就ok. 等比数列通过奇偶性的判断处理下就行了。 A^B=p1^(a1*B)*p2^(a2*B)*...*pn^(an*B); #include #include #include #include ...
  • imutzcy
  • imutzcy
  • 2014-01-01 15:04
  • 541
    个人资料
    • 访问:79897次
    • 积分:3561
    • 等级:
    • 排名:第10994名
    • 原创:277篇
    • 转载:6篇
    • 译文:0篇
    • 评论:17条
    最新评论