筛法和欧拉函数

最新推荐文章于 2022-02-20 19:13:24 发布

klc1

最新推荐文章于 2022-02-20 19:13:24 发布

阅读量348

点赞数 1

分类专栏： ----------数学---------- 数论算法

本文链接：https://blog.csdn.net/k_l_c_/article/details/52012711

版权

----------数学---------- 同时被 3 个专栏收录

25 篇文章 0 订阅

订阅专栏

数论

12 篇文章 0 订阅

订阅专栏

算法

6 篇文章 0 订阅

订阅专栏

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<queue>
using namespace std;
#define M 100010
bool visit[M];
int p[M];
void iii()
{
    memset(visit,true,sizeof(visit));
    visit[1]=0;
    for(int i=2; i<M; ++i)//如果M*M>int的话类型改为long long
    {
        if(visit[i]==true)
        {
            p[np++]=i;
            for(int j=i*i; j<M; j+=i)
            {
                visit[j]=false;
            }
        }
    }
}
//v[i]==1 是素数

//线性筛法

欧拉函数打表
const int MAX=3000001;
long long e[MAX];
void geta()
{
    e[1]=0;
    for(int i=2;i<MAX;++i)
        e[i]=i;
    for(int i=2;i<MAX;++i)
        if(e[i]==i)
            for(int j=i;j<MAX;j+=i)
                e[j]=e[j]/i*(i-1);

    for(int i=2;i<MAX;++i)//去掉这句话,e[i]表示i的欧拉值
        e[i]+=e[i-1];
}
//e[i]表示从1到i每个数欧拉值的总和

int euler(int n)//返回一个数的欧拉值
{
    int res=n,a=n;
    for(int i=2; i*i<=a; i++)
    {
        if(a%i==0)
        {
            res=res/i*(i-1);
            while(a%i==0) a/=i;
        }
    }
    if(a>1) res=res/a*(a-1);
    return res;
}

#include<iostream>//区间筛法
#include<cstdio>
#include<cmath>
#include<cstring>
#include<queue>
using namespace std;
#define M 65538
#define LL long long
const int maxn=1000001;//区间长度最大值
bool visit[M];
bool is[maxn];
int ans[maxn];
LL p[M];
int cnt=0;
void iii()//找2~sqrt(INF)的素数,65537是素数且大于sqrt(2^31-1);
{
    memset(visit,true,sizeof(visit));
    visit[1]=0;
    for(long long i=2; i<M; ++i)
    {
        if(visit[i]==true)
        {
            p[cnt++]=i;
            for(long long  j=i*i; j<M; j+=i)
            {
                visit[j]=false;
            }
        }
    }
}
//区间筛法,找到大于等于n,小于等于m的素数
for(int i=0; i<maxn; i++)
    is[i]=1;
int d=sqrt(m+0.5);
for (LL i = 0; i < cnt && p[i]<= d; i++)
{

    int s = n/ p[i] + (n% p[i] > 0);
    if (s == 1) s = 2;
    for (LL j = s; j*p[i] <= m; j++)
        if (j*p[i] >= n)
            is[j*p[i]-n] = 0;
}