幸运素数

跳皮筋 我第一,马兰开花二十一;

二五六 二五七,二八二九三十一;

三五六 三五七,三八三九四十一;

童谣《马兰花》作为建国以来最大的数学难题,里面数字的含义困扰了全国人民几十年,并且一直未解

难道这些数字真的只是为了押韵生搬硬凑的数字吗?

并!不!是!

经过wzy的多年研究,他发现了里面潜藏的惊天大秘密(并没有):

28=256,而257是离256最近的素数

28是第二个完美数,29是离28最近的素数

31是第四个幸运素数

仅仅这一句话中,就出现了二进制、素数、完美数、幸运数!

这哪里是童谣?这分明就是中国版的达芬奇密码!

今天wzy就来考考你《马兰花》里出现的数学知识:

幸运数是经由类似埃拉托斯特尼筛法的演算法后留下的整数集合,是在1955年波兰数学家乌拉姆提出。

由一组由1开始的数列为例:

1, 2, 3, 4, 5, 6, 7, 8, 9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25, …
先将所有偶数删去,只留下奇数:

1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23, 25,…
然后把数列中的第 2 个数字(设该数字为 x )的倍数对应的数删除,即把所有第 nx,x∈Z+ 个数删除,例如上述例子中,第 2 数字是 3 ,所以删去所有第 3n 个数:

1, 3, 7, 9, 13, 15, 19, 21, 25,…
新数列的第 3 项(每次都加上 1 )为 7 ,因此将新数列的第 7n 个数删除:

1, 3, 7, 9, 13, 15, 21, 25,…
若一直重复上述的步骤,最后剩下的数就是幸运数

(以上内容来自维基百科幸运数)

我们将既是幸运数又是素数的数叫做幸运素数

现在给你一个数N,请判断N是否是一个幸运素数

输入格式

第一行一个数T,代表有T个数(1≤T≤2×105)

第1∼T行,每行一个正整数N(1≤N≤2×105)

输出格式

对于每个输入的数N,如果N是幸运素数,输出Yes,否则输出No

样例

input
3
1
2
3
output
No
No
Yes

解题思路

思路很简单,就是判断一个数既是素数也是幸运数

#pragma GCC optimize(2)
#include<bits/stdc++.h>
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=2e5+100;
typedef long long ll;
int as[maxn],r[maxn+1],vis[maxn+2];//这一点最好加上个数字,不然有时候会溢出 
vector<int> ve;
void ss()
{
	vis[0]=vis[1]=1;
	for(int i=2;i*i<maxn;i++)
	{
		if(!vis[i])
		{
			for(int j=i*i;j<maxn;j+=i)
			vis[j]=1;//素数打表 
		}
	}
}
void luck()
{
	ve.push_back(0);
	for(int i=1;i<=maxn;i+=2)
	{
		ve.push_back(i);
	}
	int len=ve.size();
	int cnt=2,res=ve[2];	
	while(res<=len)
	{
		for(int i=1;i*res<=len;i++)
		{
			ve[i*res]=0;
		}
		for(auto it=ve.begin()+1;it!=ve.end();)
		{
			if(*it==0)
				ve.erase(it);
			else
				it++;
		}    //删除是零的数,也就是 
		len=ve.size();
		cnt++;
		res=ve[cnt];
	}
	for(int i=0;i<ve.size();i++)
		r[ve[i]]=1;
}

int main()
{
	ios::sync_with_stdio(false);
	int t;
	luck();
	ss();
	cin>>t;
	for(int i=1;i<=t;i++)
	{
		int x;
		cin>>x;
		if(!vis[x]&&r[x]) 
		cout<<"Yes"<<endl;
		else cout<<"No"<<"\n";
	}
	return 0;
}
  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值