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

歡迎訪問 生活随笔!

生活随笔

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

编程问答

莫比乌斯反演模版

發布時間:2023/12/31 编程问答 36 豆豆
生活随笔 收集整理的這篇文章主要介紹了 莫比乌斯反演模版 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
//--莫比烏斯反演函數--// //說明:利用線性素數篩選順便求了個mu //注釋部分為求從區間[1,b]和區間[1,d]中取兩個數,互質對數O(n^0.5) //復雜度:O(n) int mu[N]; //int sum[N];void mobus() {bool mark[N];int prime[N];int pcnt=0;memset(mark,0,sizeof(mark));mu[1] = 1;for(int i=2;i<N;i++){if(mark[i] == 0){prime[pcnt++] = i;mu[i] = -1;}for(int j=0;j<pcnt && i*prime[j]<N;j++){int tmp = i*prime[j];mark[tmp] = 1;if( i%prime[j] == 0 ){mu[tmp] = 0;break;}mu[tmp] = mu[i]*-1;}} // for(int i=1;i<N;i++) // sum[i] += sum[i-1]+mu[i]; }//long long gaobili(int b,int d) //{ // if(b<=0||d<=0) return 0; // int m = min(b,d); // long long ans = 0; // while(m>=1) // { // int tb = b/( b/m +1 )+1; // int td = d/( d/m +1 )+1; // //前進的最大位置 // int tm = max(tb,td); // ans += (long long)(sum[m]-sum[tm-1])*(b/m)*(d/m); // m = tm-1; // } // return ans; //}

?

代碼量超少的求一些特別情況的mobus,復雜度O( nlog(n) )

for (int i = 1; i <= n; i++)g[i] = f[i];for (int i = 1; i <= n; i++)for (int j = i + i; j <= n; j += i)g[j] -= g[i];

?

總結

以上是生活随笔為你收集整理的莫比乌斯反演模版的全部內容,希望文章能夠幫你解決所遇到的問題。

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