莫比烏斯反演

適用題目特征

常見和式

\sum_{i=x}^{n}\sum_{j=y}^{m}[\gcd(i,j)=k]\qquad (1\leqslant T,x,y,n,m,k\leqslant 5\times 10^4)

常見和式

\sum_{i=1}^n \operatorname{lcm}(i,n)\quad \text{s.t.}\ 1\leqslant T\leqslant 3\times 10^5,1\leqslant n\leqslant 10^6

常見和式

\sum_{i=1}^n\sum_{j=1}^m\operatorname{lcm}(i,j)\qquad (n,m\leqslant 10^7)

常見和式

\sum_{i=1}^n\sum_{j=1}^md(i\cdot j)\\ \left(d(n)=\sum_{i \mid n}1\right) n,m,T\leq5\times10^4
其中 d(n) 表示 n 的約數(shù)個數(shù)。

常見和式

\sum_{i=1}^n\sum_{j=1}^ni\cdot j\cdot \gcd(i,j)\bmod p
n\leq10^{10},5\times10^8\leq p\leq1.1\times10^9,p \text{是質(zhì)數(shù)}

原理

例題

Luogu P2522
代碼如下

/*
Luogu P2522
*/
#define method_1
#ifdef method_1
/*

*/
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<set>
#include<map>
#include<queue>
#include<stack>
#include<vector>
#include<cstring>
#include<cstdlib>
#include<iomanip>
#include<ctime>
#include<string>
#include<bitset>
#define D(x) cout<<#x<<" = "<<x<<"  "
#define E cout<<endl
#define rep(i,s,t) for(int i=(s);i<=(t);i++)
#define rep0(i,s,t) for(int i=(s);i<(t);i++)
#define rrep(i,t,s) for(int i=(t);i>=(s);i--)
using namespace std;
typedef long long ll;
typedef pair<int,int>pii;
const int maxn=5e4+5;
const int INF=0x3f3f3f3f;
int T,a,b,c,d,k;
int mu[maxn],p[maxn],v[maxn],m;
void pre(){
    mu[1]=1;
    rep(i,2,maxn-5){
        if(!v[i]) p[++m]=i,mu[i]=-1;
        for(int j=1;j<=m&&i*p[j]<=maxn-5;j++){
            v[i*p[j]]=1;
            if(i%p[j]==0) {mu[i*p[j]]=0;break;}
            mu[i*p[j]]=-mu[i];
        }
    }
    rep(i,1,maxn-5) mu[i]+=mu[i-1];
}
int calc(int n,int m){
    int res=0;
    for(int i=1,j;i<=min(n,m);i=j+1){
        j=min(n/(n/i),m/(m/i));
        res+=(mu[j]-mu[i-1])*(n/i)*(m/i);
    }
    return res;
}
void solve(){
    int res=calc(b/k,d/k)-calc(b/k,(c-1)/k)-calc((a-1)/k,d/k)+calc((a-1)/k,(c-1)/k);
    cout<<res<<endl;
}
int main() {
//  ios::sync_with_stdio(false);
//  freopen("Problem b.in","r",stdin);
    cin>>T;
    pre();
    while(T--) cin>>a>>b>>c>>d>>k,solve();
    return 0;
}
#endif
#ifdef method_2
/*

*/

#endif
#ifdef method_3
/*

*/

#endif
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時請結(jié)合常識與多方信息審慎甄別。
平臺聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點(diǎn),簡書系信息發(fā)布平臺,僅提供信息存儲服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

  • LUOGU 3455DescriptionFGD正在破解一段密碼,他需要回答很多類似的問題:對于給定的整數(shù)和,有多...
    蘇子旃閱讀 802評論 0 0
  • 提示:別用莫比烏斯反演公式,會炸的 只需要記?。?[gcd(i,j)=1]=\sum_{d|gcd(i,j)}\m...
    An_Account閱讀 2,041評論 0 0
  • 前言 高級數(shù)據(jù)結(jié)構(gòu)難點(diǎn)很多,而且小編接近一年沒有碰過代碼了,簡書一天能發(fā)布的文章數(shù)目有限??,所以今天決定爆肝一個晚...
    gzr666閱讀 2,793評論 2 4
  • 積性函數(shù) 定義 一個數(shù)論函數(shù)f,,有,那么稱f為積性函數(shù)。一個數(shù)論函數(shù)f,對于,有 ,那么稱f為完全積性函數(shù)其中數(shù)...
    dachengquan閱讀 853評論 0 1
  • 久違的晴天,家長會。 家長大會開好到教室時,離放學(xué)已經(jīng)沒多少時間了。班主任說已經(jīng)安排了三個家長分享經(jīng)驗(yàn)。 放學(xué)鈴聲...
    飄雪兒5閱讀 7,810評論 16 22

友情鏈接更多精彩內(nèi)容