与Min_25筛有关的一些模板

最新推荐文章于 2022-07-01 17:32:12 发布
最新推荐文章于 2022-07-01 17:32:12 发布
阅读量241 收藏 1
点赞数
分类专栏： 数论
版权声明：本文为博主原创文章，遵循 CC 4.0 BY 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_45483201/article/details/108761930
版权

模板

求 ∑ i = 1 n f ( i ) , f ( p k ) = p k × ( p k − 1 ) \sum \limits_{i = 1} ^{n} f(i), f(p ^ k) = p ^ k \times(p ^ k - 1) i=1∑n​f(i),f(pk)=pk×(pk−1)，最后对   m o d   1 e 9 + 7 \bmod 1e9 + 7 mod1e9+7，这个函数是个积性函数。

/*
  Author : lifehappy
*/
#pragma GCC optimize(2)
#pragma GCC optimize(3)
#include <bits/stdc++.h>

#define mp make_pair
#define pb push_back
#define endl '\n'
#define mid (l + r >> 1)
#define lson rt << 1, l, mid
#define rson rt << 1 | 1, mid + 1, r
#define ls rt << 1
#define rs rt << 1 | 1

using namespace std;

typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> pii;

const double pi = acos(-1.0);
const double eps = 1e-7;
const int inf = 0x3f3f3f3f;

inline ll read() {
    ll f = 1, x = 0;
    char c = getchar();
    while(c < '0' || c > '9') {
        if(c == '-')    f = -1;
        c = getchar();
    }
    while(c >= '0' && c <= '9') {
        x = (x << 1) + (x << 3) + (c ^ 48);
        c = getchar();
    }
    return f * x;
}

const int N = 1e6 + 10, mod = 1e9 + 7, inv2 = 500000004, inv6 = 166666668;

namespace MIN_25 {
    int prime[N], id1[N], id2[N], cnt, m, T;

    ll g1[N], g2[N], sum1[N], sum2[N], a[N], n;

    bool st[N];

    int ID(ll x) {
        return x <= T ? id1[x] : id2[n / x];
    }

    ll calc2(ll x) {
        x %= mod;
        return (x * (x + 1) % mod * inv2 % mod - 1 + mod) % mod;
    }
    
    ll calc1(ll x) {
        x %= mod;
        return (x * (x + 1) % mod * (2 * x + 1) % mod * inv6 % mod - 1 + mod) % mod;
    }

    ll f(ll x) {
        x %= mod;
        return x * (x - 1) % mod;
    }

    void init() {
        T = sqrt(n + 0.5);
        for(int i = 2; i <= T; i++) {
            if(!st[i]) {
                prime[++cnt] = i;
                sum1[cnt] = (sum1[cnt - 1] + 1ll * i * i) % mod;
                sum2[cnt] = (sum2[cnt - 1] + i) % mod;
            }
            for(int j = 1; j <= cnt && 1ll * i * prime[j] <= T; j++) {
                st[i * prime[j]] = 1;
                if(i % prime[j] == 0) {
                    break;
                }
            }
        }
        for(ll l = 1, r; l <= n; l = r + 1) {
            r = n / (n / l);
            a[++m] = n / l;
            if(a[m] <= T) id1[a[m]] = m;
            else id2[n / a[m]] = m;
            g1[m] = calc1(a[m]);
            g2[m] = calc2(a[m]);
        }
        for(int j = 1; j <= cnt; j++) {
            for(int i = 1; i <= m && 1ll * prime[j] * prime[j] <= a[i]; i++) {
                g1[i] = ((g1[i] - 1ll * prime[j] * prime[j] % mod * (g1[ID(a[i] / prime[j])] - sum1[j - 1]) % mod) % mod + mod) % mod;
                g2[i] = ((g2[i] - 1ll * prime[j] * (g2[ID(a[i] / prime[j])] - sum2[j - 1]) % mod) % mod + mod) % mod;
            }
        }
    }

    ll solve(ll n, int m) {
        if(n < prime[m]) return 0;
        ll ans = ((g1[ID(n)] - sum1[m - 1] - g2[ID(n)] + sum2[m - 1] % mod) + mod) % mod;
        for(int j = m; j <= cnt && prime[j] * prime[j] <= n; j++) {
            for(ll i = prime[j]; i * prime[j] <= n; i *= prime[j]) {
                ans = (ans + f(i) * solve(n / i, j + 1) % mod + f(i * prime[j])) % mod;
            }
        }
        return ans;
    }

    ll solve(ll x) {
        n = x;
        init();
        return solve(n, 1) + 1;
    }
}

int main() {
    // freopen("in.txt", "r", stdin);
    // freopen("out.txt", "w", stdout);
    // ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
    cout << MIN_25::solve(read()) << endl;
    return 0;
}

求 ∑ i = 1 n i [ i ∈ p r i m e s ] \sum\limits_{i = 1} ^{n} i[i \in primes] i=1∑n​i[i∈primes]

即是求 f ( p ) = p f(p) = p f(p)=p，直接求它的 g g g函数即可。

/*
  Author : lifehappy
*/
#pragma GCC optimize(2)
#pragma GCC optimize(3)
#include <bits/stdc++.h>

#define mp make_pair
#define pb push_back
#define endl '\n'
#define mid (l + r >> 1)
#define lson rt << 1, l, mid
#define rson rt << 1 | 1, mid + 1, r
#define ls rt << 1
#define rs rt << 1 | 1

using namespace std;

typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> pii;

const double pi = acos(-1.0);
const double eps = 1e-7;
const int inf = 0x3f3f3f3f;

inline ll read() {
    ll f = 1, x = 0;
    char c = getchar();
    while(c < '0' || c > '9') {
        if(c == '-')    f = -1;
        c = getchar();
    }
    while(c >= '0' && c <= '9') {
        x = (x << 1) + (x << 3) + (c ^ 48);
        c = getchar();
    }
    return f * x;
}

const int N = 1e6 + 10, mod = 1e9 + 7, inv2 = 500000004, inv6 = 166666668;

int prime[N], id1[N], id2[N], cnt, m, T;

ll g[N], sum[N], a[N], n;

bool st[N];

int ID(ll x) {
    return x <= T ? id1[x] : id2[n / x];
}

ll calc(ll x) {
    return x * (x + 1) / 2 - 1;
}

void init() {
    T = sqrt(n + 0.5);
    for(int i = 2; i <= T; i++) {
        if(!st[i]) {
            prime[++cnt] = i;
            sum[cnt] = sum[cnt - 1] + i;
        }
        for(int j = 1; j <= cnt && 1ll * i * prime[j] <= T; j++) {
            st[i * prime[j]] = 1;
            if(i % prime[j] == 0) {
                break;
            }
        }
    }
    for(ll l = 1, r; l <= n; l = r + 1) {
        r = n / (n / l);
        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 j = 1; j <= cnt; j++) {
        for(int i = 1; i <= m && 1ll * prime[j] * prime[j] <= a[i]; i++) {
            g[i] = g[i] - 1ll * prime[j] * (g[ID(a[i] / prime[j])] - sum[j - 1]);
        }
    }
}

ll solve(ll x) {
    n = x;
    init();
    return g[ID(x)];
}

int main() {
    // freopen("in.txt", "r", stdin);
    // freopen("out.txt", "w", stdout);
    // ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
    cout << solve(read()) << endl;
    return 0;
}

#6053. 简单的函数

f ( 1 ) = 1 f ( p c ) = p ⊕ c f ( a b ) = f ( a ) ( b )   a , b 互 质 求 ∑ i = 1 n f ( i ) f(1) = 1\\ f(p ^ c) = p \oplus c\\ f(ab) = f(a)(b)\ a, b互质\\ 求\sum_{i = 1} ^{n} f(i)\\ f(1)=1f(pc)=p⊕cf(ab)=f(a)(b) a,b互质求i=1∑n​f(i)

/*
  Author : lifehappy
*/
#pragma GCC optimize(2)
#pragma GCC optimize(3)
#include <bits/stdc++.h>

#define mp make_pair
#define pb push_back
#define endl '\n'
#define mid (l + r >> 1)
#define lson rt << 1, l, mid
#define rson rt << 1 | 1, mid + 1, r
#define ls rt << 1
#define rs rt << 1 | 1

using namespace std;

typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> pii;

const double pi = acos(-1.0);
const double eps = 1e-7;
const int inf = 0x3f3f3f3f;

inline ll read() {
    ll f = 1, x = 0;
    char c = getchar();
    while(c < '0' || c > '9') {
        if(c == '-')    f = -1;
        c = getchar();
    }
    while(c >= '0' && c <= '9') {
        x = (x << 1) + (x << 3) + (c ^ 48);
        c = getchar();
    }
    return f * x;
}

const int N = 1e6 + 10, mod = 1e9 + 7, inv2 = 500000004, inv6 = 166666668;

namespace MIN_25 {
    int prime[N], id1[N], id2[N], cnt, m, T;

    ll g1[N], g2[N], sum1[N], sum2[N], a[N], n;

    bool st[N];

    int ID(ll x) {
        return x <= T ? id1[x] : id2[n / x];
    }

    ll calc(ll x) {
        x %= mod;
        return (x * (x + 1) % mod * inv2 % mod - 1 + mod) % mod;
    }

    void init() {
        T = sqrt(n + 0.5);
        for(int i = 2; i <= T; i++) {
            if(!st[i]) {
                prime[++cnt] = i;
                sum1[cnt] = (sum1[cnt - 1] + i) % mod;
                sum2[cnt] = (sum2[cnt - 1] + 1) % mod;
            }
            for(int j = 1; j <= cnt && 1ll * i * prime[j] <= T; j++) {
                st[i * prime[j]] = 1;
                if(i % prime[j] == 0) {
                    break;
                }
            }
        }
        for(ll l = 1, r; l <= n; l = r + 1) {
            r = n / (n / l);
            a[++m] = n / l;
            if(a[m] <= T) id1[a[m]] = m;
            else id2[n / a[m]] = m;
            g1[m] = calc(a[m]);
            g2[m] = (a[m] - 1) % mod;
        }
        for(int j = 1; j <= cnt; j++) {
            for(int i = 1; i <= m && 1ll * prime[j] * prime[j] <= a[i]; i++) {
                g1[i] = ((g1[i] - 1ll * prime[j] * (g1[ID(a[i] / prime[j])] - sum1[j - 1]) % mod) % mod + mod) % mod;
                g2[i] = ((g2[i] - 1ll * (g2[ID(a[i] / prime[j])] - sum2[j - 1]) % mod) % mod + mod) % mod;
            }
        }
    }

    ll solve(ll n, int m) {
        if(n < prime[m]) return 0;
        ll ans = ((g1[ID(n)] - sum1[m - 1] - g2[ID(n)] + sum2[m - 1]) % mod + mod) % mod;
        for(int j = m; j <= cnt && 1ll * prime[j] * prime[j] <= n; j++) {
            for(ll i = prime[j], c = 1; i * prime[j] <= n; i *= prime[j], c++) {
                ans = (ans + 1ll * (prime[j] ^ c) * solve(n / i, j + 1) % mod + 1ll * (prime[j] ^ (c + 1))) % mod;
            }
        }
        return ans;
    }

    ll solve(ll x) {
        n = x;
        init();
        return (solve(x, 1) + 1 + (x > 1) * 2) % mod;
    }
}

int main() {
    // freopen("in.txt", "r", stdin);
    // freopen("out.txt", "w", stdout);
    // ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
    cout << MIN_25::solve(read()) << endl;
    return 0;
}

#6235. 区间素数个数

/*
  Author : lifehappy
*/
#pragma GCC optimize(2)
#pragma GCC optimize(3)
#include <bits/stdc++.h>

#define mp make_pair
#define pb push_back
#define endl '\n'
#define mid (l + r >> 1)
#define lson rt << 1, l, mid
#define rson rt << 1 | 1, mid + 1, r
#define ls rt << 1
#define rs rt << 1 | 1

using namespace std;

typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> pii;

const double pi = acos(-1.0);
const double eps = 1e-7;
const int inf = 0x3f3f3f3f;

inline ll read() {
    ll f = 1, x = 0;
    char c = getchar();
    while(c < '0' || c > '9') {
        if(c == '-')    f = -1;
        c = getchar();
    }
    while(c >= '0' && c <= '9') {
        x = (x << 1) + (x << 3) + (c ^ 48);
        c = getchar();
    }
    return f * x;
}

const int N = 1e6 + 10, mod = 1e9 + 7, inv2 = 500000004, inv6 = 166666668;

namespace MIN_25 {
    int prime[N], id1[N], id2[N], m, cnt, T;

    ll a[N], g[N], sum[N], n;

    bool st[N];

    int ID(ll x) {
        return x <= T ? id1[x] : id2[n / x];
    }

    void init() {
        T = sqrt(n + 0.5);
        for(int i = 2; i <= T; i++) {
            if(!st[i]) {
                prime[++cnt] = i;
                sum[cnt] = sum[cnt - 1] + 1;
            }
            for(int j = 1; j <= cnt && 1ll * i * prime[j] <= T; j++) {
                st[i * prime[j]] = 1;
                if(i % prime[j] == 0) {
                    break;
                }
            }
        }
        for(ll l = 1, r; l <= n; l = r + 1) {
            r = n / (n / l);
            a[++m] = n / l;
            if(a[m] <= T) id1[a[m]] = m;
            else id2[n / a[m]] = m;
            g[m] = a[m] - 1;
        }
        for(int j = 1; j <= cnt; j++) {
            for(int i = 1; i <= m && 1ll * prime[j] * prime[j] <= a[i]; i++) {
                g[i] = g[i] - 1ll * (g[ID(a[i] / prime[j])] - sum[j - 1]);
            }
        }
    }

    ll solve(ll x) {
        n = x;
        init();
        return g[ID(n)];
    }
}

int main() {
    // freopen("in.txt", "r", stdin);
    // freopen("out.txt", "w", stdout);
    // ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
    cout << MIN_25::solve(read()) << endl;
    return 0;
}
_lifehappy_
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
专栏目录
Min-25筛超实用型模板
Shu_B的博客
08-11 957
前言 原本以为会了杜教筛就能走遍积性函数题的天下都不怕的了，结果5月省赛时遇到了专门卡杜教筛的题，只好滚回去学了…结果省赛之后一直在忙这忙那的，然后一拖拖到了暑假的现在才有时间来补Min-25筛的内容。 复杂度 O(n34log⁡n)O(\frac{n^\frac{3}{4}}{\log n})O(lognn43​​)，而且据说常数非常优秀 证明？那是什么，能吃吗? 相关链接 Min-25筛原理 ...
模板 - min25筛
weixin_30882895的博客
07-30 209
好像在某些情况下杜教筛会遇到瓶颈，先看着。暑假学一些和队友交错的知识的同时开这个大坑。 2019/7/30 求一个前缀和 $\sum\limits_{i=1}^n f(i) $ ，其中 \(f(x)\) 是积性函数，且 \(f(p^k)\) 是一个关于 \(p\) 的低次多项式。 #include<cstdio> #include<algorithm> #incl...
Min_25筛 模板
Harris-H的博客
09-25 376
Min_25筛 模板 题目传送门 现在还不是完全明白，意会一下。 #include<bits/stdc++.h> using namespace std; #define il inline typedef long long ll; const int N=2e5+5,mod=1e9+7; bool isp[N];//判断是否为质数 int cnt,m,id1[N],id2[N]; //cnt质数的个数,m待处理的数个数 //id1[]储存<=sq的m下标, id2[]储存>=s
【模板】Min_25筛
Brute♂force
09-13 300
题目 思路 我曾以为线性筛就够了，直到我看到n<=1e10 我曾以为杜教筛就够了，直到我发现根本不知道卷什么函数 于是，我就来学min25筛了 我曾以为min25筛就够了…… min25筛其实有些像DP，我们要做的就是用预处理的东西转移求出前缀和 这里推荐https://www.cnblogs.com/zhoushuyu/p/9187319.html 讲的很透彻 代码 #include<bits/stdc++.h> #define ll long long using namespace
Min_25筛 代码及模板
KXL5180的博客
01-19 577
LOJ简单的函数 以这道题为例子，讲一下min_25筛如何用代码实现。 首先min_25是一种亚线性筛，可以处理1e9以上的数据，老版min_25筛的复杂度为O(n0.75logn)O(\frac{n^{0.75}}{logn})O(lognn0.75​)，新版已经优化到O(n23)O(n^{\frac{2}{3}})O(n32​)，复杂度，我也不会证明。 据说常数比较大，确实是跑得快，但不过差距...
ACM巨全模板 .pdf
10-07
1.循环排列 (与环有关的排列组合) 计算几何： 1.三角形 (求面积)) 2.多边形 3.三点求圆心和半径 4.扫描线 (矩形覆盖求面积) (矩形覆盖求周长) 5.凸包 (平面上最远点对) 6.求凸多边形的直径 7.求凸多边形的宽度 8.求...
常用算法模板库(C++)
03-08
### 常用算法模板库(C++) #### 目录 1. **排序算法** 1. **冒泡排序** 2. **选择排序** 3. **插入排序** 4. **快速排序** 5. **哈希排序** 2. **数学问题** 1. **求最大公约数最小公倍数** 2. **求素数** ...
数论模板精心整编
07-10
根据提供的文档信息，我们可以整理出一系列与数论相关的知识点，并对每个部分进行详细的解析和解释。 ### 1. 负数取模 在计算机科学中，取模运算通常用于获取除法运算的余数。当涉及到负数时，取模操作可能会有...
表面活性剂模板剂应用常识及实验经验.doc
10-02
表面活性剂作为模板剂在材料科学中的应用广泛，特别是在合成介孔材料领域。模板剂的选择和使用方式直接影响最终产品的结构和性能。以下是基于文档内容的详细解释： 1. **模板剂的溶解**： 高分子类的表面活性剂，...
P5325 【模板】Min_25筛（min25筛）
FLY_WHITE的博客
05-24 541
P5325 【模板】Min_25筛（min25筛） 题目大意 定义积性函数f(x)，且f(pk)=pk(pk−1)f(p^k)=p^k(p^k-1)f(pk)=pk(pk−1)（pp是一个质数），求 ∑i=1nf(x)​ \sum_{i=1}^n f(x)​ i=1∑n​f(x)​ 对10^9+7取模。 解题思路 min25筛版题https://blog.csdn.net/baiyifeifei/...
【洛谷P5325】 【模板】Min_25筛
Coldfresh的博客
05-31 484
n的范围是小于101010^{10}1010 抛弃ees和洲阁筛，拥抱min_25 代码： #include<cstdio> #include<cstring> #include<iostream> #include<cmath> #include<vector> #include<algorithm> #define IN...
min_25筛模板及应用 / loj6053 / loj6235
m0_54646258的博客
07-01 206
讲述算法的实现以及应用场景
【数论】Min_25筛
05-16 4485
听说这玩意玩爆洲阁筛？？？ Min_25筛可以解决一类积性函数求和问题， 筛质数 假设我们现在要对n以内的质数求和， 先线筛出小于n−−√n\sqrt n的所有质数，设第i个为pipip_i，共有p0p0p0个质数，pspsps为p的前缀和， 设函数：S(x,j)=∑xi=2i∗[i为质数&amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;nbsp;或&amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;
LOJ #6053. 简单的函数 （min25筛裸题）
Ark
01-15 594
题面 n≤1010n\le10^{10}n≤1010 题解 裸的min25筛，要学习min25筛的戳这里 我主要用的是下面那种写法（比较好写） CODE #include <bits/stdc++.h> using namespace std; #define int long long const int N = 200005; const int mod = 1e9 + 7; ...
模板 - Min_25 筛法求素数和
繁凡さん的博客
10-07 1161
不能算是Min_25筛，该模板只是在较短时间内求1e11~1e12内的素数和 #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 ? i.
Min_25筛学习小记
beginend
03-24 1755
前言 Min_25筛的时间复杂度和洲阁筛（我不会）一样，都是O(n34logn)O(n34logn)O(\frac{n^{\frac{3}{4}}}{logn})。但Min_25筛的时空常数，代码复杂度据说比洲阁筛要优秀，甚至可以替代掉洲阁筛，于是我这个学不懂洲阁筛的菜鸡就赶紧跑来学习了一发。 我看的是zzq的博客 目标 设F(x)F(x)F(x)是一个积性函数，现在要求∑i=1nF(...
min25筛学习笔记&模板详解
FLY_WHITE的博客
05-22 3182
min25筛算法推导&模板详解 问题引入 引入函数σ0(n)\sigma _0(n)σ0​(n)为n的正因数的数量，求S(n,k)=∑i=1nσ0(ik)S(n,k)=\sum_{i=1}^n\sigma_0(i^k)S(n,k)=∑i=1n​σ0​(ik),n,k≤1010n,k\le 10^{10}n,k≤1010 简化分析与引入 对于σ0(n)\sigma_0(n)σ0​(n)函数，...
【学习笔记】Min25筛
cz_xuyixuan的博客
08-09 7397
【定理简介】 Min25Min25Min25筛是一种能够求解积性函数f(x)f(x)f(x)的前缀和∑Ni=1f(i)∑i=1Nf(i)\sum_{i=1}^{N}f(i)的筛法，其前提条件为∑Ni=1[i&amp;nbsp;is&amp;nbsp;a&amp;nbsp;prime]∗f(i)∑i=1N[i&amp;nbsp;is&amp;nbsp;a&amp;nbsp;prime]∗f(i)\sum_{i=1}^{N}[i\ i...
% 获取运动目标模板 template = gray_frame(min_row:max_row, min_col:max_col); % 缩小模板 resized_template = imresize(template, 0.5);min_row:max_row, min_col:max_col怎么获得
最新发布
05-26
这里的 `min_row`, `max_row`, `min_col`, `max_col` 可以通过目标物体的位置信息获得。假设你使用了一个物体检测器来检测运动物体，那么物体检测器会返回物体的位置信息，包括物体在图像中的左上角和右下角坐标 `...
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值