POJ_2478_Farey Sequence
来源:互联网 发布:ubuntu vi复制命令 编辑:程序博客网 时间:2024/06/11 21:13
Farey Sequence
Time Limit: 1000MS Memory Limit: 65536KTotal Submissions: 11067 Accepted: 4267
Description
The Farey Sequence Fn for any integer n with n >= 2 is the set of irreducible rational numbers a/b with 0 < a < b <= n and gcd(a,b) = 1 arranged in increasing order. The first few are
F2 = {1/2}
F3 = {1/3, 1/2, 2/3}
F4 = {1/4, 1/3, 1/2, 2/3, 3/4}
F5 = {1/5, 1/4, 1/3, 2/5, 1/2, 3/5, 2/3, 3/4, 4/5}
You task is to calculate the number of terms in the Farey sequence Fn.
F2 = {1/2}
F3 = {1/3, 1/2, 2/3}
F4 = {1/4, 1/3, 1/2, 2/3, 3/4}
F5 = {1/5, 1/4, 1/3, 2/5, 1/2, 3/5, 2/3, 3/4, 4/5}
You task is to calculate the number of terms in the Farey sequence Fn.
Input
There are several test cases. Each test case has only one line, which contains a positive integer n (2 <= n <= 106). There are no blank lines between cases. A line with a single 0 terminates the input.
Output
For each test case, you should output one line, which contains N(n) ---- the number of terms in the Farey sequence Fn.
Sample Input
23450
Sample Output
1359
Source
POJ Contest,Author:Mathematica@ZSU
题型:数论
题意:求n以内所有互质的数的对数
分析:
先预处理出10^7以内的所有的欧拉函数值φ( i ),然后累加从2到n的phi即可。(题目中将讲n<10^6,开了123456返了个RE,改成1234567就AC了,无解。。。)
代码:
#include<iostream>#include<cstring>#include<cstdio>#include<cmath>#define LL long longusing namespace std;const int maxn=1234567;int phi[maxn];void getphi(){ memset(phi,0,sizeof(phi)); phi[1]=1; for(int i=2;i<maxn;i++){ if(!phi[i]){ for(int j=i;j<maxn;j+=i){ if(!phi[j]) phi[j]=j; phi[j]=phi[j]/i*(i-1); } } }}int main(){ getphi(); int n; while(scanf("%d",&n)&&n){ LL sum=0; for(int i=2;i<=n;i++){ sum+=(LL)phi[i]; } printf("%I64d\n",sum); } return 0;}
- POJ_2478_Farey Sequence
- sequence
- SEQUENCE
- Sequence
- sequence
- sequence
- Sequence
- sequence
- sequence
- sequence
- Sequence
- Sequence
- Sequence
- Sequence
- sequence
- sequence
- sequence
- Sequence
- SecureCRT 自动登陆设置方法
- 在RedHat5.5安装oracle10g中遇到错误
- MD5密码阴谋
- 黑马程序员---读取文本文档的内容
- QSemaphore
- POJ_2478_Farey Sequence
- hdu2767 Proving Equivalences 强连通(缩点染色)
- hdu 1297 Children’s Queue(高精度加法+情况分析+打表)
- 计算机的开机启动过程
- 开源项目代码阅读小技巧
- NYOJ 22 素数求和问题 2013年8月20日
- poj 3335 Rotating Scoreboard(半平面交)
- _declspec(naked) 使用(裸函数)
- 在c文档中C2143问题出现的一种方式及解决方法