博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
HDU 3501 Calculation 2
阅读量:6321 次
发布时间:2019-06-22

本文共 1525 字,大约阅读时间需要 5 分钟。

题意:

给定正整数n,求所有小于n的且与n为非互质数的和。输入对于每一个测试用例,都有一个包含正整数n(1≤n≤1000000000)的行。包含单个0的行将遵循最后一个测试用例。输出对于每个测试用例,您应该在一行中打印求和模1000000007。

思路:

考虑互质的数的一个性质,n与m互质,那么m-n与m也一定互质。

好,我们再来考虑这个题。

让求与n为非互质数的数的和,那么我们是不是就可以先把所有数的和求出来然后再减去与n互质的数的和???!!

好,我们就这样干。 有人又要问了,为什么要这样啊,不是更麻烦吗??!!

好像不是这样的,因为我们可以用欧拉函数求与n互质的数的个数啊!!这样我们就不用暴力枚举了啊。。。

(好像知道与n互质的数的个数也没有用啊。。。)

(⊙o⊙)…      我们考虑数学里的小技巧。再求n内所有数的和时,我们是不是可以用等差数列前n项和公式??(公式是什么。。。)

等差数列求和公式                  s=n*(n-1)/2

好,这样我们就轻松求出n以内所有数的和了。 我们接下来再求与n互质的数的和。(怎么求??)

我们继续看上面给出的性质,一共有phi(n)个与n互质的数,而且这phi(n)个数一定存在m、n-m的关系,那么这两个数两两相加不就是n吗??我们可以凑出phi(n)/2对n,这样这phi(n)个数的和不就是n*phi(n)/2吗??!!

好了,这样我们就完工了。。。。

代码:

#include
#include
#include
#include
#include
#define mod 1000000007using namespace std;int t,n;long long ans;int read(){ int x=0,f=1; char ch=getchar(); while(ch<'0'||ch>'9'){
if(ch=='-') f=-1; ch=getchar();} while(ch>='0'&&ch<='9') {x=x*10+ch-'0'; ch=getchar();} return x*f;}int get_phi(int x){ int sum=x; if(x%2==0) { while(x%2==0) x/=2; sum/=2; } for(int i=3;i*i<=x;i+=2) { if(x%i==0) { while(x%i==0) x/=i; sum=sum/i*(i-1); } } if(x>1) sum=sum/x*(x-1); return sum;}int main(){ while(1) { n=read(); if(n==0) break; ans=(long long )n*(n-1)/2; ans-=(long long )n*get_phi(n)/2; printf("%lld\n",ans%mod); } return 0;}

转载于:https://www.cnblogs.com/z360/p/7305052.html

你可能感兴趣的文章
java okhttp发送post请求
查看>>
laravel查找某个类拥有的方法:
查看>>
Handler相关
查看>>
关于Java 枚举类型的自定义属性
查看>>
Entity Framework Code First属性映射约定
查看>>
解决zabbix“ZBX_NOTSUPPORTED: Timeout while executing a shell script”报错
查看>>
IntelliJ IDEA2018.1、2017.3激活
查看>>
MSSQL 调用C#程序集 实现C#字符串到字符的转化
查看>>
【微信】1.微信小程序开发--入门
查看>>
DELL平板如何安装WIN10系统 -PE启动问题
查看>>
Mac上Vim的配置文件及插件
查看>>
js进阶 9-12 如何将数组的信息添加到下拉列表
查看>>
js 进阶 10 js选择器大全
查看>>
Linux 下MongoDb的安装
查看>>
Amazon RDS多区域高可用测试
查看>>
为什么分布式一定要有Redis?
查看>>
微服务之分布式跟踪系统(springboot+pinpoint)
查看>>
FAT AP v200R005 配置二层透明模式(web&命令行,开局)
查看>>
01-股票预测-买卖急切度
查看>>
系列文章--一步一步学Silverlight2
查看>>