二次剩余

定义

$x^2\equiv a\pmod{p}$

若有解$a$为二次剩余,否则为二次非剩余

性质

二次剩余 乘 二次剩余 = 二次剩余, 显然

二次剩余 乘 二次非剩余 = 二次非剩余, 显然

二次非剩余 乘 二次非剩余 = 二次剩余, 暂时无法证明,后续证明


$1\cdots p-1$中有$\frac{p-1}{2}$个二次剩余 和 $\frac{p-1}{2}$个二次非剩余

二次剩余 $1^2,2^2,\cdots,(\frac{p-1}{2})^2$, 显然包含所有,且都是二次剩余,且两两不等

延伸

勒让德符号(legendre symbol)

$xy\equiv a\pmod p$

显然 对于 给定x,只有唯一$y\equiv ax^{-1}\pmod{p}$


如果a二次非剩余,所以每个x配对的都不是x自己,所以是p-1个有序的不等的配对$(i,ai^{-1})$,

而 每个(x,y)对应也有一个(y,x), 所以是$\frac{p-1}{2}$对 前后可以颠倒的

把所有(x<y)的xy全部乘起来 $(p-1)!\equiv a^{\frac{p-1}{2}}\pmod{p}$


a是二次剩余, x是解,p-x也是解,显然只有这两个是解

所以$\frac{p-1}{2}-1$对 前后不等可交换的,和$(x,x),(p-x,p-x)$

同样 成对的只选(x<y)的全部乘起来,$(p-1)!\equiv a^{\frac{p-1}{2}-1}x(p-x)\equiv -a^{\frac{p-1}{2}-1}x^2\equiv-a^{\frac{p-1}{2}}\pmod{p}$


由上面的性质 引入 勒让德符号(legendre symbol): $\left(\frac{a}{p}\right)$

$=1$则$a$是$\pmod{p}$二次剩余

$=-1$则$a$是$\pmod{p}$二次非剩余


所以$(p-1)!\equiv -\left(\frac{a}{p}\right)a^{\frac{p-1}{2}}\pmod{p}$,即

推论

威尔逊定理$(p-1)!\equiv -1\pmod{p}$

$\left(\frac{a}{p}\right)\equiv a^{\frac{p-1}{2}}\pmod{p}$, 因此可以快速知道一个数是否是$\pmod{p}$的二次剩余,本章节核心

由此 我们知道

$\left(\frac{a}{p}\right)\left(\frac{b}{p}\right)\equiv a^{\frac{p-1}{2}}b^{\frac{p-1}{2}} \equiv (ab)^{\frac{p-1}{2}}\equiv \left(\frac{ab}{p}\right) \pmod{p}$, 也就是勒让德符号的乘法分配率,也证明了上面的 两个二次非剩余的乘积是二次剩余

$\left(\frac{-1}{p}\right)\equiv (-1)^{\frac{p-1}{2}}\pmod{p} = (-1)^{\frac{p-1}{2}}$

所以$p=4n+3$则$-1$不是二次剩余,如果$p=4n+1$则$-1$是二次剩余


$\gcd(a,b)=1$则$a^2+b^2$的素因子都是$4n+1$型

反证: $p=4n+3,p|a^2+b^2$,

$a^2\not\equiv 0\pmod{p}$

$b^2\not\equiv 0\pmod{p}$

$a^2+b^2\equiv 0\pmod{p}$

$-a^2\equiv b^2 \pmod{p}$

$-1 \equiv b^2(a^{-1})^2 \pmod{p}$

$-1 \equiv (ba^{-1})^2 \pmod{p}$, 这说明 -1是二次剩余

$\left(\frac{-1}{p}\right) = (-1)^{\frac{p-1}{2}} = (-1)^{2n+1} = -1$, 这说明-1是二次非剩余,矛盾

高斯引理

例子: $\left(\frac{2}{p}\right)=2^{\frac{p-1}{2}}\mod p$

$a=1,2,\cdots,\frac{p-1}{2}$

$b=2a=2,4,\cdots,p-1$

$2^{\frac{p-1}{2}}(\frac{p-1}{2})!=\prod b=(\frac{p-1}{2})!(-1)^{\mathrm{count}(\mathrm{even} > \frac{p-1}{2})}$, 对于 $> \frac{p-1}{2}$的用$p-(p-v)$替换

$2^{\frac{p-1}{2}}=(-1)^{\mathrm{count}(\mathrm{even} > \frac{p-1}{2})}$

$p=4n+3$是$\left(\frac{2}{p}\right)=(-1)^{n+1}$

$p=4n+1$是$\left(\frac{2}{p}\right)=(-1)^n$

再细分

$p=8n+1,\left(\frac{2}{p}\right)=1$

$p=8n+3,\left(\frac{2}{p}\right)=-1$

$p=8n+5,\left(\frac{2}{p}\right)=-1$

$p=8n+7,\left(\frac{2}{p}\right)=1$

所以$\left(\frac{2}{p}\right)=(-1)^{\frac{p^2-1}{8}}$


一般化:

$\left(\frac{m}{p}\right)$

$a=1,2,\cdots,\frac{p-1}{2}$

$b=ma=1m,2m,\cdots,\frac{p-1}{2}m$, 两两不同余

$\mu_{m,p}=$ 数列$b$中$\mod p$后$> \frac{p-1}{2}$的个数

$\left(\frac{m}{p}\right)=(-1)^{\mu_{m,p}}$

二次互反律

二次也就是上面二次剩余的二次,也就是legendre符号要判定的内容

$p,q$是两个不同奇素数则$\displaystyle \left(\frac{p}{q}\right)\left(\frac{q}{p}\right)=(-1)^{\frac{p-1}{2}\frac{q-1}{2}}$

好处,计算时,可以辗转 反复颠倒 从而快速算出p,q都很大时的 结果


埃森斯坦的证明法

$\displaystyle \frac{p^2-1}{8}q = \sum_{i=1}^{\frac{p-1}{2}} iq$

$\displaystyle =\sum_{i=1}^{\frac{p-1}{2}}([\frac{iq}{p}]p+r_i)$, 写成余数和商

$\displaystyle =\sum_{i=1}^{\frac{p-1}{2}} [\frac{iq}{p}]p+\sum_{i=1,r_i\le \frac{p-1}{2}}^{\frac{p-1}{2}} r_i+\sum_{i=1,r_i > \frac{p-1}{2}}^{\frac{p-1}{2}} r_i$, 根据$r_i$大小拆分

$\displaystyle =\sum_{i=1}^{\frac{p-1}{2}} [\frac{iq}{p}]p+\sum_{i=1,r_i\le \frac{p-1}{2}}^{\frac{p-1}{2}} r_i+\sum_{i=1,r_i > \frac{p-1}{2}}^{\frac{p-1}{2}} (p-(p-r_i))$, 把$r_i$转化成为$[1,\frac{p-1}{2}]$

$\displaystyle =\sum_{i=1}^{\frac{p-1}{2}} [\frac{iq}{p}]p+\sum_{i=1,r_i\le \frac{p-1}{2}}^{\frac{p-1}{2}} r_i+\sum_{i=1,r_i > \frac{p-1}{2}}^{\frac{p-1}{2}} p-\sum_{i=1,r_i > \frac{p-1}{2}}^{\frac{p-1}{2}} (p-r_i)$, 把$r_i$转化成为$[1,\frac{p-1}{2}]$

$\displaystyle =\sum_{i=1}^{\frac{p-1}{2}} [\frac{iq}{p}]p+\sum_{i=1,r_i\le \frac{p-1}{2}}^{\frac{p-1}{2}} r_i+\sum_{i=1,r_i > \frac{p-1}{2}}^{\frac{p-1}{2}} (p-r_i)+\sum_{i=1,r_i > \frac{p-1}{2}}^{\frac{p-1}{2}} p-2\sum_{i=1,r_i > \frac{p-1}{2}}^{\frac{p-1}{2}} (p-r_i)$, 拆出2倍

$\displaystyle =\sum_{i=1}^{\frac{p-1}{2}} [\frac{iq}{p}]p+\sum_{i=1}^{\frac{p-1}{2}} i+p\mu_{q,p}-2\sum_{i=1,r_i > \frac{p-1}{2}}^{\frac{p-1}{2}} (p-r_i)$, 处理(2+3),4项

$\displaystyle =\sum_{i=1}^{\frac{p-1}{2}} [\frac{iq}{p}]p+\frac{p^2-1}{8}+p\mu_{q,p}-2\sum_{i=1,r_i > \frac{p-1}{2}}^{\frac{p-1}{2}} (p-r_i)$, 处理第2项

$\displaystyle \sum_{i=1}^{\frac{p-1}{2}} [\frac{iq}{p}]p+p\mu_{q,p}\equiv \frac{p^2-1}{8}(q-1) +2\sum_{i=1,r_i > \frac{p-1}{2}}^{\frac{p-1}{2}} (p-r_i) \equiv 0\pmod{2}$, 奇偶性

$\displaystyle \sum_{i=1}^{\frac{p-1}{2}} [\frac{iq}{p}]+\mu_{q,p}\equiv 0\pmod{2}$, 奇偶性

$\displaystyle \left(\frac{q}{p}\right)=(-1)^{\mu_{q,p}}=(-1)^{\sum_{i=1}^{\frac{p-1}{2}} [\frac{iq}{p}]}$ 利用高斯引理

$\displaystyle \left(\frac{q}{p}\right)\left(\frac{p}{q}\right)=(-1)^{\sum_{i=1}^{\frac{p-1}{2}} [\frac{iq}{p}]+\sum_{i=1}^{\frac{q-1}{2}} [\frac{ip}{q}]}$ 对称带入

$\displaystyle \sum_{i=1}^{\frac{p-1}{2}} [\frac{iq}{p}]+\sum_{i=1}^{\frac{q-1}{2}} [\frac{ip}{q}]$ 在二维平面中表示的是 $y=\frac{q}{p}x$ 这条线下方和上方,与$x=1,y=1,x=\frac{p-1}{2},y=\frac{q-1}{2}$ 围成三角形边界和内部的格点之和 =$\frac{p-1}{2}\frac{q-1}{2}$

这里 发现中间用到了奇偶性,毕竟是(-1)的次方,所以会有一点不那么精确得到的感觉

展望

三次互反律,四次互反律,eisenstein互反律,Artin互反律,Langlands纲领

https://zhuanlan.zhihu.com/p/601704520

https://zhuanlan.zhihu.com/p/646169244

lemma 2

n非零整数, p奇素数,$p\not{|}n$,$(x,y)=1$

$p|x^2+ny^2 \leftrightarrow \left(\frac{-n}{p}\right)=1$

证明:

必要性: 这个形式的质因子 满足 -n是p的二次剩余

$x^2+ny^2\equiv 0\pmod{p}$

$p\not{|}y$

$(xy^{-1})^2+n\equiv 0\pmod{p}$

$(xy^{-1})^2\equiv -n\pmod{p}$

充分性: -n是p的二次剩余的p可以表示成$x^2+ny^2$从而多个素因子相乘依然是$x^2+ny^2$

$x^2\equiv -n\pmod{p}$

取$x=x,y=1$即可

延伸 $(x,y)=1,x^2+ny^2$的质因数

$p | x^2+y^2, (x,y)=1 \leftrightarrow \left(\frac{-1}{p}\right)=1 \leftrightarrow p\equiv 1\pmod{4}$

$p | x^2+2y^2, (x,y)=1 \leftrightarrow \left(\frac{-2}{p}\right)=1 \leftrightarrow p\equiv 1\pmod{4}$

$p | x^2+3y^2, (x,y)=1 \leftrightarrow \left(\frac{-3}{p}\right)=1 \leftrightarrow p\equiv 1\pmod{3}$,

$p | x^2+7y^2, (x,y)=1 \leftrightarrow \left(\frac{-7}{p}\right)=1 \leftrightarrow p\equiv 1,2,4\pmod{7}$,

注: 这里上面p可以直接是$n$,因为x取n的倍数,而y不是时,$(x,y)=1$,

Refs

What’s the “best” proof of quadratic reciprocity?

https://atcoder.jp/contests/abc340

G - Leaf Color

n点树,点有颜色,问又多少个 导出子图的“端点”都是同色 mod 998244353

n 2e5

2s

1024mb

我的思路

每个颜色单独考虑

那实际上是选择同色的点,要让任何三点不共线,没想到怎么统计

閱讀全文 »

https://codeforces.com/contest/1930

E. 2..3…4…. Wonderful! Wonderful!(3s,256mb)

a[i]=i

for k = 1..(n-1)/2:

输出 原序列 任意次操作(>=0) 后能得到的不同的序列的方案数 mod 998244353

每次操作选择2k+1个数(不需要连续),删掉被选中的里的左边k个和右边k个

n 1e6

2e3组测试

我的思路

容易想到 操作w次,那么就要删除2kw个

w \in [0,n/k]

所以sum(w) = O(n log n)

所以如果对于给定的(k,w)能够 常数时间或均摊常数时间求出来就好了

先不考虑时间限制

考虑被删除的 [k-1]个随便选,中间的[2kw-(2k-1)]个不能全部连续,[k-1]个随便选

必要性: 显然每个实际方案最后一次删除一定保证了最左k和最右侧k之间有未被删除的

充分性: 倒着恢复,如果满足,则可以让最后一次删除的一侧恢复到所有删除中的中间部分,归纳得证

这样就是容斥了 binom(n,2kw)- 固定中间连续的方案

for 连续的一坨的 起点是i

$\sum_{i=k}^{n-(2kw-(k-1))+1} \binom{i-1}{k-1}\binom{n-(i+(2kw-(2k-1))-1)}{k-1}$

令$a=k-1$ 会发现是$\sum \binom{a+t}{a}\binom{b-t}{a}$

的形式

$=[x^{n-2kw}] (\sum_{i=0}^{\infty} \binom{a+i}{i}x^i)^2$

或者

$=\frac{1}{(a!)^2}[x^{n-2kw}] (\sum_{i=0}^{\infty} \frac{(a+i)!}{i!}x^i)^2$

然后这个的话, 对于给定的a也就是给定的k来说

可以一次性计算,这样子不同的w只需要取值一下就好了

但 如果 for k {ntt, for w} 的话, ntt的部分的总代价还是 (n^2logn)

不知道 这样能变成 什么函数 从而快速求得系数

閱讀全文 »

https://atcoder.jp/contests/abc339

F - Product Equality

n个数字a[n], 问有多少个有序对(i,j,k)满足

A[i]*A[j]=A[k]

n 1000

$a[i] \in [0,10^{10}]$

2s

1024mb

我的思路

光是 枚举 (i,j) 那么也是 n^2了,那么乘法和对比就要足够快, 然而最快的就算ntt,也是 (长度 log(长度)) 会tle

一个想法 是 通过减少比较范围,但如果正好 一半长度 [len/2] 一半长度[len]

第二个想法是 选取一些 质数p,然后通过类似hash的想法来完成

a[i] * a[j] == a[k]

那么 (a[i] % p)*(a[j] % p) % p== a[k] %p

没有任何确定性算法

然后就是 十进制压8位,但这样也只是降低到 int[125] x int[125]

n的个数没有减少

閱讀全文 »

Kirchhoff’s theorem(基尔霍夫定律)

无自环,可重边无向图$G$的生成树个数 = G的Laplacian矩阵(基尔霍夫矩阵)的$n-1$阶余子式的行列式

閱讀全文 »
0%