码迷,mamicode.com
首页 > 其他好文 > 详细

1240 莫比乌斯函数

时间:2018-11-10 18:05:51      阅读:208      评论:0      收藏:0      [点我收藏+]

标签:action   using   icon   fas   amp   lock   void   clu   线性筛   

1240 莫比乌斯函数

 
莫比乌斯函数,由德国数学家和天文学家莫比乌斯提出。梅滕斯(Mertens)首先使用μ(n)(miu(n))作为莫比乌斯函数的记号。(据说,高斯(Gauss)比莫比乌斯早三十年就曾考虑过这个函数)。
 
具体定义如下:
如果一个数包含平方因子,那么miu(n) = 0。例如:miu(4), miu(12), miu(18) = 0。
如果一个数不包含平方因子,并且有k个不同的质因子,那么miu(n) = (-1)^k。例如:miu(2), miu(3), miu(30) = -1,miu(1), miu(6), miu(10) = 1。
给出一个数n, 计算miu(n)。
 

输入

输入包括一个数n,(2 <= n <= 10^9)

输出

输出miu(n)。

输入样例

5

输出样例

-1


莫比乌斯函数模板题

 1 //莫比乌斯函数 
 2 
 3 #include <bits/stdc++.h>
 4 #define N 1000005
 5 using namespace std;
 6 
 7 bool check[N];
 8 int prime[N];
 9 int mu[N];
10 
11 void moblus(){ //莫比乌斯函数  线性筛
12     mu[1] = 1;
13     int tot = 0;
14     for(int i = 2; i <= N; ++i){
15         if(!check[i]){
16             prime[tot++] = i;
17             mu[i] = -1;
18         }
19         for(int j = 0; j < tot; ++j){
20             if(i*prime[j] > N) break;
21 
22             check[i*prime[j]] = true;
23 
24             if( i%prime[j] ){
25                 mu[i*prime[j]] = -mu[i];
26             }else{
27                 mu[i*prime[j]] = 0;
28                 break;
29             }
30         }
31     }
32 }
33 
34 
35 int miu(int n){ // 单点计算
36     int ans = 0;
37     for(int i = 2; i*i <= n; ++i){
38         if(n%i)
39             continue;
40         int cnt = 0;
41         ans ++;
42         while(n%i == 0){
43             n /= i;
44             cnt ++;
45         }
46         if(cnt >= 2){
47             return 0;
48         }
49     }
50     if(n != 1)
51         ans ++;
52     return ans%2 ? -1 : 1;
53 }
54 
55 int main(){
56     int n;
57     cin>>n;
58     // moblus();
59     // for(int i = 1; i <= n; i++){
60     //     cout<<mu[i]<<" ";
61     // }
62     // cout<<endl;
63     cout<<miu(n)<<endl;
64     return 0;
65 }

 




1240 莫比乌斯函数

标签:action   using   icon   fas   amp   lock   void   clu   线性筛   

原文地址:https://www.cnblogs.com/zllwxm123/p/9939846.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!