GCD BZOJ2818 [省队互测] 数学
题目描述
给定整数N,求1<=x,y<=N且Gcd(x,y)为素数的数对(x,y)有多少对.
输入输出格式
输入格式:
一个整数N
输出格式:
答案
输入输出样例
4
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 [省队互测] 数学相关推荐
- P6271 [湖北省队互测2014]一个人的数论(莫比乌斯反演,拉格朗日插值)
整理的算法模板合集: ACM模板 点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划 P6271 [湖北省队互测2014]一个人的数论(莫比乌斯反演,拉格朗日插值) Problem Sol ...
- 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 ...
- [bzoj3670][2014湖北省队互测week2]似乎在梦中见过的样子
Description 已知一个字符串S,求它有多少个形如A+B+A的子串(len(A)>=k,len(B)>=1 ). Input 第一行一个字符串,第二行一个数 k. Output 仅 ...
- 2017.10.24队内互测——压轴出场的互测终曲|(*_-)
出题人: Sherlock, Frank, WWQ, MurasameKatana 终于到了我们组出题啦. 题面都是我自己写的2333 Problem 1 :令咒 题目来源:http://codevs ...
- 2017.9.16队内互测——老年组Day1
2017.9.16队内互测--老年组Day1 出题人:feather,MeiCo,Summer,Black Problem 1: 对于100%的数据,n<=1000 截图比较奇怪-还请见谅 简单 ...
- 2108889队2021年数学建模美赛C题花絮视频!
2108889队2021年数学建模美赛C题花絮视频
- UOJ#191. 【集训队互测2016】Unknown
UOJ#191. [集训队互测2016]Unknown 题目描述 Solution 二进制分组. 每一个组内维护一个斜率单调减的凸包. 因为有删点,避免出现反复横跳产生的爆炸复杂度,需要等到同一深度的 ...
- [2015国家集训队互测]口胡
比赛链接 http://uoj.ac/contest/11 口胡题解 A.[集训队互测2015]Robot 直接果断打暴力了...这个暴力很好写,我就不废话了 B.[集训队互测2015]Marketi ...
- 【LOJ3077】「2019 集训队互测 Day 4」绝目编诗
[题目链接] 点击打开链接 [思路要点] 不难发现各个边双连通分量可以分开处理,桥边可以直接删除. 可以证明,对于每一个边双连通分量,当 M−NM-NM−N 超过 O(N)O(\sqrt{N})O(N ...
最新文章
- Windows 7 搭建 Mobile 6 真机调试开发环境
- Hibernate的延迟检索和立即检索
- LED调光,PFM即pulse frequence modulation
- Redis 如何保持和MySQL数据一致【一】
- 牛客网 [编程题]餐馆
- 我的docker随笔25:一个测试用的镜像制作过程
- 软件测试为什么需要学习Linux的知识?Linux学到什么程度?-log5
- 开发平台怎么选?来看看专业人士怎么说
- EMNLP'21 | 检索式对话情感回复
- 狂神说ElasitcSearch笔记(一)
- ORACLE 查询表空间使用情况
- 如何在photoshop中等比例缩放一张图
- php公众号解决投票刷票,微信公众号投票活动如何防止刷票现象
- 金融市场数据平台的架构设计之道
- 杜骡的前端面试题(大全)
- GESD 离群值检验——理论与 Python 实现
- java 中的网络编程(Socket、TCP三次握手四次挥手、TCP/UDP/URL)
- 华为ensp防火墙nat64案例配置
- 服务器性能计时器如何关闭,如何在游戏关闭时继续倒数计时器?
- 《Unity 游戏设计与实现 南梦宫一线程序员的开发实例》学习笔记01