题目描述

for i=1 to n

for j=1 to n

 sum+=gcd(i,j)

给出n求sum. gcd(x,y)表示x,y的最大公约数.

输入输出格式

输入格式:

n

输出格式:

sum

输入输出样例

输入样例#1:

2
输出样例#1:

5

说明

数据范围

30% n<=3000

60% 7000<=n<=7100

100% n<=100000


题目的意思大概是这样的

O(n2)枚举当然是不行的啦。

考虑枚举k,求gcd为k的“数对”的个数。

而可以证明gcd为k的“数对”的个数为

利用容斥把gcd为2k,3k,4k的“数对”的个数减去就好啦?

注意k要从大到小枚举。

 #include<cstdio>
#include<cstring>
#include<iostream>
using namespace std; typedef long long ll; const int maxn=; int n;
ll dp[maxn],ans=; int main(){
scanf("%d",&n);
for(int i=n;i>;i--){
dp[i]=1ll*(n/i)*(n/i);
for(int j=(i<<);j<=n;j+=i)
dp[i]-=dp[j];
ans+=dp[i]*i;
}
printf("%lld\n",ans);
return ;
}

最新文章

  1. 关于JavaScript中的delete操作
  2. iOS开发之多线程技术(NSThread、OperationQueue、GCD)
  3. HDU 3999 二叉排序树
  4. 修改wamp默认网站目录
  5. lua 操作中文字符串之截取和长度竖排显示
  6. How can I pretty-print JSON in python?
  7. C++之路进阶——bzoj1823(满汉全席)
  8. 课堂所讲整理:super和转型(修改版)
  9. ios的runtime为什么可以动态添加方法
  10. CC++刚開始学习的人编程教程(9) Windows8.1安装VS2013并捆绑QT与编程助手
  11. crawler_x-requested-with 请求头
  12. tableviewcell滑动显示多个按钮UITableViewRowAction(转载)
  13. oracle删除数据后表空间仍过大问题解决方法
  14. bam文件softclip , hardclip ,markduplicate的探究
  15. VxWorks中的中断应用设计要点
  16. Linux之centos系统常用命令总结
  17. Android: 待机时如何让程序继续运行 extends Service
  18. java多线程快速入门(九)
  19. 【iCore4 双核心板_ARM】例程六:IWDG看门狗实验——复位ARM
  20. C#中的事件(event)处理机制

热门文章

  1. Java中equals和hashcode的区别?
  2. 使用fastClick.js所产生的一些问题
  3. Center OS 7安装 Apollo
  4. map 结构体
  5. Delphi locate函数
  6. cnn知识点汇总
  7. csp-s模拟测试95
  8. (转)使用Apache的ab工具进行压力测试
  9. python3 线程 threading.Thread GIL性能详解(2.3)
  10. (转)python资料汇总(建议收藏)零基础必看