BSGS       const int MOD=76543; ll hs[MOD],head[MOD],nxt[MOD],id[MOD],top; void insert(ll x,int y){     ll k=x%MOD;     hs[++top]=x,id[top]=y,nxt[top]=head[k],head[k]=top; } ll find(int x){     ll k=x%MOD;     for(int i=head[k];i!=-1;i=nxt[i])     if(hs[i]==x)return id[i];     return -1; } ll bsgs(l...

整数  #include<cstdio> #include<iostream> #include<algorithm> #include<string> #include<cmath> #include<cstring> using namespace std; const int N=105; int a[N][N]; int x[N]; bool free_x[N]; inline int gcd(int a,int b){return b?gcd(b,a%b):a;} inline int lcm(int a,int...

int phi[maxn+5],prime[maxn+5],cnt; bool notp[maxn+5]; void getphi() {     phi[1]=1,cnt=0;     for(int i=2;i<=maxn;i++)     {         if(!notp[i])         {             prime[++cnt]=i;             phi[i]=i-1;         }         for(int j=1;j<=cnt&&i*prime[j]<=maxn;j++)         {      ...

递推     inv[1]=1;for(int i=2;i<MAX;i++)inv[i]=(p-p/i)*inv[p%i]%p;  2.阶乘逆元  inv[0]=inv[1]=inv_fac[0]=fac[0]=1; rep(i,2,maxn) inv[i]=inv[mod%i]*(mod-mod/i)%mod; rep(i,1,maxn) fac[i]=fac[i-1]*i%mod; rep(i,1,maxn) inv_fac[i]=inv_fac[i-1]*inv[i]%mod;  3.素数分解  #include<bits/stdc++.h> using namespa...

//互质 void ex_gcd(int a, int b, int &x, int &y){     if(!b){         x = 1,y = 0;         return ;     }     ex_gcd(b, a%b, y, x);y -= (a / b) * x; } LL CRT(int remid[], int MOD[]) {     LL M = 1;     LL ans = 0;     for(int i = 0; i < n; i++)M *= MOD[i];     for(int i = 0; i < n; i++) ...