HNOI2019 DAY1 题解

最新推荐文章于 2022-10-11 17:49:17 发布
最新推荐文章于 2022-10-11 17:49:17 发布
阅读量555 收藏 1
点赞数 2
分类专栏： 比赛
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/WAautomaton/article/details/89302539
版权

鱼

题目链接

题解：
容易发现等腰三角形个数似乎是 O ( n 2 ) O(n^2) O(n2) 的？于是可以计算出任意两个点作为底边，向两个方向的等腰三角形个数。之后显然枚举中心点，然后把其它点极角排序，算出所有等腰三角形的中垂线，然后再枚举一个作为鱼尾，二分+前缀和查询即可。复杂度 O ( n 2 l o g n ) O(n^2logn) O(n2logn)。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int, int> P;

const double PI = 3.14159265358979324, EPS = 1E-10;
const int MAXN = 1005, MAXM = 1000005;
struct Point { int x, y, id; ll d; double ang; } po[MAXN][MAXN];
bool cmp(const Point &a, const Point &b) {
	return a.d == b.d ? a.ang < b.ang : a.d < b.d;
}
int xx[MAXN], yy[MAXN], app[MAXN][MAXN][2], sum[MAXM], n;
struct Node { int s; double ang; } arr[MAXM];
double angs[MAXM];
double range(double x) {
	return x > PI ? x - PI * 2 : (x - EPS < -PI ? x + PI * 2 : x);
}
int main() {
	scanf("%d", &n);
	for (int i = 1; i <= n; i++) scanf("%d%d", xx + i, yy + i);
	ll res = 0;
	for (int i = 1; i <= n; i++) {
		int cnt = 0;
		for (int j = 1; j <= n; j++) if (i != j) {
			po[i][++cnt].id = j;
			po[i][cnt].x = xx[j] - xx[i];
			po[i][cnt].y = yy[j] - yy[i];
			po[i][cnt].d = (ll)po[i][cnt].x * po[i][cnt].x + (ll)po[i][cnt].y * po[i][cnt].y;
			po[i][cnt].ang = range(atan2(po[i][cnt].y, po[i][cnt].x));
		}
		sort(po[i] + 1, po[i] + n, cmp);
		for (int j = 1; j < n; j++) {
			int k = j;
			while (k < n - 1 && po[i][k + 1].d == po[i][j].d) ++k;
			for (int l = j; l <= k; l++)
			for (int o = l + 1; o <= k; o++) {
				int a = po[i][l].id, b = po[i][o].id;
				if (b < a) swap(a, b);
				ll c = (ll)(xx[i] - xx[a]) * (yy[b] - yy[a]) - (ll)(yy[i] - yy[a]) * (xx[b] - xx[a]);
				if (c == 0) continue;
				++app[a][b][c > 0];
			}
			j = k;
		}
	}
	for (int i = 1; i <= n; i++) {
		int cnt = 0;
		for (int j = 1; j < n; j++) {
			int k = j;
			while (k < n - 1 && po[i][k + 1].d == po[i][j].d) ++k;
			for (int l = j; l <= k; l++)
			for (int o = l + 1; o <= k; o++) {
				int a = po[i][l].id, b = po[i][o].id;
				if (b < a) swap(a, b);
				ll c = (ll)(xx[i] - xx[a]) * (yy[b] - yy[a]) - (ll)(yy[i] - yy[a]) * (xx[b] - xx[a]);
				if (c == 0) continue;
				if (po[i][o].ang - po[i][l].ang - EPS > PI)
					arr[++cnt].ang = range((po[i][l].ang + po[i][o].ang) * 0.5 + PI);
				else arr[++cnt].ang = (po[i][o].ang + po[i][l].ang) * 0.5;
				arr[cnt].s = app[a][b][c < 0];
			}
			j = k;
		}
		sort(arr + 1, arr + cnt + 1, [](const Node &a, const Node &b) { return a.ang < b.ang; });
		for (int l = 1; l <= cnt; l++) angs[l] = arr[l].ang;
		for (int l = 1; l <= cnt; l++) sum[l] = sum[l - 1] + arr[l].s;
		for (int j = 1; j < n; j++) {
			int k = j;
			while (k < n - 1 && po[i][k + 1].d == po[i][j].d) ++k;
			for (int l = j; l <= k; l++)
			for (int o = l + 1; o <= k; o++) {
				double a, b;
				if (po[i][o].ang - po[i][l].ang - EPS > PI)
					a = range(po[i][o].ang - PI * 0.5), b = range(po[i][l].ang + PI * 0.5);
				else if (po[i][o].ang - po[i][l].ang + EPS < PI)
					a = range(po[i][o].ang + PI * 0.5), b = range(po[i][l].ang - PI * 0.5);
				else continue;
				if (a > b) swap(a, b);
				if (b - a > PI) {
					int c = lower_bound(angs + 1, angs + 1 + cnt, a - EPS) - angs - 1;
					int d = lower_bound(angs + 1, angs + 1 + cnt, b + EPS) - angs - 1;
					res += sum[cnt] - sum[d] + sum[c];
				} else {
					int c = lower_bound(angs + 1, angs + 1 + cnt, a + EPS) - angs - 1;
					int d = lower_bound(angs + 1, angs + 1 + cnt, b - EPS) - angs - 1;
					res += sum[d] - sum[c];
				}
			}
			j = k;
		}
		//printf("%d %lld\n", i, res);
	}
	printf("%lld\n", res << 2);
	return 0;
}

JOJO

题目链接

题解:
没有2操作的话应该是比较好做的吧。注意到相邻插入的字符都不一样，因此我们对二元组做KMP。跳 f a i l fail fail 时，我们需要找到一个和当前完全一样的二元组才行，否则由于后面插入的字符不一样，就一定不能接下去。特殊的，当跳到第一个二元组时，如果字符一样且当前二元组长度大于第一个二元组长度，那么 f a i l fail fail 可以指向第一个二元组。

考虑如何统计答案。在跳 f a i l fail fail 的过程中，每跳一次可能会在当前二元组上多匹配一段区间，直接加上贡献即可。于是此部分可以 O ( n ) O(n) O(n) 解决。

但是这个复杂度是均摊的，可持久化的时候就不对了（但听说暴力直接过了？）。考虑优化跳 f a i l fail fail 的过程，即对于每个二元组，记 f ( i , j , k ) f(i,j,k) f(i,j,k) 表示从 i − 1 i-1 i−1 位置跳二元组 ( j , k ) (j,k) (j,k) 会到达的节点。更新的话直接从 f ( f a i l i ) f(fail_i) f(faili​) 复制一遍，并且把 f ( i , j , k ) f(i,j,k) f(i,j,k) 指向 i i i。计算答案的话，当前点对后面产生的贡献实际上是一个等差数列，区间打 s e t set set 标记即可。

于是离线建树，然后 d f s dfs dfs 的时候拿可持久化线段树维护就行了，复杂度 O ( n l o g n ) O(nlogn) O(nlogn)。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int, int> P;

const int MAXN = 100005, MAXT = 10000005;
struct Node { Node *ls, *rs; ll sum; int tag, id, l, r; }
	nd[MAXT], *rt[MAXN][26], *tot = nd;
struct Edge { int to, next; } edge[MAXN];
ll ans[MAXN];
int head[MAXN], id[MAXN], val[MAXN], len[MAXN], n, m, turn, ecnt;
int pre[MAXN], mx[MAXN][26], lqs[MAXN];
char opt[5];
Node *newnode(int l, int r, int t) {
	Node *p = tot++; p->sum = (ll)(r - l + 1) * t;
	p->l = l, p->r = r, p->tag = t;
	return p;
}
Node *copy(Node *d, int l, int r) {
	Node *p = tot++;
	if (d != nullptr) return &(*p = *d);
	p->l = l, p->r = r; return p;
}
void pushdown(Node *d) {
	if (!d->tag) return;
	int mid = (d->l + d->r) >> 1;
	d->ls = newnode(d->l, mid, d->tag);
	d->rs = newnode(mid + 1, d->r, d->tag);
	d->tag = 0;
}
void modify(int r, int val, int id, Node *&d, int a = 1, int b = m) {
	d = copy(d, a, b);
	if (d->l > r) return;
	if (d->r < r) { d->sum = (ll)(d->r - d->l + 1) * (d->tag = val); return; }
	if (d->l == d->r) { d->sum = d->tag = val, d->id = id; return; }
	pushdown(d);
	int mid = (a + b) >> 1;
	modify(r, val, id, d->ls, a, mid);
	modify(r, val, id, d->rs, mid + 1, b);
	d->sum = d->ls->sum + d->rs->sum;
}
void query(int r, ll &ans, int &id, Node *d) {
	if (d == nullptr || d->l > r) return;
	if (d->r < r) { ans += d->sum; return; }
	if (d->l == d->r) { ans += d->sum, id = d->id; return; }
	pushdown(d);
	query(r, ans, id, d->ls);
	query(r, ans, id, d->rs);
}
ll get_sum(int x) { return (ll)x * (x + 1) >> 1; }
void addedge(int u, int v) {
	edge[++ecnt] = (Edge) { v, head[u] };
	head[u] = ecnt;
}
void dfs(int u, int dep) {
	int fail = 0;
	lqs[dep] = val[u];
	if (dep > 1) {
		ans[u] += get_sum(min(len[u], mx[dep][val[u]]));
		query(len[u], ans[u], fail, rt[dep][val[u]]);
		if (!fail && len[u] > pre[1] && lqs[1] == lqs[dep])
			ans[u] += (ll)pre[1] * max(0, len[u] - mx[dep][val[u]]), fail = 1;
	} else if (dep == 1) ans[u] += get_sum(len[u] - 1);
	else fail = -1;
	mx[dep][val[u]] = max(mx[dep][val[u]], len[u]);
	if (dep > 0) modify(len[u], pre[dep - 1], dep, rt[dep][val[u]]);
	for (int i = head[u]; i; i = edge[i].next) {
		for (int j = 0; j < 26; j++) {
			mx[dep + 1][j] = mx[fail + 1][j];
			rt[dep + 1][j] = rt[fail + 1][j];
		}
		int v = edge[i].to;
		pre[dep + 1] = pre[dep] + len[v];
		ans[v] = ans[u];
		dfs(v, dep + 1);
	}
}
int main() {
	scanf("%d", &n);
	for (int i = 1; i <= n; i++) {
		int a, b; scanf("%d%d", &a, &b);
		if (a == 1) {
			addedge(id[i - 1], id[i] = ++turn);
			scanf("%s", opt);
			len[turn] = b, val[turn] = opt[0] - 'a';
			m = max(m, b);
		} else id[i] = id[b];
	}
	dfs(0, 0);
	for (int i = 1; i <= n; i++) printf("%lld\n", ans[id[i]] % 998244353);
	return 0;
}

多边形

题目链接

题解：
凸多边形的三角剖分？赶紧转二叉树……具体的最小左旋建对偶图即可。然后会发现，题目中的“旋转”操作实际上就是二叉树的旋转。考虑令靠近 1 − n 1-n 1−n 那条边的三角形为根，对于每个三角形，除了父节点外还有两条边，令边上点编号较小的为左儿子，较大的为右儿子。

不难发现，题目显然让我们旋转到所有边都从 n n n 出发为止，否则可以证明一定能够继续旋转。也就是说，我们要旋转到二叉树变成一条除叶节点外只有右子树的链。考虑最优策略最多每次能让一个点进入最右边的那条链，因此答案就是 n − 2 − n-2- n−2−最右边链的长度。显然，一个父节点必须在子节点之前进入链，已经在链上的节点可以不考虑，也就是不计入子树 s i z e size size。于是就是一个很经典的计数 d p dp dp 了， f [ i ] f[i] f[i] 表示 i i i 的子树中所有不在链上的节点，组成的不同拓扑序个数。于是大概就是儿子的 f f f 乘起来再乘上一个组合数。

最后考虑每次旋转后如何动态维护这个 d p dp dp，动态 d p dp dp （雾）？容易发现，根节点的 d p dp dp 值实际上就是所有节点把儿子插起来的组合数的乘积。于是如果当前旋转没有把任何点旋转到链上，那么直接除掉那两个变化点的贡献，再乘上他们的新贡献即可。否则的话，第一问的答案减1，然后从它开始到根节点所有点的 s i z e size size 减1，于是我们还要维护一个新的 d p dp dp，表示从根到某个点上的 s i z e size size 全部减1，答案会乘上多少。然后就可以在 O ( n l o g n ) O(nlogn) O(nlogn) 的时间内解决了。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int, int> P;

const int MAXN = 100005, MAXM = 400005, MOD = 1000000007;
set<int> edge[MAXN]; map<P, int> mp;
int vis[MAXM], fr[MAXM], to[MAXM], id[MAXN][3], n, W, ecnt;
void addedge(int u, int v) {
	edge[u].insert(v);
	mp[P(u, v)] = ++ecnt;
	fr[ecnt] = u, to[ecnt] = v;
}
vector<int> tr[MAXN];
int sz[MAXN], par[MAXN], num[MAXN], son[MAXN][2], ans1;
ll fac[MAXN], rev[MAXN], f[MAXN], g[MAXN], ans2 = 1;
ll C(int a, int b) { return fac[a] * rev[b] % MOD * rev[a - b] % MOD; }
ll modpow(ll a, int b) {
	ll res = 1;
	for (; b; b >>= 1) {
		if (b & 1) res = res * a % MOD;
		a = a * a % MOD;
	}
	return res;
}
void dfs(int u, int fa) {
	num[u] = max(id[u][0], max(id[u][1], id[u][2]));
	int t = 0;
	for (int v : tr[u]) {
		if (v == fa) continue;
		dfs(v, u);
		sz[u] += sz[v];
		son[u][t++] = v;
	}
	f[u] = C(sz[u], sz[son[u][0]]);
	ans2 = ans2 * f[u] % MOD;
	if (t == 2 && num[son[u][0]] > num[son[u][1]]) swap(son[u][0], son[u][1]);
	if (t == 1 && num[son[u][0]] == num[u]) swap(son[u][0], son[u][1]);
	if (num[u] == n) ++ans1;
	else ++sz[u];
}
void dfs2(int u, int fa) {
	for (int v : tr[u]) {
		if (v == fa) continue;
		g[v] = g[u] * sz[v] % MOD * modpow(sz[u], MOD - 2) % MOD;
		dfs2(v, u);
	}
}
void print(int a, ll b) {
	if (!W) printf("%d\n", a);
	else printf("%d %lld\n", a, b);
}
int main() {
	scanf("%d%d", &W, &n);
	assert(n > 3);
	for (int i = 1; i <= n - 3; i++) {
		int x, y; scanf("%d%d", &x, &y);
		addedge(x, y), addedge(y, x);
	}
	for (int i = 1; i <= n; i++) {
		int x = i, y = i % n + 1;
		addedge(x, y), addedge(y, x);
	}
	int tot = 0, rt, bes;
	for (int i = 1; i <= ecnt; i++) if (!vis[i]) {
		vis[i] = ++tot;
		edge[fr[i]].erase(to[i]);
		int lst = fr[i], cnt = 0;
		id[tot][cnt++] = fr[i];
		for (int u = to[i]; u != fr[i];) {
			if (cnt < 3) id[tot][cnt] = u;
			++cnt;
			int x = lst;
			if (lst == *edge[u].rbegin())
				lst = u, u = *edge[lst].begin();
			else lst = u, u = *edge[lst].upper_bound(x);
			edge[lst].erase(u), vis[mp[P(lst, u)]] = tot;
		}
		if (cnt >= 3) bes = tot;
	}
	for (int i = 1; i <= ecnt; i += 2) {
		if (vis[i + 1] == bes || vis[i] == bes) {
			if (fr[i] == n && to[i] == 1) rt = vis[i] == bes ? vis[i + 1] : vis[i];
			continue;
		}
		tr[vis[i]].push_back(vis[i + 1]);
		tr[vis[i + 1]].push_back(vis[i]);
	}
	for (int i = fac[0] = 1; i <= n; i++) fac[i] = fac[i - 1] * i % MOD;
	rev[n] = modpow(fac[n], MOD - 2);
	for (int i = n; i > 0; i--) rev[i - 1] = rev[i] * i % MOD;
	dfs(rt, 0);
	ans1 = tot - 1 - ans1, g[rt] = 1;
	dfs2(rt, 0);
	print(ans1, ans2);
	int m; scanf("%d", &m);
	while (m--) {
		int a, b; scanf("%d%d", &a, &b);
		int id1 = vis[mp[P(a, b)]], id2 = vis[mp[P(b, a)]];
		if (par[id2] != id1) swap(id1, id2);
		if (num[id1] < n) {
			print(ans1, ans2 * C(sz[id1] - 1, sz[son[id2][0]]) % MOD *
				C(sz[son[id2][1]] + sz[son[id1][1]], sz[son[id2][1]]) % MOD * modpow(f[id1] * f[id2] % MOD, MOD - 2) % MOD);
		} else print(ans1 - 1, ans2 * g[id2] % MOD);
	}
	return 0;
}
WAautomaton
关注
专栏目录
HNOI2019游记
Yu Gi Oh!
04-08 955
前言 本来计算着退役稳了根本就不想写什么游记啥的，但是既然老天并没有让我退役，那就还是写一些记录一下吧。 Day0 省选前一天随便打了一些板子和复习然后就过了。 Day1 开场看题。 T1计算几何?? 感觉有点不好。 T2字符串题，50分可以直接写暴力。 T3多边形?这又是个什么鬼题。 感觉今年省选画风不对啊， 稍微分析一下 T3 最后的终态显然是所有边都连在了 n 号点上，那么和 n 相连的边把...
题解-[HNOI2004]宠物收养场
KonnyWen的博客
02-28 222
题解-[HNOI2004]宠物收养场 Introduction\texttt{Introduction}Introduction 蒟蒻初学平衡树，见此题无妙解，乃记之。 Description\texttt{Description}Description [HNOI2004]宠物收养场 有 nnn 个宠物或主人会来，每个来者都有性格 aaa。如果当前主人有剩余，宠物 xxx 就会领走一个和它性...
[HNOI2019]JOJO
weixin_34032792的博客
05-13 163
不错的好题 有一个比较弱化的问题： [题集]串 已知小串的子串问题，例题2 本题： 首先单纯快速跳kmp就要用到AC自动机的fail了（单串AC自动机） 所谓回撤，就是一个可持久化， 建出时间树，直接dfs，回来的时候撤销这一路的操作。 问题就是往S后插入x个c，考虑新加入的fail之和 我们把(x,c)尽量看成一个整体 （请感性理解这句话23333） 不妨称一...
【HNOI 2019】校园旅行
洛水·锦依卫's Blog
04-09 340
【HNOI 2019】校园旅行 Problem Description 某学校的每个建筑都有一个独特的编号。一天你在校园里无聊，决定在校园内随意地漫步。 你已经在校园里呆过一段时间，对校园内每个建筑的编号非常熟悉，于是你情不自禁的把周围每个建筑的编号都记了下来——但其实你没有真的记下来，而是把每个建筑的编号除以 \(2\) 取余数得到 \(0\) ...
HNOI2019
weixin_30682415的博客
04-09 152
省选总结 day0 写了下平时不经常写的模板，像什么\(LCT\),圆方树，\(exlucas\)之类的，但是一个都没考。 day1 提前十几分钟进了考场，可以提前动电脑，赶紧把\(vimrc\)打了。八点钟解压之后先看了下三道题，三道题题面都很长，第一道题\(fish\)看了半天感觉是个计算几何，先放到一边。 \(T2\)也看了半天，发现就是要对每个前缀求最长的\(border\)。感觉还挺有想...
【hnoi2007】
cjoilmd的专栏
03-10 3131
盾哥说hnoi年份靠前的题目价值不大，所以考虑07年选做一些题后转战poi吧。       hnoi2007简直是奇葩的一年，day1明显难于day2（day1就算算上的第四题（那个时候还没有cdq的论文），也只有两道可做题，而day2有三道可做题）              day1       海盗分宝：       。。。。。。       求神牛解释题意。。。。。。
算法竞赛进阶指南题解目录
m0_73035684的博客
10-11 774
题解目录
2019.4.summary
LMB_001的博客
07-04 1064
2019.4.1 BZOJ1061: [Noi2008]志愿者招募 真心有点难QAQ https://www.byvoid.com/zhs/blog/noi-2008-employee 看void爷的博客看了三遍才懂。这种网络流建模看不出来啊 首先可以列出一些数相加减的式子 然后可以考虑差分，得到一些式子和为0的情况 这不就很像每一个点流量平衡对吧 正的是流进来的，负的是流出去的 常数是关于s，t...
cf757G. Can Bash Save the Day?
miaom的博客
03-02 968
简述题意：给一棵n个节点的数和一个n的数组（n的排列），支持两个操作：1 询问a[l]~a[r]的所有点到某个点的距离，2 交换a[x]与a[x+1]。 一开始想到了这道题bzoj4012: [HNOI2015]开店，大概看了之后发现是开店的加强版。我想到像开店一样使用主席树处理，显然操作2很容易。然而 他的时间是一个log。大概看了题解，感觉套路~   毕竟开店的题解大多是log^2。。
HNOI 2019 简要题解
wxh010910
04-15 1150
fish 枚举 ADADAD ， BCBCBC 的限制相当于 ADADAD 是 BCBCBC 中垂线且 BCBCBC 中点在 ADADAD 上，预处理所有中垂线，排序之后二分即可；EFEFEF 的限制可以双指针扫一下。 #include <bits/stdc++.h> using namespace std; struct point { int x, y; point(...
「HNOI2019」校园旅行
qq_29361257的博客
11-10 227
题目 某学校的每个建筑都有一个独特的编号。一天你在校园里无聊，决定在校园内随意地漫步。 你已经在校园里呆过一段时间，对校园内每个建筑的编号非常熟悉，于是你情不自禁的把周围每个建筑的编号都记了下来——但其实你没有真的记下来，而是把每个建筑的编号除以 222 取余数得到 000 或 111，作为该建筑的标记，多个建筑物的标记连在一起形成一个 010101 串。 你对这个串很感兴趣，尤其是对于这个串是回文串的情况，于是你决定研究这个问题。 学校可以看成一张图，建筑是图中的顶点，而某些顶点之间存在无向边。对于每个顶
HNOI2019无盘可翻记及简要题解
Rayment的博客
04-14 740
Day 0 感觉高一和高二考试时都抱有完全不同的心态，高一的时候只是想着去玩，甚至都没有好好复习。有时候想如果自己还是高一就好了，就还有机会，也不至于没有退路。 我已经做好了面对各种结果的心理准备了，自己应该能接受的 祝自己以及各位好运吧 Day 1 到考点发现大家都已经进去了，赶紧跟上去 开考前发现系统时间是错的，真坑。。调了下时间，然后把配置打了 怎么只有一个PDF？？？我拿错包了？ 今年连大...
HNOI2019 游记
weixin_30847939的博客
04-08 121
链接：HNOI2019 游记 转载于:https://www.cnblogs.com/water-mi/p/10668617.html
HNOI2019 头铁记爆肝记
litble的成(tui)长(fei)史
04-14 1244
Day0 省选前的模拟考试连环爆炸，说是爆炸，其实也不算，因为出来的分基本上都和预估分差不多，至于为啥考那么点分，原因就一条——我菜。我还真没犯多少傻逼错误，真的就是单纯地做不出来而已。 结果弄得心态挺爆炸的。 不过一想，我去年也打过Ag了，至于集训队真是想都不敢想的事情，于是我就只有两种结果： 若是进了省队，那很大可能不会有新收获，要多浪费学文化课的时间，然后还要大夏天地清早起床坐公交车跑到雅礼...
【LOJ3055】「HNOI2019」JOJO
cz_xuyixuan的博客
04-23 843
【题目链接】 点击打开链接 【思路要点】 考虑将字符串缩成若干个相同字符的段 (cnt,char)(cnt,char)(cnt,char) ，相邻的字符不相同。 考虑一个与某一个前缀匹配的后缀，若该后缀横跨了 x (x≥2)x\ (x\geq2)x (x≥2) 段，那么其对应的第 222 段至第 x−1x-1x−1 段的 (cnt,char)(cnt,char)(c...
HNOI2019选讲
imalyd的博客
05-09 264
D1T1 鱼 传送门 首先O(N6)O(N^6)O(N6)的做法是显然的：直接枚举即可。但是显然枚举量太大了。 这类题有一种思路，就是考虑哪些点是比较关键的，固定比较关键的点。 这题中，显然D和A是最关键的。因此固定D，把其他点绕D极角排序，并按顺序枚举A。 B、C的方案数： 显然BC⊥\perp⊥AD，且BC中点在AD上。 因此事先把每条线段的方向、中点hash并存下来，查找时lower_bo...
[HNOI2019]白兔之舞
wyy603的博客
07-08 239
题解 我们枚举 yyy 坐标序列，并使用组合数计算 xxx 坐标序列的个数，所以我们要对 ∀0≤t<k\forall 0\le t<k∀0≤t<k 将 ∑m mod k=t∑a0=x,a1,a2,⋯ ,am=y(Lm)∏i=1m−1w(ai,ai+1) \sum_{m\bmod k=t} \sum_{a_0=x,a_1,a_2,\cdots,a_m=y} \binom{L}{m} \prod_{i=1}^{m-1} w(a_i,a_{i+1}) mmodk=t∑​a0​=x,a1​,a2​
P3722 [AH2017/HNOI2017] 影魔
最新发布
07-28
根据引用\[1\]和引用\[2\]的描述，题目中的影魔拥有n个灵魂，每个灵魂有一个战斗力ki。对于任意一对灵魂对i,j (i)，如果不存在ks (i)大于ki或者kj，则会为影魔提供p1的攻击力。另一种情况是，如果存在一个位置k，...
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值