下述理论主要参考书目： 电子版pdf:算法竞赛进阶指南(p133-150)

首先补充两个前置知识：约数和最大约数有关性质，可选择性看~

有关素数筛的算法 埃式$nlog(n)$&欧拉素数筛$o(n)$&大区间素数筛

约数

最大公约数

1.gcd(a,b)=gcd(a,a-b)=gcd(b,a-b)=d;

2.欧几里得算法： b!=0; gcd(a,b)=gcd(b,a%b)=d;

(ps:其实这里有个一般性的结论：gcd(m,$n$)=gcd(m,$k\ast m+r_1$),n=$k\ast m+r_1$,k为整数,0<=$r_1$<m)

1.证明：d|a 且 d|b 那么d| (a-b) 得证

2.证明：a<b 显然易得
a>=b: 令 a=$q\ast b$+r;
由于d|a d|$(q\ast b)$
那么易得 d|(a-$q\ast b$) 即得d|r
又r=a mod b ,所以gcd(a,b)=gcd(b,a%mod);

欧拉函数

φ(N)=N*(1 - $\frac{1}{p_1}$)(1 - $\frac{1}{p_2}$)(1 - $\frac{1}{p_3}$)…(1 - $\frac{1}{p_m}$)

证明：

具体用到了容斥原理：
φ(N)=N - ($\frac{N}{p_1}$ + $\frac{N}{p_2}$ + $\frac{N}{p_3}$ + $\frac{N}{p_4}$- $\frac{N}{p_1{p}_2}$ - $\frac{N}{p_2{p}_3}$-$\frac{N}{p_1{p}_3}$ - $\frac{N}{p_1{p}_4}$ - $\frac{N}{p_2{p}_4}$ - $\frac{N}{p_3{p}_4}$ + $\frac{N}{p_1{p}_2{p}_3}$ + $\frac{N}{p_1{p}_2{p}_4}$…-…+…-…$\frac{N}{p_1{p}_2{p}_3......p_m}$)
(ps:分式这小东西麻烦啦啊aaa…)
上面的式子因式分解就得到了:
即是： φ(N)=N*(1 - $\frac{1}{p_1}$)(1 - $\frac{1}{p_2}$)(1 - $\frac{1}{p_3}$)…(1 - $\frac{1}{p_m}$)

欧拉函数的一些性质

2.性质2我们知道了欧拉函数是积性函数(PS:积性函数：对于任意互质的整数a和b有性质f(ab)=f(a)f(b)的数论函数)

3.p 为质数，若p|n 且 $p^2$|n，则φ(n)=$\phi (\frac{n}{p})\ast p$;

4.p 为质数，若p|n 且 $p^2$不n整除，则φ(n)=$\phi (\frac{n}{p})\ast (p-1)$;

5.${\sum }_{d|n}\phi (d)=n$

证明：

1.由于gcd(n,x)=gcd(n,n-x)=1,所以与n互质的数是"成对"的出现。进而n+(n-x)=n. 就得出了
2.性质2就很显然了，根据欧拉函数的公式 φ(N). a,b互质，所以分解质因数没有一样的，根据欧拉函数就可以将 φ(a) φ(b)合并得出了φ(N)。
所以欧拉函数是积性函数
3.用到了 φ(N)=N $\ast$ (1 - $\frac{1}{p_1}$)(1 - $\frac{1}{p_2}$)(1 -$\frac{1}{p_3}$)…(1 - $\frac{1}{p_m}$)
$p^2$|n主要告诉咱们N分解后p的质因数个数不少于2. 所以$\frac{N}{p}$分解后的质因数还有p.
φ($N$)和φ($\frac{N}{p}$)的欧拉函数表达式中(1 - $\frac{1}{p_i}$)全都一样，只是一个是N和$\frac{N}{p}$。
所以就得出了φ(N)=$\phi (\frac{N}{p})\ast p$;
4.是同理3的证明方法。
5.性质5：$f(n)={\sum }_{d|n}\phi (d)=n$ 证明比较难了，首先得证明$f(n)$是积性函数. $f(n)$积性函数证明：其实有点显然了, 若n,m互质，出发点仍然是质因数分解考虑. 举个栗子:n=6, m=35, n,m互质.
n中的d：1,2,3. m中的d：1,5,7. 然后我们知道欧拉函数 φ(N)是个积性函数(见性质2)。
$(1,2,3)\ast (1,5,7)$的分解得出集合d全是能整除$n\ast m$，不多不少。 ps:$(1,2,3)\ast (1,5,7)$ 得到的d集合(1,5,7,2,10,14,3,15,21).
所以就知道了嘛: n,m互质 那么:$f(n\ast m)=f(n)f(m)$ 得证：$f(n)$是积性函数

ps:φ(1)=1

同余

定义：若整数a和整数b除以正整数m的余数相等，则称a,b模m同余。记为a $\equiv$ b(mod m)。
eg：10 $\equiv$ 24(mod 7) 同余数3.

同余的基本性质

在写性质之前，先得有种这样的思路。
a=$p\ast m$+$r_1$
b=$q\ast m$+$r_2$
如果a$\equiv$b(mod m)，那么$r_1$=$r_2$。

性质一: 如果a$\equiv$b(mod m),c$\equiv$d(mod m).那么($a\pm c$)$\equiv$($b\pm d$) (mod m)。

(ps:从$r_1$ 和$r_2$这方面思考)。

性质二:如果a$\equiv$b(mod m),c$\equiv$d(mod m).那么($a\ast c$)$\equiv$($b\ast d$) (mod m)。

证明: 首先从已知得到a和b的余数是$r_1$，c和d的余数是$r_2$.
a=$p\ast m$+$r_1$, c=$q\ast m$+$r_2$
然后呢$a\ast c$=$p\ast q\ast m$ +$(r_1+r_2)\ast m$ +$r_1\ast r_2$.
对于$b\ast d$同理的，那些($p\ast q\ast m$ +$(r_1+r_2)\ast m$)能被整除的不用管它,然后后面的尾巴都是$r_1\ast r_2$.所以得以证明。

性质三: 如果ac$\equiv$bc(mod m), 且 gcd(c,m)=1(c和m互质) , 则有a$\equiv$b(mod m)。

这条性质很重要，反复品味。
证明： 对于ac$\equiv$bc(mod m)有ac=$p\ast m$+$r_1$,bc=$q\ast m$+$r_1$。
所以ac-$p\ast m$=bc-$q\ast m$.
得到：$(a-b)\ast c$=$(p-q)\ast m$。
再可得到$(a-b)\ast c$中能被m整除。
然后根据算术基本定理：n=$p_1^{c_1}$$p_2^{c_2}$$p_3^{c_3}$$p_4^{c_4}$…$p_m^{c_m}$
由于gcd(c,m)=1, m不在c中，那么m只能在$a-b$中。(即是$a-b$=$k\ast m$)
用符号表示：m|(a-b).
a=$p\ast m$+$r_1$
b=$q\ast m$+$r_2$
a-b=$(p-q)\ast m$+$(r_1-r_2)$
能被m整除，所以$r_1$=$r_2$
这就说明了a和b的余数都是$r_1$。 所以a$\equiv$b(mod m)。

性质四：若ac≡bd,c≡d(mod m)，且gcd(c,m)=1，得到a≡b(mod m)。

(ps:其实性质三是性质四的特殊情况，c=d )
证明: 首先c与d的余数相等，记为$r_1$,
又因为gcd(c,m)=1;
所以$r_1$=1
所以我们得出了gcd(d,m)=1.
a=$p\ast m$+$r_2$
c=$q\ast m$+$1$
b=$k\ast m$+$r_3$
ac=($p\ast q\ast m$) + ($p\ast m$) + ($q\ast m\ast r_1$) + $r_2$
bd=…+$r_3$
由于ac与bd余数相等,所以$r_2=r_3$。
所以a≡b(mod m)。

同余类与剩余系

gcd($a\ast b$,m)=1; gcd($a\ast b$,m)=gcd($(a\ast b)$%m,m)=1. (ps：见最大公约数性质二)
所以$(a\ast b)$%m 与 m互质。 所以得出了简化剩余系关于模m乘法封闭。
(ps:这个关系怎么理解呢？对于所有的a，gcd(a,m)=1，那么gcd(a%m,m)=1，即是所有a都能映射到m的简化剩余系中)

费马小定理 & 欧拉定理

1.费马小定理：若p是质数，则对于任意整数a，有$a^p\equiv a$(mod p)。
2.欧拉定理：若正整数a,n互质，则$a^{\phi (n)}\equiv 1$(mod n).
ps:我个人一般会写成$a^{\phi (p)}\equiv 1$(mod p).
证明：
设n的简化剩余系为{$\overline{a_1}$,$\overline{a_2}$,$\overline{a_3}$,$\cdots$$\overline{a_{\phi (n)}}$}。对于$\forall a_i,a_j$,若$a\ast a_i\equiv a\ast a_j$(mod n)。
根据同余性质三的结论(如果ac$\equiv$bc(mod m), 且 gcd(c,m)=1(c和m互质) , 则有a$\equiv$b(mod m)。) 可知道$a_i\equiv a_j$ (mod n)又因为$a_i,a_j\in (1,n)$,所以$a_i\ne a_j$。
这个证明有什么用呢？意味着若$a_i\ne a_j$,则$(a\ast a_i)$%n和$(a\ast a_j)$%n不等。
那么$a\ast a_i$等一系列是一一对应的。
因此，集合{$\overline{a_1}$,$\overline{a_2}$,$\overline{a_3}$,$\cdots$$\overline{a_{\phi (n)}}$}与集合{$\overline{a{a}_1}$,$\overline{a{a}_2}$,$\overline{a{a}_3}$,$\cdots$$\overline{a{a}_{\phi (n)}}$}都能表示n的简化剩余系。
所以我们得出了这样一个东东…
$a^{\phi (n)}{a}_1{a}_2{a}_3\cdots a_{\phi (n)}\equiv (a{a}_1)(a{a}_2)(a{a}_3)\cdots a_{\phi (n)}$
$(a{a}_1)(a{a}_2)(a{a}_3)\cdots a_{\phi (n)}\equiv a_1{a}_2{a}_3\cdots a_{\phi (n)}$(mod n)
$a{a}_1=p_1\ast m$+$r_i$,$a{a}_2=p_2\ast m$+$r_j$,$a{a}_3=p_3\ast m$+$r_k$,然后像这样所有的$a{a}_s$的余数$r_i,r_j,\cdots r_k$组成的集合就是n的简化剩余系(ps:用到了上面证明的两个集合一一对应的性质，我相信您能看到这里对同余的理解应该很清楚了)就得出了上面的等式。
结合两个等式$a^{\phi (n)}{a}_1{a}_2{a}_3\cdots a_{\phi (n)}\equiv a_1{a}_2{a}_3\cdots a_{\phi (n)}$(mod n)
$a_1{a}_2{a}_3\cdots a_{\phi (n)}$与n互质(根据同余性质三)
得出了$a^{\phi (n)}\equiv 1$(mod n)
在上述中我们的前提是a与n是互质的。故我们证明了欧拉定理。
如果n是质数，那么φ(n)=n-1。 在欧拉定理的基础上左右乘以a，得到$a^n\equiv a$(mod n)
即证明了费马小定理

趁热打铁来一道关于BSBG的算法题(依据费马下定理)

Discrete Logging POJ - 2417

#include<iostream>
#include<cstdio>
#include<cmath>
#include<map>
#define ll long long
#define rep(i,a,b) for(int i=a;i<=b;i++)
ll gcd(ll a,ll b){ return b? gcd(b,a%b):a;}
//const int N=2e5+10;
const ll mod=1e9+7;
using namespace std;
ll read(){
    ll s = 0, f = 1; char ch = getchar();
    while(!isdigit(ch)){
        if(ch == '-') f = -1;
        ch = getchar();
    }
    while(isdigit(ch)) s = (s << 3) + (s << 1) + (ch ^ 48), ch = getchar();
    return s * f;
}
ll a,b,p;
ll qpow(ll x,ll y){
	ll ans=1;
	while(y){
		if(y&1){
			ans=(ans*x)%p;
		}
			x=(x*x)%p;
			y>>=1;
	}
	return ans;
}
ll x;
map<ll,int> mp;
void solve(){
//	cout<<qpow(2,10)<<endl;
    while(scanf("%lld %lld %lld",&p,&a,&b)!=EOF){
         a%=p;
         if(a==0 and b==0) {
             printf("%d\n",1);
             return ;
         }
         if(a==0)
         {  
         printf("no solution\n");
             return ;
         }
    	 //注意一个小细节，就是在快速幂中是p哦，别写成了mod啦，一个小小的坑~
    	  ll m=ceil(sqrt(p));
    	  mp.clear();
    	  ll ans=b%p;
    	  //枚举区间i:1-m 
    	  for(int i=1;i<=m;i++){
				ans=ans*a%p;
				mp[ans]=i;
//				cout<<ans<<" "<<(i+1)<<endl;
		  }
		  ans=1;
		  int flag=0;
		  ll  t=qpow(a,m);
//		  ll top=p/m+1;
		  for(int i=1;i<=m;i++){
			  ans=ans*t%p;
		  	if(mp[ans]){
//		  		cout<<(i*m-mp[ans]+p)%p<<endl;
                printf("%lld\n",(i*m-mp[ans]+p)%p);
		  		flag=1;
		  		break;
			  }
		  }
		  if(!flag) printf("no solution\n");
//		   cout<<"no solution"<<endl;
	}
}
int main (){
//  freopen("in.txt","r",stdin);
//  freopen("out.txt","w",stdout);
  solve();
  return 0;
}

欧拉定理的推论

最小正整数$x_0$

欧拉的指数循环节

其他理论

1. i在区间[x,g(x)]内，$\lfloor k/i\rfloor$ 的值都相等。

说明：

求出i的一段连续区间使得：$\lfloor k/i\rfloor$ 是个定值。 eg： k=33，i=10; $\lfloor 33/10\rfloor$=3; 但是呢？ $\lfloor 33/8\rfloor$=3; 满足 $\lfloor 33/9\rfloor$=3; 满足 $\lfloor 33/11\rfloor$=3; 满足 $\lfloor 33/12\rfloor$=2; 不满足 所以他的一段连续的区间是：[8,11]。 i在[8,11]内， $\lfloor k/i\rfloor$是个定值。

证明：

由于g(x)>=x,所以i在这个区间内是个定值。

来看一到题：余数之和

#include<bits/stdc++.h>
#define ll long long
#define rep(i,a,b) for(int i=a;i<=b;i++)
ll gcd(ll a,ll b){ return b? gcd(b,a%b):a;}
const int N=2e5+10;
const ll mod=1e9+7;
ll read(){
    ll s = 0, f = 1; char ch = getchar();
    while(!isdigit(ch)){
        if(ch == '-') f = -1;
        ch = getchar();
    }
    while(isdigit(ch)) s = (s << 3) + (s << 1) + (ch ^ 48), ch = getchar();
    return s * f;
}
using namespace std;
ll n,k;
ll sum=0;
void solve(){
    n=read();k=read();
	sum=n*k;
	ll x=1;
	ll y;
	ll cnt;
	while(x<=n){
		if(k/x !=0)
		 y=k/(k/x);
		 else y=n;
		 cnt=y-x+1;
		if(y<=n)
		sum-= (k/x)*cnt*(x+y)/2;
		else {
			sum-=(k/x)*(n-x+1)*(x+n)/2;
		}
		x=y+1;
	}
	cout<< sum<<endl;
}
int main (){
//   freopen("in.txt","r",stdin);
//   freopen("out.txt","w",stdout);
  solve();
//   while(1);
  return 0;
}