码迷,mamicode.com
首页 > 其他好文 > 详细

poj1845 Sumdiv

时间:2018-12-17 02:18:57      阅读:195      评论:0      收藏:0      [点我收藏+]

标签:first   back   分解   sed   play   src   for   tor   turn   

将a质因数分解 S=p1^c1*b+p2^c2*b+...+pn^cn*b

根据乘法分配率可知答案=(p1^0+p1^1...+p1^c1*b)*(p2^0+p2^1...+p2^c2*b)*...*(pn^0+pn^1...+pn^cn*b)

1.如何快速求等比数列呢?我知道通项公式!S=(p^(c+1)-1)/(p-1)% 9901   除法部分乘逆元就是了pow(p-1,9901-2)

2.我不会逆元!   那就分治!  大概思想就是求一半数,通过乘某个数来求另一半,毕竟等比数列很特殊。

法二:

技术分享图片
#include<cstdio>
#include<iostream>
#include<vector>
using namespace std;
const int MOD=9901;
typedef pair<int,int> P;
vector<P>v;
int power(int a,int b){
    int ret=1;
    for(;b;b>>=1){
        if(b&1)ret=ret*a%MOD;
        a=a*a%MOD;
    }
    return ret;
}
int sum(int p,int c){
    if(c==0)return 1;
    if(c%2)return (1+power(p,(c+1)/2))*sum(p,(c-1)/2)%MOD;
    else return ((1+power(p,c/2+1))*sum(p,(c-1)/2)+power(p,c/2))%MOD;
}
int main(){
    int a,b;
    scanf("%d%d",&a,&b);
    int ta=a;
    for(int i=2;i*i<=ta;++i){
        int cnt=0;
        while(a%i==0)a/=i,++cnt;
        v.push_back(P(i,cnt*b));
    }
    if(a!=1)v.push_back(P(a%MOD,b));
    int ans=1;
    for(int i=0;i<v.size();++i){
        cout<<v[i].first<<" "<<v[i].second<<endl;
        printf("%d\n",sum(v[i].first%MOD,v[i].second)%MOD);
        ans=ans*sum(v[i].first%MOD,v[i].second)%MOD;
    }
    printf("%d",ans);
    return 0;
}
View Code

 

poj1845 Sumdiv

标签:first   back   分解   sed   play   src   for   tor   turn   

原文地址:https://www.cnblogs.com/Dream-Runner/p/10128791.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!