一.进位计数制
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