FLOYD算法模板练习题，+FLOODFILL，没什么可说

#include <cstdio>
#include <iostream>
#include <vector>
#include <algorithm>
#include <math.h>
#include <cstring>
using namespace std;

const int MAXN = 165;
const double INF = 1e20;
int n, m, id = 0, field[MAXN], vis[MAXN];
double tot_ans, dis[MAXN][MAXN], md[MAXN], maxdis[MAXN];
struct node{
	double x, y;
}a[MAXN];

double calc(node x, node y){
	double disx = abs(x.x - y.x);
	double disy = abs(x.y - y.y);
	return sqrt(disx*disx + disy*disy);
}

void dfs(int u, int answer){
	vis[u] = 1; field[u] = answer;
	for (int i = 1; i <= n; i++)
		if (!vis[i] && dis[i][u] != INF)
			dfs(i, answer);
}

int main(){
	freopen("cowtour.in", "r", stdin);
	freopen("cowtour.out", "w", stdout);
	cin >> n;
	for (int i = 1; i <= n; i++)
		cin >> a[i].x >> a[i].y;
	memset(dis, 0x3f, sizeof(dis));
	for (int i = 1; i <= n; i++)
		for (int j = 1; j <= n; j++){
			char c; cin >> c;
			if (c == '1' || i == j)
				dis[i][j] = calc(a[i], a[j]);
			else dis[i][j] = INF;
		}
	
	for (int i = 1; i <= n; i++)
		if (!vis[i]) dfs(i, ++id);
	
	for (int k = 1; k <= n; k++)
		for (int i = 1; i <= n; i++)
			for (int j = 1; j <= n; j++)
				if (dis[i][j] > dis[i][k] + dis[k][j])
					dis[i][j] = dis[i][k] + dis[k][j];
	
	for (int i = 1; i <= n; i++){
		for (int j = 1; j <= n; j++)
			if (dis[i][j] < INF)
				md[i] = max(md[i], dis[i][j]);
		maxdis[field[i]] = max(maxdis[field[i]], md[i]);
	}
	
	tot_ans = INF;
	for (int i = 1; i <= n; i++)
		for (int j = i+1; j <= n; j++)
			if(field[i] != field[j]){
				double maxd = max(max(maxdis[field[i]],maxdis[field[j]]),md[i]+md[j]+calc(a[i],a[j]));
				tot_ans = min(tot_ans, maxd);
			}
	
	printf("%.6f\n", tot_ans);
	
	return 0;
}

首先直接按题意模拟一下，发现“所有质数都要被选上”这个条件很烦，因为选上一个卡牌后有很多质数会受到影响，非常不好做。换言之，题中的限制很强。

考虑正难则反，钦定一些质数使得它们 恰好 没有被选上，其它质数都被选上。然后会发现 恰好 这个条件还是很强，不好搞，直接上容斥，把 恰好

弱化成 一部分 被选上，另一部分不用管。这样限制条件就松了很多。

回到题目，$n\le 10^6$ 这个条件是唬人的，真正互不相同的数字只有 $2000$ 个，拿桶记一下数 $n$ 就可以扔了。

对于 $s_i\le 30$ 的部分分，显然可以状压预处理然后容斥计算。这为我们提供了一些思路。

[NOI2015] 寿司晚宴 这道题给出了一个相当经典的处理手段：一个数至多有一个 $\gt \sqrt n$ 的质因子，对于 $\le \sqrt n$ 的质因子状压，对于 $\gt \sqrt n$ 的质因子单独考虑。为了方便，称这类质因子为“大质因子”。

本题 $V=2\times 10^3$，发现 $\le \sqrt V$ 的质数只有 $2\sim 43$ 这 14 个，进一步地，$43\times 47\gt 2000$，所以 43 也不用参与状压，进一步压缩了状态。

然后统计，令 $f(i,j)$ 表示前 13 个质数没有被选的状态为 $i$，当前大质因子为 $j$ 时的方案数。预处理是简单的。

然后我们考虑容斥，容易发现容斥系数为 $(-1)^{|S|}$。

计算答案并不难。对于一个状态 $i$，枚举 $j\in [43,2000],j\in \mathbb P$，如果 $j$ 在询问集合中，那么有 $2^{f(i,j)}-1$ 的贡献（筛去空集），反之有 $2^{f(i,j)}$ 的贡献。

把询问集合存成 std::vector，每次就不用枚举 $j$，直接扫一遍 std::vector 然后统计答案，和容斥系数乘起来加到答案里即可。实现可以参考代码。

时间复杂度 $\mathcal O(\sum c_i\times 2^{13})$。

合格的签到题。

首先发现式子是一个 0-1 分数规划的模型，所以先二分答案 $mid$，男生 $i$ 和女生 $j$ 配对的价值为 $a_{i,j}-b_{i,j}\times mid$。

进一步地，发现男生和女生的关系形成一张二分图，并且男女双双配对，配对有权值，相当于二分图带权最大完美匹配。

可以使用 KM 算法求解，但我不会，学了一个费用流打了一下。核心思想就是把 EK 算法的 BFS 找增广路改为 SPFA 找增广路，顺便求一下最长路。

时间复杂度 $\mathcal O(n^4)$，在 COGS 的机器上可能有点难跑，但这只是理论上界，而且这可是二分图，SPFA 和 EK 都很难卡满，足以通过本题。

easy tea，不过有点小作弊：想完以后看了 Asm.Def 学长的评论才确定自己的思路是对的，并不是很自信。

考虑转化题意，求最多的说真话的人。

注意到如果有 $a_i$ 个人比 $i$ 低，有 $b_i$ 个人比 $i$ 高，说明 $i$ 的成绩排名 $\in [a_i+1,n-b_i]$。

将区间视作一条线段，然后把 $[l,r]$ 相同的线段合并，给 $[l,r]$ 这条线段赋权。本题就转化为了最大线段带权覆盖问题。$\rm dp + std::map$ 可以 $\mathcal O(n\log n)$ 地解决这个问题。

小陷阱：如果有超过 $r-l+1$ 个人的成绩排名 $\in [l,r]$，那么 $[l,r]$ 的权值为 $r-l+1$，$\rm dp$ 的时候取 $\rm min$ 判断一下就好了。

首先，01 背包是 NP 的，直接做不行，要关注题中的小性质。

发现 $C_i\le 300$，于是对 $C_i$ 分组，同组内价值降序排序，设为 $V_1\sim V_m$，那么如果选上 $V_k$，则必然选上 $V_1\sim V_{k-1}$，挺显然的。把它看作一个前缀和的形式。

然后发现 dp 方程可以各个同余类分开计算。

进一步的，一个同余类中的 dp 转移具有决策单调性。

感性理解一下，如果决策不具有单调性，那么会取的 $k$ 会更大，而我们对元素进行了降序排序，$k$ 更大是不优的。

决策单调性分治计算即可。

时间复杂度 $\mathcal O(kc\log k)$。

Editorial of srf.

不知道状压可不可行阿，先说一个我的错误状压思路：

设 $f(i,S)$ 表示前 $i$ 个点的叶子节点状态为 $S$，最大深度在点 $i$ 处取得的方案数。然后会发现这个东西没办法转移。

继续瞎搞，考虑枚举深度 $d$，划分为若干层，进行状压 dp，但是仍然没啥好方法转移。也可能是我菜。

树上深度计数考虑组合数。计算总数，最后除以 $(n-1)!$。

设 $f(i,j)$ 表示 $i$ 个点的树，深度为 $j$ 的方案数。

因为点 $2$ 一定接在 $1$ 上，所以我们枚举点 $2$ 的子树大小及深度，然后合并其它子树和点 $2$ 的子树。这样显然不会算重。

这个过程可以保证方案是合法的，只关心点编号之间的大小关系。用组合数辅助计算。

取整可以用 long double 或者 __int128 存一下。

时间复杂度 $\mathcal O(n^4)$。