【算法4.4】快速幂（完结）

 

一、875 快速幂

活动 - AcWing 

求a^b%p的值

#include <iostream>
#include <cstdio>
using namespace std;

typedef long long LL;

int qmi(int a,int b,int p)
{
    LL res=1;
    while(b)
    {
        if(b&1) res=(LL)res*a%p;
        b>>=1;
        a=(LL)a*a%p;
    }
    return res;
}

int main()
{
    int n;
    scanf("%d",&n);
    while(n--)
    {
        int a,b,p;
        scanf("%d %d %d",&a,&b,&p);
        printf("%d\n",qmi(a,b,p));
    }
    return 0;
}

二、876 快速幂求逆元

活动 - AcWing

#include <iostream>
#include <cstdio>
using namespace std;

typedef long long LL;

int qmi(int a,int b,int p)
{
    LL res=1;
    while(b)
    {
        if(b&1) res=(LL)res*a%p;
        b>>=1;
        a=(LL)a*a%p;
    }
    return res;
}

int main()
{
    int n;
    scanf("%d",&n);
    while(n--)
    {
        int a,p;
        scanf("%d %d",&a,&p);
        int res=qmi(a,p-2,p);
        if(a%p) printf("%d\n",res);
        else puts("impossible");
    }
    return 0;
}