数学公式

转自：http://blog.csdn.net/hanhai768/article/details/38018991

（1）   整数的唯一分解定理：

      任意正整数都有且只有一种方式写出其素因子的乘积表达式。

      A=(p1^k1)*(p2^k2)*(p3^k3)*....*(pn^kn)   其中pi均为素数

（2）   约数和公式：

对于已经分解的整数A=(p1^k1)*(p2^k2)*(p3^k3)*....*(pn^kn)

有A的所有因子之和为

    S = (1+p1+p1^2+p1^3+...p1^k1) * (1+p2+p2^2+p2^3+….p2^k2) * (1+p3+ p3^3+…+ p3^k3) * .... * (1+pn+pn^2+pn^3+...pn^kn)

（3）   同余模公式：

(a+b)%m=(a%m+b%m)%m

(a*b)%m=(a%m*b%m)%m

（4）用递归二分求等比数列1+pi+pi^2+pi^3+...+pi^n：

一、若n为奇数,一共有偶数项，则：
      1 + p + p^2 + p^3 +...+ p^n

      = (1+p^(n/2+1)) + p * (1+p^(n/2+1)) +...+ p^(n/2) * (1+p^(n/2+1))   

     这里将其分成两半，提后一半的公因式p^(n/2+1)
      = (1 + p + p^2 +...+ p^(n/2)) * (1 + p^(n/2+1))

 

二、若n为偶数,一共有奇数项,则:
      1 + p + p^2 + p^3 +...+ p^n

      = (1+p^(n/2+1)) + p * (1+p^(n/2+1)) +...+ p^(n/2-1) * (1+p^(n/2+1)) + p^(n/2)

     可知n-1为偶数，对前偶数项进行上述操作再加上最后一项p^n
      = (1 + p + p^2 +...+ p^(n/2-1)) * (1+p^(n/2)) + p^n;

（5）反复平方法计算幂次式p^n

   这是本题关键所在，求n次幂方法的好坏，决定了本题是否TLE。

   以p=2，n=8为例

   常规是通过连乘法求幂，即2^8=2*2*2*2*2*2*2*2

   这样做的要做8次乘法

 

   而反复平方法则不同，

   定义幂sq=1，再检查n是否大于0，

While，循环过程若发现n为奇数，则把此时的p值乘到sq

{

   n=8>0 ，把p自乘一次， p=p*p=4     ，n取半 n=4

   n=4>0 ，再把p自乘一次， p=p*p=16   ，n取半 n=2

n=2>0 ，再把p自乘一次， p=p*p=256  ，n取半 n=1，sq=sq*p

n=1>0 ，再把p自乘一次， p=p*p=256^2  ，n取半 n=0，弹出循环

}

则sq=256就是所求，显然反复平方法只做了3次乘法

（6）费马小定理

费马小定理是数论中的一个重要定理，其内容为： 假如p是质数，且(a,p)=1，那么 a^(p-1) ≡1（mod p）。即：假如a是整数，p是质数，且a,p互质，那么a的(p-1)次方除以p的余数恒等于1。

（7）组合数快速幂求模

[cpp]  view plain  copy

1. ll c[M];  
2. ll pow_mod(ll a,ll b)  
3. {  
4.     ll s=1;  
5.     while(b)  
6.     {  
7.         if(b&1)s=s*a%mod;  
8.         a=a*a%mod;  
9.         b=b>>1;  
10.     }  
11.     return s;  
12. }  
13. c[0]=1;  
14.  for(ll k=1;k<=m;k++)  
15.  {  
16.         if(m-k<k)c[k]=c[m-k];  
17.          else c[k]=c[k-1]*(m-k+1)%mod*pow_mod(k,mod-2)%mod;  
18.  }  

（8）SG函数！！

直接上模板！！！

[cpp]  view plain  copy

1. int sg[N];  
2. bool hash[N];  
3. void sg_solve(int *s,int t,int N)   //N求解范围 S[]数组是可以每次取的值，t是s的长度。  
4. {  
5.     int i,j;  
6.     memset(sg,0,sizeof(sg));  
7.     for(i=1;i<=N;i++)  
8.     {  
9.         memset(hash,0,sizeof(hash));  
10.         for(j=0;j<t;j++)  
11.             if(i - s[j] >= 0)  
12.                 hash[sg[i-s[j]]] = 1;  
13.         for(j=0;j<=N;j++)  
14.             if(!hash[j])  
15.                 break;  
16.         sg[i] = j;  
17.     }  
18. }  

转载声明！！！！！内容基本由各处复制整理得！！