位运算详解

首先我们需要清楚位运算的一些基本知识。

与 a&b
或 a|b
非 ~a
异或 a^b

int 32位

1:0000000...1 
2:0000000...10
3:0000000...11
补码(表示计算机中的负数): ~a+1是a的补码即(-a)

memset(a,0x3f,sizeof a)//把这个a所占的全部字节初始化为0x3f(十六进制)

0x3f3f3f3f*2 < 0xffffffff //乘2的话不会溢出

1 -> 10 -> 100

1<<n => 2^n
n>>x => n/2^x

这是位运算的一些基本知识。通过位运算来进行加减乘除，能够很大的降低时间复杂度。

有了这些知识点之后我们可以用快速幂来做一道题，如下题所示。

这题需要求出a的b次方，有一个问题就是数据范围是a,b,p<=10e9，如果用for循环去遍历相乘的话时间复杂度最坏是10e9，而计算机一般一秒能运算5.0*10e8，所以o(n)的时间复杂度肯定是不行的，这里我们需要用到快速幂。

举个例子:

3 7 6

我们知道要求3^7的值需要for循环7次，结果是2187

而如果使用快速幂的话是log级别的时间复杂度

因为7的二进制代码是111,我们可以通过二进制代码来一个个的与3相乘。

即：

3^1=3

3^2=9

3^4=81，

最后我们把这三个数相乘的话就能得到3^7的值是2187，实际上我们执行循环的次数只有三次，这题让求mod 6之后的数据我们可以在每次相乘得到结果后%6这样不会影响最后结果，而且还能防止溢出，因为最后的数据会很大，所以必须不断乘不断模，要注意的一点是这个两个数相乘一次的时候可能会溢出，所以我转换为了long long型的数据防止溢出，注意下方代码注释。

数据比较小还看不出来什么，

如果是3 10000 6的话可以自行试一下，会发现只有二十次的相乘。

c++版本:

#include <iostream>

using namespace std;

int main()
{
 int a,b,p;
 int res=1%p;
     while(b)
    {
        if(b&1)res=res*1ll*a%p;
        a=a*1ll*a%p;
        b>>=1;
    }
    cout<<res;    
    return 0;
}

c语言版本:

#include <stdio.h>

int main()
{
 int a,b,p;
scanf("%d%d%d",&a,&b,&p);
 int res=1%p;//结果值初始化为1,mod p的原因是要处理123456789 0 1这种特殊情况
     while(b)
    {
        if(b&1)res=res*1ll*a%p;//b的个位是1的话就相乘，乘1ll的话会把res*a这个表达式的值转换为long long int型的数据
        a=a*1ll*a%p;//把下一位准备好，要注意%p
        b>>=1;//个位去掉
    }
    printf("%d",res);   
    return 0;
}

这里说说为什么要a=a*1ll%p，举个例子

3 5 2

5:101

3^1=3

3^4=27

1.res=1*3%2=1

a=3*3%2=1

2.a=1*1%2=1

3.res=1*1%2=1

a=1

从这里我们可以发现res保存的是每次相乘取模之后的结果，对a再进行一次取模运算的话不会影响结果，但是如果不取模的话a会溢出的。（这里因为题意是取模运算，如果只是普通的相乘的话不要Mod啥。直接乘就行。)

下面我们来看第二道题:

这题的数据范围是a,b,p<=10e8，也就最大是64位的整数，两个64位的整数相乘能得到一个128位的整数，如果用高精度算法的话肯定能做出来，但是这题有一个Mod p,说明我们可以通过快速幂来去算出来。

举个例子

3 5 7

5二进制:101

3*1=3

3*4=12

12+3=15

 

c语言版本:

#include <stdio.h>

int main()
{
    long long int a,b,p,res=0;
    scanf("%lld%lld%lld",&a,&b,&p);
    while(b)
    {
        if(b&1)res=(res+a)%p;
        a=2*a*1ll%p;
        b>>=1;
    }
    printf("%lld",res);
    return 0;
}

以上就是快速幂了，下面还要说几个地方。

/*

1.用异或来实现配偶
0,1
2,3
4,5
6,7

0^1=1 1^0=1
2^1=3 3^1=2
4^1=5 5^1=4

e[index]
e[index^1]
(在图论里挺有用的)

2.lowbit运算

假如有一个二进制数，我们想就求出来最右边的1到末尾的值，
即
x=1110001000 =>1000
~x+1=0001111000
(~x+1)&x=1000
所以lowbit运算可以这么总结
int lowbit(int n)
{
return (-x)&x;
}
3.
n>>j&1 判断二进制第j位是否为1(二进制)
i-(1<<j)>>k&1 把i减去2的j次方之后判断第k位是否为1
M=1<<20 2的20次方
state|1<<j  把state的第j位置1


*/