杜教筛学习报告(随训练更新)

最新推荐文章于 2023-05-10 13:48:04 发布
最新推荐文章于 2023-05-10 13:48:04 发布
阅读量155 收藏 1
点赞数
分类专栏： 数论 集训或补题
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_41791981/article/details/98030489
版权

目录

关于杜教筛的简述

训练记录：

51nod1244 莫比乌斯函数之和

51nod1239 欧拉函数之和

bzoj3944 sum

hdu5608 Function

 

 

关于杜教筛的简述

看了skywalkert的博客大概明白了。

author: skywalkert 
original article: http://blog.csdn.net/skywalkert/article/details/50500009 
last update time : 2017-04-01

大体思路是如果能通过狄利克雷卷积构造一个更好计算前缀和的函数，且用于卷积的另一个函数也易计算，则可以简化计算过程。

构造函数后前n^ \frac{2}{3} 项的前缀和预处理出来，剩下的部分通过前面构造的函数进行递归。大概构造成：f(n) = g(i) - \sum_{i = 2}^{n}f(\frac{n}{i})

这种形式，f(n)是待求前缀和，g(n)是构造的函数。然后进行递归。

特别注意！！！由于亚线性筛的复杂度低于O(n)所以一般数据范围是1e8-1e11，取模！！！一定注意取模！！！

前置技能：莫比乌斯反演，常用积性函数，线性筛

训练记录：

51nod1244 莫比乌斯函数之和

题意：求莫比乌斯函数的前缀和。

思路：由于数据范围是1e10，显然不能采用线性做法，采取杜教筛，构造一个函数满足f(n) = \sum_{i = 1}^{n}g(i) - \sum_{i = 2}^{n}f(\frac{n}{i})，注意到莫比乌斯函数的性质[n==1]=\sum_{d|n}^{n}\mu(d),所以有1=\sum_{i = 1}^{n}[n==1]=\sum_{i=1}^{n}\sum_{d|n}^{n}\mu(d)=\sum_{i=1}^{n}\sum_{d=1}^{\lfloor \frac{n}{i} \rfloor}\mu(d),于是乎，我们构造g(n)=1.

然后式子就成了f(n)=g(n)-\sum_{i=2}^{n}f(\frac{n}{i})=1-\sum_{i=2}^{n}f(\frac{n}{i}),然后预处理一下4e6范围的莫比乌斯函数的前缀和递归就完了，中间用map记录一下已计算的项减低复杂度。

复杂度：O(n^\frac{2}{3})

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 4641600;
const int inf = 0x3f3f3f3f;
const int mod = 1e6 + 3;
//const double pi = acos(-1.0);
const double eps = 1e-6;
const int inv = 5e8 + 4;
int gcd(int a,int b){return b ? gcd(b,a % b) : a;}
ll n,m;
int mu[N],prime[N],cnt;
bool vis[N];
inline ll read()
{
    ll x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0'; ch=getchar();}
    return x*f;
}
void Mobius()
{
    cnt = 0;
    mu[1] = 1;mu[0] = 0;
    for(int i = 2; i < N; ++i){
        if(!vis[i]){
            prime[++cnt] = i;
            mu[i] = -1;
        }
        for(int j = 1; j <= cnt; ++j){
            if(i * prime[j] >= N) break;
            vis[i * prime[j]] = true;
            if(i % prime[j]== 0){
                mu[i * prime[j]] = 0;
                break;
            }
            mu[i * prime[j]] = -mu[i];
        }
        mu[i] += mu[i - 1];
    }
}
map<ll ,int>mp;
int work(ll jq)
{
    if(jq < N) return mu[jq];
    if(mp[jq]) return mp[jq];
    int res = 1;ll last;
    for(ll i = 2;i <= jq;i = last + 1){
        last = jq / (jq / i);
        res -= (last - i + 1) * work(jq / i);
    }
    mp[jq] = res;
    return res;
}
int main()
{
    Mobius();//cout << pow((ll)10000000000,2.0/ 3);
    n = read(),m = read();
    printf("%d\n",work(m) - work(n - 1));
    return 0;
}

51nod1239 欧拉函数之和

题意：求欧拉函数前缀和

思路：注意到欧拉函数的性质n=\sum_{n|d}^{n}\phi(d)=\sum_{d = 1}^{\lfloor \frac{n}{d} \rfloor}\phi(d),所以有\sum_{i=1}^{n}i=\sum_{i=1}^{n}\sum_{d=1}^{\lfloor \frac{n}{\frac{i}{d}} \rfloor }\phi(d)=\sum_{i=1}^{n}\phi(\lfloor \frac{n}{i} \rfloor ),所以构造g(n)=\sum_{i=1}^{n}i,所以就变成了求解f(n)=g(n)-\sum_{i=2}^{n}\phi(\lfloor \frac{n}{i} \rfloor ).

复杂度：O(n^\frac{2}{3})

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 4641600;
const int inf = 0x3f3f3f3f;
const int mod = 1e9 + 7;
//const double pi = acos(-1.0);
const double eps = 1e-6;
const int inv = 5e8 + 4;
int gcd(int a,int b){return b ? gcd(b,a % b) : a;}
ll phi[N];
int prime[N],cnt;
bool isprime[N];
ll n;
inline ll add(ll x)
{
    return x >= mod ? x - mod : x;
}
inline ll sub(ll x)
{
    return x < 0 ? x + mod : x;
}
void get_phi()
{
   phi[1] = 1;
   for(int i = 2;i < N;++i){
       if(!isprime[i]){
             prime[++cnt] = i;
             phi[i] = i - 1;
        }
       for(int j = 1;j <= cnt && i * prime[j] < N;++j)
       {
          isprime[i * prime[j]] = 1;
          if(i % prime[j] == 0){
             phi[i * prime[j]] = phi[i] * prime[j];
             break;
          }
          else phi[i * prime[j]] = phi[i] * (prime[j] - 1);
       }
   }
   for(int i = 2;i < N;++i) phi[i] = add(phi[i] + phi[i - 1]);
}

map<ll ,ll>mp;
ll work(ll jq)
{
    if(jq < N) return phi[jq];
    if(mp[jq]) return mp[jq];
    ll res = jq % mod * ((jq + 1) % mod) % mod * inv % mod,last;
    for(ll i = 2;i <= jq;i = last + 1){
        last = jq / (jq / i);
        res = sub(res - (last - i + 1) * work(jq / i) % mod);
    }
    mp[jq] = res;
    return res;
}
int main()
{
    get_phi();
    scanf("%lld",&n);
    printf("%lld\n",work(n));
    return 0;
}

bzoj3944 sum

题意：求欧拉前缀和和莫比乌斯函数前缀和

思路：同上，注意开结构体避免栈溢出！！！

#pragma comment(linker, “/STACK:1024000000,1024000000”
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 5e6 + 7;
const int inf = 0x3f3f3f3f;
const int mod = 1e9 + 7;
//const double pi = acos(-1.0);
const double eps = 1e-6;
const int inv = 5e8 + 4;
int gcd(int a,int b){return b ? gcd(b,a % b) : a;}
ll phi[N];
int prime[N],cnt,mu[N];
bool isprime[N];
ll n;
void get_phi()
{
   phi[1] = 1;mu[1] = 1;mu[0] = 0;
   for(int i = 2;i < N - 2;++i){
       if(!isprime[i]){
             prime[++cnt] = i;
             phi[i] = i - 1;
             mu[i] = -1;
        }
       for(int j = 1;j <= cnt && i * prime[j] < N - 2;++j){
          isprime[i * prime[j]] = 1;
          if(i % prime[j] == 0){
             phi[i * prime[j]] = phi[i] * prime[j];
             mu[i * prime[j]] = 0;
             break;
          }
          phi[i * prime[j]] = phi[i] * (prime[j] - 1);
          mu[i * prime[j]] = -mu[i];
       }
       mu[i] += mu[i - 1];
   }
   for(int i = 2;i < N - 2;++i) phi[i] = phi[i] + phi[i - 1];
}
struct node{
    ll p;
    int m;
};
map<int ,ll>mp,mp1;
node work(ll jq)
{
    if(jq < N - 3) return (node){phi[jq],mu[jq]};
    if(mp[jq]) return (node){mp[jq],mp1[jq]};
    node res = (node){(jq & 1 ? jq : jq / 2) * (jq & 1 ? (jq + 1) / 2 : jq + 1),1},tmp;
    ll last;
    for(ll i = 2;i <= jq;i = last + 1){
        last = jq / (jq / i),tmp = work(jq / i);
        res.p -= (last - i + 1) * tmp.p;
        res.m -= (last - i + 1) * tmp.m;
    }
    if(!mp[jq]) mp[jq] = res.p;
    if(!mp1[jq]) mp1[jq] = res.m;
    return res;
}
int main()
{
    get_phi();
    int t;
    scanf("%d",&t);
    while(t--){
        scanf("%lld",&n);
        if(!n){
            printf("0 0\n");
            continue;
        }
        node r = work(n);
        printf("%lld %d\n",r.p,r.m);
    }
    return 0;
}

hdu5608 Function

题意：已知f(x)满足n^2-3n+2=\sum_{d|n}^{n}f(d),F(x)=\sum_{i=1}^{n}f(n)求F(x).

思路：看都整出想到莫比乌斯反演。令f(n) = \sum_{d|n}^{n}g(d)*\mu(\frac{n}{d}),所以有F(n)=\sum_{i=1}^{n}f(n)=\sum_{i=1}^{n}\sum_{d=1}^{\lfloor \frac{n}{d} \rfloor }g(i)*\mu(d),由前面莫比乌斯推的经验，我们直接构造G(n)=\sum_{i=1}^{n}g(x)*1,是不是突然清晰明了呢，这题就解决了。

//#pragma comment(linker, “/STACK:1024000000,1024000000”
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 1e6 + 7;
const int inf = 0x3f3f3f3f;
const int mod = 1e9 + 7;
const int inv2 = 5e8 + 4;
//const double pi = acos(-1.0);
const double eps = 1e-6;
const int inv = 333333336;
int gcd(int a,int b){return b ? gcd(b,a % b) : a;}
int prime[N],cnt,mu[N];
ll g[N];
bool isprime[N];
ll n;
inline ll gao(ll x)
{
    return (x - 1) * (x - 2) % mod * x % mod * inv % mod;
}
void get_phi()
{
   mu[1] = 1;
   for(int i = 2;i < N;++i){
       if(!isprime[i]){
             prime[++cnt] = i;
             mu[i] = -1;
        }
       for(int j = 1;j <= cnt && i * prime[j] < N;++j){
          isprime[i * prime[j]] = 1;
          if(i % prime[j] == 0){
             mu[i * prime[j]] = 0;
             break;
          }
          mu[i * prime[j]] = -mu[i];
       }
       //mu[i] += mu[i - 1];
   }
   for(int i = 1;i < N;++i)
       for(int j = i;j < N;j += i)
           g[j] = ((g[j] + (ll)mu[j / i] * (i - 1) % mod * (i - 2) % mod) % mod + mod) % mod;
    for(int i = 2;i < N;++i) g[i] = (g[i] + g[i - 1]) % mod;
}
map<ll ,ll>mp;
ll solve(ll pos)
{
    if(pos < N) return g[pos];
    if(mp[pos]) return mp[pos];
    ll res = gao(pos),last;
    for(ll i = 2;i <= pos;i = last + 1){
        last = pos / (pos / i);
        res = ((res - (last - i + 1) * solve(pos / i) % mod) % mod + mod) % mod;
    }
    mp[pos] = res;
    return res;
}
int main()
{
    get_phi();
    int t;
    scanf("%d",&t);
    while(t--){
        scanf("%lld",&n);
        printf("%lld\n",solve(n));
    }
    return 0;
}

 

xiange_hu
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
把ICDAR2019-LSVT原始数据集转为REC任务识别数据
Andrwin的博客
10-09 6444
ICDAR2019-LSVT数据集有两个： ①Baidu Research Open-Access Dataset - Introductionhttps://ai.baidu.com/broad/introduction?dataset=lsvt 是包含定位和标记的数据类型 ②飞桨AI Studio - 人工智能学习与实训社区集开放数据、开源算法、免费算力三位一体，为开发者提供高效学习和开发环境、高价值高奖金竞赛项目，支撑高校老师轻松实现AI教学，并助力开发者学习交流，加速落地AI业务场景ht.
杜教筛的某道练习题
weixin_30305735的博客
12-11 112
from here. \(Description\) 求\[\sum_{i=1}^n\varphi(i)*i\] \(Solution\) 把\(f(x)=\sum_{i=1}^n\varphi(i)*i\)与\(g(x)=x\)做狄利克雷卷积并求前缀和，会得到：\[\sum_{i=1}^n\sum_{d|i}f(d)g(\frac id)=\sum_{i=1}^n\sum_{d|i}\varph...
杜教筛学习
auu55357的博客
03-22 88
转载于:https://www.cnblogs.com/zzyer/p/8626741.html
杜教筛详解+入门题
weixin_30820077的博客
07-31 391
首先看看大佬的博客 下面是推理过程，然而我们只需要记结论（h，g，f 都是假定的函数，S是f的前缀和函数）： 例题1：N=1e10 ，求下式 思路：S(i)是f的前缀和，我们想找到一个g使得h的前缀和好求（h=f*g）。 先将h的式子写出来： 再把f套成u，我们会发现，前面那一部分是可以套莫比乌斯的性质的： 那么把g带成1函数，f*g就等于元函数，而元函数...
深度学习-VGG19网络识别手写汉字（基于HWDB1.1trn_gnt）
xingxingzhilong的专栏
05-10 2243
这段代码呢，使用了迭代器，能够减少内存的使用，不会一下次吃掉所有内存（不明白迭代器的小伙伴，自行百度学习）。图片按照上面六个顺序，进行处理，当然了，我这里还不够全面，比如还可以添加，检测图片是否摆正，是否需要翻转等等。，图片容易受到外在因素的干扰，比如有噪点，有色彩，有缺失等等一系列的问题。我罗列了一下我的处理过程。好的，做完这些之后，我们接下来就到正餐了，训练需要的样本都准备好了，接下来，我们要把我们的人工网络搭建起来。说明下，这两个字不是测试集和训练集当中的，是找两个人写的，识别的精准度还是很高的。
【长更】一句话题解（组队训练的俄罗斯题、oj、camp）
热门推荐
Four's
04-15 2万+
还是太长了，第二次分裂。。
NOI前训练记录
weixin_30614109的博客
05-21 127
从今年年初开始就没怎么碰oi，学了三个月文化课，然后就是省选到处跑(浪)，进了省队后又去北京打(浪)了波(七)铁(天)，回家后又生了一星期病，感觉自己noi凉透了。。。 ctsc因为运气的原因有人放弃D3自己才拿到au，apio什么牌都没拿到，自己这么菜的主要原因可能还是没怎么做题吧，五个月不碰oi，已经连dfs都不会了，还有两个月就noi了，再不训练又要为JS丢脸了，自己拿不到au却占了个省队名...
日常（更新至2019.8.14）
afwrri2311的博客
02-19 295
八月 Camp完之后就开始HDU多校自闭之旅了 第一场连签到题都是网络流这种级别的orz 因为这几场比赛 强行学习了线性基、划分树、支配树和一些奇奇怪怪的可持久化数据结构 算是颇有收获了 还有两场，坚持！ 说起网络流，我发现自己最小费用最大流多路增广的部分还没实现 以及最近才开始看的带上下界网络流也是个坑 之前还说想写一篇Blog 现在看来得推迟了 本月除了学了一大堆算法与数据结构外 最...
tf 文字识别 多个文字 data 数据流
阿豪
06-14 1999
原文链接: tf 文字识别 多个文字 data 数据流 ...
杜云海ARM学习报告
01-19
《杜云海ARM学习报告》是一份由资深技术专家杜云海编写的关于ARM架构学习的宝贵资料，特别适合ARM技术的初学者。ARM（Advanced RISC Machines）是一种广泛应用于移动设备、嵌入式系统以及服务器等领域的处理器架构，...
算法竞赛专题解析--杜教筛1
08-04
杜教筛是一种高效计算积性函数的算法，尤其在处理较大规模数据时具有显著优势。在算法竞赛中，掌握杜教筛能帮助参赛者解决涉及数论和组合优化的问题。本解析将深入介绍积性函数及其相关概念，并通过实例演示杜教筛的...
ACM算法竞赛中的杜教筛的C++代码实现
最新发布
02-12
杜教筛是一种用于解决数论问题的算法。它主要用于计算在给定区间内数的质因数个数之和。该算法的基本思想是结合了区间筛和积性函数的性质，在一定范围内高效地计算出积性函数的前缀和。 具体来说，杜教筛的步骤包括...
ARM学习报告
11-29
在本文中，我们将深入探讨杜云海的ARM学习报告，这份报告涵盖了三个关键主题：ARM映象文件及其执行机理、GNU的映象构建过程以及BIOS的源码分析。ARM架构是现代嵌入式系统和移动设备中广泛使用的微处理器架构，理解其...
Miller_Rabin算法详解
xiange的自闭路
08-22 1万+
目录 基本引理： 1，费马定理： 2，二次探测定理： 作用： 证明： 代码实现： 目录 基本引理： 1，费马定理： 2，二次探测定理： 基本引理： 1，费马定理： 费马定理的证明链接 2，二次探测定理： 二次探测定理的证明链接 作用： 有效的检测大整数是否为素数。 证明： 由费马定理，可以排除大部分非素数的情况（满足费马定理是素数的必要条件），给出一个奇素数n...
二次探测定理（数论，理解Miller_Rabin算法所需要引理）
xiange的自闭路
08-22 4865
最近学习的过程中接触了随机算法，Miller_Rabin算法，关于大数判断是否为素数的随机算法，然后就学习的数论知识。 定理很简单，如果p为一个素数，则的解为，. 证明过程如下： 由p为一个素数可以推出。...
逆元详解（加扩展欧几里得和费马小定理的证明）
xiange的自闭路
08-08 3247
最近，wyb小朋友老是不好好搞他的数据结构，跑过来问我数学，没办法，所以我决定每天发一篇数论的博客，骗骗流量（以后wyb有不会的就看我博客，哈哈哈）先从基础的更起吧。 逆元： 我第一次接触逆元是在离散数学的代数系统中，对于一种运算满足（为该运算的单位）则称是的逆元。 逆元在算法中的运用： 逆元在算法中主要是为了整数的除法取模，显然除法是不能直接取模的。但是我们可以转化一下，因为乘法是可以直...
卡特兰数详解
xiange的自闭路
08-29 2194
目录 关于卡特兰数： 计算公式 一般性质 代码实现： 关于卡特兰数： 卡特兰数是一种经典的组合数，经常出现在各种计算中，其前几项为 :     1, 2, 5, 14, 42, 132, 429, 1430, 4862, 16796, 58786…… 计算公式 卡特兰数一般的计算公式：  另类递推公式：C(n)=C(n-1)*((4*n-2)/(n+1)); 一般性质 Cn的另...
欧拉降幂（欧拉-费马定理）
xiange的自闭路
07-30 1492
map&lt;ll,ll&gt; map_phi; ll get_phi(ll n) { ll res = n,a = n; for(ll i = 2;i * i &lt;= n;++i){ if(a % i == 0){ res = res / i * (i - 1); while(a % i == 0) a...
杜教筛详解：积性函数与欧拉函数在算法竞赛中的应用
算法竞赛专题解析--杜教筛1深入探讨了算法竞赛中常见的积性函数及其在欧拉函数中的应用。首先，章节1介绍了积性函数的概念，包括其定义、基本问题和常见类型，如欧拉函数，这是一种在数学中用于衡量正整数因子个数的...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值