ET算法【位运算与二分】

一、位运算

P1068 - 我们需要0 - ETOJ (eriktse.com)

P1083 - Mex and Xor - ETOJ (eriktse.com) 思维题

1.记一下异或的性质。

①结合律：(a^b)^c=a^(b^c)

②交换律：a^b=b^a

没有分配律

a^a=0        a^0=a

若a^b=c，则有a^c=b和b^c=a.

2.a+b=2(a&b)+(a^b)    a&b进位  a^b加法

1.Problem - 1300C - Codeforces

发现只要该位上有两个以上的1，最终就会变为0

找到第一位只有一个数的二进制位为1的数，放到第一个，剩下的顺序输出

#include<bits/stdc++.h>
using namespace std;
#define qio ios::sync_with_stdio(0), cin.tie(0), cout.tie(0);
typedef long long ll;
typedef double db;
void solve();
#define int ll
signed main() {
	qio
	int T = 1;
//	cin >> T;
	while (T--)solve();
}

const int N = 1e5 + 10;
int a[N];
void solve() {
	int n; cin >> n;
	for (int i = 1 ; i <= n; i++) {
		cin >> a[i];
	}
	for(int i = 30; i >= 0; i--) {
		int t = 1 << i;
		int cnt = 0, pos = 1;
		for (int j = 1; j <= n; j++) {
			if (a[j] & t) {
				cnt++;
				pos = j;
			}
		}
		if (cnt == 1) {
			swap(a[pos], a[1]);
			break;
		}
	}
	for(int i = 1; i <= n; i++) {
		cout << a[i] << " ";
	}	
}

二、二分查找

1.模板

P1088 - 查找 - ETOJ (eriktse.com)

二分套此模板：

if (a[mid] < x)此处可根据需要修改，此时搜索到的是目标元素中任意一个的位置

if (a[mid] <= x)此时搜索到的是相同目标元素中最左侧的一个，输出r

if (a[mid] >= x)此时搜索到的是相同目标元素中最右侧的一个，输出l

带等号即为遇到相同元素时仍然推进l或者r所以能推到端点处，推进哪个最终就输出哪个

#include<bits/stdc++.h>
using namespace std;
#define qio ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
const int N = 2e5 + 10;
int a[N];
int main() {
	qio
	int n, q;
	cin >> n >> q;//n为数组元素个数，q为询问次数
	for (int i = 1; i <= n; i++) cin >> a[i];	
	while (q--) {
		int l = 0, r = n;
		int x;
		cin >> x;
		while (l + 1 != r) {
			int mid = (l + r) >> 1;
			if (a[mid] < x)l = mid;
			else r = mid;
		}
		if (a[r] == x) cout << r << ' ';
		else cout << "-1" << ' ';
	}
}

---

2.思维提升：

P1048 - 进击的奶牛 - ETOJ (eriktse.com)

寻找变量之间的关系，建立平面直角坐标系帮助理解。此题中，纵坐标为放置的奶牛数，横坐标为最大的最近距离，因此搜索的是目标奶牛数量c对应的横坐标的最大值。抽象理解为奶牛数量数组4 4 4 4 3 3 3 2 2 2 2 1 1 1 1 1 1 1中，找到最右端c=3对应的位置 

#include<bits/stdc++.h>
using namespace std;
#define qio ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
using ll = long long;
int a[N], n, c;
const int N = 2e5 + 10;
int check(int mid) {
	int res = 0, pre = -1e9;
	for (int i = 1; i <= n; i++) {
		if (a[i] - pre >= mid) {
			res++;
			pre = a[i];
		}
	}
	return res;
}
int main() {
	qio
    cin >> n >> c;
	for (int i = 1; i <= n; i++)cin >> a[i];
	sort(a + 1, a + n + 1);
	int l = 0, r = 1e9;
	while (l + 1 != r) {
		int mid = (l + r) >> 1;
		if (check(mid) >= c)l = mid;
		else r = mid;
	}
	cout << l << '\n';
}

---

STL库查找函数：

lower_bound(a,a+n,k)         返回第一个大于等于k的元素在a到a+n中出现的地址

upper_bound(a,a+n,k)        返回第一个大于k的元素在a到a+n中出现的地址

C++ 折半查找之 lower_bound 和 upper_bound_RYW2022的博客-CSDN博客

3.D - Pond (atcoder.jp)

求第几大的数——最小的中位数，二分假设答案为x，比x大的变为1，小的变为0，二维前缀和统计个数

#include<bits/stdc++.h>
using namespace std;
#define qio ios::sync_with_stdio(0), cin.tie(0), cout.tie(0);
typedef long long ll;
#define int ll

const int N = 800 + 10;
int mp[N][N];
int n, k;
int pre[N][N];

bool isok(int mid, int x) {
	for (int i = 1; i <= n; i++) {
		for (int j = 1; j <= n; j++) {
			pre[i][j] = mp[i][j] > mid;
			pre[i][j] = pre[i][j] + pre[i - 1][j] + pre[i][j - 1] - pre[i - 1][j - 1];
		}
	}
	for (int i = k; i <= n; i++) {
		for (int j = k; j <= n; j++) {
			int res = pre[i][j] - pre[i - k][j] - pre[i][j - k] + pre[i - k][j - k];
			if (res <= x) return 1;
		}
	}
	return 0;
}

int main() {
    qio
	cin >> n >> k;
	int x = (k * k) / 2;
	int l = -1, r = 0;
	for (int i = 1; i <= n; i++) {
		for (int j = 1; j <= n; j++) {
			cin >> mp[i][j];
			r = max(r,mp[i][j]);
		}
	}
	
	while (l + 1 != r) {
		int mid = (l + r) >> 1;
		if (isok(mid, x))r = mid;
		else l = mid;
	}
	cout << r;
}