求模幂

**
Problem Description

找出一种有效的方法来计算 a^n( mod m) 的值,其中a, n为非负整数,m为正整数, 且a<1018。

Input Description

可能会有多组,每组输入满足题设的三个整数a, n 和 m.

Output Description

对于每组输入,输出a^n( mod m) 的值。

Sample Input

2 3 5
0 456 36
12345678901234567 123456789012345 9876543210123456

**

#include<stdio.h>
typedef long long ll;
ll multi(ll a,ll b,ll m) //快速幂//
{
    ll exp=a%m,res=0;
    while(b)
    {
        if(b&1)
        {
            res+=exp;
            if(res>=m)res-=m;
        }
        exp<<=1;
        if(exp>m)exp-=m;
        b>>=1;
    }
    return res;
}
ll lg(ll a,ll n,ll m)
{
    if(n==0)return 1;
    ll ans=lg(a,n>>1,m);
    ans=multi(ans,ans,m);
    if(n&1)ans=multi(ans,a,m);
    return ans;
}
int main()
{
    ll a,n,m;
    while(~scanf("%lld %lld %lld",&a,&n,&m))
    {
        printf("%lld\n",lg(a,n,m));
    }
 
    return 0;
}

  • 2
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
好的,我来为您解答。以下是 C 语言的代码实现: ```c #include <stdio.h> // 欧几里得算法 int gcd(int a, int b) { int r; while (b) { r = a % b; a = b; b = r; } return a; } // 扩展欧几里得算法 int exgcd(int a, int b, int *x, int *y) { int d, t; if (b == 0) { *x = 1; *y = 0; return a; } d = exgcd(b, a % b, x, y); t = *x; *x = *y; *y = t - a / b * (*y); return d; } // 平方-乘算法求模 int modpow(int x, int y, int m) { int res = 1; while (y > 0) { if (y & 1) { res = (res * x) % m; } x = (x * x) % m; y >>= 1; } return res; } int main() { int p, q, n, phi, e, d, x, y, m, c, d2; // 输入两个质数 p 和 q printf("请输入两个质数 p 和 q:\n"); scanf("%d%d", &p, &q); // 计算 n 和 phi(n) n = p * q; phi = (p - 1) * (q - 1); // 输入加密密钥 e printf("请输入加密密钥 e(需满足 1 < e < %d 且 gcd(e, %d) = 1):\n", phi, phi); scanf("%d", &e); // 检查 e 是否符合要求 if (e <= 1 || e >= phi || gcd(e, phi) != 1) { printf("加密密钥 e 不符合要求!\n"); return 1; } // 计算解密密钥 d exgcd(e, phi, &x, &y); d = x < 0 ? x + phi : x; // 输出加密密钥和解密密钥 printf("加密密钥 e = %d,解密密钥 d = %d\n", e, d); // 输入要加密的明文 m printf("请输入要加密的明文 m(需满足 0 <= m < %d):\n", n); scanf("%d", &m); // 加密明文 m c = modpow(m, e, n); // 输出密文 c printf("密文 c = %d\n", c); // 解密密文 c d2 = modpow(c, d, n); // 输出解密后的明文 printf("解密后的明文 m' = %d\n", d2); return 0; } ``` 以上代码实现了要求的三个功能: 1. 利用欧几里得算法设定加密密钥 e; 2. 采用扩展欧几里得算法由加密密钥 e,求出解密密钥 d; 3. 利用平方-乘算法进行加密和解密,也即求模。 具体的实现过程在代码注释中有详细说明。希望能对您有所帮助!

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值