數(shù)學(xué)day2-數(shù)論maths in acm程序設(shè)計(jì)競(jìng)賽中的

1、費(fèi)馬小定理與素p互質(zhì)，則aap 1mod1偽素?cái)?shù)：n是費(fèi)馬小定理與素p互質(zhì)，則aap 1mod1偽素?cái)?shù)：n是一個(gè)正整數(shù)，如果存在和n互素的正數(shù)a滿足1modn，n是基于a的偽- Miller-Rabbin測(cè)試：不斷的選取不超過(guò)n-1的不同基b，計(jì)算1mod n是否成立，如果每次- 成立，則n是素輔助以二次探測(cè)算法，可以避免將偽素?cái)?shù)判成如果p是素?cái)?shù)，x21modp解只能是x1或xp- TwiceDetect(longlonga,longlong b,longlongt= TwiceDetect(longlonga,longlong b,longlongt= long longx,while(b&

2、 1)=b=1; y = x = while(t-)erMod(a,b,y= MultiMod(x,x,if(y=1&x!=1&x!=k-1) return 0;x =returny !=1 ?0 :bool MillerRabin(long long longlong bool MillerRabin(long long longlong for (i = 0;i MAX; tmp =rand() % (n -1) + if(TwiceDetect(tmp, n-1,n)!=1) return (i = erMod(a, b, k) =a b%k, b, k) =a * b % P = p1

3、 a1 * p2 P = p1 a1 * p2 a2* * pn 約數(shù)和: 1991 Happy O(n + n /2 + n / 3+ + n / n) HOJ2807Patting給定n(n = 10 5)個(gè)數(shù)，每個(gè)數(shù)范圍1 106，對(duì)于每個(gè)數(shù)ai，輸出除ai外所有能夠整除aiHOJ2807Patting給定n(n = 10 5)個(gè)數(shù)，每個(gè)數(shù)范圍1 106，對(duì)于每個(gè)數(shù)ai，輸出除ai外所有能夠整除ai的數(shù)的for i = 2; i * i N; if for i = 2; i * i N; if (tagi) forj =i;j*j -a m a mod m (a a % m = a a

4、/m * (a + b) % m a % m+ b % m mod (a - b) % m a % m- b % m mod (a *b) % m (a % m) * (b % m) mod(a / b) % m (a % m) / (b % m) mod 如何計(jì)算a / b % 定理：如果(b, m) = 1, 1如何計(jì)算a / b % 定理：如果(b, m) = 1, 1 mod 題目中通常bm且m是一個(gè)素?cái)?shù)，滿足如果b*b -11modm,那么稱(chēng)b1為b在模 定b在模m下存在逆元的充要條件是(b, m) = b -1 b -1 * bb -1 b -1 * b b (m)-1) mod

5、這樣，如果(b, m)互素，求解a / b和a*b -1m就缺點(diǎn)：有很強(qiáng)的限制條件(b, m) = 如何計(jì)算a / b % 如果(b, m) != 1, 設(shè)a / b如何計(jì)算a / b % 如果(b, m) != 1, 設(shè)a / b % m有：a / b = km + 則：a = kbm + 先計(jì)算a % bm = br, 然后再除以b如何計(jì)算a / b % 每個(gè)素因子pi的個(gè)數(shù)ki (ki0)如何計(jì)算a / b % 每個(gè)素因子pi的個(gè)數(shù)ki (ki0)，之后對(duì)每個(gè)素因子分別計(jì)算pi kim，最優(yōu)點(diǎn)：適合計(jì)算a! / b! % 缺點(diǎn)：a、b HOJ 2342、1528、1953 (POJ 10

6、91、POJ 1619、POJ 1845、POJ 2478、 2342、1528、1953 (POJ 1091、POJ 1619、POJ 1845、POJ 2478、 POJ2480、POJ2603、POJ2649、POJ2773、 POJ 2992、POJ 3292C(n, ans = if (k n k = nC(n, ans = if (k n k = n for (i = 1;i = k; ans = ans * (n i + 1) / 溢出，最好用long 一般n i N nnN N nnNjn-i=1ji會(huì)法語(yǔ)的有C人，blah blah，問(wèn)至少會(huì)一種語(yǔ)言會(huì)法語(yǔ)的有C人，blah

7、blah，問(wèn)至少會(huì)一種語(yǔ)言2, 3, 4, 1 2, 4, 3, 12, 3, 4, 1 2, 4, 3, 1 | Fnn pa1 a2 1n pa1 a2 1|Ai 2k| Ai AjAk HOJ給定n個(gè)數(shù)(n10)，求m以?xún)?nèi)有多少個(gè)數(shù)能至m = HOJ給定n個(gè)數(shù)(n10)，求m以?xún)?nèi)有多少個(gè)數(shù)能至m = 10, n= 2:3, 3, 4, 6, 89求解前n項(xiàng)(n 求解前n項(xiàng)(n 求解前n項(xiàng)(n 求解前n項(xiàng)(n 的定義所要求的f(n) = Ak-1 * f(n - 1) + Ak-2 * f(n - 2) + . + A0 * f(n-k)，其中f(0).f(k-1)的初構(gòu)造k * k的矩陣

8、其中定義所要求的f(n) = Ak-1 * f(n - 1) + Ak-2 * f(n - 2) + . + A0 * f(n-k)，其中f(0).f(k-1)的初構(gòu)造k * k的矩陣其中A =(Ak-1 Ak-2 . A1)，I然后構(gòu)造一個(gè)k*1列向量這樣，M*b之后b0的值就是f(k)，以此類(lèi)推， M n * b然后構(gòu)造一個(gè)k*1列向量這樣，M*b之后b0的值就是f(k)，以此類(lèi)推， M n * b之后b0的值就是f(k-1+n)，算法復(fù) 雜度O(k 3 * logn)。求和式：A0 + A1A2對(duì)應(yīng)圖論中的路| A I求和式：A0 + A1A2對(duì)應(yīng)圖論中的路| A I| IHOJ2471

9、 Learning = 109)H:HOJ2471 Learning = 109)H:利用輸入構(gòu)造26 * 26的鄰接矩陣題目所求為A + A2 + + A(L-、在在1. All FromeveryN- one move to a P-Fromeveryitions are P- ition,there1. All FromeveryN- one move to a P-Fromeveryitions are P- ition,thereisition,everymoveistog(x)isthesmallestnon-g(x)isthesmallestnon-egerfoundamongtheSprague-Grundyvaluesofthe followers of x.HOJ 2756Who Will Be TheA vs. B Start HOJ 2756Who