题目链接:

Farey Sequence

Time Limit: 1000MS   Memory Limit: 65536K
Total Submissions: 14230   Accepted: 5624

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.

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

2
3
4
5
0

Sample Output

1
3
5
9
题意:问满足a/b with 0 < a < b <= n and gcd(a,b) = 1,的数对有多少个;
思路:dp[i]=dp[i-1]+n的欧拉函数;
AC代码:
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int N=1e6+;
long long dp[N],a[N];
int get_a()
{
memset(a,,sizeof(a));
for(int i=;i<N;i++)
{
if(!a[i])
{
for(int j=i;j<N;j+=i)
{
if(!a[j])a[j]=j;
a[j]=a[j]/i*(i-);
}
}
} }
int fun()
{
get_a();
dp[]=;
dp[]=;
for(int i=;i<N;i++)
{
dp[i]=dp[i-]+a[i];
}
}
int main()
{
int n;
fun();
while()
{
scanf("%d",&n);
if(!n)break;
cout<<dp[n]<<"\n";
}
return ;
}

最新文章

  1. iOS,多媒体,地图相关
  2. 【解决】同一url的http请求所获取的结果总是相同
  3. python基础——实例属性和类属性
  4. If A wants to use B
  5. 关于favicon.ico的使用
  6. 关于iframe调用父页面元素操作
  7. js本地存储解决方案(localStorage与userData)
  8. LED驅動芯片對LED壽命的影響
  9. MapReduce的C#实现及单元测试(试验)
  10. PHP自练项目中个人中心创建,修改,验证(服务器端和客户端验证)
  11. 重识 CSS
  12. JavaWeb(二)会话管理之细说cookie与session
  13. ORM规约变更经典案例---mysql军规
  14. VUE 简单属性操作
  15. 《连载 | 物联网框架ServerSuperIO教程》- 16.集成OPC Server,及使用步骤。附:3.3 发布与版本更新说明。
  16. [Swift]LeetCode713. 乘积小于K的子数组 | Subarray Product Less Than K
  17. JAVA常用加密解密算法Encryption and decryption
  18. WordPress &lt;= 4.6 命令执行漏洞(PHPMailer)复现分析
  19. Linux命令大全总结
  20. es6/ts for in/ for of

热门文章

  1. Lumen开发:结合Redis实现消息队列(3)
  2. Android Studio报Error:Execution failed for task &amp;#39;:Companion:preDexDebug&amp;#39;.
  3. 【bzoj2226】[Spoj 5971] LCMSum 欧拉函数
  4. 九度OJ 1203:IP地址 (字符串处理)
  5. maven filters 和 resource
  6. PAT 1060. 爱丁顿数(25)
  7. struts2 获取表单数据封装到list和map集合
  8. 教你使用SQL查询(1-12)
  9. mysql设置指定ip远程访问连接的方法
  10. 保护眼睛,win7家庭版如何修改窗口的背景颜色