当前位置: 首页 > news >正文

算法基础-快速幂

(a + b) % p = (a % p + b % p) % p

(a - b) % p = (a % p - b % p) % p

(a * b) % p = (a % p * b % p) % p

(a ^ b) % p = ((a % p)^b) % p

快速幂 

3^10
10 = 1010 = 2^2 + 2^3
3^10 = 3 * 2^2 + 3 * 2^3 = 3^4 + 3^8

a = a * a => a   a^2   a^4   a^8
所以当 b & 1 == 1 时, res = res *a;
 

long long qmi(long long a,int b,int p)
{
    long long res=1;
    while(b)
    {
        if(b&1) res = res *a %p;
        b>>=1;
        //更新a,a依次为a^{2^0},a^{2^1},a^{2^2},....,a^{2^logb}
        a=a*a%p;
    }
    return res;
}

快速幂求逆元

求某个数的逆元,实际上就是求一个数x,使得(a * x) % p == 1
即原本的 a%p != 1,使得:ax % p == 1

当b、m互质时,b 的倍数 a 在模m情况下也有相应逆元,否则不存在

public class Main {
    public static void main(String[] args) throws Exception {
        Scanner in = new Scanner(System.in);
        int n = in.nextInt();
        for(int i = 0; i < n; i ++) {
            long a = in.nextLong();
            long p = in.nextLong();
            long res = quick(a, p - 2,p);
            if(a % p == 0)
                System.out.println("impossible");
            else
                System.out.println(res);
        }
    }
    private static long quick(long a, long b, long p) {
        long res = 1;
        while(b != 0) {
            if(b % 2 == 1)
                res = res * a % p;
            b = b >> 1;
            a = a * a % p;
        }
        return res;
    }
}

相关文章:

  • 怎么看待日本IT营业这个岗位呢?
  • Redis常用操作及springboot整合redis
  • 内网nginx配置Https证书
  • Mac 电脑 git credential osxkeychain问题之一
  • mfc140u.dll文件错误的相关修复方法,4种方法修复mfc140u.dll
  • Vim笔记
  • LeetCode 3. 无重复字符的最长子串
  • uni-app环境搭建
  • 关于若依flowable的安装
  • OpenAI发布全新o1 AI模型具备推理能力
  • MariaDB 导入报错: Error at line 1: Unknown command ‘\-‘. 版本冲突的解决方法
  • 什么是AIGC?有哪些免费工具?
  • openvpn两个路由器之间呢,需要同一个网段么
  • 嵌入式硬件基础知识
  • 瑞芯微rv1126 Linux 系统,修改系统时区,包有效方法
  • 用Python爬虫制作一个简易翻译器
  • UDP 53端口主要用于DNS(Domain Name System,域名系统)的域名解析服务
  • PFC理论基础与Matlab仿真模型学习笔记(1)--PFC电路概述
  • ORACLE 导出/导入表空间
  • 2024年9月12日(k8s环境及测试 常用命令)
  • 小耳朵等来了春天:公益义诊筛查专家走进安徽安庆
  • 甘肃发布外卖食品安全违法行为典型案例:一商家用鸭肉冒充牛肉被罚
  • 山西省委常委李金科添新职
  • 新疆交通运输厅厅长西尔艾力·外力履新吐鲁番市市长候选人
  • 体坛联播|巴萨4比3打服皇马,利物浦2比2战平阿森纳
  • 中美经贸高层会谈11日在日内瓦将继续进行