模板 - Min_25 筛法求素数和

整理的算法模板合集： ACM模板

#include <bits/stdc++.h>
#define ll long long
using namespace std;
const int N = 1000010;
int prime[N], id1[N], id2[N], flag[N], ncnt, m;
ll g[N], sum[N], a[N], T;
ll n;
int ID(ll x) {return x <= T ? id1[x] : id2[n / x];
}
ll calc(ll x) {return x * (x + 1) / 2 - 1;
}
ll f(ll x) {return x;
}
ll init(ll n){T = sqrt(n + 0.5);for (int i = 2; i <= T; i++) {if (!flag[i]) prime[++ncnt] = i, sum[ncnt] = sum[ncnt - 1] + i;for (int j = 1; j <= ncnt && i * prime[j] <= T; j++) {flag[i * prime[j]] = 1;if (i % prime[j] == 0) break;}}for (ll l = 1; l <= n; l = n / (n / l) + 1) {a[++m] = n / l;if (a[m] <= T) id1[a[m]] = m; else id2[n / a[m]] = m;g[m] = calc(a[m]);}for (int i = 1; i <= ncnt; i++)for (int j = 1; j <= m && (ll)prime[i] * prime[i] <= a[j]; j++)g[j] = g[j] - (ll)prime[i] * (g[ID(a[j] / prime[i])] - sum[i - 1]);
}
ll solve(ll x){if(x<=1){return x;}return n=x,init(n),g[ID(n)];
}
int main() {int t = 1;while(t-- ) {memset(g,0,sizeof(g));memset(a,0,sizeof(a));memset(sum,0,sizeof(sum));memset(prime,0,sizeof(prime));memset(id1,0,sizeof(id1));memset(id2,0,sizeof(id2));memset(flag,0,sizeof(flag));ncnt=m=0;scanf("%lld", &n);printf("%lld\n", solve(n));}
}

模板 - Min_25 筛法求素数和相关推荐

普通素数筛法求素数二次筛法求素数 MillerRabin素数测试【模板】
素数和合数共同的性质: 1.a > 1是合数,当且仅当a = b * c,其中1 < b < a,1 < c < a. 2.合数必有素数因子. 3.如果d > 1, ...
一般筛法和快速线性筛法求素数求素数的一点总结
素数总是一个比较常涉及到的内容,掌握求素数的方法是一项基本功. 基本原则就是题目如果只需要判断少量数字是否为素数,直接枚举因子2 ..N^(0.5) ,看看能否整除N. 如果需要判断的次数较多,则先用 ...
线性筛法求素数c语言,[算法]素数筛法(埃氏筛法线性筛法)
一.素数筛的定义给定一个整数n,求出[1,n]之间的所有质数(素数),这样的问题为素数筛(素数的筛选问题). 二.埃氏筛法(Eratosthenes筛法) 埃氏筛法又叫做Eratosthenes筛法 ...
一般筛法求素数+快速线性筛法求素数
一般筛法求素数+快速线性筛法求素数标签: 正则表达式算法优化扩展c 2010-08-22 01:28 28738人阅读评论(8) 收藏举报分类: 算法学习资料(5) 版权声明:本文为博主原 ...
CodeForce 236B Easy Number Challenge（筛法求素数 + 整数因式分解）
题目链接:http://codeforces.com/problemset/problem/236/B Easy Number Challenge time limit per test 2 seco ...
1042: 筛法求素数
1042: 筛法求素数 Time Limit: 1 Sec Memory Limit: 128 MB Submit: 1387 Solved: 918 [Submit][Status][Web B ...
素数筛法求素数（线性时间）
摘自:http://blog.csdn.net/once_hnu/article/details/6302283 逆向思维的典型应用! 1)普通方法判断素数:O(n*sqrt(n)) 2)筛法求素数: ...
数论判断素数：普通素数判别线性筛二次筛法求素数米勒拉宾素数检验
普通的素数判断法当我们要判断一个数字是否是素数的时候,往往会直接看这个数字模1到这个数字的根号,看有没有等于零的,从而判断这个数字是不是素数,这样做的时间复杂度为O(sqrt(n)) bool is ...
筛法求素数 (20分)
题干: 素数是仅仅能被它本身和1整除的任何整数.筛法求素数是一种查找素数的方法.它的算法如下: 1.创建一个数组,并将所有元素初始化为1(真).具有素数下标的数组元素将保持为1,而其它数组元素最终将被 ...

模板 - Min_25 筛法求素数和

目录

模板 - Min_25 筛法求素数和相关推荐

最新文章

热门文章