【数学一本通第一章】Discrete Logging [POJ 2417]

最新推荐文章于 2021-01-02 20:43:56 发布

blng

最新推荐文章于 2021-01-02 20:43:56 发布

阅读量103

点赞数

分类专栏：数学一本通文章标签：数论

本文链接：https://blog.csdn.net/secretzyl/article/details/93757725

版权

数学一本通专栏收录该内容

15 篇文章 0 订阅

订阅专栏

题目：

题目链接：Discrete Logging
题解：
一个裸的BSGS，很好的练习题，也认识了新的写hash的方法。
不是很会就看一下这个：BSGS：大佬博客

代码：

//#include<bits/stdc++.h>
#include<stdio.h>
#include<cmath>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<map>
#include<set>//最烦这么长的库o(╥﹏╥)o了
#define LL long long
using namespace std;
LL A,B,C,ans,D,y,x,m,t;
class hash 
{
	public:
    hash(){memset(a,0xff,sizeof(a));}
    int locate(LL x)
	{
        LL l=x%MOD;
        while(a[l]!=x&&a[l]!=-1) l=(l+1)%MOD;
        return l;
    }
    void insert(LL x,LL va)
	{
        LL l=locate(x);
        if(a[l]==-1) a[l]=x,v[l]=va;
    }
    LL get(LL x) 
	{
       LL l=locate(x);
        return a[l]==x?v[l]:-1;
    }
    void clear(){memset(a,0xff,sizeof(a));}
	private:
    static const LL MOD=100007;
    LL a[MOD+100],v[MOD+100];
}H;
LL exgcd(LL a,LL b,LL &x ,LL &y)
{
	LL t,s;
	if(!b){x=1,y=0; return a;}
	s=exgcd(b,a%b,x,y);
	t=x,x=y,y=t-a/b*y;
	return s;
}
int main()
{
	while(scanf("%lld%lld%lld",&C,&A,&B)!=EOF)
	{
		H.clear();
		m=ceil(sqrt(double(C))); t=1;
		for(int i=0;i<m;i++) H.insert(t,i),t=t*A%C;
		D=1; ans=-1;
		for(int i=0;i<m;i++)
		{
			exgcd(D,C,x,y);//exgcd求逆元
			x=((x*B)%C+C)%C;///B*x=B*D^(-i*m)
			y=H.get(x);//匹配即可
			if(y!=-1){ans=i*m+y; break;}
			D=(D*t)%C;///D=t^i,(t=A^m)
		}
		if(ans==-1) printf("no solution\n");
		else printf("%lld\n",ans);
	}	
	return 0;	
}