筛法思想的题目

这道题目比较经典,或者说这种思想比较经典。

这种筛法的思想。

我们正着想对于每一个 n 、 n − 1 、 n − 2 、 . . . 、 2 、 1 n、 n-1、n-2、...、2、1 n、n−1、n−2、...、2、1都分解一遍质因数显然是来不及的时间复杂度达到 O ( n n ) O(n \sqrt{n}) O(nn )

我们考虑对于每一个1e6以内的质因数的个数

跑了一下程序是 78498 78498 78498个

素数定理告诉我们不超过x的素数近似有 x l n x \frac{x}{lnx} lnxx个

对于每一个质因子我们看一下那些数有一个这个质因子,那些数有两个,那些数有3个这个每次是平方的增长很快

时间复杂度: O ( n l o g x ∗ l o g n ) O(\frac{n}{logx} * logn) O(logxn∗logn)近似 O ( n ) O(n) O(n)

cpp 复制代码
#include <bits/stdc++.h> 
#define int long long
#define rep(i,a,b) for(int i = (a); i <= (b); ++i)
#define fep(i,a,b) for(int i = (a); i >= (b); --i)
#define pii pair<int, int>
#define pll pair<long long, long long>
#define ll long long
#define db double
#define endl '\n'
#define x first
#define y second
#define pb push_back

using namespace std;

const int N=1e6+10,mod=1e9+7;
int vis[N];
vector<int>p;
void solve()
{
	int n;cin>>n;
	auto get=[&](int x){
		rep(i,2,x){
			if(!vis[i])	p.pb(i);
			for(int j=0;p[j]*i<=x;++j){
				vis[i*p[j]]=1;
				if(i%p[j]==0)	break;
			}
		}
	};
	get(1e6);
	cout<<p.size()<<endl;
	int res=0;
	vector<int>ans(n+1);
	rep(i,0,p.size()-1){
		for(int j=p[i];j<=n;j*=p[i]){
			ans[i]+=n/j;
			res+=n/j;
		}
	}
	rep(i,0,p.size()-1){
		cout<<p[i]<<' '<<ans[i]<<endl;
	}
	return;
}

signed main(){
	ios::sync_with_stdio(false); cin.tie(0); cout.tie(0);
//   	freopen("1.in", "r", stdin);
  	int _;
//	cin>>_;
//	while(_--)
	solve();
	return 0;
}
相关推荐
whoarethenext3 分钟前
使用 C++ 实现 MFCC 特征提取与说话人识别系统
开发语言·c++·语音识别·mfcc
R-G-B4 分钟前
【MFC】Combobox下拉框中4个选项,运行后点击下拉框选项不能全部展示出来,只能显示2个选项,需要垂直滚动条滚动显示其余选项
c++·mfc
写个博客34 分钟前
暑假算法日记第二天
算法
ChaITSimpleLove1 小时前
.NET9 实现排序算法(MergeSortTest 和 QuickSortTest)性能测试
算法·排序算法·.net·benchmarkdotnet·datadog.trace
CVer儿1 小时前
svd分解求旋转平移矩阵
线性代数·算法·矩阵
Owen_Q1 小时前
Denso Create Programming Contest 2025(AtCoder Beginner Contest 413)
开发语言·算法·职场和发展
视觉人机器视觉2 小时前
Visual Studio2022和C++opencv的配置保姆级教程
c++·opencv·visual studio
liulilittle2 小时前
C++ i386/AMD64平台汇编指令对齐长度获取实现
c语言·开发语言·汇编·c++
Wilber的技术分享2 小时前
【机器学习实战笔记 14】集成学习:XGBoost算法(一) 原理简介与快速应用
人工智能·笔记·算法·随机森林·机器学习·集成学习·xgboost
巴里巴气2 小时前
selenium基础知识 和 模拟登录selenium版本
爬虫·python·selenium·爬虫模拟登录