//把n的约数的莫比乌斯函数值使用map的形式返回
map<int,int> moebius(int n){
map<int,int>res;
vector<int>primes;
//枚举n的质因数
for(int i=2;i*i<=n;i++){
if(n%i==0){
primes.push_back(i);
while(n%i==0) n/=i;
}
}
if(n!=1) primes.push_back(n);
int m=primes.size();
for(int i=0;i<(1<<m);i++){
int mu=1,d=1;
for(int j=0;j<m;j++){
if(i>>j & 1){
mu*=-1;
d*=primes[j];
}
}
res[d]=mu;
}
}
const int MOD=10009;
int n;
void solve(){
int res=0;
map<int,int>mu=moebius(n);
for(map<int,int>::iterator it=mu.begin();it!=mu.end();it++){
res+=it->second*mod_pow(26,n/it->first,MOD);
res=(res%MOD+MOD)%MOD;
}
printf("%d\n",res);
}
2019-05-21 莫比乌斯反演函数
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
推荐阅读更多精彩内容
- 题目：http://www.lydsy.com/JudgeOnline/problem.php?id=2005
- Description神犇YY虐完数论后给kAc出了一题:给定求且为质数的有多少对。kAc不想做，并把这道题扔给了...
