题目大意：令F(i)为i的约数和，多次询问对于1
n,m
首先如果不考虑a的限制 令g(i)为1
那么显然有
利用线性筛处理出F(i) 那么答案显然是
治好了我多年的公式恐惧症。。。
现在我们只需要求出的前缀和 这个问题就能在O(√n)的时间内出解
枚举每一个i 枚举i的倍数 暴力即可求出这个函数 然后处理前缀和即可 复杂度是O(nlogn)的
那么现在有了a的...
分类：
编程语言 时间：
2014-12-22 11:20:34
阅读次数：
221
http://www.lydsy.com/JudgeOnline/problem.php?id=1101无限膜拜数论和分块orz首先莫比乌斯函数的一些性质可以看《初等数论》或《具体数学》或贾志鹏的《线性筛法和积性函数》我写一些笔记啥的吧。。首先莫比乌斯函数的定义及一些性质（免去证明）：$$\mu (...
分类：
其他好文 时间：
2014-11-28 15:57:22
阅读次数：
351
题目大意：求1~n的排列能组成多少种小根堆
考虑一个1~i的排列所构成的堆，l为左儿子大小，r为右儿子的大小
那么1一定是堆顶 左儿子和右儿子分别是一个堆 显然如果选出l个数给左儿子 那么左儿子的方案数显然是f[l]，右儿子的方案数为f[r]
于是有f[i]=C(i-1,l)*f[l]*f[r]
于是我们线性筛处理出阶乘和阶乘的逆元 代入即可得到WA
原因是这题n可以大于p 此时要用到L...
分类：
其他好文 时间：
2014-11-21 16:21:39
阅读次数：
156
题目：输出n！中素数因数的个数。
分析：数论。这里使用欧拉筛法计算素数，在计算过程中求解即可。
传统筛法是利用每个素数，筛掉自己的整数倍；
欧拉筛法是利用当前计算出的所有素数，乘以当前数字筛数；
所以每个前驱的素椅子个数一定比当前数的素因子个数少一个。
说明：又一次用了“线性筛法”。
#include
#incl...
分类：
其他好文 时间：
2014-11-17 14:04:17
阅读次数：
138
题目大意：求n以内所有数的约数个数和
100W，n√n别想了
线性筛可以处理，对于每个数记录最小质因数的次数
令factoral[i]为i的因数个数 cnt[i]为i的最小质因数的次数
若i为质数 则factoral[i]=2 cnt[i]=1
若i%prime[j]!=0 则factoral[prime[j]*i]=factorial[i]*2 cnt[prime[j]*i]=1
若...
分类：
其他好文 时间：
2014-11-15 18:53:52
阅读次数：
125
一、线性筛法众所周知。。。线性筛就是在O(n)的时间里找出所有素数的方法code：void get_prime(int N){ int i, j, k; memset(Flag, sizeof(Flag), 0); for (i = 2; i N) break; ...
分类：
其他好文 时间：
2014-11-01 13:25:26
阅读次数：
152
题目大意：给出一个等差数列，问这个等差数列的第n个素数是什么。
思路：这题主要考如何筛素数，线性筛。详见代码。
CODE：
#include
#include
#include
#include
#define MAX 1000010
using namespace std;
int prime[MAX],primes;
bool notp[MAX];
int a...
分类：
其他好文 时间：
2014-10-31 12:02:10
阅读次数：
199
时间复杂度为O(n)线性筛选n以内的素数， 关键在于if i mod p[j]=0 then break;保证每个数最多被筛选一次。
分类：
其他好文 时间：
2014-10-28 17:21:45
阅读次数：
277
题目大意：给两个数，第一个数的范文是10^100，第二个数10^6，第一个数是两个质数的乘积，问有没有不超过第二个数的数是第一个树的因子。
思路：10^6中只有7w+个素数，只要挨个判定能不能整除即可。然后这个题素数必须线性，不然就T。还有高精度一开始我压了4位，之后就wa，调了很长时间发现判断整除的过程中爆int了，换成long long就是T，最后十分生气，直接压到了7位，果断A了，...
分类：
其他好文 时间：
2014-10-23 14:34:07
阅读次数：
212
题目大意：给定n和m，求Σ(1
i和j的限制不同，传统的线性筛法失效了，这里我们考虑容斥原理
令f[x]为GCD(i,j)=x的数对(i,j)的个数，这个不是很好求
我们令g[x]为存在公因数=x的数对(i,j)的个数(注意不是最大公因数！)，显然有g[x]=(n/x)*(m/x)
但是这些数对中有一些的最大公因数为2d,3d,4d，我们要把他们减掉
于是最终f[x]=(n/x)*(m/...
分类：
其他好文 时间：
2014-10-09 15:46:48
阅读次数：
366
