2021年hznu寒假集训第十二天 数论基础

2021年hznu寒假集训第十二天

常用运算符

基本常识

基本定理

常用算法

埃及筛

从整数中筛去合数，留下素数
考虑给定一个素数x 那么很明显所有x的倍数都不是素数，我们只需要枚举一个系数i,将所有的xi标记为
合数即可
复杂度参考前面的常识~O（nloglogn）

for(int i=2;i<maxn;++i)
{
	for(int j=i;j<maxn;j+=i)
	{
		notPrime[j]=1;
	}
}

一点点小优化，可以发现对于所有 j<i^2
一定在之前已经被筛掉了，所以可以从i^2开始筛

for(int i=2;i<maxn;++i)
{
	for(int j=i*i;j<maxn;j+=i)
	{
		notPrime[j]=1;
	}
}

线性筛

for(int i=2;i<maxn;++i)
{
	if(!notPrime[i])
	{
		prime[tot++] = i;
	}
	for(int j=0;j<tot&&i*prime[j]<maxn;++j)
	{
		notPrime[i*prime[j]]=1;
		if(i%prime[j]==0)
			break;
	}
}

用线性筛求积性函数

for(int i=2;i<maxn;++i)
{
	if(!notPrime[i])
	{
		prime[tot++]=i;
		f[i]=func();
		g[i]=f[i];
	}
	for(int j=0;j<tot&&i*prime[j]<maxn;++j)
	{
		notPrime[i*prime[j]]=1;
		if(i%prime[j]==0)
		{
			int ret = 1,x = i;
			while(x%prime[j]==0)
			{
				ret*=prime[j];
				x/=prime[j];
			}
			f[i*prime[j]]=f[i]/g[i]*f[ret];
			g[i*prime[j]]=f[ret];
			break;
		}
		f[i*prime[j]]=f[i]*f[prime[j]];
		g[i*prime[j]]=f[prime[j]];
	}
}

一般来说这些函数都会有一些方便计算的性质，中间求ret的部分可以被极大的简化

快速幂

int quickpow(int a,int n,int mod)
{
	int ans = 1;
	while(n)
	{
		if(n&1)
			ans=ans*a%mod;
		a=a*a%mod;
		n>>=1;
	}
	return ans;
}

矩阵快速幂

辗转相除法/欧几里得算法

int gcd(int a,int b)
{
	return b ? gcd(b, a % b) : a;
}

扩展欧几里得

int exgcd(int a,int b,int& x,int&y)
{
	if(b==0)
	{
		x=1,y=0;
		return a;
	}
	int d = exgcd(b,a%b,x,y);
	int tmp=x;
	x=y;
	y=tmp-y*(a/b);
	return d;
}

更相减损术

第一步：任意给定两个正整数；判断它们是否都是偶数。若是，则用2约简；若不是则执行第二步。
第二步：以较大的数减较小的数，接着把所得的差与较小的数比较，并以大数减小数。继续这个操作，直到所得的减数和差相等为止。
则第一步中约掉的若干个2的积与第二步中等数的乘积就是所求的最大公约数。
其中所说的“等数”，就是公约数。求“等数”的办法是“更相减损”法。

例1、用更相减损术求98与63的最大公约数。

解：由于63不是偶数，把98和63以大数减小数，并辗转相减：
98-63=35
63-35=28
35-28=7
28-7=21
21-7=14
14-7=7
所以，98和63的最大公约数等于7。
例2、用更相减损术求260和104的最大公约数。

解：由于260和104均为偶数，首先用2约简得到130和52，再用2约简得到65和26。
此时65是奇数而26不是奇数，故把65和26辗转相减：
65-26=39
39-26=13
26-13=13
所以，260与104的最大公约数等于13乘以第一步中约掉的两个2，即1322=52。

中国剩余定理/CRT

扩展中国剩余定理