0X01 位运算

目录

A.快速幂取模

题目链接

题解

代码

B.64位乘法取模

题目链接

题解

C.最短Hamilton路径

题目链接

题解

代码

---

A.快速幂取模

题目链接

　　快速幂取模

题解

　　把b看作2的次方倍的拆分，然后计算。

　　注意点：如果求1的b次方去mod 1，快速幂开始时res赋值应该为(1%p)，而不是1，因为如果p也等于1会造成结果出错。

代码

ll pow(ll a,ll b,ll p) {
	ll res = 1%p;
	for(;b;b>>=1) {
		if(b&1) res = res*a%p; //如果b在这一位二进制为1，说明需要乘上a的对应次方
		a = a*a%p; 
	}
	return res;
}

---

B.64位乘法取模

题目链接

　　64位乘法取模

题解

　　方法1：依据之前的快速幂的思想，(a*b)%p 可以把b进行二进制拆分，然后每次mod p，这样求出来的就是最终结果。

ll mul(ll a,ll b,ll p) {
	ll ans = 0;
	for(;b;b>>=1) {
		if(b&1)	ans = (ans+a)%p;
		a = a*2%p;
	}	
	return ans;
}

　　方法2：(a*b)%p = ((a%p)*(b%p))%p ,即使这样了a%p和b%p依然有可能是很大的数，它们的乘积会溢出，所以我们用

(a%p)*(b%p)-[(a%p)*(b%p)/p]*p , 其中(a%p)*(b%p)/p这一步用long double来保存（有效位数在18-19）位，符合条件，然后用long long来舍去。

ll mul(ll a,ll b,ll p) {
	a %= p;
	b %= p;
	ll c = (long double)a*b/p;
	cout<<c<<endl;
	ll ans = a*b-c*p;
	if(ans>=p)	ans -= p;
	else if(ans<0)	ans += p;
	return ans;
}

---

C.最短Hamilton路径

题目链接

　　最短Hamilton路径

题解

　　典型的TSP问题（编号从0-n-1），在规模小的情况下采用状压DP求解。

　　用f[i][j]来表示当前的节点为j，已经便利过的状态是i（的二进制）。那么起点f[1][0] = 0,终点就是f[(1<<n)-1][n-1]；

　　递推式子即为 f[i][j] = min(f[i][j],f[i^(1<<j)][k]+w[k][j]) , 其中 i^(1<<j) 表示的是j这个点之前没访问过，所以为0（用异或的形式来表示i的前一个状态），k表示 i^(1<<j) 所表示的位置中任何一个为1的位置k。

代码

int hamilton(int n) {
	memset(f,inf,sizeof(f));
	f[1][0] = 0; //遍历状态为0001,当前点为0，所需花费为0 
	for(int i=1;i<(1<<n);i++) {
		for(int j=0;j<n;j++) {
			if((i>>j)&1) { //如果j这个点是当前正在遍历的点 
				for(int k=0;k<n;k++)
					if(((i^(1<<j))>>k)&1) //上一刻所处的位置是 i xor (1<<j) 中任意一个是1的数位k 
						f[i][j] = min(f[i][j],f[i^(1<<j)][k]+w[k][j]); //上一刻状态i的第j位必然为0 
			}
		}
	}
	return f[(1<<n)-1][n-1];
}