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

tyvjP1346 MMT数

时间:2017-11-07 00:10:48      阅读:118      评论:0      收藏:0      [点我收藏+]

标签:个数   using   描述   ade   red   ret   研究   /tmp   main   

tyvjP1346 MMT数

描述

FF博士最近在研究MMT数(莫明堂数-_-)。
如果对于一个数n,存在gcd(n,x)<>1并且n mod x<>0 那么x叫做n的MMT数
显然这样的数可以有无限个。
FF博士现在想知道在所有小于n的正整数里面有多少个n的MMT数

输入格式

仅一行一个数,为n

输出格式

 所有小于n的正整数里面有多少个n的MMT数

测试样例1

输入

10

输出

3

备注

样例解释: 3个数分别是 4 6 8 
gcd(n,x)的意思是求n和x的最大公约数

对于50%的数据 n<=1000000
对于100%的数据n<=maxlongint

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<0||ch>9){if(ch=-)f=-1;ch=getchar();}
    while(ch>=0&&ch<=9){x=x*10+ch-0;ch=getchar();}
    return x*f;
}
int n,phi,tot,ans,num[500010],sum=1;
int gcd(int a,int b){
    if(!b) return a;
    return gcd(b,a%b);
}
int main(){
    phi=n=read(); int tmp=n;
    for(int i=2;i*1ll*i<=tmp;++i){
        if(tmp%i==0){
            ++tot; phi=phi/i*(i-1);
            while(tmp%i==0){ ++num[tot]; tmp/=i;}
        }
        if(tmp==1) break;
    }
    if(tmp>1){phi=phi/tmp*(tmp-1); num[++tot]=1;}
    for(int i=1;i<=tot;++i)
    sum=sum*(num[i]+1);
    printf("%d\n",n-phi-sum+1);
    return 0;
}
    

 

tyvjP1346 MMT数

标签:个数   using   描述   ade   red   ret   研究   /tmp   main   

原文地址:http://www.cnblogs.com/huihao/p/7795579.html

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