题目:https://loj.ac/problem/6053

参考博客:http://www.cnblogs.com/zhoushuyu/p/9187319.html

算 id 也可以不存下来,因为 \( \left \lfloor \frac{i}{n} \right \rfloor \) 的取值是连续的,当 \( i \leqslant \sqrt{n} \) 时取值就是 \( i \);

而 \( i > \sqrt{n} \) 时,因为 \( i \) 越大,\( \left \lfloor \frac{i}{n} \right \rfloor \) 越小直到变成 1,所以编号可以直接数到;

注意 int 类型的 \( \sqrt{n} \) 是下取整的,所以判断条件最好写成 \( x > \sqrt{n} \) 而非 \( x \geqslant \sqrt{n} \);

因为 \( n \) 是 1e10 的,所以各种地方要注意不能爆 long long;

质数的 \( f \) 应该也可以直接筛吧(就直接令所有 \( f[i] = i-1 \)),不过质数和 - 质数个数也很方便。

代码如下:

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;
int const xn=1e6+,mod=1e9+;
int m,f[xn],h[xn],pri[xn],cnt,sqr,ps[xn];
ll n,w[xn]; bool vis[xn];
ll pw(ll a,ll b){a%=mod; b%=(mod-); ll ret=; for(;b;b>>=,a=a*a%mod)if(b&)ret=ret*a%mod; return ret;}
int upt(ll x){while(x>=mod)x-=mod; while(x<)x+=mod; return x;}
void init()
{
for(int i=;i<=sqr;i++)
{
if(!vis[i])pri[++cnt]=i,ps[cnt]=upt(ps[cnt-]+i);
for(int j=;j<=cnt&&(ll)i*pri[j]<=sqr;j++)
{
vis[i*pri[j]]=;
if(i%pri[j]==)break;
}
}
}
int Id(ll x)
{
if(x>sqr)return n/x;//>!
else return m-x+;
}
int S(ll x,int y)
{
if(x<=||pri[y]>x)return ;
int k=Id(x);
int ret=upt((ll)f[k]-h[k]-(ps[y-]-y+));
if(y==)ret=upt(ret+);
for(int i=y;i<=cnt&&(ll)pri[i]*pri[i]<=x;i++)
{
ll p0=pri[i],p1=(ll)pri[i]*pri[i];//ll
for(int k=;p1<=x;k++,p0=p1,p1*=pri[i])
ret=upt(ret+(ll)(pri[i]^k)*S(x/p0,i+)%mod+(pri[i]^(k+)));
}
return ret;
}
int main()
{
scanf("%lld",&n); sqr=sqrt(n); int inv2=pw(,mod-); init();
for(ll i=,j;i<=n;i=j+)
{
w[++m]=n/i; j=n/(n/i);
f[m]=(ll)upt(w[m]+)*upt(w[m]-)%mod*inv2%mod;//upt!
f[m]=upt(f[m]);//筛质数和
h[m]=upt(w[m]-);//筛质数个数
}
for(int j=;j<=cnt;j++)//j=1
for(int i=;i<=m&&w[i]>=(ll)pri[j]*pri[j];i++)
{
int k=Id(w[i]/pri[j]);
f[i]=upt(f[i]-(ll)pri[j]*(f[k]-ps[j-])%mod);
h[i]=upt(h[i]-h[k]+j-);
}
printf("%d\n",upt(S(n,)+));
return ;
}

loj 6053 简单的函数 —— min_25筛的更多相关文章

  1. LOJ 6053 简单的函数——min_25筛

    题目:https://loj.ac/problem/6053 min_25筛:https://www.cnblogs.com/cjyyb/p/9185093.html 这里把计算 s( n , j ) ...

  2. LOJ.6053.简单的函数(Min_25筛)

    题目链接 Min_25筛见这里: https://www.cnblogs.com/cjyyb/p/9185093.html https://www.cnblogs.com/zhoushuyu/p/91 ...

  3. loj#6053. 简单的函数(Min_25筛)

    传送门 题解 \(Min\_25\)筛有毒啊--肝了一个下午才看懂是个什么东西-- \(zsy\)巨巨强无敌-- //minamoto #include<bits/stdc++.h> #d ...

  4. 简单的函数——Min_25筛

    %%yyb %%zsy 就是实现一下Min-25筛 筛积性函数的操作 首先要得到 $G(M,j)=\sum_{t=j}^{cnt} \sum_{e=1}^{p_t^{e+1}<=M} [\phi ...

  5. LOJ #6053. 简单的函数

    $Min$_$25$筛模版题 为什么泥萌常数都那么小啊$ QAQ$ 传送门:Here 题意: $ f(1)=1$$ f(p^c)=p⊕c(p 为质数,⊕ 表示异或)$$ f(ab)=f(a)f(b)( ...

  6. LOJ.6235.区间素数个数(Min_25筛)

    题目链接 \(Description\) 给定\(n\),求\(1\sim n\)中的素数个数. \(2\leq n\leq10^{11}\). \(Solution\) Min_25筛.只需要求出\ ...

  7. LOJ #6202. 叶氏筛法(min_25 筛)

    题意 求 \([L, R]\) 之间的素数之和 . \(L≤10^{10},2×10^{10} \le R \le 10^{11}\) 题解 一个有点裸的 min_25筛 ? 现在我只会筛素数的前缀和 ...

  8. min_25筛

    min_25筛 用来干啥? 考虑一个积性函数\(F(x)\),用来快速计算前缀和\[\sum_{i=1}^nF(i)\] 当然,这个积性函数要满足\(F(x),x\in Prime\)可以用多项式表示 ...

  9. 「算法笔记」Min_25 筛

    戳 这里(加了密码).虽然写的可能还算清楚,但还是不公开了吧 QwQ. 真的想看的 私信可能会考虑给密码 qwq.就放个板子: //LOJ 6053 简单的函数 f(p^c)=p xor c #inc ...

随机推荐

  1. php编写的抽奖程序中奖概率算法

    本文给大家分享的是php中奖概率算法,可用于刮刮卡,大转盘等抽奖算法.用法很简单,代码里有详细注释说明,一看就懂,有需要的小伙伴参考下吧. 我们先完成后台PHP的流程,PHP的主要工作是负责配置奖项及 ...

  2. HBase 协处理器编程详解,第二部分:客户端代码编写

    实现 Client 端代码 HBase 提供了客户端 Java 包 org.apache.hadoop.hbase.client.coprocessor.它提供以下三种方法来调用协处理器提供的服务: ...

  3. Navicat Premium 10/12——破解激活

    Navicat Premium 12官方Windows64位百度云 链接:https://pan.baidu.com/s/1hGmDljszQsUoi194CYdfmA 密码:1xff 官方下载链接 ...

  4. Pandas时间差(Timedelta)

    时间差(Timedelta)是时间上的差异,以不同的单位来表示.例如:日,小时,分钟,秒.它们可以是正值,也可以是负值.可以使用各种参数创建Timedelta对象,如下所示 - 字符串 通过传递字符串 ...

  5. npm和git代理

    npm 删除代理设置:npm config delete proxynpm config delete https-proxynpm 设置代理:npm config set proxy http:// ...

  6. 编辑器上的vim插件,常用到的快捷件

    [vim]常用到的插件命令 [Esc] 退出编辑(输入,插入)模式,回到命令模式中(常用)i, I进入插入模式(Insert mode):i 为『从目前光标所在处插入』, I 为『在目前所在行的第一个 ...

  7. javascript异步编程的几种方法

    目前工作中用的比较多的异步模式编程有如下几种方法 一 回调函数 这是异步编程最基本的方法,假设有两个函数f1和f2,后者等待前者的执行结果 f1(); f2(); 如果f1是一个很耗时的任务,可以考虑 ...

  8. Java Maven项目使用CXF插件生成WebService代理

    CXF生成代理类插件名称:cxf-codegen-plugin 实现功能: 指定代理类所在的包 生成soapheader 操作: eclipse中run as -> maven build -& ...

  9. Idea_00_资源贴

    一.精选资料 tengj/IntelliJ-IDEA-Tutorial IntelliJ IDEA 使用教程-极客学院 二.参考资料 eclipse&Myeclipse&Intelli ...

  10. 关于linux 编程

    emacs 编辑器 gcc/g++ 编译器 gdb 调试工具 valgrind 内存泄露检查 doxygen 文档组织工具