题目描述

给定整数N,求1<=x,y<=N且Gcd(x,y)为素数的数对(x,y)有多少对.

输入输出格式

输入格式:

一个整数N

输出格式:

答案

输入输出样例

输入样例#1: 复制

4

输出样例#1: 复制

4

说明

对于样例(2,2),(2,4),(3,3),(4,2)

1<=N<=10^7

来源:bzoj2818

本题数据为洛谷自造数据,使用CYaRon耗时5分钟完成数据制作。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstdlib>
#include<cstring>
#include<string>
#include<cmath>
#include<map>
#include<set>
#include<vector>
#include<queue>
#include<bitset>
#include<ctime>
#include<time.h>
#include<deque>
#include<stack>
#include<functional>
#include<sstream>
//#include<cctype>
//#pragma GCC optimize(2)
using namespace std;
#define maxn 2000005
#define inf 0x7fffffff
//#define INF 1e18
#define rdint(x) scanf("%d",&x)
#define rdllt(x) scanf("%lld",&x)
#define rdult(x) scanf("%lu",&x)
#define rdlf(x) scanf("%lf",&x)
#define rdstr(x) scanf("%s",x)
#define mclr(x,a) memset((x),a,sizeof(x))
typedef long long  ll;
typedef unsigned long long ull;
typedef unsigned int U;
#define ms(x) memset((x),0,sizeof(x))
const long long int mod = 1e9 + 7;
#define Mod 1000000000
#define sq(x) (x)*(x)
#define eps 1e-5
typedef pair<int, int> pii;
#define pi acos(-1.0)
//const int N = 1005;
#define REP(i,n) for(int i=0;i<(n);i++)
typedef pair<int, int> pii;inline int rd() {int x = 0;char c = getchar();bool f = false;while (!isdigit(c)) {if (c == '-') f = true;c = getchar();}while (isdigit(c)) {x = (x << 1) + (x << 3) + (c ^ 48);c = getchar();}return f ? -x : x;
}ll gcd(ll a, ll b) {return b == 0 ? a : gcd(b, a%b);
}
int sqr(int x) { return x * x; }/*ll ans;
ll exgcd(ll a, ll b, ll &x, ll &y) {if (!b) {x = 1; y = 0; return a;}ans = exgcd(b, a%b, x, y);ll t = x; x = y; y = t - a / b * y;return ans;
}
*/bool vis[10000002];
int phi[10000002];
ll sum[10000002 >> 1];
int pri[10000002];
int tot;
int N;
void init() {phi[1] = 1;for (int i = 2; i <= 10000000; i++) {if (!vis[i]) {pri[++tot] = i; phi[i] = i - 1;}for (int j = 1; j <= tot && i*pri[j] <= 10000000; j++) {vis[i*pri[j]] = true;phi[i*pri[j]] = phi[i] * phi[pri[j]];if (i%pri[j] == 0) {phi[i*pri[j]] = phi[i] * pri[j];break;}}}
}int main()
{// ios::sync_with_stdio(0);N = rd();init();for (int i = 1; i <= 10000000 / 2; i++)sum[i] = 1ll*sum[i - 1] + 1ll*phi[i];ll ans = 0;for (int i = 1; i <= tot; i++) {if (pri[i] > N)break;ans += 1ll * 2 * sum[N / pri[i]] - 1;}cout << (ll)ans << endl;return 0;
}

转载于:https://www.cnblogs.com/zxyqzy/p/10363171.html

GCD BZOJ2818 [省队互测] 数学相关推荐

  1. P6271 [湖北省队互测2014]一个人的数论(莫比乌斯反演,拉格朗日插值)

    整理的算法模板合集: ACM模板 点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划 P6271 [湖北省队互测2014]一个人的数论(莫比乌斯反演,拉格朗日插值) Problem Sol ...

  2. P6271 [湖北省队互测2014]一个人的数论(莫比乌斯反演 + 伯努利数)

    P6271 [湖北省队互测2014]一个人的数论 ∑i=1nim[gcd⁡(i,n)=1]∑d∣nμ(d)dm∑i=1ndim由伯努利数可知∑i=0nim=1m+1∑i=0mCm+1iBi(n+1)m ...

  3. [bzoj3670][2014湖北省队互测week2]似乎在梦中见过的样子

    Description 已知一个字符串S,求它有多少个形如A+B+A的子串(len(A)>=k,len(B)>=1 ). Input 第一行一个字符串,第二行一个数 k. Output 仅 ...

  4. 2017.10.24队内互测——压轴出场的互测终曲|(*_-)

    出题人: Sherlock, Frank, WWQ, MurasameKatana 终于到了我们组出题啦. 题面都是我自己写的2333 Problem 1 :令咒 题目来源:http://codevs ...

  5. 2017.9.16队内互测——老年组Day1

    2017.9.16队内互测--老年组Day1 出题人:feather,MeiCo,Summer,Black Problem 1: 对于100%的数据,n<=1000 截图比较奇怪-还请见谅 简单 ...

  6. 2108889队2021年数学建模美赛C题花絮视频!

    2108889队2021年数学建模美赛C题花絮视频

  7. UOJ#191. 【集训队互测2016】Unknown

    UOJ#191. [集训队互测2016]Unknown 题目描述 Solution 二进制分组. 每一个组内维护一个斜率单调减的凸包. 因为有删点,避免出现反复横跳产生的爆炸复杂度,需要等到同一深度的 ...

  8. [2015国家集训队互测]口胡

    比赛链接 http://uoj.ac/contest/11 口胡题解 A.[集训队互测2015]Robot 直接果断打暴力了...这个暴力很好写,我就不废话了 B.[集训队互测2015]Marketi ...

  9. 【LOJ3077】「2019 集训队互测 Day 4」绝目编诗

    [题目链接] 点击打开链接 [思路要点] 不难发现各个边双连通分量可以分开处理,桥边可以直接删除. 可以证明,对于每一个边双连通分量,当 M−NM-NM−N 超过 O(N)O(\sqrt{N})O(N ...

最新文章

  1. Windows 7 搭建 Mobile 6 真机调试开发环境
  2. Hibernate的延迟检索和立即检索
  3. LED调光,PFM即pulse frequence modulation
  4. Redis 如何保持和MySQL数据一致【一】
  5. 牛客网 [编程题]餐馆
  6. 我的docker随笔25:一个测试用的镜像制作过程
  7. 软件测试为什么需要学习Linux的知识?Linux学到什么程度?-log5
  8. 开发平台怎么选?来看看专业人士怎么说
  9. EMNLP'21 | 检索式对话情感回复
  10. 狂神说ElasitcSearch笔记(一)
  11. ORACLE 查询表空间使用情况
  12. 如何在photoshop中等比例缩放一张图
  13. php公众号解决投票刷票,微信公众号投票活动如何防止刷票现象
  14. 金融市场数据平台的架构设计之道
  15. 杜骡的前端面试题(大全)
  16. GESD 离群值检验——理论与 Python 实现
  17. java 中的网络编程(Socket、TCP三次握手四次挥手、TCP/UDP/URL)
  18. 华为ensp防火墙nat64案例配置
  19. 服务器性能计时器如何关闭,如何在游戏关闭时继续倒数计时器?
  20. 《Unity 游戏设计与实现 南梦宫一线程序员的开发实例》学习笔记01

热门文章

  1. DPM2007轻松恢复Exchange邮件,DPM2007系列之三
  2. 用hmmlearn学习隐马尔科夫模型HMM
  3. 2015年度打榜之数据备份硬件
  4. PHP(十二)文件操作
  5. 浅谈java中的四个核心概念【转】
  6. overflow:hidden;zoom:1 理解转__
  7. PingingLab传世经典系列《CCNA完全配置宝典》-2.5 RIPv2基本配置
  8. Luhn(卢恩)算法,检测信用卡号的合法性
  9. Facebook和Google,做不了恋人做仇人
  10. 【太阳软件】动态换ip软件:反爬虫对抗对数据分析的影响