LOJ6053简单的函数

https://loj.ac/p/6053
min_25筛模板题，但是要注意质数点有2比较特殊，需要在y==0的时候特判质数包含2的情况。

#include<bits/stdc++.h>
#define LL long long
using namespace std;
inline LL read()
{char x='\0';LL fh=1,sum=0;for(x=getchar();x<'0'||x>'9';x=getchar())if(x=='-')fh=-1;for(;x>='0'&&x<='9';x=getchar())sum=sum*10+x-'0';return fh*sum;
}
const LL N=1000009;
const LL mod=1e9+7;
LL n,sq;
LL pri[N],pcnt,sp[N];
bool vis[N];
inline void init(int mx)
{vis[0]=vis[1]=true;for(LL i=2;i<=mx;i++){if(!vis[i])pri[++pcnt]=i,sp[pcnt]=(sp[pcnt-1]+i)%mod;for(int j=1;j<=pcnt&&i*pri[j]<=mx;j++){vis[i*pri[j]]=true;if(i%pri[j]==0)break;}}
}
LL w[N],tot,ind1[N],ind2[N],g1[N],g0[N];
inline LL S(LL x,int y)
{if(x<=1||pri[y]>=x)return 0;LL k=(x<=sq)?ind1[x]:ind2[n/x];LL ans=g1[k]-g0[k]-(sp[y]-y);if(y==0)ans+=2;for(int i=y+1;i<=pcnt&&pri[i]*pri[i]<=x;i++){for(LL j=1,pe=pri[i];pe<=x;j++,pe*=pri[i]){ans=(ans+(pri[i]^j)*(S(x/pe,i)+(j!=1))%mod)%mod;}}return ans;
}
int main()
{n=read();sq=sqrt(n);init(sq);for(LL lp=1,rp=0;lp<=n;lp=rp+1){rp=n/(n/lp);w[++tot]=n/lp;if(w[tot]<=sq)ind1[w[tot]]=tot;else ind2[n/w[tot]]=tot;LL now=w[tot]%mod;g1[tot]=(now*(now+1)/2%mod-1+mod)%mod;g0[tot]=(now-1+mod)%mod;}
//  cout<<g1[1]<<endl;//for(int i=1;i<=pcnt;i++){for(int j=1;j<=tot&&pri[i]*pri[i]<=w[j];j++){LL now=w[j]/pri[i];LL k=(now<=sq)?ind1[now]:ind2[n/now];g1[j]=(g1[j]-pri[i]*(g1[k]-sp[i-1])%mod+mod)%mod;g0[j]=(g0[j]-(g0[k]-(i-1))+mod)%mod;}}
//  cout<<g1[1]<<' '<<g0[1]<<endl;//printf("%lld",(S(n,0)+1)%mod);return 0;
}

LOJ6053简单的函数（min_25筛）相关推荐

LOJ.6053.简单的函数(Min_25筛)
题目链接 Min_25筛见这里: https://www.cnblogs.com/cjyyb/p/9185093.html https://www.cnblogs.com/zhoushuyu/p/91 ...
loj #6053 简单的函数 min_25筛
\(\color{#0066ff}{ 题目描述 }\) 某一天,你发现了一个神奇的函数\(f(x)\),它满足很多神奇的性质: \(f(1)=1\) \(f(p^c)=p \oplus c\) (\( ...
Min_25筛（LibreOJ #6053: 简单的函数）
一个简单的问题:求 (i为质数) 设小于的质数有t个,其中表示第i个质数再设 (i为质数或者i的最小质因子不小于) 那么有其中 , 就是答案因为每次S(x, y)都是 ...
Min_25 筛小结
Min_25 筛这个东西,完全理解花了我很长的时间,所以写点东西来记录一些自己的理解. 它能做什么对于某个数论函数 \(f\),如果满足以下几个条件,那么它就可以用 Min_25 筛来快速求出这个函 ...
[复习]莫比乌斯反演,杜教筛,min_25筛
[复习]莫比乌斯反演,杜教筛,min_25筛莫比乌斯反演做题的时候的常用形式: \[\begin{aligned}g(n)&=\sum_{n|d}f(d)\\f(n)&=\sum_ ...
Min_25筛学习笔记
引入问题:求一个积性函数\(f(i)\)的前缀和 \[ \sum_{i=1}^nf(i) \] 其中\(f(i)\)满足对质数\(p\)有\(f(p)\)是关于\(p\)的低次多项式.\(f(p^c) ...
Min_25筛学习Tip+链接
前言机房里差不多都会Min_25筛了,我也赶紧补一波坑v_v 参考: txc的Min_25筛学习笔记 yx的Min_25筛学习笔记由于前面两位dalao的标题都是笔记,所以我这里就是小记了,因为这 ...
【UOJ#188】Sanrd（min_25筛）
[UOJ#188]Sanrd(min_25筛) 题面 UOJ 题解今天菊开讲的题目.(千古神犇陈菊开,扑通扑通跪下来) 题目要求的就是所有数的次大质因子的和. 这个部分和\(min\_25\)筛中枚 ...
与Min_25筛有关的一些模板
模板求∑i=1nf(i),f(pk)=pk×(pk−1)\sum \limits_{i = 1} ^{n} f(i), f(p ^ k) = p ^ k \times(p ^ k - 1)i=1∑n ...

LOJ6053简单的函数（min_25筛）

LOJ6053简单的函数

LOJ6053简单的函数（min_25筛）相关推荐

最新文章

热门文章