文章目录
数论之神(待理解)
扩展BSGS
- 首先一般BSGS处理的是a,p互质的情况,因为在这种情况下适用欧拉定理,所以解的范围比较小可以直接分块,但是虽然此时仍然有扩展欧拉定理,但是它不优。
- 我们可以利用裴蜀定理的性质将其转换为a,p互质的情况,那么就是每次提取(a,p),然后有解的必要条件就是(a,p)|b或b=1,x=0,所以发现三个数都是(a,p)的倍数,所以我们可以同时除掉,然后就变为了一个规模更小的问题,迭代不超过 O ( l o g n ) O(logn) O(logn)次必然就有(a,p)=1,此时我们就可以直接使用BSGS处理了
- 注意一些b=1特判,还有迭代过程中出现B=C也需要特判此时解就是k,因为之后BSGS的结果必然大于k,过程中(a,p)不整除b的时候直接输出no solution
#include<bits/stdc++.h>
#define LL long long
#define I inline int
#define V inline void
#define FOR(i,a,b) for(register int i=a,end##i=b;i<=end##i;++i)
#define REP(i,a,b) for(register int i=a,end##i=b;i>=end##i;--i)
#define go(i,a) for(register i=hed[a];i;i=e[i].pre)
using namespace std;
inline int read()
{
char x='\0';
int fh=1,sum=0;
for(x=getchar();x<'0'||x>'9';x=getchar())if(x=='-')fh=-1;
for(;x>='0'&&x<='9';x=getchar())sum=sum*10+x-'0';
return fh*sum;
}
I gcd(int a,int b){return (b==0)?a:gcd(b,a%b);}
I ksm(int a,int b,int mod)
{
int sum=1;
while(b)
{
if(b&1)sum=1LL*sum*a%mod;
b>>=1,a=1LL*a*a%mod;
}
return sum;
}
map<int,int>mp;
I exBSGS(int a,int b,int p)
{
int d=gcd(a,p),c=1,k=0;
if(b==1)return 0;
while(d!=1)
{
if(b%d!=0)return -1;
p/=d,b/=d,c=1LL*c*(a/d)%p,k++;
d=gcd(a,p);
if(b==c)return k;
}
// cout<<"k a b p "<<k<<' '<<a<<' '<<b<<' '<<p<<endl;//
int m=ceil(sqrt(p));
mp.clear();
for(int i=0,t=b;i<m;i++,t=1LL*t*a%p)mp[t]=i;
for(int i=1,t=ksm(a,m,p),s=1LL*c*t%p;i<=m;i++,s=1LL*s*t%p)
if(mp.find(s)!=mp.end())return m*i-mp[s]+k;
return -1;
}
int main()
{
int A=0,B=0,P=0;
while(scanf("%d%d%d",&A,&P,&B)!=EOF)
{
if(!A&&!B&&!P)return 0;
int ans=exBSGS(A%P,B%P,P);
if(ans<0)printf("No Solution\n");
else printf("%d\n",ans);
}
return 0;
}
/*
5 58 33
2 4 3
0 0 0
*/
P3846 [TJOI2007] 可爱的质数/【模板】BSGS
主要是注意几个细节:
- 首先是ceil(sqrt(mod))
- 然后就是b%=mod,这一步很关键
- map清空
- 枚举的时候要在map里存b*a^i的值
- 特判a=0的情况,b=0则return 1否则return -1
- 最后没有找到解要return -1
#4128. Matrix
矩阵BSGS,现在要求的就是对一个矩阵求对数,那么我们的做法还是BSGS,然后分块处理即可,只不过有一些代码细节需要注意。
- 首先就是map的比较运算符需要自己定义,并且还要满足我们的要求,必须满足偏序关系,实际上我们可以把矩阵哈希一下也是可以的。但是这里有一个细节就是两个矩阵相等时必须要return false只有这样才是合法的运算。
- 重载中括号的时候不能加friend
- 出现expected initialzer…的时候就是忘记加分号了。
#3122. [Sdoi2013]随机数生成器
这道题已知 x i = a x i − 1 + b x_i=ax_{i-1}+b xi=axi−1+b给出 x 1 , a , b , x n x_1,a,b,x_{n} x1,a,b,xn求解一个最小的n
首先这是一个线性变换的形式,所以我们可以想到利用矩阵求解,所以就是一个类似于4128的问题了。
另外这道题还有一个更简洁的方法,就是利用数列知识,将其化为等比数列,这样直接BSGS即可。
【BZOJ5104】Fib数列(BSGS,二次剩余)
对于一个给出fn求解n的问题,首先我们利用其通项公式,可以表示成一个 1 5 ( A n − B n ) \frac{1}{\sqrt{5}}(A^n-B^n) 51(An−Bn)其中 A = 1 + 5 2 , B = 1 − 5 2 A=\frac{1+\sqrt{5}}{2},B=\frac{1-\sqrt{5}}{2} A=21+5,B=21−5然后可以发现这个东西实际上 A = − 1 B A=-\frac{1}{B} A=−B1所以改写这个式子 1 5 ( 1 ( − B ) n − B n ) \frac{1}{\sqrt{5}}(\frac{1}{(-B)^n}-B^n) 51((−B)n1−Bn)然后根据n的奇偶性进行分类讨论即可,然后可以得到一个一元二次方程,得出两个解检验即可,所以总共可能会有4个解,我们选择最小的即可。
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
#include<map>
using namespace std;
#define MOD 1000000009
int fpow(int a,int b){int s=1;while(b){if(b&1)s=1ll*s*a%MOD;a=1ll*a*a%MOD;b>>=1;}return s;}
namespace SQRT
{
bool check(int n){return fpow(n,(MOD-1)>>1)==1;}
int w;struct Num{int a,b;};
Num operator*(Num a,Num b){return (Num){(1ll*a.a*b.a+1ll*a.b*b.b%MOD*w)%MOD,(1ll*a.a*b.b+1ll*a.b*b.a)%MOD};}
Num fpow(Num a,int b){Num s=(Num){1,0};while(b){if(b&1)s=s*a;a=a*a;b>>=1;}return s;}
int Sqrt(int n)
{
if(!check(n))return -1;
int a;do a=rand()%MOD;while(check((1ll*a*a+MOD-n)%MOD));
w=(1ll*a*a-n+MOD)%MOD;
return fpow((Num){a,1},(MOD+1)>>1).a;
}
}
namespace BSGS
{
map<int,int> M;
const int m=sqrt(MOD);
int BSGS(int x,int y)
{
M.clear();
for(int i=0,t=1;i<m;++i,t=1ll*t*x%MOD)M[1ll*y*t%MOD]=i;
for(int i=1,p=fpow(x,m),t=p;i<=m;++i,t=1ll*t*p%MOD)
if(M.find(t)!=M.end())
return i*m-M[t];
return 2e9;
}
}
int n,ans=2e9;
int main()
{
int sqrt5=SQRT::Sqrt(5);
scanf("%d",&n);
int A=MOD-1,B=(MOD-1ll*sqrt5*n%MOD)%MOD;
int d1=(1ll*B*B%MOD+MOD-4ll*A%MOD)%MOD;
int d2=(1ll*B*B%MOD+4ll*A%MOD)%MOD;
if(SQRT::check(d1))
{
int vd=SQRT::Sqrt(d1);
int x1=1ll*(0ll+MOD-B+vd)%MOD*fpow(A+A,MOD-2)%MOD;
int x2=1ll*(0ll+MOD-B+MOD-vd)%MOD*fpow(A+A,MOD-2)%MOD;
int beta=1ll*(1+MOD-sqrt5)*fpow(2,MOD-2)%MOD;
int n1=BSGS::BSGS(beta,x1);
int n2=BSGS::BSGS(beta,x2);
if(!(n1&1))ans=min(ans,n1);
if(!(n2&1))ans=min(ans,n2);
}
if(SQRT::check(d2))
{
int vd=SQRT::Sqrt(d2);
int x1=1ll*(0ll+MOD-B+vd)%MOD*fpow(A+A,MOD-2)%MOD;
int x2=1ll*(0ll+MOD-B+MOD-vd)%MOD*fpow(A+A,MOD-2)%MOD;
int beta=1ll*(1+MOD-sqrt5)*fpow(2,MOD-2)%MOD;
int n1=BSGS::BSGS(beta,x1);
int n2=BSGS::BSGS(beta,x2);
if(n1&1)ans=min(ans,n1);
if(n2&1)ans=min(ans,n2);
}
printf("%d\n",ans<MOD?ans:-1);
return 0;
}