数论常用模板部分整理

最新推荐文章于 2019-08-16 21:31:22 发布

等我学会后缀自动机

最新推荐文章于 2019-08-16 21:31:22 发布

阅读量337

点赞数

CC 4.0 BY-SA版权

分类专栏： ——————个人模板整理——————

本文链接：https://blog.youkuaiyun.com/qq_37451344/article/details/82049528

——————个人模板整理—————— 专栏收录该内容

7 篇文章

订阅专栏

博客提供了BSGS算法的相关链接，可通过该链接https://blog.youkuaiyun.com/chenxiaoran666/article/details/83450304进一步了解BSGS算法的详细内容。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#pragma comment(linker, "/STACK:102400000,102400000")
#include<bits/stdc++.h>
using namespace std;

#define debug puts("YES");
#define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++)

#define lrt int l,int r,int rt
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define ll long long
const int  maxn =3e3+5;
const int mod=1e9+7;

ll gcd(ll x,ll y){return y==0?x:gcd(y,x%y);}
ll powmod(ll x,ll y){ll t;for(t=1;y;y>>=1,x=x*x%mod) if(y&1) t=t*x%mod;return t;}

///扩展欧几里得求逆元
void exgcd(ll a,ll b,ll &d,ll &x,ll &y)
{
    if(b==0)
    {
        d=a;
        x=1;
        y=0;
        return ;
    }
    exgcd(b,a%b,d,y,x);
    y-=x*(a/b);
}
ll Inv(ll a,ll n){
    ll d,x,y;
    exgcd(a,n,d,x,y);
    return d == 1 ? (x+n)%n : -1;
}

///阶乘逆元+求阶乘组合数
ll fac[maxn],inv[maxn];
void init()
{
    fac[0]=inv[0]=1;
    for(int i=1;i<maxn;i++)   fac[i]=fac[i-1]*i%mod;
   inv[maxn-1]=powmod(fac[maxn-1],mod-2)%mod;
   for(int i=maxn-2;i>=0;i--) inv[i]=inv[i+1]*(i+1)%mod;
}
inline ll C(ll n,ll m)
{
    if(m>n) return 0;
    return fac[n]*inv[m]%mod*inv[n-m]%mod;
}

///mod为质数的情况下的特殊逆元线性求法
ll inv[maxn];
void get_inv() { inv[1]=1; for(int i=2;i<maxn;i++) inv[i]=1LL*(mod-mod/i)*inv[mod%i]%mod; }

///筛法筛欧拉函数，莫比乌斯函数，质数
int vis[maxn],miu[maxn],phi[maxn];
int prim[maxn],cnt=0;
void sieve()
{
    int N=maxn-1;
    memset(vis,0,sizeof(vis));
    memset(miu,0,sizeof(miu));
    for(int i=2;i<=N;i++;
    {
        if(!vis[i])
        {
            prim[cnt++]=i;
            miu[i]=-1;
            phi[i]=i-1;
        }
        for(int j=0;j<cnt;j++)
        {
            ll k=1LL*i*prim[j];
            if(k>N) break;
            vis[k]=1;
            if(i%prim[j])
            {
                phi[k]=phi[i]*(prim[j]-1);
                miu[k]=-miu[i];
            }
            else
            {
                phi[k]=phi[i]*prim[j];
                break;
            }
        }
    }
}

BSGS算法:https://blog.youkuaiyun.com/chenxiaoran666/article/details/83450304

map<int,int> s;//定义一个map
inline int BSGS(int x,int y,int MOD)//对于一个式子x^t=y(mod MOD)，求出t的值
{
	register int i,t=1,base,Size=ceil(sqrt(MOD));//注意此处要用ceil函数向上取整
	for(i=0;i<=Size;++i) s[1LL*t*y%MOD]=i,base=t,t=1LL*t*x%MOD;//预处理将(x^j)*y的值全部用map存下对应的j，并用base存储下x^Size
	for(t=base,i=1;i<=Size;++i,t=1LL*t*base%MOD)//枚举i，每次将t乘上x^Size 
		if(s[t]) return i*Size-s[t];//找到一个合法的i，则答案就是i*Size-j
	return 0;//无解返回0
}