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;
}

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

POJ1845-Sumdiv

转载请注明出处:優YoU  http://user.qzone.qq.com/289065406/blog/1309237394   大致题意: 求A^B的所有约数(即因子)之和,并对其取模...
  • lyy289065406
  • lyy289065406
  • 2011年07月31日 16:29
  • 5712

POJ 1845:Sumdiv 快速幂+逆元

Sumdiv Time Limit: 1000MS   Memory Limit: 30000K Total Submissions: 16466   Accepted: 4101 ...
  • u010885899
  • u010885899
  • 2015年09月17日 10:41
  • 1330

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

逆元求分数取模代码: #include #include #include #include using namespace std; const int mod=9901; int pow...
  • a601025382s
  • a601025382s
  • 2013年10月02日 12:10
  • 1468

poj(1845)

大致题意: 求A^B的所有约数(即因子)之和,并对其取模 9901再输出。   解题思路: 要求有较强 数学思维 的题 应用定理主要有三个: 要求有较强 数学思维 的题 应用定理主要有三...
  • u011519618
  • u011519618
  • 2013年08月29日 18:02
  • 347

poj 1845

const LL mod = 9901LL ; const int maxn = 10008 ; bool is[maxn] ; int ps ; int ...
  • u013491262
  • u013491262
  • 2014年12月13日 13:29
  • 345

Poj - 1845

Sumdiv Description Consider two natural numbers A and B. Let S be the sum of all natural divisors ...
  • qq_35608674
  • qq_35608674
  • 2016年07月14日 22:19
  • 41

POJ - 1845

题目链接:https://vjudge.net/problem/POJ-1845题目大意:求A^B的所有因数和解题思路:由约数和定理可以知道,A可以分解成p1^a1*p2^a2*……pk^ak,那么A...
  • Nightmare_ak
  • Nightmare_ak
  • 2017年09月20日 23:15
  • 62

POJ 1845

以前做的忘了保存了,找不到博客,才记得这个大数类的经典没保留下来,重新做一个才得,这样的公式挺难得的,谨记着…… #include using namespace std; long long ...
  • u011466175
  • u011466175
  • 2013年10月10日 13:05
  • 547

POJ 1845 Sumdiv (唯一分解定理+求等比数列前n项和)

题目大意:先求出A^B(0
  • u013961718
  • u013961718
  • 2015年09月09日 23:41
  • 521

算术基本定理(唯一分解定理)

摘自 Loi_Peacefuldogの蒟蒻博客
  • LOI_QER
  • LOI_QER
  • 2016年10月12日 20:08
  • 1430
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:POJ 1845 唯一分解定理
举报原因:
原因补充:

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