素数


这次想写一篇关于素数的这个博客,因为在上一次的中国大学生acm的比赛中,我和我的队友们就困在了哪里,然后其实我们的思路是对的,但是在素数求的过程中有很多的问题,然后最近就看了一些题目,然后呢我就最近再看一些人写得然后总结再把这些基本的代码总结一下。
搬运的大佬对于lunch的解析链接

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<bitset>
#include<cassert>
#include<cctype>
#include<cmath>
#include<cstdlib>
#include<ctime>
#include<deque>
#include<iomanip>
#include<list>
#include<map>
#include<queue>
#include<set>
#include<stack>
#include<vector>
#include<unordered_set>
#include<unordered_map>
using namespace std;
//extern "C"{void *__dso_handle=0;}
typedef long long ll;
typedef long double ld;
#define fi first
#define se second
#define pb push_back
#define mp make_pair
#define pii pair<int,int>
#define lowbit(x) x&-x

const double PI=acos(-1.0);
const double eps=1e-6;
const ll mod=1e9+7;
const int inf=0x3f3f3f3f;
const int maxn=4e4+10;
const int maxm=1e4+10;
#define ios ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
int n,m;

int prime[maxn+10];
void getprime()
{
	for(int i=2;i<=maxn;i++)
	{
		if(!prime[i]) prime[++prime[0]]=i;
		for(int j=1;j<=prime[0]&&i*prime[j]<=maxn;j++)
		{
			prime[prime[j]*i]=1;
			if(i%prime[j]==0) break;
		}
	}
}
int SG(int x)
{
	int res=0;
	if(x%2==0) res=1;
	int tmp=x,cnt=0;
	for(int i=1;i<=prime[0];i++)
	{
		if(tmp%prime[i]==0)
		{
			while(tmp%prime[i]==0)
			{
				tmp/=prime[i];
				if(prime[i]!=2) cnt++;
			}
		}
		if(tmp==1) break;
	}
	if(tmp!=1 && tmp!=2) cnt++;
	return cnt+res;
}
int main(void)
{
	getprime();
	int t;
	scanf("%d",&t);
	while(t--)
	{
		scanf("%d",&n);
		int res = 0;
		for(int i=0;i<n;++i){
			int x;
			cin>>x;
			res ^= SG(x);
		}
		if(res) printf("W\n");
		else printf("L\n");
	}
}

素数的判定

bool isprime1(int n){	
	if(n==1) return false;
    else {
    		for(int i=2;i*i<=n;i++)
    		if(n%i==0)
    		return false;
    	return true;
    	}
   }
bool isprime2(int n){
 bool yes=true;
 for(int i=2;i<=sqrt(n);i++){
  if(n%i==0){
   yes=false;
   break;
  }
 }
 return yes;
}

求素数

#define maxn 1000000
bool isprime[maxn+1];
void eratos(int n){
	int i,j;
	isprime[0]=isprime[1]=false;
	for(i=2;i<=n;++i)
	isprime[i]=true;
	for(i=2;i*i<=n;++i)
		if(isprime[i]){
		for(j=i*i;j<=n;j+=i)
		isprime[j]=false;
		}
}

求素数的个数

网络中大佬的讲解
在这里插入图片描述
在这里插入图片描述

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>

using namespace std;

int main()
{
    int n, cnt = 0;
    int prime[100001],vis[100001];
    while(~scanf("%d",&n)&&n)
    {
        memset(prime, 0, sizeof(prime));//初始化
        memset(vis, 0, sizeof(vis));
        for(int i = 2; i <= n; i++)
        {
            if(!vis[i])
                prime[cnt++] = i;
            for(int j = 0; j<cnt && i*prime[j]<=n; j++)
            {
                vis[i*prime[j]] = 1;
                if(i % prime[j] == 0) break;//不同于埃式筛的关键
            }
        }
        printf("%d\n", cnt);
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值