数论（min_25筛）

类欧几里得算法

扩展欧拉定理（对于a,p不互质的情况计算a%p）：https://blog.csdn.net/hzj1054689699/article/details/80693756

二次剩余

https://blog.csdn.net/kele52he/article/details/78897187
对于n整数a若存在整数x 满足 x^2%n=a那么称a是mod n 的二次剩余
判定：当且仅当 a^(p-1)/2=1
引理：对于奇素数p下的二次剩余n，有两个不同的x满足x^2%p=n
一个奇素数的完全剩余系下有(p-1)/2 个二次剩余。

Cipolla 算法

求解 x^2=n mod p，$O(n\log n)$。

步骤：
1.随机一个 a，使得 $(a^2-n{)}^{(p-1)/2}=-1\mathrm{m}\mathrm{o}\mathrm{d}p$，期望次数 2 次。
2.令 $\omega =\sqrt{a^2-n}$，扩域为虚根。
3.$(a+\omega {)}^{(p+1)/2}$ 即为所求。

//xy%p=xy-(xy)/p*p
inline ll qmul(ll a,ll b,ll mod)
{
    a%=mod,b%=mod;
    return (((a*b)-(ll)((ll)((long double)a/mod*b+1e-3)*mod))%mod+mod)%mod;
}

扩展BSGS

狄利克雷卷积的逆

$f\ast f^{-1}=e$。通过构造，可以 $O(n\log n)$ 构造出函数的逆。
如果积性函数 $f$ 和 $g$ 对于任意质数 $p$ 满足 $f(p)=g(p)=1$，我们就可以快速求出 $f$ 的前缀和。
设 $h=f\ast g^{-1}$，那么就有 $f(p)=g(1)h(p)+g(p)h(1)=h(p)+1=1$，推出 $h(p)=0$。
我们要求的是：
$$\sum _{i=1}^{n}f(i)=\sum _{i=1}^n\sum _{d|i}g(d)h(\frac{i}{d})=\sum _{d=1}^{n}h(d)\sum _{i=1}^{n/d}g(i)$$
由于 $h$ 只在所有质因子次数都大于 1 的情况下才不为零，这样的数可以写成 $a^2{b}^3$ 的形式，不超过 $O(\sqrt{n})$ 个。枚举 $a$，枚举 $b$，就得到了 $a^2{b}^3$ 的质因数分解，即可快速计算。

狄利克雷生成函数

$$F(s)=\sum _{i=0}^{\infty }f(i)/i^s$$
多项式卷积即为狄利克雷卷积。

定义求导算子 ${}^{\prime }$:当 $p$ 为质数时，$p^{\prime }=1$，对于整数 $a$，$b$，有 $(ab{)}^{\prime }=a^{\prime }b+a{b}^{\prime }$。

1

求 ${\sum }_{i=1}^n\gcd (i,i^{\prime })$。
??? ??? ???

min_25 筛

作用：计算积性函数前缀和。复杂度 $O(\frac{n^{3/4}}{\log n})$。
要求：$f(p)$ 为关于 $p$ 的低次多项式，$f(p^c)$ 能很快计算。

第一步

首先我们想要求所有质数的函数值。构造完全积性函数 $F(n)$，满足 $F(p)=f(p)$。记 $g(n,j)={\sum }_{i=1}^n[mind(i)>p_j]F(i)$ ，$mind(i)$ 表示 $i$ 的最小质因子。换句话说统计的是质数和最小值因子大于 $p_j$ 的数的贡献。我们的目的便是求 $g(n,|P|)$。

构造完全积性函数的目的一是使 $g(n,0)$ 好算。比如如果 $f(p)$ 是低次单项式，那么 $g(n,0)$ 就是自然数幂和；如果 $f(p)$ 是低次多项式，因为我们是求和，因此可以拆成若干个单项式再相加，转化为上一种情况。

初始，我们知道 $g(n,0)$ 的值，考虑用质数把合数的贡献筛掉。构造完全积性函数的目的之二就是方便我们能够筛掉合数的贡献。我们按照 $j$ 从小到大的顺序转移 $g$。考虑 $g(n,j-1)$ 转移到 $g(n,j)$，有哪些数变得不合法了。就是那些最小质因子刚好为 $p_j$ 的那些数。

注意到，如果 $p_j^2$ 如果大于 $x$，由于没有小于 $x$ 的合数含有 $p_j$，因此直接转移过来。
$$g(n,j)=\{\begin{array}{ll}g(n,j-1)-p_j^k(g(n/p_j,j-1)-g(p_j-1,j-1))& ，p_j^2\le n\\ g(n,j-1)& ，p_j^2>n\end{array}$$
第一维只有根号 $n$ 大小。发现求 $g$ 的过程与原函数是不是积性函数没关系（$F(x)$ 是我们构造的），因此可以求 $1$ 到 $n$ 质数个数。

第二步

考虑用这个东西求函数的前缀和。定义 $s(n,j)={\sum }_{i=1}^n[mind(i)\ge p_j]f(i)$。目标是求 $s(n,1)$，边界条件是 $s(n,|p|)$ 为所有质数的答案。因此我们递归计算 $s$。转移的时候枚举最小质因子，枚举最小质因子 $p_k$ 的次数 $e$，$s[n/p_k^e][k+1]$。注意单独考虑 $f(p_k^e)$。

$$s(n,j)=g(n,|p|)-\sum _{i=1}^{j-1}f(p_i)+\sum _{k\ge j}\sum _e(f(p_k^e)s(n/p_k^e,k+1)+f(p_k^{e+1}))$$

Loj6053 简单的函数
Notice： 1 拿出来最后单独统计，val 数组里的值从大到小，数组开到两倍根号。

#include<bits/stdc++.h>
#define ll long long
#define pb push_back
#define ld long double
using namespace std;
const int mod=1e9+7,M=200010,inv2=(mod+1)/2;
int ind1[M],ind2[M];
bool np[M];
ll sum[M],p[M],cnt,tot,h[M],g[M],val[M],n,B;
int read()
{
	int x=0;char c=getchar(),flag='+';
	while(!isdigit(c)) flag=c,c=getchar();
	while(isdigit(c)) x=x*10+c-'0',c=getchar();
	return flag=='-'?-x:x;
}
void init(int n)
{
	np[1]=1;
	for(int i=2;i<=n;i++)
	{
		if(!np[i]) p[++cnt]=i;
		for(int j=1;j<=cnt&&i*p[j]<=n;j++) 
		{
			np[i*p[j]]=1;
			if(i%p[j]==0) break;
		}
	}
	for(int i=1;i<=cnt;i++) sum[i]=(sum[i-1]+p[i])%mod;
}
ll S(ll x,int j)
{
	if(x<=1||x<p[j]) return 0;
	int t=(x<=B)?ind1[x]:ind2[n/x];
	ll ans=g[t]-sum[j-1]-h[t]+(j-1);
	if(j==1) ans+=2;
	for(int k=j;k<=cnt&&p[k]*p[k]<=x;k++)
	{
		ll t1=p[k],t2=p[k]*p[k];
		for(int e=1;t2<=x;e++,t1=t2,t2*=p[k])
			ans=(ans+(p[k]^e)*S(x/t1,k+1)+(p[k]^(e+1)))%mod;
	}
	return ans;
}
int main()
{
	cin>>n;
	init(B=sqrt(n)+1);
	for(ll l=1,r;l<=n;l=r+1)
	{
		r=n/(n/l);
		val[++tot]=n/l;
		if(n/l<=B) ind1[n/l]=tot;
		else ind2[r]=tot;
		ll t=(n/l)%mod;
		h[tot]=(t-1)%mod;
		g[tot]=(2+t)*(t-1)%mod*inv2%mod;
	}
	for(int j=1;j<=cnt;j++)
	{
		for(int i=1;i<=tot&&p[j]*p[j]<=val[i];i++)
		{
			int t=(val[i]/p[j]<=B)?ind1[val[i]/p[j]]:ind2[n/(val[i]/p[j])];
			g[i]=(g[i]-p[j]*(g[t]-sum[j-1]))%mod;
			h[i]=(h[i]-h[t]+j-1)%mod; 
		}
	}
	cout<<(S(n,1)%mod+1+mod)%mod;
	return 0;
}
/*by DT_Kang*/