2190: [SDOI2008]仪仗队

Time Limit: 10 Sec Memory Limit: 259 MB
Submit: 2112 Solved: 1337
[Submit][Status][Discuss]
Description

  作为体育委员,C君负责这次运动会仪仗队的训练。仪仗队是由学生组成的N * N的方阵,为了保证队伍在行进中整齐划一,C君会跟在仪仗队的左后方,根据其视线所及的学生人数来判断队伍是否整齐(如下图)。
  
  现在,C君希望你告诉他队伍整齐时能看到的学生人数。

Input

  共一个数N。

Output

  共一个数,即C君应看到的学生人数。

Sample Input

  4

Sample Output

  9

HINT

【数据规模和约定】   对于 100% 的数据,1 ≤ N ≤ 40000


写在前面:胡写数论题只求不挂啊
思路:


代码写的非常非常丑,如果想看好的code请务必别来这里= =
代码:

#include<bits/stdc++.h>
using namespace std;
int minp[40010],prime[40010],phi[40010],n;//数组从左到右分别表示最小质因子,素数,φ
long long ans;
main()
{scanf("%d",&n);for (int i=2;i<=n;i++){if (minp[i]==0)phi[i]=i-1,minp[i]=i,prime[++prime[0]]=i;for (int j=1;j<=prime[0];j++){if (i*prime[j]>n) break;if (minp[prime[j]*i]!=0) minp[prime[j]*i]=min(minp[prime[j]],min(prime[j],minp[i]));else minp[prime[j]*i]=min(prime[j],minp[i]);if (i*prime[j]/minp[prime[j]*i]%minp[prime[j]*i]==0)phi[i*prime[j]]=phi[i]*prime[j];else phi[i*prime[j]]=phi[i]*phi[prime[j]]; }}for (int i=2;i<n;i++) ans+=phi[i];ans=ans*2+3;printf("%lld",ans);
} 

【BZOJ2190】【codevs2296】仪仗队,数论练习之φ的线性筛相关推荐

  1. 中高级数论 [欧拉函数线性筛,二次剩余]

    欧拉函数线性筛 对于素数ppp, φ(p∗i)={p−1i=1p∗φ(i)p∣i(p−1)∗φ(i)p∤i\varphi (p*i)= \begin{cases} p-1& i=1\\ p*\ ...

  2. [CodeForces1603D] Artistic Partition(四边形不等式 + 决策单调性优化dp + 分治 + 线性筛 + 数论分块)

    problem codeforces 对于给定的正整数 l≤l\leql≤,定义 c(l,r)c(l,r)c(l,r) 为满足下列条件的正整数对 (i,j)(i,j)(i,j) 的数量: l≤i≤j≤ ...

  3. matlab狄利克雷函数,数论入门1——莫比乌斯函数,欧拉函数,狄利克雷卷积,线性筛,莫比乌斯反演,杜教筛...

    数论入门1 一个菜鸡对数论的一点点理解... 莫比乌斯函数 定义函数$\mu(n)$为: 当n有平方因子时,$\mu(n)=0$. 当n没有平方因子时,$\mu(n)=(-1)^{\omega(n)} ...

  4. 第二十九章 数论——中国剩余定理与线性同余方程组

    第二十九章 数论--中国剩余定理与线性同余方程组 一.中国剩余定理 1.作用: 2.内容: 3.证明: (1)逆元的存在性 (2)验证定理的正确性 4.代码实现: (1)步骤: (2)问题: (3)代 ...

  5. 洛谷 - P1217 [USACO1.5]回文质数 Prime Palindromes(欧拉线性筛+数论小知识)

    题目链接:点击查看 题目大意:给出一个a和b,求出区间[a,b]中所有回文质数,回文质数指的是一个数既是回文数也是质数 题目分析:首先质数好判断,我们只需要用欧拉筛预处理打好表即可,判断回文数我们也只 ...

  6. [数论]线性筛——约数个数与约数和

    参考博客 参考博客 参考博客 这个讲的挺好 预备知识点: 大于1的数n可以分解质因数: n=p1a1×p2a2×p3a3*-*pka n的约数的个数是(a1+1) * (a2+1) * (a3+1)- ...

  7. 【bzoj3309】DZY Loves Math 莫比乌斯反演+线性筛

    Description 对于正整数n,定义f(n)为n所含质因子的最大幂指数.例如f(1960)=f(2^3 * 5^1 * 7^2)=3, f(10007)=1, f(1)=0. 给定正整数a,b, ...

  8. P3383 【模板】线性筛素数

    https://www.luogu.com.cn/problem/P3383 //线性筛法 /* P3383 [模板]线性筛素数 https://www.luogu.com.cn/problem/P3 ...

  9. 素数一套:Miller-Rabin 素性检验算法Pollard-Rho算法线性筛——Upside down primesDivisions

    部分目录 Solved 94 / 304 K Gym 100753K Upside down primes 高效判断素数 快速幂取模 继续 Miller-Rabin 素性检验算法 Unsolved 6 ...

最新文章

  1. 聚类小分子数据集(基于RDKit的Python脚本)
  2. 刘永鑫:20分钟讲解微生物组数据分析与可视化实战
  3. Java-ReentrantLock初探
  4. Xpath语法学习记录
  5. “假冒hao123”“北大青鸟”被黑 钓鱼挂马两不误
  6. Linux 硬盘管理
  7. Flutter进阶—自定义主题风格
  8. php 正则 修饰符,php 正则修饰符
  9. J2SE、J2EE与J2ME
  10. c语言简单系统实例,简单的c语言编程实例
  11. Win7 远程桌面限制IP
  12. 20210319_23期_心跳检测_Task02_数据探索性分析
  13. avr单片机c语言计算log,AVR单片机定时器初值计算公式与方法 - 单片机定时器初值计算公式(51单片机和AVR单片机的初值计算三种方法)...
  14. 多家软件厂商卷入360与腾讯之争
  15. windows“运行”自定义“命令”的实现
  16. 华硕笔记本电脑一键重装系统详细教程
  17. 帧同步优化难点及解决方案
  18. 如何上传自己的项目到Maven中央仓库
  19. 长生生物事件的反思:质量是生命
  20. spotlight on oracle

热门文章

  1. 从React 编程到好莱坞
  2. 多元算力加持,华为云鲲鹏大数据服务公测上线
  3. 三伏天里小试牛刀andriod 开发 #华为云·寻找黑马程序员#
  4. 某学校有计算机主机1300台,网络互联技术与实践第14章:构建基于静态路由的多层网络.ppt...
  5. c语言程序2048_C语言2048小游戏演示和说明
  6. mongodb php update,MongoDB文档的更新(php代码实例)
  7. mysql二分法查找亿行_算法——二分法查找(binarySearch)
  8. router vue 页签文字_vue-router实现tab标签页(单页面)详解
  9. libevent c++高并发网络编程_【多线程高并发编程】Callable源码分析
  10. 华为p20Android怎么解开,华为P20如何获得root权限来解决自启动手机应用程序的问题...