标签：P4454 int 题解 ll CQOI2018 res include ak BSGS

题面比较长,但其实就是一道板子题

根据题意, \(g^a\equiv A(mod~p),g^b\equiv B(mod~p)\),可以直接通过BSGS求出 \(a,b\) 的值,然后直接用快速幂算出来即可

关于BSGS,可以看看我的博客

代码

#include<iostream>
#include<cmath>
#include<unordered_map>
using namespace std;
typedef long long ll;

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

int BSGS(int a, int b, int p)
{
    if (1 % p == b % p)
        return 0;
    unordered_map<int, int> hash;
    int k = sqrt(p) + 1;
    for (int i = 0, j = b % p; i < k; i++, j = (ll)j * a % p)
        hash[j] = i;
    int ak = ksm(a, k, p);
    for (int i = 1, j = ak; i <= k; i++, j = (ll)j * ak % p)
        if (hash.count(j))
            return i * k - hash[j];
    return -1;
}

int main()
{
    int g, p;
    scanf("%d%d", &g, &p);
    int n;
    scanf("%d", &n);
    while (n--)
    {
        int A, B;
        scanf("%d%d", &A, &B);
        int a = BSGS(g, A, p);
        int b = BSGS(g, B, p);
        printf("%d\n", ksm(g, (ll)a * b % (p - 1), p));
    }
    return 0;
}

标签：P4454,int,题解,ll,CQOI2018,res,include,ak,BSGS
来源： https://www.cnblogs.com/A2484337545/p/14695479.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。