码迷,mamicode.com
首页 > 其他好文 > 详细

数论:入门到入土 Three

时间:2021-03-08 13:23:25      阅读:0      评论:0      收藏:0      [点我收藏+]

标签:std   war   namespace   break   amp   数论   cin   clu   mes   

超短文警告!

埃氏筛:

#include <iostream>
#include <cstring>
#include <stdio.h>
#define HRiver2 return
#define Warma 0
#define ll long long
#define maxn 114514
using namespace std;
bool isprime[maxn];
ll prime[maxn];
int As(int n)//埃氏筛 
{
	memset(isprime,true,sizeof(isprime));
	isprime[1]=false;
	for(int i=1;i*i<=n;i++)
	{
		if(isprime[i])
		{
			int j=i*i;
			while(j<=n)
			{
				isprime[j]=false;
				j+=i;
			}
		}
	}
}

int main()
{
	int n;
	cin>>n;
	As(n);
	for(int i=1;i<=n;i++)
	{
		if(isprime[i])
		cout<<i<<‘ ‘;
	}
	HRiver2 Warma;
}

欧拉筛:

#include <iostream>
#include <cstring>
#include <stdio.h>
#define HRiver2 return
#define Warma 0
#define ll long long
#define maxn 114514
using namespace std;
bool isprime[maxn];
ll prime[maxn];
ll sum;
int Es(int n)//欧拉筛 
{
	memset(isprime,true,sizeof(isprime));
	isprime[1]=false;
	for(int i=1;i<=n;i++)
	{
		if(isprime[i])prime[++sum]=i;
		for(int j=1;j<=sum&&i*prime[j]<=n;j++)
		{
			isprime[i*prime[j]]=false;
			if(i%prime[j]==0) break;
		}
	 } 
}
int main()
{
	int n;
	cin>>n;
	Es(n);
	for(int i=1;i<=n;i++)
	{
		if(isprime[i])
		cout<<i<<‘ ‘;
	}
	HRiver2 Warma;
}

数论:入门到入土 Three

标签:std   war   namespace   break   amp   数论   cin   clu   mes   

原文地址:https://www.cnblogs.com/HRiver2/p/HR2noteM3.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!