数论笔记

ps:以算法进阶指南为学习资料，主要用于复习背

质数：

埃氏筛

代码

for (int i=2;i<=sqrt(n);i++)
  if (!vis[i])
    for (int j=2;j<=n/i;j++)
	  vis[i*j]=true;

线性筛

代码

for (int i=2;i<=n;i++)
{
	if (!v[i]) {p[++k]=i;v[i]=i;}
	for (int j=1;j<=k;j++)
	{
		 if (p[j]>v[i]||p[j]>n/i) break;
		 v[i*p[j]]=p[j];
	}
}

约数

唯一分解定理相关(含互质内容）

定理：对于大于1的正整数N可以分解为
$N=p_1^{c_1}{p}_2^{c_2}...p_n^{c_n}(p_i为质数)$
推论：
N的正约数个数为：${\prod }_{i=1}^n{c}_i+1$
N的所有正约数之和为：${\prod }_{i=1}^m({\sum }_{j=0}^{c_i}(p_i{)}^j)$

欧拉函数：

与正整数N互质的数的个数成为欧拉函数，记为$\phi$(N)
$\phi (N)$=$N\ast {\prod }_{质数p|N}1-\frac{1}{p}$

线形求欧拉函数：

代码

	for (int i=2;i<=n;i++)
	{
        if (!p[i]) f[i]=i-1;
        else if (i%(p[i]*p[i])==0) f[i]=f[i/p[i]]*p[i];
        else f[i]=f[i/p[i]]*(p[i]-1);
        ans+=(f[i])*2; 
    }

公约数

欧几里得算法：

用于求最大公约数
定理：
gcd(a,b)=gcd(b,a mod b)
代码：

int gcd(int a,int b)
{
   if (b==0) return a;
   return gcd(b,a%b);
}

更相减损术略

同余

扩展欧几里得（exgcd）

用于求ax+by=c，线性同余方程

裴蜀定理：

对于任意整数a,b,存在一对整数 x,y,满足ax+by=gcd(a,b);
证明：

欧几里得算法的最后，会得到gcd(a,0)，即最大公约数为a。 当x=1是等式显然成立(a1+0y=a);
假设对于a,b当a=b,b=a mod b时等式成立 有一组解x1,y1 即bx1+(a mod b)y1=gcd(b,a mod b)
因为 a mod b=a-$\lfloor a/b\rfloor$b 所以bx1+(a mod b)y1=bx1+(a-$\lfloor a/b\rfloor$b)y1 将等号后的式子展开得 bx1+(a-b$\lfloor a/b\rfloor$)y1=bx1+ay1-b$\lfloor a/b\rfloor$y1=ay1+b(x1-$\lfloor a/b\rfloor$y1) 所以可以对于a,b可以得到一组解x=y1,y=x1-$\lfloor a/b\rfloor$y1 a=b,b=a
mod b即为欧几里得算法的过程 有因为最后gcd(a,0)时有一组整数解，所以根据数学归纳法一定会有解。 据此就可以求出一组特解
对于更一般的方程ax+by=c有解仅当c是gcd(a,b)的约数时（我懒得写证明了） 设ax+by=gcd(a,b)的有解为x1,y1
则对于方程ax+by=c 有解x=x1/(gcd(a,b)/c),y=y1/(gcd(a,b)/c);

代码:

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

对于求x为最小正整数时的解（我不太会证）

int exgcd(int a,int b,int &x,int &y)
{
	if (b==0) {x=1;y=0;return a;}
    int  Q=exgcd(b,a%b,x,y);	
	int  x1=x;
	x=y;y=x1-(a/b)*y;
	return Q;
}
int main()
{
	cin>>a>>b>>c;
	int ans=exgcd(a,b,x,y);
	if (c%ans!=0) return 0;
	x=x*(c/ans);y=y*(c/ans);
	int d=b/ans;
	if (d<0) d=-d;
	x=x%d;
	if (x<0) x+=d;
	cout<<x;
}

费马小定理：

若p是质数，对于任意整数a ，$a^p\equiv a(modp)$

欧拉定理：

若正整数a,n互质 则$a^{\phi (n)}\equiv 1(modn)$

欧拉定理推论：

若正整数a,n互质，则对于任意正整数b,有$a^b\equiv a^{bmod\phi (n)}(modn)$

求解乘法逆元：

emm;
当你要求a/b%p的时候，
a/b%p貌似不等于a%p/b%p
所以你需要求得一个数x
使a/b$\equiv$a*x(mod p)
当 p为质数时，$x=b^{p-2}$就可以啦
不然啊，解方程呗- -

中国剩余定理：

设$m_1$$m_2$…$m_n$是两两互质的整数，m=${\prod }_{i=1}^n{m}_i$,$M_i$=m/$m_i$,$t_i$是线性同余方程$m_i$$t_i$ $\equiv$ 1(mod p)的一个解。对于任意的n个正整数$a_1$,$a_2$…$a_n$,方程组
$$\{\begin{array}{l}x\equiv a_1(mod{m}_1)\\ x\equiv a_2(mod{m}_2)\\ ...\\ x\equiv a_n(mod{m}_n)\end{array}$$
有整数解,解为${\sum }_{i=0}^n{a}_i{m}_i{t}_i$

求解高次同余方程(BSGS=北上广深=拔山盖世)

ps:不熟
求解$a^x$$\equiv$b(mod p)
蒟蒻只会a,p互质的
代码

long long m=(long long)ceil(sqrt(p));
long long ss=1;
h.clear();
for (int i=0;i<m;i++,ss=ss*a%p)
  h[b*ss%p]=i;
long long s=1;
for (int i=0;i<=m;i++)
{
	if (h.find(s)!=h.end()&&i*m-h[s]>=0) {cout<<i*m-h[s]<<endl;break;}
	s=s*ss%p;
}

矩阵乘法

这个看缘分了啦

#include <bits/stdc++.h>
using namespace std;
long long n,a[2][2],f[2],P=10000;
void aaa()
{
    long long c[2][2];
    memset(c,0,sizeof(c));
    for (int i=0;i<2;i++)
      for (int j=0;j<2;j++)
        for (int k=0;k<2;k++)
          c[i][j]=(c[i][j]+a[i][k]*a[k][j]%P)%P;
    memcpy(a,c,sizeof(c));
}
void aaaa()
{
   long long c[2];
   memset(c,0,sizeof(c));  
   for (int j=0;j<2;j++)
     for (int k=0;k<2;k++)
       c[j]=(c[j]+f[k]*a[k][j]%P)%P;
   memcpy(f,c,sizeof(c));
}
int main()
{
    cin>>n;
    for (;n!=-1;)
    {
        if (n==0) {cout<<0<<endl;cin>>n;continue;}
        f[0]=1;f[1]=0; a[0][0]=1; a[0][1]=1;a[1][0]=1;a[1][1]=0;
        n=n-1;
        for (;n>0;n=n>>1)
        {
            if (n&1) aaaa();
            aaa();

        }
        cout<<f[0]<<endl;
        cin>>n;

    }
    return 0;
}

斐波那契以表敬意

高斯消元

球形空间产生器：https://www.acwing.com/problem/content/209/
正常消元

#include <bits/stdc++.h>
using namespace std;
int n;
double a[22][22],b[22],c[22][22]; 
int main()
{
	freopen("test.in","r",stdin);
	freopen("test.out","w",stdout); 
    cin>>n;
    for (int i=1;i<=n+1;i++)
      for (int j=1;j<=n;j++)
        cin>>a[i][j];
    for (int i=1;i<=n;i++)
      for (int j=1;j<=n;j++) {
        c[i][j]=2*(a[i][j]-a[i+1][j]);
        b[i]+=a[i][j]*a[i][j]-a[i+1][j]*a[i+1][j];
      }
    for (int i=1;i<=n;i++)
      for (int j=i;j<=n;j++) {
        if (fabs(c[i][j])<1e-8) {
           for (int k=1;k<=n;k++) swap(c[i][k],c[j][k]);
           swap(b[i],b[j]);
           break;
        }
        for (int j=1;j<=n;j++) {
            if (j==i) continue;
            double m=c[j][i]/c[i][i];
            for (int k=i;k<=n;k++)
              c[j][k]-=c[i][k]*m;
            b[j]-=b[i]*m;
        }
      }
    for (int i=1;i<=n;i++)
      printf("%.3lf ",b[i]/c[i][i]);
    return 0;
}

开关问题：https://www.acwing.com/problem/content/210/
xor只好感性理解

#include <bits/stdc++.h>
using namespace std;
int T,n,a[30];
int J,x,y;
int ans;
int main() {
	freopen("switch.in","r",stdin);
	freopen("switch.out","w",stdout); 
    scanf("%d",&T);
    for (int i1=1;i1<=T;i1++)
    {
    	scanf("%d",&n);
    	for (int i=1;i<=n;i++)
    	  scanf("%d",&a[i]),a[i]|=1<<i;
    	for (int i=1;i<=n;i++)
    	  scanf("%d",&J),a[i]^=J;
    	for (;;) {
		  scanf("%d%d",&x,&y);
		  if (!x) break;
		  a[y]|=1<<x;
	    }
    	ans=1;
    	for (int i=1;i<=n;i++) {
          for (int j=i+1;j<=n;j++)
		    if (a[j]>a[i]) swap(a[j],a[i]);
		  if (a[i]==0) {ans=1<<(n-i+1);break;}
		  if (a[i]==1) {ans=0;break;}
		  for (int k=n;k;k--)
		    if (a[i]>>k&1) {
		      for (int j=1;j<=n;j++)
			    if (i!=j&&(a[j]>>k&1)) a[j]^=a[i];
			    //a[j]^=a[i],a[i]的这一位和a[j]的这一位都为一,，需要将a[j]
				//的这一位变为0,相同为0，异或运算,那取反可以吗？
				//取反前面就变成1了,并且这个消元似乎和a[i]没有什么关系
				//(总之不对吧)，又或者整个题就是一个异或？？ 
			  break;
	 	    }
	    }
	    if (ans==0) printf("Oh,it's impossible~!!\n");
	      else printf("%d\n",ans);
    }
}    	

线性基：

我只能感性理解
装备购买：https://www.acwing.com/problem/content/211/

#include <bits/stdc++.h>
using namespace std;
int n,m,c[510],ans=0,ans_=0;
long double a[510][510];
int main()
{
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++)
      for (int j=1;j<=m;j++)
        scanf("%Lf",&a[i][j]);
    for (int i=1;i<=n;i++)
      scanf("%d",&c[i]);
    for (int i=1;i<=m;i++)
    {
    	int Min=1e9;
    	int r=0;
    	for (int j=i;j<=n;j++)
		  if (fabs(a[j][i])>1e-8&&c[j]<Min) Min=c[j],r=j;
		if (r==0) continue;
		swap(a[i],a[r]),swap(c[r],c[i]);
		ans+=c[i];ans_++;
		for (int j=1;j<=n;j++)
		{
			if (j==i) continue;
			long double g=a[j][i]/a[i][i];
			for (int k=i;k<=m;k++)
			  a[j][k]-=a[i][k]*g;
	    }
    }
    printf("%d %d",ans_,ans);
    return 0;
}

xor的
XOR：https://www.acwing.com/problem/content/212/

#include <bits/stdc++.h>
using namespace std;
int T;
unsigned long a[10010];
int main()
{
    cin>>T;
	for (int i1=1;i1<=T;i1++) {
	  int n,Q;
	  unsigned long t;
	  bool fla=false;
	  cin>>n;
	  t=n;
	  for (int i=1;i<=n;i++)
		cin>>a[i];
	  for (int i=1;i<=n;i++) { 
	    for (int j=i+1;j<=n;j++)
		  if (a[j]>a[i]) swap(a[j],a[i]);
		if (a[i]==0) {fla=true;t=i-1;break;}
		for (int k=63;k>=0;k--)
		  if (a[i]>>k&1) {
		  	for (int j=1;j<=n;j++)
		  	  if (i!=j&&(a[j]>>k&1)) a[j]^=a[i];
		  	break;
	      }
	  }
	  cin>>Q;
	  cout<<"Case #"<<i1<<':'<<endl;
	  for (int i=1;i<=Q;i++) {
	  	unsigned long k;
	  	unsigned long ans=0;
	  	cin>>k;
	  	if (fla) k--;//能得到0
		if (k>=(1LL<<t)) {cout<<-1<<endl;continue;}
		for (int j=t-1;j>=0;j--)
		  if (k>>j&1) ans^=a[t-j];
		cout<<ans<<endl;
      }	 
    }
    return 0;
}

组合计数

多重集的排列数

$\frac{n!}{n_1!n_2!...n_k!}$

多重集的组合数

k种数，选r个(r<=$n_i$)$n_i$表示第i种数的个数
$\left(\genfrac{}{}{0ex}{}{k-1}{k+r-1}\right)$

Lucas定理

若p是质数，则对于任意整数1$\le$m$\le$n
$\left(\genfrac{}{}{0ex}{}{m}{n}\right)$$\equiv$$\left(\genfrac{}{}{0ex}{}{mmodp}{nmodp}\right)$*$\left(\genfrac{}{}{0ex}{}{m/p}{n/p}\right)$ (mod p)

Catalan 数列

$Ca{t}_n$=$\frac{\left(\genfrac{}{}{0ex}{}{n}{2n}\right)}{n+1}$

莫比乌斯反演

定义：

莫比乌斯函数：
$$\mu (n)=\{\begin{array}{ll}1,& n=1\\ (-1{)}^k,& n=p_1{p}_2{p}_3...p_k(p_i为互不相同的质数)\\ 0,& othervalues\end{array}$$

莫比乌斯反演定理：

形式1：

$$F(n)=\sum _{d|n}f(d)\Rightarrow f(n)=\sum _{d|n}\mu (d)F(\frac{n}{d})$$

形式2：

$$F(n)=\sum _{n|d}f(d)\Rightarrow f(n)=\sum _{n|d}\mu (\frac{d}{n})F(d)$$

变形1：

$$f(i)=\sum _{d=1}^{\lfloor \frac{n}{i}\rfloor }g(d\ast i)\Rightarrow g(i)=\sum _{d=1}^{\lfloor \frac{n}{i}\rfloor }\mu (d)f(d\ast i)$$

变形2：

$$f(i)=\sum _{i|d,d\le n}g(d)\Rightarrow g(i)=\sum _{i|d,d\le n}\mu (\frac{d}{i})f(d)$$

性质：

1.对于任意正整数n有
$$\sum _{d|n}\mu (d)=\{\begin{array}{ll}1,& n=1\\ 0,& n>1\end{array}$$
2.对于任意正整数n有
$$\sum _{d|n}\frac{\mu (d)}{d}=\frac{\phi (n)}{n}$$

概率和期望

这个东西似乎总是和DP搞一起，