比赛-thh学长的训练赛 (16 Aug, 2018)
1.) 欧拉回路
打表找规律,打了斯特林,打了组合数,找不出,弃疗。我太蠢了,这题和前面这些东西没什么关系啊。考虑 \(n\) 个点时的完全图 \(n \cdot (n - 1) / 2\) 条边,每条边选与不选两种决策,所以可以得到神奇的数字 \(2^{n(n-1)/2}\) ,而答案就是 \(2^{(n-1)(n-2)/2}\) 。具体证的话,对于 \(n\) 个点的图,拿出一个点,把 \(n-1\) 个点胡乱连边,方案数 \(2^{(n-1)(n-2)/2}\) 。然后让单独拿出来的点与其他点连边。如果某点度数是奇度就把单独拿出来的点与这个点连一条边,否则不连。这样的话,连边方案是唯一的。可以得到:A.) 除了单独拿出来的点外,其他点的度数在连边后为偶数。同时:B.) 根据无向图性质,所有点度数之和为偶数。由 A 和 B 可以得到,单独拿出来的点的度数也为偶数。因为每个点度数都是偶数,所以这样构造的 \(n\) 个点的图含有欧拉回路。并且这样显然会不遗漏地讨论完所有方案。呼……
#include <cstdio>
typedef long long ll;
const ll MOD = 998244353ll;
ll mul(ll a, ll b)
{
ll t = 0;
while (b) {
if (b & 1)
t = (t + a) % MOD;
b >>= 1;
a = (a << 1) % MOD;
}
return t;
}
ll mont(ll a, ll b)
{
ll t = 1;
a %= MOD;
while (b) {
if (b & 1)
t = mul(t, a);
b >>= 1;
a = mul(a, a);
}
return t;
}
int main()
{
ll n;
scanf("%lld", &n);
if (n == 1) {
printf("1\n");
return 0;
}
ll ans;
if (n & 1)
ans = mont(mont(2, (n-1)/2), n-2);
else
ans = mont(mont(2, (n-2)/2), n-1);
printf("%lld\n", ans);
return 0;
}
2.) 乌鸦坐飞机
ljh 大佬考场上想的线段树优化建图,多么窒息的操作……我非常幸运地一眼看出考倍增,然后玄学写错爆 0 ,考完重写了一次就 A 了,暂时不知道错哪。
推一下 dp 方程 \(f_i = min\{f_j + w_i\}\) ,要求 \(j\) 在 \(i\) 到根的路径上。发现是树上极值查询,所以就是倍增啊。
#include <stdio.h>
#include <stack>
#include <cstring>
#include <vector>
#include <cmath>
#include <ctype.h>
using namespace std;
typedef long long ll;
const int _N = 320000;
const int _lgN = 21;
vector<int> G[_N];
int dad[_N][_lgN], D[_N], W[_N], depth[_N], LG[_N], N, lgN;
ll f[_N][_lgN], ans[_N];
int moveup(int p, int s)
{
for (int i = 0; i <= lgN; ++i)
if ((s >> i) & 1) p = dad[p][i];
return p;
}
void dfs(int p, int d)
{
if (p != 1) {
depth[p] = depth[d] + 1;
dad[p][0] = d;
f[p][0] = ans[d];
for (int i = 1; i <= lgN; ++i) {
dad[p][i] = dad[dad[p][i-1]][i-1];
f[p][i] = min(f[p][i-1], f[dad[p][i-1]][i-1]);
}
int x = min(D[p], depth[p]);
int tmp = moveup(p, x-(1<<LG[x]));
ans[p] = min(f[tmp][LG[x]], f[p][LG[x]]) + W[p];
}
for (int i = G[p].size()-1; i >= 0; --i) {
int v = G[p][i];
if (v != d) dfs(v, p);
}
return;
}
int main()
{
scanf("%d", &N);
for (int x, y, i = 1; i < N; ++i) {
scanf("%d%d", &x, &y);
G[x].push_back(y);
G[y].push_back(x);
}
for (int i = 1; i <= N; ++i)
scanf("%d%d", &W[i], &D[i]);
lgN = log2(N + 0.5) + 1;
LG[0] = -1;
for (int i = 1; i <= N; ++i)
LG[i] = LG[i >> 1] + 1;
dfs(1, 0);
int T;
scanf("%d", &T);
while (T--) {
int x;
scanf("%d", &x);
printf("%lld\n", ans[x]);
}
return 0;
}
3.) 四边形上树
真是妙啊。先考虑暴力,存在四边形的充要条件是任意三边和大于第四边。枚举四边形的最长边 \(d\) ,只要有 \(a+b+c > d\) 就可以了。直接把两点路径上的点排序,然后贪心地找到比枚举的 \(d\) 小的,最大的三个数,判断是否可以作为边长构成四边形。然后就是很精妙的操作了。推理可得,要构造出连续很多个数,还要使他们总是满足 \(a+b+c \leq d\) 的话,当这些数的个数为 40 时,最大的数已经超过 int 上限了。所以对于路径长度 \(>40\) 的询问直接输出 \(Yes\) ,否则暴力排序并判断就好了。
#include <stdio.h>
#include <vector>
#include <algorithm>
const int _N = 120000;
const int X = 40;
using namespace std;
#define SC(a, b) (static_cast<a>(b))
typedef long long ll;
vector<int> G[_N];
int dad[_N], A[_N], B[_N], Bcnt, N;
bool mk[_N];
void dfs(int p)
{
for (int i = G[p].size()-1; i >= 0; --i) {
int v = G[p][i];
if (v != dad[p]) {
dad[v] = p;
dfs(v);
}
}
return;
}
int main()
{
int T;
scanf("%d", &N);
for (int i = 1; i <= N; ++i)
scanf("%d", &A[i]);
for (int x, y, i = 1; i < N; ++i) {
scanf("%d%d", &x, &y);
G[x].push_back(y), G[y].push_back(x);
}
scanf("%d", &T);
dfs(1);
while (T--) {
int ins, x, y;
scanf("%d%d%d", &ins, &x, &y);
if (ins == 1) {
A[x] = y;
} else {
int a, i;
for (a = x, i = 1; i <= X; ++i)
mk[a] = 1, a = dad[a];
for (a = y, i = 1; i <= X && !mk[a]; ++i)
a = dad[a];
if (!mk[a] || !a) {
printf("Yes\n");
for (a = x, i = 1; i <= X; ++i)
mk[a] = 0, a = dad[a];
continue;
}
int lca = a;
Bcnt = 0;
for (a = x; a != lca; a = dad[a])
B[++Bcnt] = A[a];
for (a = y; a != lca; a = dad[a])
B[++Bcnt] = A[a];
B[++Bcnt] = A[lca];
if (Bcnt < 4) {
printf("No\n");
} else {
sort(B+1, B+1+Bcnt);
for (i = 4; i <= Bcnt; ++i)
if (SC(ll, B[i-3])+B[i-2]+B[i-1] > SC(ll, B[i])) break;
if (i == Bcnt+1)
printf("No\n");
else
printf("Yes\n");
}
for (a = x, i = 1; i <= X; ++i)
mk[a] = 0, a = dad[a];
}
}
return 0;
}