POJ 1845 唯一分解定理

原创 2016年08月29日 14:58:08
#include <iostream>
using namespace std;
typedef unsigned long long llong;
const int mod=9901;
int p[4000],len=0,top[32][2],tlen;
bool s[7072]={true,true};
void split(llong n)
{
    tlen=0;
    int i;
    for(i=0;i<len;i++)
        {
            if(n%p[i]==0)
                {
                    top[tlen][0]=0;
                    top[tlen][1]=p[i];
                    while(n%p[i]==0)
                        {
                            top[tlen][0]++;
                            n/=p[i];
                        }
                    tlen++;
                    if(n==1||(n<=7071&&!s[n]))break;
                }
        }
    if(n!=1){top[tlen][0]=1;top[tlen++][1]=n;}
}
llong pow_mod(llong a,llong b,llong c)
{
    if(!b) return 1;
    return pow_mod(a*a%c,b/2,c)*((b&1)?a:1)%c;
}
llong sum(llong a,llong n)
{
   if(n==1) return 1;
   llong k=sum(a,n/2)%mod;
   if(n%2==0) return (k+(pow_mod(a,n/2,mod)*k)%mod)%mod;
   else return (pow_mod(a,n/2,mod)+(k+(pow_mod(a,n/2+1,mod)*k)%mod)%mod)%mod;
}
int main()
{
    llong a,b,ans;
    int i,j;
    for(i=2;i*i<=7071;i++)if(!s[i])for(j=i*i;j<=7071;j+=i)if(!s[j])s[j]=true;
    for(i=2;i<=7071;i++)if(!s[i])p[len++]=i;
    while(cin>>a>>b)
        {
            ans=1;
            split(a);
            for(i=0;i<tlen;i++)
                {
                    ans*=sum(top[i][1],top[i][0]*b+1);
                    ans%=9901;
                }
            cout<<ans<<endl;
        }
    return 0;
}

版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

poj 1845Sumdiv(唯一分解定理&&约数和公式&&二分求等比数列和&&反复平方法计算p^n幂~~~好多定理啊)

Sumdiv Time Limit: 1000MS   Memory Limit: 30000K Total Submissions: 19354   Accept...

poj1845(唯一分解定理,等比数列求和,约数个数公式)

/* translation: 输入两个数A,B。求A^B%9901的值 solution: 唯一分解定理, 首先将a进行质因数分解,存储在factor二维数组当中。其形式可以描述为: p1^...

poj3421(多重集排列,唯一分解定理)

/* translation: 给出一个数,求能够形成的x数链最长是多少?(x数链是指以1开始,以n结尾,其中任意相邻的两个数,前面的数都能整除后面) solution: 多重集排列计数,唯一分解...

poj3904 Sky Code —— 唯一分解定理 + 容斥原理 + 组合

题目链接:http://poj.org/problem?id=3904 Sky Code Time Limit: 1000MS   Memory Limit: 65536K Tot...

poj1845Sumdiv+约数和定理

Description Consider two natural numbers A and B. Let S be the sum of all natural divisors of A^B. ...

poj-1845 Sumdiv (逆元+费马小定理+因子和)

分析:与hdu 1452很类似,就不详细解释了,详细解释; #include #include #include typedef long long LL; using namespace std ;...

POJ 1845 Sumdiv 逆元 费马小定理 Trick

原题见POJ 1845给A,B,求ABA^B的所有因子的和模9901。首先把A分解成素数的幂的形式AB=ps1B1ps2B2...pskBkA^B=p_1^{s_1B}p_2^{s_2B}...p_k...

(POJ - 1845)Sumdiv(分治法/费马小定理+快速幂)

Consider two natural numbers A and B. Let S be the sum of all natural divisors of A^B. Determine S m...

poj1845 质因数分解+等比数列求和

题意:求A^B约数之和。 分析:首先可以讲A分解质因数: A=p1^a1+p2^a2+p3^a3+……+pn^an 则A^B=p1^(a1*B)+p2^(a2*B)……+pn^(...
  • Frods
  • Frods
  • 2016年09月25日 21:55
  • 147

poj 1845 Sumdiv(二分递归求等比数列+素因子分解)

题目:http://poj.org/problem?id=1845 Sumdiv Time Limit: 1000MS   Memory Limit: 30000K Total Sub...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:POJ 1845 唯一分解定理
举报原因:
原因补充:

(最多只允许输入30个字)