SPOJ GCDEX (数论)

最新推荐文章于 2019-03-28 16:04:00 发布
最新推荐文章于 2019-03-28 16:04:00 发布
阅读量4.5k 收藏
点赞数 2
分类专栏： ACM_数学类
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/ACM_cxlove/article/details/11834983
版权
博客介绍了如何解决SPOJ的GCDEX问题，该问题涉及到数论中的sigma函数与最大公约数。由于直接枚举求解会导致超时，作者提出了利用sigma(gcd(i, n)) = sigma(d * phi[n / d])公式进行预处理phi数组，并以gcd为d进行枚举，实现O(n)预处理和O(1)查询的方法。" 109173632,10184517,0805合金贴片电阻参数解析与应用,"['电子工程', '电阻技术', '元器件', '电路设计', '功率电子']
摘要由CSDN通过智能技术生成

转载请注明出处，谢谢http://blog.csdn.net/ACM_cxlove?viewmode=contents    by---cxlove

题意：求sigma (gcd (i , j))  1 <= i < j <= n

和POJ 2480类似，如果枚举j，求的话，还是会TLE的。。。

考虑sigma(gcd (i , n)) = sigma (d * phi[n / d]) d | n。

做法同样是先预处理出phi，然后枚举gcd = d。

O(n)预处理，O(1)查询

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <set>
#include <map>
#include <string>
#include <vector>
#include <
最低0.47元/天 解锁文章
ACM_cxlove
关注
专栏目录
GCD - Extreme (II) （数论，欧拉函数）
Crazy
05-09 1186
题目来源：https://vjudge.net/problem/UVA-11426 【题意】 求在1-n之间所有任意两个数的的最大公因数的和。 【思路】 做完这道题，依旧感觉头脑懵懵的，看了好多篇博客方才看懂了这道题的解法，因为之前做过几道数论题，关于素数的比较多，这到题因为最大公因数就自然而然的想到了素数筛法，一直在找所谓的规律。最后终于崩溃了。。。 便搜了博客慰藉我不安的心灵，但是很多
spoj GCDEX - GCD Extreme
weixin_30955341的博客
10-12 112
传送门：http://www.spoj.com/problems/GCDEX/ 思路：令g(n)=ΣΣgcd(i,j)(i<=n,j<=i)，再令f(n)=Σgcd(i,n) 很显然ans=g(n)-(n+1)*n/2 那么g(n)=Σf(i)(i<=n) f(n)在上一题已证为积性函数，线筛求出，就可以在O(n)时间打出表来。 为了方便先设mindiv[i]为...
数论 - Funny scales(SPOJ - SCALE)
weixin_34326429的博客
04-25 91
Funny scales  Problem's Link  ---------------------------------------------------------------------------- Mean:  给定两个数n和x，有一个天平，初始时左盘为x，你需要从以下集合中选一些数字来放到两个盘中，使得两个盘相等（note:每个数字只...
spoj gcdex
a5199519的博客
04-15 122
题解： 首先我们设gcd(i,j)=k 所以我们就要求对于所有k的方案总数 可以线性帅选欧拉函数 然后算法一：枚举k，O(NT) 算法二：考虑到我们只要n/k的整数部分 容易证明是sqrt(n)级别的 所以就可以在O(Tsqrt(n))的时间内解决 但是要考虑卡常数 代码： #include<bits/stdc++.h> typedef long l...
【SPOJ-GCDEX】GCD Extreme【欧拉函数】【线性筛】
BraketBN
02-24 570
题意： 求∑(1 白书上原题...注意卡常就可以过。 #include typedef unsigned long long ULL; typedef unsigned int uint; const int maxn = 1000001; int cnt; int prime[maxn], phi[maxn]; ULL f[maxn]; bool isnotpri
SPOJ - GCDEX GCD Extreme
ACVevtor的博客
08-21 373
题目：给你一个N，求sigma(gcd(i,j))，1 思路：枚举gcd，利用欧拉函数计算 代码： #pragma comment(linker, "/STACK:1024000000,1024000000") #include #include #include #include #include #include #include #include #include #include #
SPOJ.rar_SPOJ
09-19
《SPOJ在线判题平台：越南版解题攻略》 SPOJ，全称Sphere Online Judge，是一个全球知名的在线编程竞赛平台，旨在提供一个环境让程序员可以测试、提交并评估他们的算法解决方案。这个名为"SPOJ.rar_SPOJ"的压缩包...
SPOJ-2.zip_SPOJ_SPOJ-TRIKA.cpp_online judge_sPOJ-Solution_spoj2
09-24
《SPOJ在线判题平台的解题策略与实例解析》 SPOJ（Sphere Online Judge）是一个全球知名的在线编程竞赛平台，它为程序员提供了一个检验自己算法能力的场所。这个名为"SPOJ-2.zip"的压缩包包含了在SPOJ上解决不同...
SPOJ.zip_SPOJ_sgu_them
09-24
《SPOJ问题解决策略：动态规划的应用》 在编程竞赛和算法研究中，SPOJ（Sphere Online Judge）是一个著名的在线平台，它提供了一系列挑战性的编程问题供参赛者解决。"SPOJ.zip_SPOJ_sgu_them"这个压缩包文件就包含...
【SPOJ-GCDEX】GCD Extreme（欧拉函数）
weixin_30770495的博客
03-28 140
题目： SPOJ-GCDEX (洛谷 Remote Judge) 分析： 求： \[\sum_{i=1}^{n}\sum_{j=i+1}^{n}gcd(i,j)\] 这道题给同届新生讲过，由于种种原因只讲了 \(O(n)\) 预处理欧拉函数 \(O(n)\) 查询的暴力做法，顺带提了一句 “这题能根号查询” 被教练嘴了 QAQ 。以及小恐龙给我说有 \(O(n\log n)\) 预处理 \(O(1...
SPOJ - DIV15(数论)
稻草人
07-29 771
题意：给一个数字字符串，要求给出能被15整除的最大的数。 思路：15是3和5的倍数，那么能够被15整除的的数，肯定也是3和5共同的整倍数。 而3和5的倍数，有很有特点，3的倍数要求每位上的数字之和是3的整倍数，而5的倍数要求个位数上是0或者5. 综合其特点，我们就可以找到思路了。。。 如果数字中有即没有0有没有5的话，则该数字不符合条件， 如果每位上的数字之和sum能被3整除，则该数字在
SPOJ：EXPOR - OR（数论）
junior19的博客
05-23 263
EXPOR - OR #goc-s01e01 Given an array of N integers A1, A2, A3…AN. If you randomly choose two indexes i ,j such that 1 ≤ i i | Aj? Input First line contains an integer T, the numbe
SPOJ GCD Extreme(推导+筛法)
Miracle_ma的专栏
09-05 575
T次询问，求∑n−1i=1∑nj=i+1gcd(i,j)T次询问，求\sum_{i=1}^{n-1}\sum_{j=i+1}^{n}gcd(i,j) T有2W，反正反演是T了，虽然感觉O(Tsqrt(n))能过T有2W，反正反演是T了，虽然感觉O(Tsqrt(n))能过 如果固定j，那么就是求T(j)=∑j−1i=1gcd(i,j)如果固定j，那么就是求T(j)=\sum_{i=1}^{j-1}
SPOJ LCMSUM&GCDEX
my coding life
09-26 1299
爱神博客上的两道数论题。SPOJ做的苦不堪言。。。灰常难受的说，各种卡时。。。。。然后因为一个小错误，WA了很长时间。。。。PS:代码能力已经真心不忍直视了，导致长沙网络赛没有过前70。。。。。。。不多说 先是LCMSUM这道题 f(n) = sigma{lcm(i,n)} =sigma{ (i*n)/gcd(i,n)} 分母都不一样，所以要去分母，将n与gcd约分后，可以看到，对于gcd(
SPOJ FACVSPOW(数论+二分)
十分残念的博客
02-25 613
传送门：http://www.spoj.com/problems/FACVSPOW/en/ 题意：给出一个数a，要求找到最小的n使得n!>a^n 题解：先将不等式两边取对数得到ln(1)+ln(2)+...+ln(n)>nln(a) 可以发现a #include using namespace std; const int MX = 3e6 + 5; double sum[MX]; voi
Spoj数论专场解题报告
beat anyone if necessary
05-29 6169
Spoj数论专题http://acm.hust.edu.cn:8080/judge/contest/view.action?cid=8035#overview  Problem A SPOJ ARCTAN Use of Function Arctan Problem B SPOJ CUTSQRS Cutting off
UVA 11426 GCD - Extreme (II) （欧拉函数）
AI蜗牛之家
05-16 1439
题目链接：http://acm.hust.edu.cn/vjudge/contest/view.action?cid=70017#problem/O 假设a、b(a     假设和b互质的数有n个，也就是n对(?,b)（?和b互质），那么在i、j循环到?、b时结果会增加n，循环到(2*?,2*b)时结果就会增加2*n...当i、j循环到k*?、k*b时结果就会增加k*n。那么我们不
gcd以及ex_gcd的总结
AC_Dream
05-22 2022
gcd()---表示最大公约数，常用方法是欧几里得算法 ex_gcd()---表示扩展欧几里得算法 定义1：a和b是两个不全为0的整数，称a与b的公因子中最大的为a和b的最大公约数，用gcd(a,b)来表示。 定义2：a和b是两个非0的整数，称a与b的公倍数中最小的为a和b的最小公倍数，用lcm(a,b)来表示。 最小公倍数与最大公约数之间的性质：     1、若a|m,b|m,那么lc
洛谷SPOJ如何注册
最新发布
07-27
2. 在网站顶部导航栏中找到“题库”选项，将鼠标悬停在上面，然后选择“SPOJ”。 3. 在SPOJ页面上，您会看到一个提示，要求您注册SPOJ账号并进行绑定。点击提示中的链接，将会跳转到SPOJ注册页面。 4. 在SPOJ注册...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值