日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

欧拉筛法的应用

發布時間:2023/12/3 编程问答 24 豆豆
生活随笔 收集整理的這篇文章主要介紹了 欧拉筛法的应用 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

[數論]-----歐拉篩法的應用

文章目錄

  • 1.求1~n之間的所有質數
  • 2.求1~n之間所有自然數的歐拉函數φ(x)
  • 3.求1~n之間的每個數的因子個數
    • 詳細推導:
    • 代碼:
  • 4.求1~n之間每個數的因數和
    • 詳細的推導:
    • 代碼:
  • 篩法求莫比烏斯函數

1.求1~n之間的所有質數

歐拉篩法核心思想:每個合數只被自己的最小質因子篩一次

2.求1~n之間所有自然數的歐拉函數φ(x)

3.求1~n之間的每個數的因子個數

詳細推導:

代碼:

void Euler(){memset(isprime,1,sizeof isprime); d[1]=1;for(int i=2;i<=n;i++){if(isprime[i]){prime[++tot]=i;d[i]=2;num[i]=1;}for(int j=1;j<=tot&&i*prime[j]<=n;j++){isprime[i*prime[j]]=0;if(i%prime[j]==0){d[i*prime[j]]=d[i]/(num[i]+1)*(num[i]+2);num[i*prime[j]]=num[i]+1; break;}else{d[i*prime[j]]=d[i]*2;num[i*prime[j]]=1;}}} }

4.求1~n之間每個數的因數和

詳細的推導:

代碼:

O(n)

void Euler() {ans[1]=1;for(LL i=2;i<=r;i++){if(!vis[i]){prime[++cnt]=i;s[i]=psum[i]=i+1;//i是質數 vis[i]=1;}for(LL j=1;j<=cnt&&prime[j]<=r/i;j++){vis[i*prime[j]]=1;if(i%prime[j]==0){psum[i*prime[j]]=psum[i]*prime[j]+1;s[i*prime[j]]=s[i]/psum[i]*psum[i*prime[j]];break;}else{psum[i*prime[j]]=prime[j]+1;s[i*prime[j]]=s[i]*s[prime[j]];}}} }

篩法求莫比烏斯函數

總結

以上是生活随笔為你收集整理的欧拉筛法的应用的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。