#include using namespace std; int P[40005]={1,1},phi[40005]; vector prime; void getphi(int n){ for(int i=2;i>n; getphi(n-1); return 0; } ...
分类：
其他好文 时间：
2017-09-18 13:24:33
阅读次数：
162
很多地方要用到素数，而能很快的写出代码筛出素数是很不错的我就单独写一个线性筛的代码和证明。 ...
分类：
其他好文 时间：
2017-09-17 17:30:35
阅读次数：
127
线性筛+莫比乌斯反演 盗波图 来自candy?大神 反演很重要的一条公式就是[gcd(i,j)==1]= 线性筛怎么推呢？ 我们分4个步骤，1.先推出f[1]，2.推出f[p],p是一个质数，3.由于线性筛筛的是积性函数，那么当gcd(i,p[j])==1的时候，f[i*p[j]]=f[i]*f[p ...
分类：
其他好文 时间：
2017-09-15 13:50:51
阅读次数：
157
题目描述 求 输入 第一行包含一个正整数T，表示有T组测试数据。接下来T<=10^5行，每行给出一个正整数N,N<=10^6。 输出 包含T行，依次给出对应的答案。 样例输入 7 1 10 100 1000 10000 100000 1000000 样例输出 1 2127 18446224 1830 ...
分类：
其他好文 时间：
2017-09-13 21:16:29
阅读次数：
152
蒟蒻要开始打数论模板了orz 线性筛都忘了怎么打，我太弱啦！ ...
分类：
其他好文 时间：
2017-09-13 20:14:52
阅读次数：
127
【BZOJ2694】Lcm Description 对于任意的>1的n gcd(a, b)不是n^2的倍数也就是说gcd(a, b)没有一个因子的次数>=2 对于任意的>1的n gcd(a, b)不是n^2的倍数也就是说gcd(a, b)没有一个因子的次数>=2 Input 一个正整数T表示数据组数 ...
分类：
其他好文 时间：
2017-09-13 19:15:41
阅读次数：
148
素数总是一个比较常涉及到的内容，掌握求素数的方法是一项基本功。 基本原则就是题目如果只需要判断少量数字是否为素数，直接枚举因子2 。。N^(0.5) ，看看能否整除N。 如果需要判断的次数较多，则先用下面介绍的办法预处理。 一般的线性筛法 首先先介绍一般的线性筛法求素数 void make_prim ...
分类：
其他好文 时间：
2017-09-13 17:09:24
阅读次数：
179
题目描述 对于正整数x，定义f(x)为x所含质因子的最大幂指数。例如f(1960)=f(2^3 * 5^1 * 7^2)=3, f(10007)=1, f(1)=0。给定正整数n,m，求$\sum\limits_{i=1}^n\sum\limits_{j=1}^mf(\gcd(i,j))$ 输入 第 ...
分类：
其他好文 时间：
2017-09-13 15:11:08
阅读次数：
195
#include using namespace std; vector prime; bool is_p[10000001]={1,1}; void getprime(int n){ for(int i=2;i>N>>M; getprime(N); while(M--){ int a; cin>>... ...
分类：
其他好文 时间：
2017-09-12 15:44:13
阅读次数：
129
思路是hdu6134的简化版，只需要在外面套上一个枚举素数就行了。 http://www.cnblogs.com/autsky-jadek/p/7491730.html ...
分类：
其他好文 时间：
2017-09-08 01:23:01
阅读次数：
164
