数论全集

Posted 2022-08-10 shxnb666

欧几里得算法、拓展欧几里得算法

欧几里得算法：$gcd(a,b)=gcd(b,a\% b)$

int gcd(int a,int b) 
  if (b==0) return a;
  return gcd(b,a%b);

快速欧几里得算法（更相减损术）：$gcd(a,b)=gcd(b,a-b)$

#define ctz __builtin_ctzll
LL gcd(LL a,LL b) 
  if(!a) return b;
  if(!b) return a;
  int t=ctz(a|b);
  a>>=ctz(a);
  do 
    b>>=ctz(b) ;
    if(a>b) swap(a,b);
    b-=a;
   while(b);
  return a<<t;

拓展欧几里得算法：解不定方程$ax+by=gcd(a,b)$

算$gcd(a,b)$时，
有$ax+by=gcd(a,b)$ $(1)$
算$gcd(b,a\%b)$时，
有$bx‘+(a\%b)y‘=gcd(b,a\%b)$ $(2)$
设$\lfloor \fracab \rfloor=k$，有$a \% b = a-k\times b$
带入$(2)$，有$ay‘+b(x‘-ky‘)=gcd(a,b)$，
因为对于任意的$a,b$恒成立，所以有
$
\left\\beginmatrix
y=x‘-ky‘
\\ x=y‘
\endmatrix\right.
$
先求出$gcd$，每次求出当前的$x‘,y‘$，然后更新上一步的。

int exgcd(int a,int b,int& x,int& y)
  int ans;
  if (b==0) 
    x=1; y=0;
    return a;
  
  ans=exgcd(b,a%b,x,y);
  int tx=x;
  x=y; y=tx-a/b*y;
  return ans;

整除与互质的性质

互质

当$gcd(a,b)=1$时，称a、b互质（素）
性质：

1、已知$gcd(a,c)=1$，若$a|bc$，则$a|b$；若$a|b$且$c|b$，则$ac|b$
2、p为素数，若$p|ab$，则$p|a$或$p|b$
3、$gcd(a,b)\times lcm(a,b)=ab$
4、$gcd(a,b)=gcd(a,b-ac)=gcd(a-bc,b)$
5、存在整数$x$、$y$，使得$ax+by=(a,b)$
6、$m \times gcd(a,b)=gcd(ma,mb)$
7、若$a|m$，$b|m$，则$lcm(a,b)|m$
8、$m\times lcm(a,b)=lcm(ma,mb)$

整除

设$a,b$为整数，$a≠0$，若有一整数$q$，使得$b=aq$，则称$a$是$b$的因数，$b$为$a$的倍数；并称$a$整除$b$，记为$a|b$；若$a$不能整除$b$，则记为$a|b$。
性质：

1、若$c|b$且$b|a$，则$c|a$
2、若$c|a$且$d|b$，则$cd|ab$
3、若$c|a$且$c|b$，则$c|(ka+nb)$。
4、若$ma|mb$，则$a|b$
5、若$a>0$，$b>0$，且$b|a$，则$b \le a$
6、若$n$为正整数，则$(a-b)|(a^n-b^n)$；若$n$为奇数，则$(a+b)|(a^n+b^n)$；若$n$为偶数，则$(a+b)|(a^n-b^n)$

同余

设$m$是正整数，若$m|(a-b)$，称$a$和$b$对模$m$同余，记作$a≡b \pmod m$
性质：

1、$a≡a \pmod m$
2、若$a≡b \pmod m$，则$b≡a \pmod m$
3、若$a≡b \pmod m$，$b≡c \pmod m$，则$a≡c \pmod m$
4、若$a≡b \pmod m$，$c≡d \pmod m$，则$a±c≡b±d \pmod m$，$ac≡bd \pmod m$
5、若$(m,n)=1$，$a≡b \pmod m$，$a≡b \pmod n$，则$a≡b \pmod mn$
6、若$a≡b \pmod m$，$a^n≡b^n \pmod m$
7、若$n|m$，$a≡b \pmod m$，则$a≡b \pmod n$
8、若$ac≡bc \pmod m$，$(c,m)=d$，则$a≡b \pmod m/d$

 

整数分块

引理：对于$\lfloor \fracni \rfloor$最多只有$2\sqrtn$种取值。

当$i \le sqrtn$时，每个$i$一种，最多只有$\sqrtn$种。
当$i > \sqrtn$时，$\lfloor \fracni \rfloor < \sqrtn$，每个整的取值一种，最多只有$\sqrtn$种。

找出这$2\sqrtn$种取值所代表的区间。

设目前区间的左端点为$l$，我们只需找出$r$是多少即可。
有$\lfloor \fracnl \rfloor = \lfloor \fracnr \rfloor$。
设$\lfloor \fracnl \rfloor = k$，即求$kr \le n$时，最大的$r$。
由“除法的定义”得$r=n/k$。

当对应单变量时：

for (ri l=1,r;l<=n;l=r+1) r=n/(n/l);

注意：不可能$r$在一次循环前小于$n$，但是在循环后大于$n$的情况，因为$n/l$最小为$1$，所以不用和$n$取$min$。

当对应双变量时（多个$2$的常数）

for(ri l=1,r;l<=min(n,m);l=r+1) r=min(n/(n/l),m/(m/l));