1268.可见点数

时间限制:1000MS内存限制:65536KB

题目描述

ZPS经过长期的努力争取,终于成为了0901班的领操员,他要带领0901班参加广播操比赛。现在0901班的队伍可以看作是一个n*n的点阵,每个人都站在格点上。现在作为领操员的ZPS站(0,0)点,他想知道如果0901班的队伍站齐了,他能看到多少个人的脸(假设每个人的身高相同,体积相同)。

输入

一个正整数n。

输出

ZPS能看到多少个人的脸(当然他是看不到自己的脸的)。

输入样例复制

3

输出样例复制

5

说明

Data Constraint 40%的数据,n<=1500。 100%的数据,n<=100000。

题解:直观看求斜率不一样的直线个数,易得出(x,y)与(0,0)得到一根新的直线的要求是x,y互质,因为x,y不互质,必然存在(x div k,y div k)在之前出现过  k=gcd(x,y)。问题转化为线性筛求欧拉函数
欧拉函数性质:1. φ(x*y)=φ(x)*φ(y)  (x,y互质)

   2.φ(x)=x-1                     (x为质数)
                         3.φ(x*y)=φ(x)*y              (x mod y=0) 

constmaxn=100000;varphi:array[1..maxn]of int64;ans,f,prime:array[0..maxn]of int64;n,i:longint;procedure init;vari,pp,j:longint;beginphi[1]:=1;  //phi(1)=1for i:=2 to maxn dobeginif f[i]=0 thenbegininc(prime[0]);prime[prime[0]]:=i;phi[i]:=i-1;end;for j:=1 to prime[0] dobeginpp:=i*prime[j];if pp>maxn then break;f[pp]:=1;if i mod prime[j]=0 then phi[pp]:=phi[i]*prime[j] //(x mod y=0)else phi[pp]:=phi[i]*phi[prime[j]]; end;end;ans[1]:=0;ans[2]:=3;for i:=3 to n doans[i]:=ans[i-1]+phi[i-1]*2;end;beginreadln(n);init;writeln(ans[n]);end.

【NOIP2018复习】可见点数【数论】相关推荐

  1. szu 寒训复习day #4数论入门详解[修改]用Mardown再改了一下

    整数的取余运算 定义:带余除法 设a,ba,ba,b是整数,且b>0b>0b>0,则存在非负整数 q,rq,rq,r ,使得 a=bq+ra=bq+ra=bq+r 且 0≤r< ...

  2. NOIP2018复习计划

    这几天考了很多场试,发现高级一点的算法题大多都有套路,但是在满分的道路上仍然十分难... 但是这种题暴力分好拿,你只要学过就能拿不少分,只不过部分分到满分是一个质变的事情 但是难度稍低一点的题就要完全 ...

  3. 【复习】---【noip2009 普及】细胞问题 (1)

    马上又要有一年的NOIP----NOIP2018 作为一个去年考挂的退役选手,在这最后的时间复习复习.临阵磨枪吧..... 学的东西都忘了好多.....毕竟一年都没有在看过信竞了 不废话看题↓ 描述 ...

  4. 【原创】【个人向】CSP-S 2019 爆炸退役记 (已完成)

    占坑. 说在前面 Day -84 蒟蒻就该有蒟蒻的亚子. 我有什么本事 我有什么进路 我自己清楚得很 Noip2016只参加了普及组. Noip2017Day2发高烧没去. Noip2018没有好好准 ...

  5. 中山纪念中学2015年暑假集训总结

    看着阳光在地上洒下的影子,在纪中的校石上矗立的黑影,我踏进了纪中这所校园.       阳光正盛,正值回忆.       Day1被一天一套题吓傻.突然明白了纪中OI强盛的原因. #include & ...

  6. CQOI2016省选游记

    4.9-4.10.重庆邮电大学.每天5小时,3道题.今年完成更新的CQTSC. 第二次省选.依稀记得第一次省选,重庆最后一次5T5H心态赛,每道题都做不来的绝望和无助.网络流是啥?可持久化数据结构?莫 ...

  7. 举步维艰——回顾CSP2020

    举步维艰--回顾CSP2020 复出 7.20 上一年 CSP-J 水过一等奖之后真丢人 ,转身跳进文化课,结果中考还是炸裂,自主招生也没有过,只有综能保底进了 SMMS . 7.25~7.30 参加 ...

  8. 2020暑假集训 心得

    第一学期 2020年称得上是最奇幻的一年,由于疫情的爆发,第一学期是在家里度过的.在家里诱惑实在是太多了,很少时间能够静下心来学习,让我感到最舒服的时候只有在半夜凌晨的时候,一个人坐在电脑桌前思考问题 ...

  9. 数论分块专题复习(余数求和+模积和+Ice Rain+The Fool)

    文章目录 前提知识复习 T1:余数求和 title solution code T2:Ice Rain title solution code T3:The Fool title solution c ...

  10. NOIP2018初赛提高组复习提纲(By HGOI LJC)

    Download:https://pan.baidu.com/s/16khhFf_0RsUjJLETreb20w (PDF) https://pan.baidu.com/s/1BVZqLs3q1clZ ...

最新文章

  1. java矩形碰撞检测_旋转矩形的Java碰撞检测?
  2. 用python的turtle画炫酷的图-使用Python的turtle画炫酷图形
  3. 2017 ACM-ICPC 亚洲区(西安赛区)网络赛
  4. 如何查看云服务器的系统版本,如何查看云服务器系统版本
  5. OpenCV TSP(TravelingSalesman Problem)的实例(附完整代码)
  6. Jquery ValidateEngine表单验证
  7. 想入门平均月薪2.58w人工智能领域?看看BAT的工程师在学什么
  8. python生成随机整数_Python中,我想得到1~255随机整数,用rand.randint(1,255),但不想让里面170这个数出现,怎么写呢?...
  9. 什么样的摄影题材更能吸引大众的目光?
  10. 写个脚本快速启动前后端
  11. redis启动以及开机自启动
  12. 1Password 7:Mac的密码管理工具
  13. 项目管理办公室(PMO) 的不同定位
  14. Hierachy Viewer无法在真机上使用的解决方法
  15. python语言程序设计袁方答案_北邮“爱课堂”教学平台
  16. 解决电脑安装mathtype后microsoft office 2016的.ppt文件打不开问题
  17. IllegalArgumentException
  18. php里面substr函数,PHP substr()函数
  19. java实现actionlistener_Java实现了ActionListener问题
  20. 利用vmware虚拟机实验从SUSE SLES12sp3 升级到SLES15sp1的过程

热门文章

  1. 5464: Star Arrangements
  2. 成人学士学位英语单词(史上最全)
  3. php 读取图片bgr,快速解决cv2.imread()读取图像为BGR的问题
  4. 民营经济挑战未来发展
  5. python网格交易法详解_期货交易技巧:等分网格交易法详解
  6. 老兵新传 Visual Basic核心编程及通用模块开发pdf
  7. 《动手学深度学习》(PyTorch版)代码注释 - 51 【Style_transfer】
  8. html5比较热门的新标签,HTML5增加的几个新的标签
  9. Chuck语言学习笔记——0.前言:我为什么要学习这门语言
  10. 分区表中GLO字段对信息收集的影响