[BZOJ2820]YY的GCD 试题描述 神犇YY虐完数论后给傻×kAc出了一题给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对kAc这种 傻×必然不会了,于是向你来请教……多组输入 输入 第一行一个整数T 表述数据组数接下来T行,每行两个正整数,表示N, M 输出 T行,每行一个整数表示第i组数据的结果 输入示例 输出示例 数据规模及约定 T = 10000N, M <= 10000000 题解 设 易知 ……式1 根据莫比…

BZOJ2820 YY的GCD Description 神犇YY虐完数论后给傻×kAc出了一题给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对kAc这种傻×必然不会了,于是向你来请教……多组输入 Input 第一行一个整数T 表述数据组数接下来T行,每行两个正整数,表示N, M Output T行,每行一个整数表示第i组数据的结果 Sample Input 2 10 10 100 100 Sample Output 30 2791 H…

/** 题目:BZOJ2820 YY的GCD 链接:http://www.cogs.pro/cogs/problem/problem.php?pid=2165 题意:神犇YY虐完数论后给傻×kAc出了一题 给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对 kAc这种傻×必然不会了,于是向你来请教…… T = 10000 N, M <= 10000000 思路: f(n)表示gcd==n的对数. g(n)表示gcd的n的倍数的对数…

不知道为什么不想写总结,只是(因为是用别人的权限号交的所以)屯一个代码 #include<iostream> #include<cstdio> #include<algorithm> #include<cstring> #include<cmath> using namespace std; #define ll long long const ll inf=(ll)2e9; ll rd(){ll z=,mk=; char ch=getchar(…

题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=2820 题意:多次询问,求1<=x<=N, 1<=y<=M且gcd(x,y)为质数有多少对. 首先,    由于这里是多次询问,并且数据很大,显然不能直接求解,需要做如下处理.. 整数的除法是满足结合律的,然后我们设T=p*d,有: 注意到后面部分是可以预处理出来的,那么整个ans就可以用分块处理来求了,设 那么有,考虑当p|x时,根据莫比菲斯mu(x)的性质,px除以其它非…

Description 求有多少对(x,y)的gcd为素数,x<=n,y<=m.n,m<=1e7,T<=1e4. Solution 因为题目要求gcd为素数的,那么我们就只考虑素数mu的贡献就行了 对于p,对于k*p的贡献是mu[k] 然后加上整除分块优化就行了 p可以筛完素数处理,处理复杂度为O(n/log*log)正好为O(n) Code #include<cstdio> #include<algorithm> #include<cstring&g…

Description 神犇YY虐完数论后给傻×kAc出了一题给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对kAc这种 傻×必然不会了,于是向你来请教……多组输入 Input 第一行一个整数T 表述数据组数接下来T行,每行两个正整数,表示N, M Output T行,每行一个整数表示第i组数据的结果 Sample Input 2 10 10 100 100 Sample Output 30 2791 HINT T = 10000 N…

题解 题意 题目链接 Sol 反演套路题.. 不多说了,就是先枚举一个质数,再枚举一个约数然后反演一下. 最后可以化成这样子 \[\sum_{i = 1}^n \frac{n}{k} \frac{n}{k} \sum_{p \in P, p | k} \mu(\frac{K}{p})\] 然后后面的那一坨可以暴力预处理..复杂度不清楚,但是显然严格小于调和级数,所以也没啥大问题. /* */ #include<bits/stdc++.h> #define LL long long //#def…

这题跟2818一样的,只不过数据水一点,可以用多一个log的办法水过去…… 原题意思是求以下式子:$Ans=\sum\limits_{isprime(p)}\sum\limits_{i=1}^{a}\sum\limits_{i=1}^{b}[gcd(i,j)==p]$首先把p拿下来,得到$Ans=\sum\limits_{isprime(p)}\sum\limits_{i=1}^{a/p}\sum\limits_{i=1}^{b/p}[gcd(i,j)==1]$然后就跟1101一样了,我就复制下…

数论什么的全都忘光了吧QAQ 做了几道简单的题练习一下. bzoj1101: [POI2007]Zap 求有多少对数满足 gcd(x,y)=d, 1<=x<=a, 1<=y<=b 首先那个d肯定是要除下去的, 变成 gcd(x,y) = 1, 1<=x<=a, 1<=y<=b这样一个问题. 这样就变成了最经典的莫比乌斯反演, 设f(d)为 gcd(x,y) = d, 1<=x<=a, 1<=y<=b时满足要求的x,y的个数,F(d)…