标签:函数 数论 LL Day2 long QBXT Afternoon 欧拉 define
数论函数相关
首先是说了欧拉筛(线性筛)
先贴钟长者的Code
1~n 所有质数找出来
not_prime[i] 代表 i 有没有被筛过
for (int a=2;a<=n;a++)
{
if (not_prime[a] == false) plist[++pcnt] = a;
for (int b=1;b<=pcnt;b++)
{
int x = a * plist[b];
if (x>n) break;
not_prime[x] = true;
if (a % plist[b] == 0) break;
}
}
其实我在很久之前是整理过怎么写欧拉筛的,不过那个时候memset直接全部赋值true会炸,所以现在在这里放上我新的写法
当然,在写之前我先说一些有关数论函数的事情
一些性质
首先对于数论函数\(f(n)\),这个n一定是正整数,也就是说\(n\in Z^*\)
欧拉函数和莫比乌斯函数就是很经典的积性函数
这里推荐一篇博文,然后下面的代码就是利用欧拉筛求欧拉函数phi和莫比乌斯函数的Code
#include <bits/stdc++.h>
#define Heriko return
#define Deltana 0
#define S signed
#define U unsigned
#define LL long long
#define R register
#define I inline
#define D double
#define LD long double
#define mst(a,b) memset(a,b,sizeof(a))
#define ON std::ios::sync_with_stdio(false)
using namespace std;
I void fr(LL &x)
{
LL f=1;char c=getchar();
x=0;
while(c<'0'||c>'9')
{
if(c=='-') f=-1;
c=getchar();
}
while(c>='0'&&c<='9')
{
x=(x<<3)+(x<<1)+c-'0';
c=getchar();
}
x*=f;
}
bool np[1000000];
LL prime[1000000],tot,phi[100000],mu[100000];//phi是欧拉函数,mu是莫比乌斯函数,这俩都是积性函数
I void EP(LL x)
{
mst(np,false);tot=0;
for(R LL i=2;i<=n;i++)
{
if(!np[i]) prime[++tot]=i,phi[i]=i-1,mu[i]=-1;
for(R LL j=1;j<=tot && i*prime[j]<=x;j++)
{
LL t=i*prime[j];
np[t]=true;
if(t==0)
{
phi[t]=phi[i]*prime[j];
mu[t]=0;
break;
}
else phi[t]=phi[i]*phi[prime[j]],mu[t]=mu[i]*mu[prime[j]];
}
}
}
S main()
{
Heriko Deltana;
}
再往下就是一些组合数的知识,今天先咕了
标签:函数,数论,LL,Day2,long,QBXT,Afternoon,欧拉,define 来源: https://www.cnblogs.com/HRiver2/p/HR2note22.html
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。