poj-1845 Sumdiv nyoj - 928 小M的因子和

题意：求解A^B的因子和 mod 9901

先求解素因子，然后二分求解等比数列

#include<cstdio>
#include<cmath>
typedef long long LL;
const LL mod = 9901 ;
LL pow(LL a,LL b)
{LL res=1;while(b){if(b&1) res=(res*a)%mod;a=(a*a)%mod;b>>=1;}return res;
}
LL fun(LL a,LL b)
{LL s,t;if(b==0) return 0;if(b==1) return a%mod;s=fun(a,b/2)%mod;if(b&1){t=pow(a,b/2+1)%mod;return (s*(t+1)+t)%mod;}else{t=pow(a,b/2)%mod;return (s*(t+1))%mod;}
}
int main()
{LL a,b,c,i,k;while(~scanf("%lld%lld",&a,&b)){if(a == 0) {printf("0\n");continue;}LL sum=1;for(i=2;i<=sqrt(a);i++){k=0;if(a%i==0){while(a%i==0){a/=i;k++;}sum=sum*(fun(i,k*b)+1)%mod;}}if(a>1) sum=sum*(fun(a,b)+1)%mod;printf("%lld\n",sum);}
}

poj-1845 Sumdiv nyoj - 928 小M的因子和相关推荐

c语言的幂乘积表达式,POJ 1845 Sumdiv [素数分解快速幂取模二分求和等比数列]
大致题意: 求A^B的所有约数(即因子)之和,并对其取模 9901再输出. 解题基础: 1) 整数的唯一分解定理: 任意正整数都有且只有一种方式写出其素因子的乘积表达式. ,其中为素数 2) 约数和 ...
NYOJ 928 小M的因子和（数论）
小M的因子和时间限制:1000 ms | 内存限制:65535 KB 难度:2 描述小M在上课时有些得意忘形,老师想出道题目难住他.小M听说是求因子和,还是非常得意,但是看完题目是求A的B次方 ...
poj 1845 Sumdiv (算数基本定理+逆元）
输入a和b,求a^b的所有因子之和. #include <iostream> #define ll long longusing namespace std;const int mod=9 ...
POJ 1845 Sumdiv 【逆元】
题意:求A^B的所有因子之和很容易知道,先把分解得到,那么得到,那么的所有因子和的表达式如下第一种做法是分治求等比数列的和用递归二分求等比数列1+pi+pi^2+pi^3+...+pi^n: ...
Linux等比数列脚本求和,POJ 1845 (约数和+二分等比数列求和)
题目大意:A^B的所有约数和,mod 9901. 解题思路: ①整数唯一分解定理: 一个整数A一定能被分成:A=(P1^K1)*(P2^K2)*(P3^K3).....*(Pn^Kn)的形式.其中Pn ...
【POJ 1845】 Sumdiv (整数唯分+约数和公式+二分等比数列前n项和+同余)
[POJ 1845] Sumdiv 用的东西挺全最主要通过这个题学了约数和公式跟二分求等比数列前n项和另一种小优化的整数拆分整数的唯一分解定理: 随意正整数都有且仅仅有一种方式写出其素因子的乘 ...
Sumdiv POJ - 1845
Sumdiv POJ - 1845 题意: 求ABA^BAB的所有约数之和mod 9901(1<=A,B<=5e7) 题解: 我们先将A分解质因子,表示为:p1c1∗p2c2∗...... ...
【简●解】POJ 1845 【Sumdiv】
POJ 1845 [Sumdiv] [题目大意] 给定\(A\)和\(B\),求\(A^B\)的所有约数之和,对\(9901\)取模. (对于全部数据,\(0<= A <= B <= ...
POJ 1845 逆元 / 分治
题意传送门 POJ 1845 题解分解质因数 A=p1e1p2e2-pnenA=p_1^{e_1}p_2^{e_2}\dots p_n^{e_n}A=p1e1p2e2-pnen,则约数 ...

poj-1845 Sumdiv nyoj - 928 小M的因子和

poj-1845 Sumdiv nyoj - 928 小M的因子和相关推荐

最新文章

热门文章