The Euler function

      Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
      Total Submission(s): 4507    Accepted Submission(s): 1872
 
 
Problem Description
The Euler function phi is an important kind of function in number theory, (n) represents the amount of the numbers which are smaller than n and coprime to n, and this function has a lot of beautiful characteristics. Here comes a very easy question: suppose you are given a, b, try to calculate (a)+ (a+1)+....+ (b)
 

Input

There are several test cases. Each line has two integers a, b (2<a<b<3000000).
 
Output
Output the result of (a)+ (a+1)+....+ (b)
 
Sample Input
3 100
 
Sample Outpu3042
 

题目大意:

给出n,m,你需要求出从n到m所有欧拉函数的值。

思路:

先看代码吧。。。。

代码:

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#define N 3000010
using namespace std;
int n,m;
long long ans,phi[N];
void get_phi()
{
    phi[]=;
    ;i<=N;++i)
    {
        ) phi[i]=i;
        ;
    }
    phi[]+=phi[];
    ;i<=N;i+=)
    {
        if(phi[i]==i)
          for(int j=i;j<=N;j+=i)
            phi[j]=phi[j]/i*(i-);
        phi[i]+=phi[i-],phi[i+]+=phi[i];
    }
}
int main()
{
    get_phi();
    while(~scanf("%d%d",&n,&m))
    {
        ans=phi[m]-phi[n-];
        printf("%lld\n",ans);
    }
    ;
}

哈哈,看不明白对吧

显然,我也看不明白。。。。

好了,大体说说吧。 我们还是用欧拉函数来解这道题(废话(#‵′)凸)

但是我们显然不能直接暴力枚举(为什么??)(T成狗啊!!!!)可能您写的好看说不定就过了是吧(虽然不大可能吧。。。。因为有多组数据啊(无奈   )),反正我是T成狗了。。。。

既然这样不行,那我们就找个方便一点的做法。这是有的大佬就想到了,要不我们打个欧拉函数表吧(说白了就是打个表,用个数组把每个欧拉函数值都存起来)(蒟蒻头一回听说这个东西,长见识),不过好像还是T。。。。因为我们有多组数据,对于每一组数据我们都要o(n)查询

我们再找个方便一点的做法,要不我们用一个前缀和吧。(欸,好像可行。。。)这样我们就只需要处理一次,再对于每一组询问o(1)查询就好了

当然又有人要问了,怎么处理??    这个。。。。蒟蒻表示也不太懂,但我还是尽量说明白吧(若果有明显错误的地方,请不要听我瞎bibi)

首先我们先把数据处理到N(为什么?? 因为这样我们可以只处理一次,减小时间复杂度啊。 可是这样我们不就多处理了很多吗??那时间复杂度不是更高了吗??  ORZ,因为我们是有多组数据啊,我们不知道他最大的m是几啊,省得我们再对于每一组询问都要处理一次了),对于循环到的数,我们先判断她是奇数还是偶数,如果是偶数的话,那么就先把他的欧拉函数值附成i/2(为什么??  你看啊,2的欧拉函数值是不是1,4的欧拉函数值是不是2,6的欧拉函数值是不是2,8的欧拉函数值是不是4.。。。。。。这样是不是就可以推出n(为偶数)的欧拉函数值一定比n/2小  好吧,我们可以这样想,偶数一定与偶数不互质,那么我们就排除了n/2个数)如果是奇数的话我们就把他的欧拉函数值附成i,接下来我们再进行一遍循环,来筛掉与每个数不互质的数的个数。

  for(int i=3;i<=N;i+=2)//这个地方我们只循环奇数,因为对于偶数我们前面已经把与他一定不互质的数(偶数)排除了,那么如果
    {           //这个数还有与她不互质的数,那么这个最大公约数一定是奇数,不然一定早被筛掉了。。。。。
        if(phi[i]==i)//如果这个数没有被处理过
          for(int j=i;j<=N;j+=i)//跟欧拉筛一样是每个数都被他最小的因子所筛去
            phi[j]=phi[j]/i*(i-1);//欧拉函数公式
        phi[i]+=phi[i-1],phi[i+1]+=phi[i];//计算前缀和
    }

应该明白了吧。。。。

HDU——2824 The Euler function的更多相关文章

  1. hdu 2824 The Euler function(欧拉函数)

    题目链接:hdu 2824 The Euler function 题意: 让你求一段区间的欧拉函数值. 题解: 直接上板子. 推导过程: 定义:对于正整数n,φ(n)是小于或等于n的正整数中,与n互质 ...

  2. hdu 2824 The Euler function

    The Euler function Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Other ...

  3. hdu 2824 The Euler function 欧拉函数打表

    The Euler function Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Other ...

  4. HDU 2824 The Euler function --------欧拉模板

    The Euler function Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Other ...

  5. HDU - 2824 The Euler function 欧拉函数筛 模板

    HDU - 2824 题意: 求[a,b]间的欧拉函数和.这道题卡内存,只能开一个数组. 思路: ϕ(n) = n * (p-1)/p * ... 可利用线性筛法求出所有ϕ(n) . #include ...

  6. hdu 2824 The Euler function(欧拉函数)

    如果打表的话会超内存,我想到了一种方法解决这个问题.题目给出的数据时3000000,我将三百万分成300个数据,将整万的数据存储下来,计算的时候,先计算x和y之间整万的数据,然后再计算零散数据. 想法 ...

  7. HDU——T 2824 The Euler function

    http://acm.hdu.edu.cn/showproblem.php?pid=2824 Time Limit: 2000/1000 MS (Java/Others)    Memory Limi ...

  8. 欧拉函数 & 【POJ】2478 Farey Sequence & 【HDU】2824 The Euler function

    http://poj.org/problem?id=2478 http://acm.hdu.edu.cn/showproblem.php?pid=2824 欧拉函数模板裸题,有两种方法求出所有的欧拉函 ...

  9. HDU 2824.The Euler function-筛选法求欧拉函数

    欧拉函数: φ(n)=n*(1-1/p1)(1-1/p2)....(1-1/pk),其中p1.p2…pk为n的所有素因子.比如:φ(12)=12*(1-1/2)(1-1/3)=4.可以用类似求素数的筛 ...

随机推荐

  1. c语言—栈区,堆区,全局区,文字常量区,程序代码区 详解

    转:http://www.cnblogs.com/xiaowenhui/p/4669684.html 一.预备知识—程序的内存分配 一个由C/C++编译的程序占用的内存分为以下几个部分1.栈区(sta ...

  2. 二分图最大匹配(匈牙利算法) POJ 3041 Asteroids

    题目传送门 /* 题意:每次能消灭一行或一列的障碍物,要求最少的次数. 匈牙利算法:把行和列看做两个集合,当有障碍物连接时连一条边,问题转换为最小点覆盖数==二分图最大匹配数 趣味入门:http:// ...

  3. php pdo oracle

    <?php/** * Created by mestars. * User: mestars * Date: 6/13/16 * Time: 10:52 PM */header('Access- ...

  4. 喜欢Swift编程语言的人主要是初学者?

    一早一起来,朋友圈除了被苹果发布会刷屏外,还漫天散布着一条类似的招聘消息:“招聘iOS程序员,要求拥有5年的Swift开发经验,有狼性,待遇月薪20K+,专车接送.” 随后身边的朋友很快就开始调侃:& ...

  5. Sql 存储过程动态添加where条件

    )= '2,3' )= '' ) if(@bussHallId is not null) set @strWhere = @strWhere + ' and bh.ID in ('+@bussHall ...

  6. python加载不了cookirlib模块的问题

    Python 3 改成 http.cookiejar了,所以import cookielib只要改成import http.cookiejar,就可以了.

  7. 对比hive和mysql查询汇总

    由于底层的处理机制大不相同,hive和mysql在查询上还是有较大差异的! 单个表的select操作 最简单的查询 ,字段2 frome 表名 where 字段 [not]in(元素1,元素2): 例 ...

  8. bower——基本使用

    基本概念 bower可以解决项目的静态文件依赖的问题 bower是用nodejs开发的,所以要现状nodejs 安装nodejs应用程序,网上自行下载 检验是否成功安装,打开电脑cmd,执行node ...

  9. python学习笔记(6)——for...in &while

    for x in ...循环就是把每个元素代入变量x,然后执行缩进块的语句. 注意:python中的格式也起到语法规则的作用,缩进块就是一个例 求和问题,比较异同 1/sum= ,,,,,,,,,]: ...

  10. POJ_3041_Asteroids

    参考自: http://user.qzone.qq.com/289065406/blog/1299322465 解题思路: 把方阵看做一个特殊的二分图(以行列分别作为两个顶点集V1.V2,其中| V1 ...