逆元
因為mod不滿足除法,所以我們需要用到逆元來求a/b的模數。逆元就相當與(1/b)mod p
逆元有兩種求法。第一種就是用擴充套件歐幾里德來求解。
我們有同餘方程 ax
≡1(modp)
如果有解,x就算a模p意義下的逆元。
可以變成 ax
−yp=
1 用擴充套件歐幾里德來求解。
線性篩逆元
其實還有線性的方法篩出**質素**p的逆元。
我們設 1−
1≡1(
modp)
設 p=x
∗a+r
,r1<
i有我們知道 p≡
0(modp
)=>x∗
a+r≡
0(modp)
兩邊同時乘以a−
1∗r−
1 得到 x∗
r−1+
a−1≡
0(modp)
整理得 a−
1≡−x
∗r−1
(modp)
a−1≡−⌊p
a⌋∗(
pmoda)
−1(modp)
**
a[1]=1;
for(int i=2;ia[i]=-(p/i)*a[p%i];
Easy Integration(逆元 數學)
原題題面 easy integration 思路 一開始跟隊友小夥伴推了蠻久,大概推出來了,我沒化簡出來,他們化簡成功了。抱大腿 正經公式是乙個叫做貝塔函式的東西。貝塔函式 b x y 01x a 1 1 x b 1 dx b x,y int x 1 x dx b x,y 01 xa 1 1 x b...
模板 乘法逆元
適用 求某乙個數在模意義下的乘法逆元。如果a p互質,那麼有ap a就是p的倍數,所以有ap a modp ap 1 1 modp 所以只要打乙個快速冪就ok了。code include using namespace std int n,p int power int x,int k return...
模板 乘法逆元
這是一道模板題 給定n,p求1 n中所有整數在模p意義下的乘法逆元。輸入格式 一行n,p 輸出格式 n行,第i行表示i在模p意義下的逆元。輸入樣例 1 複製10 13 輸出樣例 1 複製179 108112 534 1 n 3 106,n1 n 3 106,n輸入保證 p p p 為質數。先扔個線性...