欧拉函数及欧拉打表

欧拉函数的功能:用于求小于n的与n互质数的个数

 

欧拉函数模板:

int eular(int n){
    int ret=1,i;
    for(i=2;i*i<=n;i++){
        if(n%i==0){
            n/=i;
            ret*=i-1;
            while(n%i==0){
                n/=i;
                ret*=i;
            }
        }
    }
    if(n>1)ret*=n-1;
    return ret;
}

欧拉打表:

int p[maxn];
void phi()
{
    me(p,0);
    p[1]=1;
    for(int i=2;i<maxn;i++)
    {
        if(!p[i])
            for(int j=i;j<=maxn;j+=i)
            {
                if(!p[j])p[j]=j;
                p[j]=p[j]/i*(i-1);
            }
    }
}

 

校赛时有一道题,关于欧拉函数的,当时看到心里一阵激动心想昨天刚整理的模板,哈哈,一阵快敲,交上去超时,orz,今天补上欧拉打表的模板和那道题。

 

 

 

黑客大战

时间限制: C/C++ 1000ms; Java 2000ms 内存限制: 65535KB

通过次数: 11 总提交次数: 37

 

问题描述

    gscsdlz作为上兰帝国ACM的背后dalao,管理着上兰帝国ACM实验室的所有服务器,当然也是一名不为人知的传奇黑客,一天他又在测试实验室的OJ的安全性,突然发现有人正在对网站攻击,于是gscsdlz开始疯狂反击,在经历几分钟的交手之后,gscsdlz成功地打败了对方,对方当然不肯服输了,于是从OJ的主页找到了gscsdlz的联系方式,并向他发了一封email:

 虽然你在网络攻防上胜过了我,但在算法上我还是非常自信可以打败你的,当然,如果你能解决下面这道题目,我就认输。

 求一个正整数作为分母的绝对值小于1的最简分数的个数。

 gscsdlz作为常年隐蔽在ACM的dalao,当然会做了,于是他也回了一封email:你的题目太简单了,我不仅可以做出来,我还可以求所有分母小于等于n且绝对值小于的1最简分数的个数。

    这当然涉及到黑客的知识盲区了,为了挽回颜面,他找到了身在上兰帝国的你,希望从内部获得答案。

 

输入描述

只有一个正整数n(1<=1000000)。

输出描述

输出一个数,表示所有分母小于等于n且绝对值小于1的最简分数的个数

样例输入
4
样例输出
10
来源
第三届山西省大学生程序设计大赛
提示

所有满足的最简分数有-3/4,-1/4,-2/3,-1/3,-1/2,1/2,1/3,2/3,1/4,3/4

 

很简单的一道题,很容易想到欧拉,但是1e6的数据范围,不采取点小措施容易超时。

#include <iostream>
#include <cstdio>
#include <cstring>
#define me(x,y) memset(x,y,sizeof(x))
#define sd(x) scanf("%d",&x)
#define ss(x) scanf("%s",x)
#define sf(x) scanf("%f",&x)
#define slf(x) scanf("%lf",&x)
#define slld(x) scanf("%lld",&x)
#define pd(x) printf("%d\n",x)
#define plld(x) printf("%lld\n",x)
#define ps(x) printf("%s\n",x)
#define max(x,y) (x>=y?x:y)
#define min(x,y) (x<y?x:y)
#define sum(x,y) (x+y)
#define INF 0x3f3f3f3f

using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int maxn=1e6+5;

int n;

ull sum[maxn];
int p[maxn];

void phi()
{
    me(p,0);
    p[1]=1;
    for(int i=2;i<maxn;i++)
    {
        if(!p[i])
            for(int j=i;j<=maxn;j+=i)
            {
                if(!p[j])p[j]=j;
                p[j]=p[j]/i*(i-1);
            }
    }
}

void init() {
    sum[1]=0;
    for(int i=2;i<maxn;i++)
        sum[i]=sum[i-1]+p[i];
}

int main() {
    phi();
    init();
    p[1]=0;
    sd(n);
    plld(sum[n]*2);
	return 0;
}

 

  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

水能zai舟

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值