POJ2478 Farey Sequence【快速求欧拉函数】
来源:互联网 发布:知乎 大胸妹子 编辑:程序博客网 时间:2024/05/22 22:45
题目链接:
http://poj.org/problem?id=2478
题目大意:
给你一个数n,对于0 < a < b <= n,求真分数a/b的个数
解题思路:
因为a/b为真分数,所以a和b互质。
求真分数a/b的个数。其实就是求0 < i <= n中,小于i的正整数中,
有多少个与i互质的数。累加起来就是真分数a/b的个数。
其实就是欧拉函数
因为n的规模为10^6,可用快速求欧拉函数的方法求得(类似于筛法求素数)。
根据推论:设P是素数,
若p是x的约数,则E(x*p)=E(x)*p.
若p不是x的约数,则E(x*p)=E(x)*E(p)=E(x)*(p-1).
根据筛法求素数的方法,由E(x)求得E(x*p)。
参考博文:http://www.cppblog.com/RyanWang/archive/2009/07/19/90512.html
AC代码:
#include<stdio.h>int prime[100010],phi[1000010];bool unprime[1000010];__int64 sum[1000010];void Euler(){ int i,j,k = 0; //phi[1] = 1; for(i = 2; i <= 1000000; i++) { if(!unprime[i]) { prime[k++] = i; phi[i] = i-1; } for(j = 0; j < k && prime[j]*i <= 1000000; j++) { unprime[prime[j] *i] = true; if(i % prime[j] != 0) { phi[prime[j]*i] = phi[i]*(prime[j]-1); } else { phi[prime[j]*i] = phi[i]*prime[j]; break; } } }}int main(){ int i,n; Euler(); for(i = 1; i <= 1000000; i++) sum[i] = sum[i-1] + phi[i]; while(~scanf("%d",&n) && n) { printf("%I64d\n",sum[n]); } return 0;}
1 0
- POJ2478 Farey Sequence【快速求欧拉函数】
- Farey Sequence(Poj2478)(快速求欧拉函数)
- POJ2478 Farey Sequence,欧拉函数
- poj2478 Farey Sequence (欧拉函数)
- POJ2478 Farey Sequence(欧拉函数)
- POJ2478 Farey Sequence 数论,欧拉函数
- poj2478 Farey Sequence (欧拉函数)
- POJ2478 Farey Sequence
- poj2478-Farey Sequence
- POJ2478 Farey Sequence
- POJ2478 Farey Sequence
- Farey Sequence(poj2478)
- [Poj2478]Farey Sequence
- poj2478 Farey Sequence
- poj2478 Farey Sequence
- POJ2478 Farey Sequence
- [poj2478] Farey Sequence
- POJ2478 Farey Sequence
- 堆栈
- ZipHelper zip压缩工具包
- 第八周项目3 多分段函数求值
- 【一些题】一些重要函数的实现
- 从零开始学android<Bitmap图形组件.四十七.>
- POJ2478 Farey Sequence【快速求欧拉函数】
- 黑马程序员——高新技术3——反射和JavaBean
- leetcode minimum depth(binary tree)
- hdu 4888 最大流慢板
- leetcode Rotate Image
- DocumentHelper xml工具包
- 程序员练级攻略
- 网络解析1——TCP滑动窗口机制解释1
- 系统原理分析架构-三--代理服务器简介及分类