IT数码 购物 网址 头条 软件 日历 阅读 图书馆
TxT小说阅读器
↓语音阅读,小说下载,古典文学↓
图片批量下载器
↓批量下载图片,美女图库↓
图片自动播放器
↓图片自动播放器↓
一键清除垃圾
↓轻轻一点,清除系统垃圾↓
开发: C++知识库 Java知识库 JavaScript Python PHP知识库 人工智能 区块链 大数据 移动开发 嵌入式 开发工具 数据结构与算法 开发测试 游戏开发 网络协议 系统运维
教程: HTML教程 CSS教程 JavaScript教程 Go语言教程 JQuery教程 VUE教程 VUE3教程 Bootstrap教程 SQL数据库教程 C语言教程 C++教程 Java教程 Python教程 Python3教程 C#教程
数码: 电脑 笔记本 显卡 显示器 固态硬盘 硬盘 耳机 手机 iphone vivo oppo 小米 华为 单反 装机 图拉丁
 
   -> C++知识库 -> TZOJ 5160: Yifan and LCM bitset筛素数+初等数论 -> 正文阅读

[C++知识库]TZOJ 5160: Yifan and LCM bitset筛素数+初等数论

题目描述:

请添加图片描述

题意:

求lcm(1,2,3…n)。

题解:

请添加图片描述
以上结论其实挺好证明得到。
因为一个数不是一个质数的n次方那么将这个数分解质因数,他的素因子一定包含于前面的lcm。

代码:

因为n的数特别大。普通的欧拉筛素数没法满足,所以就需要用bitset来储存bool数组。

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N=1e8+5;
const int M=6e6;
bitset<N>factor;
int prime[M],ans;
unsigned int sum[M];
unsigned int ksm(int a,int n)
{
	unsigned int ans=1;
	while(n)
	{
		if(n&1)ans*=a;
		a*=a;
		n=n>>1;
	}
	return ans;
}
void init()
{
	for(int i=2;i<N;i++)
	{
		if(!factor[i])
		{
			prime[++ans]=i;
		}
		for(int j=1;j<=ans&&i<N/prime[j];j++)
		{
			factor[i*prime[j]]=1;
			if(i%prime[j]==0)
			{
				break;
			}
		}
	}
	sum[1]=2;
	for(int i=2;i<=ans;i++)
	{
		sum[i]=sum[i-1]*prime[i];
	}
}
void solve(int n)
{
	int pos=upper_bound(prime+1,prime+ans+1,n)-prime-1;
	unsigned int ans=sum[pos];
	//cout<<pos<<" "<<ans<<endl;
	for(int i=1;i<=ans&&prime[i]*prime[i]<=n;i++)
	{
		int tmp=n,res=0;
		while(tmp/prime[i])
		{
			res++;
			tmp/=prime[i];
		}
		res--;
		ans*=ksm(prime[i],res);
	}
	printf("%u\n",ans);
}
int main()
{
	init();
	int n;
	while(~scanf("%d",&n))
	{
		solve(n);
	}
}  
  C++知识库 最新文章
【C++】友元、嵌套类、异常、RTTI、类型转换
通讯录的思路与实现(C语言)
C++PrimerPlus 第七章 函数-C++的编程模块(
Problem C: 算法9-9~9-12:平衡二叉树的基本
MSVC C++ UTF-8编程
C++进阶 多态原理
简单string类c++实现
我的年度总结
【C语言】以深厚地基筑伟岸高楼-基础篇(六
c语言常见错误合集
上一篇文章           查看所有文章
加:2021-09-08 10:29:29  更:2021-09-08 10:32:11 
 
开发: C++知识库 Java知识库 JavaScript Python PHP知识库 人工智能 区块链 大数据 移动开发 嵌入式 开发工具 数据结构与算法 开发测试 游戏开发 网络协议 系统运维
教程: HTML教程 CSS教程 JavaScript教程 Go语言教程 JQuery教程 VUE教程 VUE3教程 Bootstrap教程 SQL数据库教程 C语言教程 C++教程 Java教程 Python教程 Python3教程 C#教程
数码: 电脑 笔记本 显卡 显示器 固态硬盘 硬盘 耳机 手机 iphone vivo oppo 小米 华为 单反 装机 图拉丁

360图书馆 购物 三丰科技 阅读网 日历 万年历 2024年11日历 -2024/11/23 20:49:17-

图片自动播放器
↓图片自动播放器↓
TxT小说阅读器
↓语音阅读,小说下载,古典文学↓
一键清除垃圾
↓轻轻一点,清除系统垃圾↓
图片批量下载器
↓批量下载图片,美女图库↓
  网站联系: qq:121756557 email:121756557@qq.com  IT数码