ShanDong Multi-University Training #2

B,J签到
D：

题意： 给你一堆点的坐标，给你一个半圆的圆心和半径，半圆可以旋转，问你半圆区域内最多容纳几个点。

思路：
看数据的数量级，On²的算法完全可以过，所以直接枚举所有点，让半圆的直径边过当前点，再遍历剩下的n-1个点，剩下的点分为在直径边左边和直径边右边两种情况，统计一下两种情况的点离圆心的距离小于半径的个数，然后取一个最大值就可。思路比较简单，代码稍微难写一点。（直接贴wushi代码了）

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstdio>
using namespace std;
int x, y;
double r;
//nt check
const int N = 1e3 + 10;
//pair<int,int>dian[N];
double ge(int x1, int y1) {
	double d1 = (x1 - x) * (x1 - x), d2 = (y1 - y) * (y1 - y);
	return sqrt(d1 + d2);
}
struct node {
	int x, y;
	double l;
};
node dian[N];
int cmp(node x, node y) {
	return x.l < y.l;
}
int n;
int check(int x1, int y1) {
	double xie = 0;
	int s = 0, xia = 0;
	if (x1 - x != 0) {
		xie = double(y1 - y) / (x1 - x);
		for (int i = 0; i < n; i++) {
			if(dian[i].l>r)break;
			double d = dian[i].x - x;
			double dy = d * xie;
			double y2 = y + dy;
			if (dian[i].y >= y2) {
				s++;
			}
			if (dian[i].y <= y2) {
				xia++;
			}
		}
	} else {
		xie = 1;
		for(int i=0;i<n;i++){
			if(dian[i].l>r)break;
			if(dian[i].x>=x)s++;
			if(dian[i].x<=x)xia++;
		}
	}
	return max(s, xia);
}
int main() {
	while (cin >> x >> y >> r) {
		if(r<=0)break;
		cin >> n;
		for (int i = 0; i < n; i++) {
			scanf("%d %d", &dian[i].x, &dian[i].y);
			dian[i].l = ge(dian[i].x, dian[i].y);
		}
//		for(int i=0;i<n;i++){
//			cout<<dian[i].l<<" "<<dian[i].x<<" "<<dian[i].y<<endl;
//		}
		sort(dian, dian + n, cmp);
		int ans = 0;
		for (int i = 0; i < n; i++) {
			ans = max(check(dian[i].x, dian[i].y), ans);
		}
		cout << ans << endl;
	}
	return 0;
}

G：

题意：
求含有 k 个不同因子的最小数。1 ≤ k ≤ 103。
思路：
1.任何一个大于1的自然数N，如果不是质数，那么可以唯一分解成有限个质数的乘积。

题意要求因子个数等于k的最小n，发现以下性质：
可以猜测 pi ≤ 37，由 3 个引理可知，只需要 dfs 搜索 ci 即可。(当然是高精或者py)

K：

题意：
给你两组1-n的排列，可间断的在两个排列里选取k个组成新的排列a，b，使得b中的每个元素都是a中对应位置的倍数。
atc官方题解：
将a，b中所有有倍数关系的数子的下标用pair存起来。然后对第一维按升序排列，第一维相同时，第二维按降序排列。此时，问题已经转化为，求pair的第二维中最长的严格上升子序列。
由于数据范围较大，只能用贪心+二分的方式求最长上升子序列 贪心+二分的方式求最长上升子序列（onlogn）：

const int N = 2e5 + 7, M = 3e6 + 7;
int a[N], b[N];
pii s[M];

bool cmp(pii a, pii b)
{
	if (a.ft == b.ft) return a.sd > b.sd;
	return a.ft < b.ft;
}

int main()
{
	IOS;

	int n; cin >> n;
	for (int i = 1, x; i <= n; i++)
		cin >> x, a[x] = i;
	for (int i = 1, x; i <= n; i++)
		cin >> x, b[x] = i;

	int k = 0;
	for (int i = 1; i <= n; i++)
		for (int j = i; j <= n; j += i)
			s[k++] = { a[i], b[j] };

	sort(s, s + k, cmp);

	int ans = 0;
	vector<int>v;
	for (int i = 0; i < k; i++)
	{
		if (v.empty() || s[i].sd > v.back())
			v.push_back(s[i].sd);
		else
		{
			int idx = lower_bound(all(v), s[i].sd) - v.begin();
			v[idx] = s[i].sd;
		}
	}

	cout << sz(v) << endl;

	return 0;
	
}