码迷,mamicode.com
首页 > 其他好文 > 详细

数论总结1(基础数论)

时间:2018-01-09 13:34:57      阅读:281      评论:0      收藏:0      [点我收藏+]

标签:顺序   prim   数字   sqrt   欧拉函数   等于   ems   最小   markdown   

一.进位计数制

1.进制表示:

技术分享图片

表示b进制下的n+1位数。

2.进制转换:

进制向十进制转换:
乘以基数并展开:
技术分享图片

十进制向b进制转换:
整数部分除以基数并倒取余数。
小数部分乘以基数,并顺取整数部分
技术分享图片

3.例题:

a.天平1:

一个天平,有N个重量未知的砝码,砝码重量可由你自由确定。
砝码可任意放在天平的左右两边,但要求称出从1到M之间所有的重量,
现给出\(N\)的值,请问M最大值为多少。
答案:\([0,3^N]\)。具体见下面解析。

b.天平2:

一个天平,砝码分别为1g、3g、9g、27g、…、6561g…,
每个砝码只有一个,要称重的物品放在天平的左侧,而砝码允许放在天平的左右两侧。
已知一个物品的质量N (N≤108 ),问如何称重?
_____________________________________________________________________
解析:
就是将\(N\)转换成三进制后,
将三进制中的0、1、2三个状态转换成 0、1 、-1
具体的说,就是0和1不变,2变成-1后,其高一位加1。
例如:
\(N\)转化为三进制后为{1,2,0,1} , 那么对应的状态为:{1,-1,-1,0,1}。
即三进制中的2是没法直接称量的,因为每个重量的砝码只有一个,所以只能通过相减来处理。

c.天平3:

一个天平,砝码分别为1g、3g、9g、27g、…、6561g…,每个砝码只有一个,
要称重的物品放在天平的左侧,而砝码只允许放在天平的右侧。
将由这个系统可以称出来的重量按从小到大的顺序进行排列,
得到下列序列:1,3,4,9,10,...。问其中的第K个重量是多少?
_____________________________________________________________________
解析:仔细观察就可以发现,将K转化为二进制后,为1的即放置对应质量的三进制重量砝码即可。

二.素数和合数

1.素数判定:

枚举\(\sqrt{N}\)内的质因子即可。 时间复杂度\(O(\sqrt{N})\)

2.欧拉筛法(线性筛素数)

原理:?每个合数必有一个最小素因子,用这个因子筛掉合数。
实现代码(筛出\([0,N]\)的素数):

void Prime(int N){    
       memset(isPrime,true,sizeof(isPrime));
       isPrime[0 ] = false ;  isPrime[1] = false ;
       for (int i = 2 ; i <= N ; i++){
             if ( isPrime[i] )prime[++ total] = i ;            //把素数记录下来
             //遍历已知素数表中比i的最小素因数小的素数,并筛去合数
             for(int j = 1 ; j <= total && i * prime[j] <= N ; j++){ 
                 isPrime[ i * prime[j] ] = false ;
                 if  (!( i % prime[j]))  break;                   //找到i的最小素因数
             }
      }
}

三.唯一分解定理

对任一整数\(a>1\),有\(a= p_1^{a_1}p_2^{a_2}…p_n^{a_n}\)
几个比较重要的结论(证明都很显然)。
技术分享图片

四.欧拉函数

1.定义:

欧拉函数: \([1,n]\)中和\(n\)互素的元素个数 \(φ(n)\)
\[φ(n) = \Sigma_{i=1}^{n} [gcd(i,n) ==1]\]

2.性质:

  • <1> 如果\(n\)为某一素数\(p\),则\(φ(p)=p-1\)

  • <2> 如果\(n\)为某一素数\(p\)的幂次\(p^a\)\(φ(p^a)=(p-1)×p^{a-1}\)
    证明:\(φ(p^{a}) = φ(p^{a-1})*p = (p^a-1)*(1-\frac{1}{p})*p = (p-1)×p^{a-1}\)

  • <3> 欧拉函数是积性函数,即当\(gcd(m,n)=1\)\(f(mn)=f(m)*f(n)\)

  • <4> 欧拉函数计算公式:
    \(p_1^{a_1}×p_2^{a_2}×…×p_k^{a_k}\)为正整数n的素数乘积式,则

\[φ(n) =n×(1-\frac{1}{p_1}) ×(1-\frac{1}{p_2})×…×(1-\frac{1}{p_k})\]

3.欧拉函数线性筛

void eular(int n){
    for(int i=2;i<=n;i++){
        if (!IsPrime[i])prime[++cnt]=i,phi[i]=i-1;
        for(int j=1;j<=cnt;j++){
            if(prime[j]*i>n)break;
            Isprime[prime[j]*i]=1;
            if(i%prime[j]==0)
               {phi[i*prime[j]]=phi[i]*prime[j];  break;}       //p是x的约数
             else
               phi[i*prime[j]]=phi[i]*(prime[j]-1);         //p不是x的约数
        }
    }return;
} 

理解:
设P是素数,
若p是x的约数,则\(φ(x*p)=φ(x)*p\)
//没有加项。
若p不是x的约数,则\(φ(x*p)=φ(x)*φ(p)=φ(x)*(p-1).\)
//加了一项,即\(φ(x*p) = φ(x)*p*(1-\frac{1}{p}) = φ(x)*(p-1);\)

例题:

例一:Poj 2773Happy 2006

给出一个数字\(M\) (\(1 \leq M \leq 10^6\)), 数字\(K\) (\(1 \leq K \leq 10^8\)).求与M互质的第K个数字
_____________________________________________________________________
解法1:
注意到若(a,m) = 1,那么(a+m*k,m)=1
所以实际上与m互质的数分布是一个分块区间。
以[1,m]为一个周期循环,此法常数要卡一卡。
______________________________________________________________________
解法2:
先二分这个值n,
然后检查小于等于n且与m互质的数的个数。
先处理出m的因数,然后用\(φ(x)\)容斥即可。

例二:[SDOI2012]Longge的问题

求:\(∑gcd(i, N)\) , \((1\leq i \leq N)\) , 其中$0<N\leq 2^{32}
$ 。
解析:见博客:http://www.cnblogs.com/GuessYCB/p/8192826.html

数论总结1(基础数论)

标签:顺序   prim   数字   sqrt   欧拉函数   等于   ems   最小   markdown   

原文地址:https://www.cnblogs.com/GuessYCB/p/8250894.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!