当前位置:网站首页>线性筛求积性函数
线性筛求积性函数
2022-08-08 23:06:00 【AC__dream】
题目链接:登录—专业IT笔试面试备考平台_牛客网
样例输入:
5
1
2
3
4
5
样例输出:
1
2
2
3
2
分析:我们可以对q个询问进行暴力求解,每次询问的复杂度就是O(n^0.5),那么总的复杂度就是q*n^0.5,大概是10^8.5,估计过不了,况且如果对于别的题,询问组数再多一些,那么就铁定过不了了。
设f[n]表示n的所有正因数的个数,那么f[n]是一个积性函数,如果不知道积性函数或者不知道为什么f[n]是一个积性函数的话可以看一下:积性函数_AC__dream的博客-CSDN博客
我们可以通过线性筛直接O(n)预处理出来所有的f[n],直接对于每组询问O(1)查询即可。
过程由于在之前积性函数那篇博客中已经讲解的非常明白了,这里就不赘述了
下面是代码:
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<vector>
#include<queue>
using namespace std;
const int N=1e7+10;
bool vis[N];
int prime[N],tt,f[N];
int cnt[N];//cnt[i]记录i的最小质因子的次数
void init()
{
f[1]=1;
for(int i=2;i<N;i++)
{
if(!vis[i])
{
prime[++tt]=i;
f[i]=2;
cnt[i]=1;
}
for(int j=1;j<=tt&&i*prime[j]<N;j++)
{
vis[i*prime[j]]=true;
if(i%prime[j]==0)//prime[j]是i的最小质因子
{
cnt[i*prime[j]]=cnt[i]+1;
f[i*prime[j]]=f[i]/(1+cnt[i])*(2+cnt[i]);
break;
}
cnt[i*prime[j]]=1;
f[i*prime[j]]=f[i]*f[prime[j]];
}
}
}
int main()
{
int q;
init();
cin>>q;
while(q--)
{
int n;
scanf("%d",&n);
printf("%d\n",f[n]);
}
return 0;
}
边栏推荐
- 4399IT运维实习生面试经历
- 请问:支付宝上买基金安全吗
- stm32 利用 串口接收空闲中断 + dma 实现不定长度dma 接收
- 微信小程序错误 undefined Expecting ‘STRING‘,‘NUMBER‘,‘NULL‘,‘TRUE‘,‘FALSE‘,‘{‘,‘[‘, got ]解决方案
- JS中的预编译(AO、GO详解)
- MES对接Simba实现展讯平台 IMEI 写号与耦合测试
- 每日一R「01」跟着大佬学 Rust
- php7.4安装ssh2扩展和使用ssh链接sftp上传下载文件
- WeChat applet wx:for loop output example
- MES docks with Simba to realize IMEI number writing and coupling test of Spreadtrum platform
猜你喜欢
Virtual router redundancy protocol VRRP - double-machine hot backup
[Bug solution] ValueError: Object arrays cannot be loaded when allow_pickle=False
影响你各应用间网速的QoS你了解吗?
wps表格怎么筛选出需要的内容?wps表格筛选出需要的内容的方法
生活中无处不在的MPLS虚拟专用网
【YOLOv5】6.0环境搭建(不定时更新)
最详树莓派4B装机流程及ifconfig不到wlan0的解决办法
使用Mongoose populate实现多表关联存储与查询,内附完整代码
【Pytorch】学习笔记(一)
虚拟路由冗余协议VRRP——双机热备份基础
随机推荐
-Wl,--start-group ... -Wl,--end-group 用于解决几个库的循环依赖关系
PHP 正则给img的src添加域名
MPLS Virtual Private Network Everywhere in Life
Qt入门(四)——连续播放图片(两种定时器的运用)
DHCP's defense mechanism - DHCP Snooping (DHCP snooping)
C language library function summary2019.10.31
Kubernetes 企业如何落地
机器学习建模高级用法!构建企业级AI建模流水线
Firewall first contact
从stm32移植ucos2的代码到GD32
2022杭电多校六 1006-Maex (树形DP)
ArcPy设置全库唯一标识码
生活中无处不在的MPLS虚拟专用网
stm32使用spi1在slave 模式下 dma 读取数据
WeChat applet wx:for loop output example
微信小程序项目--订单
Xcode creates a Dylib plugin deb project
wps表格分两页断开怎么办?wps表格分两页断开的解决方法
Dynamic Host Configuration Protocol DHCP (DHCPv4)
数组去重的几种方法