UVA 11426 - GCD - Extreme (II) (数论)
2015-07-12 09:59
316 查看
UVA 11426 - GCD - Extreme (II)
题目链接题意:给定N。求∑i<=ni=1∑j<nj=1gcd(i,j)的值。
思路:lrj白书上的例题,设f(n) = gcd(1, n) + gcd(2, n) + ... + gcd(n - 1, n).这种话,就能够得到递推式S(n) = f(2) + f(3) + ... + f(n) ==> S(n) = S(n - 1) + f(n);.
这样问题变成怎样求f(n).设g(n, i),表示满足gcd(x, n) = i的个数,这样f(n) = sum{i * g(n, i)}. 那么问题又转化为怎么求g(n, i),gcd(x, n) = i满足的条件为gcd(x / i, n / i) = 1,因此仅仅要求出欧拉函数phi(n / i),就能够得到与x / i互质的个数,从而求出gcd(x , n) = i的个数,这样总体就能够求解了
代码:
#include <stdio.h> #include <string.h> const int N = 4000005; int n; long long phi , s , f ; int main() { phi[1] = 1; for (int i = 2; i < N; i++) { if (phi[i]) continue; for (int j = i; j < N; j += i) { if (!phi[j]) phi[j] = j; phi[j] = phi[j] / i * (i - 1); } } for (int i = 1; i < N; i++) { for (int j = i * 2; j < N; j += i) { f[j] += phi[j / i] * i; } } s[2] = f[2]; for (int i = 3; i < N; i++) s[i] = s[i - 1] + f[i]; while (~scanf("%d", &n) && n) { printf("%lld\n", s ); } return 0; }
相关文章推荐
- 系统上线那点事续
- Oracle-PLSQL Developer使用笔记
- 日经社説 20150712 いまこそ暴力団への資金断て
- Oracle-Oracle10 数据空间建立,导入,导出--oracle10g 删除步骤
- Tomcat配置HTTPS方式(单向)
- 迷宫游戏(二维数组的应用)
- 日经春秋 20150712
- Terminal(终端) 开发常用指令
- lol英雄联盟打开公告或者直播出现界面卡死解决方法
- UVA 11892 ENimEN (简单博弈)
- Hibernate Dialect
- Oracle-数据实现竖排打印
- 使用轻量级ORM Dapper进行增删改查
- 理解三层转发原理
- 设计模式---策略模式
- SAP HR工资配置项1---工资计算周期配置
- Python绘制CPU曲线图
- Notepad++中设置Windows、Unix、Mac三种行尾换行符格式间的转换
- macbook air 开机黑屏解决方法
- qibuild的环境配置