计蒜客12月CSP-S解题报告（线性筛+状压DP+树状数组逆序对+2-SAT）

A. 幸运数字

#include <bits/stdc++.h>
using namespace std;
const int maxn = 1e7 + 10;
int t;
bool notprime[maxn];
int cnt = 0;
int ans[maxn];
int prime[maxn];
int tmp[maxn];
//存下之前已经确定的质数 然后相乘 把相乘得到的数标记好
void init() {
	for (int i = 2; i <= maxn; i++) {//质数a
		if (!notprime[i]) {
			tmp[i] = 1;
			prime[cnt++] = i;
			for (int j = 0; j < cnt; j++) {//枚举质数b
				if (i * prime[j] > maxn) break;
				tmp[i * prime[j]] = 1;//标记：质数a*质数b满足要求
				notprime[i * prime[j]] = 1;//标记：质数a*质数b不为质数
			}
		}
		for (int j = 0; j < cnt; j++) {//板子 不解释
			if (i * prime[j] > maxn) break;
			notprime[i * prime[j]] = true;
			if (i % prime[j] == 0) break;
		}
	}
	for (int i = 1; i <= maxn; i++) ans[i] = ans[i - 1] + tmp[i];
}
int main() {
	freopen("lucky_number.in", "r", stdin);
	freopen("lucky_number.out", "w", stdout);
	cin >> t;
	init();
	while (t--) {
		int l, r;
		cin >> l >> r;
		cout << ans[r] - ans[l - 1] << endl;
	}
	return 0;
}

B. 蒜头君健身

直接上代码   解释得很全了

#include <bits/stdc++.h>
using namespace std;
const int maxn = 21, maxm = 1100000;
int a[maxn][maxn], dp[maxm];
int n;
int main() {
	cin >> n;
	for (int i = 0; i < n; i++) 
		for (int j = 0; j < n; j++) 
			cin >> a[i][j];
	for (int sta = 0; sta < (1 << n); sta++) //枚举状态 
		for (int i = 0; i < n; i++) { //枚举该状态下的项目 
			if (!(sta & (1 << i))) { //如果该项目没有被确定 
				int sum = 0;
				for (int j = 0; j < n; j++) { //枚举j直到j项目确定 
					if (sta & (1 << j)) sum += a[j][i]; //如果该项目已经确定了  那么就把i项目放到j项目之后 
				}
				dp[sta | (1 << i)] = max(dp[sta | (1 << i)], dp[sta] + sum); //状态转移：dp[sta | (1 << i)]表示在该状态的基础上 又增加了第i个项目 训练效果之和
				// dp[sta | (1 << i)]就等于自身（dp[sta | (1 << i)]）和原先状态再加上把i项目放到j项目之后训练效果变化之和（dp[sta] + sum） 
			}
		}
	cout << dp[(1 << n) - 1] << endl;
	return 0;
}
/*
考虑进行状态压缩 DP，设 dp[sta] 为已经确定的训练项目状态是 sta 时的训练效果之和的最大值。 
状态为 sta 时，如果 i 项目在当前状态下是还没有被确定顺序的项目，
那么将 i 项目安排在 当前状态下已经确定顺序的 j 项目后面，
则此时增加 i 项目可以增加的训练效果之和 sum 为：
int sum = 0; for (int j = 0; j < n; j++) { if (sta & (1 << j)) { sum += a[j][i]; } } 
那么该状态下增加 i 项目后，训练效果之和为 dp[sta | (1 << i)] = max(dp[sta | (1 << i)], dp[sta] + sum); 
那么最终的答案为 dp[(1 << n) - 1]
*/