我还是太菜了,只签了一道签到,实现了一个正解,打了点暴力和对拍。
B. Bitwise Maximization
异或线性基,模板记住了。
cpp
#include <iostream>
using namespace std;
const int N = 500010;
int a[N], p[30];
void insert(int x) {
for (int i = 29; i >= 0; --i) {
if (!(x >> i & 1)) continue;
if (!p[i]) {
p[i] = x;
break;
}
else x ^= p[i];
}
}
int main() {
ios::sync_with_stdio(0);
cin.tie(0), cout.tie(0);
int T;
cin >> T;
while (T--) {
fill(p, p + 30, 0);
int n, s = 0;
cin >> n;
for (int i = 1; i <= n; ++i) {
cin >> a[i];
s ^= a[i];
}
for (int i = 1; i <= n; ++i) {
insert(a[i] & (~s));
}
int res = 0;
for (int i = 30; i >= 0; --i) {
res = max(res, res ^ p[i]);
}
cout << (res << 1) + s << endl;
}
return 0;
}F. Fabulous Tree
是一个比较抽象的树形 DP
G. GCD Graph
打了一堆表之后发现,1e7 以内素数间隔最大只有 150,间隔以内的暴力 DP 求,以外的统计 互质 ? 1 : 2 步,求和两部分。
cpp
#include <bits/stdc++.h>
#define vi vector<int>
using namespace std;
const int N = 10000010;
typedef long long LL;
int dis[N], ne[N];
int v[N], prime[N], len;
void init() {
int n = 10000000;
for (int i = 2; i <= n; ++i) {
if (v[i] == 0) {
v[i] = i;
prime[++len] = i;
}
for (int j = 1; j <= len; ++j) {
if (prime[j] > v[i] || prime[j] > n / i) break;
v[i * prime[j]] = prime[j];
}
}
}
int getp(int n) {
int l = 1, r = len;
while (l < r) {
int mid = l + r + 1 >> 1;
if (prime[mid] <= n) l = mid;
else r = mid - 1;;
}
return prime[l];
}
vi get(int a) {
vi ret;
for (int i = 1; i * i <= a; i++) {
if (a % i == 0) {
ret.push_back(i);
if (i != a / i)
ret.push_back(a / i);
}
}
return ret;
}
int get1(int l, int r,int n) {
vi a = get(n);
sort(a.begin(), a.end());
vi cnt(a.size());
for (int i = 0; i < a.size(); i++) {
cnt[i] = r / a[i] - (l - 1) / a[i];
}
for (int i = a.size() - 1; i >= 0; i--) {
for (int j = i + 1; j < a.size(); j++) {
if (a[j] % a[i] == 0)
cnt[i] -= cnt[j];
}
}
return cnt[0];
}
int main() {
ios::sync_with_stdio(0);
cin.tie(0), cout.tie(0);
init();
int T;
cin >> T;
while (T--) {
int l, r, n;
cin >> l >> r >> n;
int p = getp(n);
if (r <= p) {
int t = get1(l, r, n);
cout << (r - l + 1) * 2 - t << '\n';
}
else if (p <= l) {
LL res = 0;
memset(dis + l, 0x3f, sizeof(int) * (n - l + 1));
dis[n] = 0;
for (int i = n - 1; i >= l; --i) {
for (int j = i + 1; j <= n; ++j) {
dis[i] = min(dis[i], __gcd(i, j) + dis[j]);
}
}
for (int i = l; i <= r; ++i) {
res += dis[i];
}
cout << res << '\n';
}
else {
LL res = 0;
memset(dis + p, 0x3f, sizeof(int) * (n - p + 1));
dis[n] = 0;
for (int i = n - 1; i > p; --i) {
for (int j = i + 1; j <= n; ++j) {
dis[i] = min(dis[i], __gcd(i, j) + dis[j]);
}
}
for (int i = p + 1; i <= r; ++i) {
res += dis[i];
}
cout << res + (p - l + 1) * 2 - get1(l, p, n) << '\n';
}
}
return 0;
}H. Hyperspace Pairing
K. Kindergarten
L. Lazy Shuffling
M. Maybe Connected
签到,尽可能连一个最大的菊花。
cpp
#include <iostream>
using namespace std;
typedef long long LL;
int main() {
ios::sync_with_stdio(0);
cin.tie(0), cout.tie(0);
int T;
cin >> T;
while (T--) {
LL n, m;
cin >> n >> m;
if (m >= n - 1) cout << (n - 1) * (n - 2) / 2 - m + (n - 1) << '\n';
else cout << (m - 1) * m / 2 << '\n';
}
return 0;
}N. Narrow to Median
cpp