Give Candies //欧拉函数模板 快速幂模板 欧拉降幂模板 (ACM-ICPC 2018 焦作赛区网络预赛)

求 2^n n<=10^100000 直接万能欧拉降幂

#include<bits/stdc++.h>
using namespace std;
#define mod 1000000007
#define LL long long
LL Euler(LL nqn){
    LL ret=nqn;
    for(LL i=2;i<=sqrt(nqn);i++)
        if(nqn%i==0){
            ret=ret/i*(i-1);//先进行除法防止溢出(ret=ret*(1-1/p(i)))
            while(nqn%i==0)
                nqn/=i;
        }
    if(nqn>1)
          ret=ret/nqn*(nqn-1);
    return ret;
}

LL quick_mod(LL aqa,LL  nqn,LL mood){
    LL ret=1;aqa=aqa%mood;
    while(nqn){
        if(nqn&1){
            ret=(ret*aqa)%mood;
        }
        aqa=(aqa*aqa)%mood;
        nqn=nqn>>1;
    }
    return ret;
}

LL string_quick_mod(LL aqa,char sqs[],LL mood){ //欧拉万能降幂
    int len=strlen(sqs);
    LL ola=Euler(mood);
    LL nqn=0;
    int i;
    for(i=0;i<len;++i){
        nqn=(nqn*10+sqs[i]-'0');
        if(nqn>=ola) {
            break;
        }
    }
    if(i==len) return quick_mod(aqa,nqn,mood);//第二种情况
    nqn=0;
    for(int i=0;i<len;++i){
        nqn=(nqn*10+sqs[i]-'0')%ola;
    }
    return quick_mod(aqa,nqn+ola,mood);  //第三种情况
}

int main(){
    ios :: sync_with_stdio(false);
    int T;char s[100005];
    cin>>T;
    LL inv=quick_mod(2,mod-2,mod);  //费马小定理
    while(T--){
        cin>>s;
        cout<<string_quick_mod(2,s,mod)*inv%mod<<endl;
    }
    return 0;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值