cheng__yu_的博客

二次剩余二次剩余勒让德符号判别式求解流程模板2019牛客第九场 B Quadratic equation二次剩余定义：对于n、p，若存在x，使得x2≡n(mod p)x^2\equiv n(mod\ p)x2≡n(mod p)，则称n为模p意义下的二次剩余。也就是说，如果n是模p的二次剩余。那么n mod p\sqrt n\ mod\ pn​ mod p，可以用二次同余方程的解 xxx 来代替n\sqrt nn​，暴力地找就是 x2%p==n

积性函数线性筛1、莫比乌斯函数 μ\muμ2、欧拉函数 φ\varphiφ3、约数个数函数d4、约数和函数σσσ5、积性函数1、莫比乌斯函数 μ\muμconst int N=6e6;int visit[N+10],prime[N+10],mu[N+10],cnt;void init(){ cnt=0,mu[1]=1; for(int i=2;i<=N;++i) { if(!visit[i]) prime[++cnt]=i,mu[i]=-1; for(int j=1;j<

欧拉降幂欧拉降幂bzoj3884. 上帝与集合的正确用法2019ICPC南昌预选赛网络赛 B. super_logCF906D. Power Tower欧拉降幂定理：当b>φ(p)\varphi(p)φ(p) 时，ab≡ab%φ(p)+φ(p)(mod  p)a^b\equiv a^{b\%\varphi(p)+\varphi(p)} (\mod p)ab≡ab%φ(p)+φ(p)(modp)细节：每次递归过程中维护的是指数，因此需要在 dfs 出口计算当前的指数，在 qpow 的过程中也通过

lucas板子int qpow(ll b,int n,int mod){ int res=1; while(n) { if(n&1) res=1ll*res*b%mod; b=1ll*b*b%mod; n>>=1; } return res;}int fac[maxn];void init(int p){ fac[0]=fac[1]=1; for(int i=2;i&lt

链接题意：给定 n、x，求 ∑a=1n∑b=1ngcd(xa−1,xb−1)\sum_{a=1}^n\sum_{b=1}^ngcd(x^a-1,x^b-1)∑a=1n​∑b=1n​gcd(xa−1,xb−1)思路：结论：(xa−1,xb−1)=x(a,b)−1(x^a-1,x^b-1)=x^{(a,b)}-1(xa−1,xb−1)=x(a,b)−1首先更相减损后，(xa−1,xb−1)=(xb(xa−b−1),xb−1)(x^a-1,x^b-1)=(x^b(x^{a-b}-1),x^b-1)(x.