BZOJ 3239 [BSGS]

最新推荐文章于 2018-03-01 21:42:56 发布

Vectorxj

最新推荐文章于 2018-03-01 21:42:56 发布

阅读量545

点赞数 2

分类专栏： BSGS 数论

本文链接：https://blog.csdn.net/Vectorxj/article/details/77855408

版权

数论同时被 2 个专栏收录

20 篇文章 0 订阅

订阅专栏

BSGS

2 篇文章 0 订阅

订阅专栏

$Description$

给定 $A,B$ ，求最小的 $x$ ，满足 $A^x \equiv B \pmod P$

$Solution$

$BSGS$ 裸题

#include <bits/stdc++.h>
using namespace std;

typedef long long ll;

int A, B, P, Ans;
map<int, int> ap;

inline int Pow(int a, int b) {
    int c = 1;
    while (b) {
        if (b & 1) c = (ll)c * a % P;
        a = (ll)a * a % P; b >>= 1;
    }
    return c;
}
inline int Inv(int a, int P) {
    return Pow(a, P - 2);
}
int Log(int a, int b, int p) {
    int m = ceil(sqrt(p)), ia = Inv(Pow(a, m), P), res;
    res = Inv(b, p);
    for (int i = 0; i < m; i++) {
        if (!ap.count(res)) ap[res] = i;
        res = (ll)res * a % P;
    }
    res = 1;
    for (int i = 0; i <= m; i++) {
        if (ap.count(res)) return ap[res] + i * m;
        res = (ll)res * ia % P;
    }
    return -1;
}

int main(void) {
    while (~scanf("%d%d%d", &P, &A, &B)) {
        ap.clear();
        Ans = Log(A, B, P);
        if (Ans == -1) puts("no solution");
        else printf("%d\n", Ans);
    }
}