题目大意:多次询问有多少个数对(x,y)满足a
首先利用容斥原理将询问分解 问题转化为求有多少个数对(x,y)满足x
这里就可以利用到莫比乌斯反演:
我们令F(d)为GCD(x,y)=d且x
f(d)为d|GCD(x,y)且x
那么显然有F(d)=(n/d)*(m/d)
但是直接套用公式还是O(n^2)级别的
考虑到(n/d)*(m/d)最多只会有2√n个商 因此我们可以枚举...
分类:
其他好文 时间:
2014-12-16 11:44:43
阅读次数:
167
http://www.lydsy.com/JudgeOnline/problem.php?id=2301和这题不是差不多的嘛~~【BZOJ】1101: [POI2007]Zap(莫比乌斯+分块)唯一不同的地方是这题有下界。。下界除以k的时候取上界,然后分块的时候因为有4个数,所以要分成4块来搞。。然...
分类:
其他好文 时间:
2014-11-29 01:25:31
阅读次数:
247
#include
#include
#include
#include
#include
#include
using namespace std;
typedef long long LL;
inline int read(){
int x = 0,f = 1; char ch = getchar();
while(ch '9'){if(ch == '-')f=...
分类:
其他好文 时间:
2014-11-04 21:15:02
阅读次数:
166
什么东西。。。搞了半天Mobius反演到底是什么还是没搞定。。。(至少会求了嘛。。。好不好)但是程序写出来了^_^,可惜意义不明T T 1 /************************************************************** 2 Problem: 2...
分类:
其他好文 时间:
2014-10-24 23:33:50
阅读次数:
272
其实就是容斥原理了代码: 1 uses math; 2 const maxn=55000; 3 var i,n,a,b,c,d,w,tot:longint; 4 ans:int64; 5 sum,mu,p:array[0..maxn] of int64; 6 procedure ge...
分类:
其他好文 时间:
2014-06-25 22:10:01
阅读次数:
299
2301: [HAOI2011]Problem bTime Limit: 50 SecMemory Limit: 256 MBSubmit: 1007Solved: 415[Submit][Status]Description对于给出的n个询问,每次求有多少个数对(x,y),满足a≤x≤b,c≤y≤...
分类:
其他好文 时间:
2014-06-18 13:44:16
阅读次数:
272
【题目描述】给你一对数a,b,你可以任意使用(a,b), (a,-b), (-a,b),
(-a,-b), (b,a), (b,-a), (-b,a),
(-b,-a)这些向量,问你能不能拼出另一个向量(x,y)。说明:这里的拼就是使得你选出的向量之和为(x,y)【输入格式】第一行数组组数t,(t ...
分类:
其他好文 时间:
2014-05-15 16:52:50
阅读次数:
260