Muscl Euler Two dimensions
標簽: dimensions Muscl Euler Two
上傳時間: 2013-12-07
上傳用戶:363186
Euler函數(shù): m = p1^r1 * p2^r2 * …… * pn^rn ai >= 1 , 1 <= i <= n Euler函數(shù): 定義:phi(m) 表示小于等于m并且與m互質(zhì)的正整數(shù)的個數(shù)。 phi(m) = p1^(r1-1)*(p1-1) * p2^(r2-1)*(p2-1) * …… * pn^(rn-1)*(pn-1) = m*(1 - 1/p1)*(1 - 1/p2)*……*(1 - 1/pn) = p1^(r1-1)*p2^(r2-1)* …… * pn^(rn-1)*phi(p1*p2*……*pn) 定理:若(a , m) = 1 則有 a^phi(m) = 1 (mod m) 即a^phi(m) - 1 整出m 在實際代碼中可以用類似素數(shù)篩法求出 for (i = 1 i < MAXN i++) phi[i] = i for (i = 2 i < MAXN i++) if (phi[i] == i) { for (j = i j < MAXN j += i) { phi[j] /= i phi[j] *= i - 1 } } 容斥原理:定義phi(p) 為比p小的與p互素的數(shù)的個數(shù) 設n的素因子有p1, p2, p3, … pk 包含p1, p2…的個數(shù)為n/p1, n/p2… 包含p1*p2, p2*p3…的個數(shù)為n/(p1*p2)… phi(n) = n - sigm_[i = 1](n/pi) + sigm_[i!=j](n/(pi*pj)) - …… +- n/(p1*p2……pk) = n*(1 - 1/p1)*(1 - 1/p2)*……*(1 - 1/pk)
上傳時間: 2014-01-10
上傳用戶:wkchong
//Euler 函數(shù)前n項和 /* phi(n) 為n的Euler原函數(shù) if( (n/p) % i == 0 ) phi(n)=phi(n/p)*i else phi(n)=phi(n/p)*(i-1) 對于約數(shù):divnum 如果i|pr[j] 那么 divnum[i*pr[j]]=divsum[i]/(e[i]+1)*(e[i]+2) //最小素因子次數(shù)加1 否則 divnum[i*pr[j]]=divnum[i]*divnum[pr[j]] //滿足積性函數(shù)條件 對于素因子的冪次 e[i] 如果i|pr[j] e[i*pr[j]]=e[i]+1 //最小素因子次數(shù)加1 否則 e[i*pr[j]]=1 //pr[j]為1次 對于本題: 1. 篩素數(shù)的時候首先會判斷i是否是素數(shù)。 根據(jù)定義,當 x 是素數(shù)時 phi[x] = x-1 因此這里我們可以直接寫上 phi[i] = i-1 2. 接著我們會看prime[j]是否是i的約數(shù) 如果是,那么根據(jù)上述推導,我們有:phi[ i * prime[j] ] = phi[i] * prime[j] 否則 phi[ i * prime[j] ] = phi[i] * (prime[j]-1) (其實這里prime[j]-1就是phi[prime[j]],利用了歐拉函數(shù)的積性) 經(jīng)過以上改良,在篩完素數(shù)后,我們就計算出了phi[]的所有值。 我們求出phi[]的前綴和 */
上傳時間: 2016-12-31
上傳用戶:gyq
【歐拉算法】 微分方程的本質(zhì)特征是方程中含有導數(shù)項,數(shù)值解法的第一步就是...歐拉(Euler)算法是數(shù)值求解中最基本、最簡單的方法,但其求解精度較低,一般不在...對于常微分方程: dy/dx=f(x,y),x∈[a,b] y(a)=y0 可以將區(qū)
上傳時間: 2014-01-09
上傳用戶:www240697738
Euler公式來源及其應用,和圖論的結(jié)合,pdf格式,
標簽: Euler
上傳時間: 2014-01-03
上傳用戶:1159797854
Euler算法,解決常微分方程的初值問題。
上傳時間: 2013-12-21
上傳用戶:www240697738
用Euler法計算常微分方程數(shù)值解,由文件輸入,以表格形式輸出
上傳時間: 2014-01-23
上傳用戶:sevenbestfei
Matching, Euler tours and the Chinese postman.pdf是Edmonds和Johnson1973年,通過使用匹配理論,首次給出中國郵遞員問題的多項式時間解法,前無古人,后無需來者。
標簽: Matching Chinese Edmonds Johnson
上傳時間: 2017-06-21
上傳用戶:jackgao
matlab code for Bernoulli s chaotic map used for encryption
標簽: encryption Bernoulli for chaotic
上傳時間: 2014-01-19
上傳用戶:yxgi5
Some ideas of Euler Algorithm
標簽: Algorithm Euler ideas Some
上傳時間: 2013-12-23
上傳用戶:問題問題