cf1371E2.Asterism，数学，同余

题目链接

题意：

给一个数组$a$，进行如下操作：
初始有一个值$x$，对于$a$的任意一个排列$P$，依次比较$a_1,a_2,...a_n$，如果$x>a_i，则x=x+1$，要求$x$大于所有的$a_i$。
令$f(x)$为满足上述条件的$a$的排列的个数。另外再给一个质数$p,p\le n$，要求$f(x)\%p\ne 0$，找到所有满足条件的$x$。

题解：

首先发现$x$的值至少为$max(a_i)-n+1$，最大值为$max(a_i)-1$，令$m=max(a_i)$，则$m-n+1\le x\le m-1$。
因为当$x\ge m$后，$f(x)=n!$，肯定被$q$整除。
所以对$x$的值进行一一遍历，然后求解即可。
那么问题就变成了对于$x$如何求$f(x)$。
首先对$a$数组排序，那么对于给定的$x$，我们首先找到最大的$i，a_i\le x$，那么对于排列$P$的第一个位置我们可以选择$P_0=a_{0}or{a}_{1}or...or{a}_i$，也就是有$i+1$种方案选择$P_0$位置上的数，对于$P_1$的选择，也是一样的，找到最大的$j,a_j\le x+1$，因为选择$P_0$已经用掉了一个数，所以选择$P_1$有$j-1$中方案，那么表达成公式就是$f(x)={\prod }_{i=x}^{x+n-1}{b}_i-(i-x)$，其中$b_i$是$a$数组中小于等于$i$的数的个数。判断$f(x)$是否满足条件只需要检查$f(x)$中每一个乘数是否包含素因子$p$就可以了。
那么数据量较小的时候，直接暴力求解即可：

#include <bits/stdc++.h>
using namespace std;

const int maxn = 2e3 + 10;

int main() {
	int n, p;
	cin >> n >> p;
	vector<int> a(n);
	for (int i = 0; i < n; i++) {
		cin >> a[i];
	}
	sort(a.begin(), a.end());
	int Max = a[n-1];
	set<int> res;
	
	for (int i = Max-n+1; i < Max; i++) {
		int x = i, v = 0, j = 0, k = 0;
		bool flag = false;
		while (k < n) {
			while (j < n && a[j] <= x) {
				v++; j++;
			}
			if (v % p == 0) {
				flag = true;
				break;
			}
			x++;
			k++;
			v = v-1;
		}
		if (!flag) {
			res.insert(i);
		}
	}
	
	cout << res.size() << endl;
	for (int x : res) {
		cout << x << " ";
	}
	cout << endl;
}

对于数据量较大的情况，需要做一点变形：
$f(x)={\prod }_{i=x}^{x+n-1}{b}_i-(i-x)={\prod }_{i=x}^{x+n-1}x-(i-b_i)$
如果乘数$x-(i-b_i)\equiv 0(modp)$，那么即$x\equiv i-b_i(modp)$.
那么问题就转换成了预处理出所以的$i-b_i$，然后依次枚举$x$，判断是否有因数与$x\%p$同余。
由于$a_i$太大，所以首先对每个$a_i=a_i-(m-n)$，这样每一个$a_i\le n$，可以顺序统计出$b_i$。
此时$x$的取值范围为$1\le x<n-1$，那么$b_i$中的$i$最大可能到$2\ast n$。
把所有的$b_i$都统计出来之后，再用一个数组统计一下余数，然后顺序遍历就可以了。

#include <bits/stdc++.h>
using namespace std;

int mod(int x, int p) {
	return (x%p + p) % p;
}

int main() {
	int n, p;
	cin >> n >> p;
	vector<int> a(n), b(2*n);
	for (int i = 0; i < n; i++) {
		cin >> a[i];
	}
	int Max = *max_element(a.begin(), a.end());
	
	for (int i = 0; i < n; i++) {
		b[max(0, a[i]-(Max-n))]++;
	}
	for (int i = 1; i < 2*n; i++) {
		b[i] += b[i-1];
	}
	
	vector<int> f(n);
	int shift = Max-n;
	for (int i = 1; i <= n; i++) {
		// x = Max-n+1，第一个x的值
		// 统计每一个乘数的余数
		// i在之前统一被减小了(Max-n)，这里要加回来
		f[mod(i+shift-b[i], p)]++;	// (i-b[i])
	}
	
	vector<int> res;
	for (int i = 1; i < n; i++) {
		// 遍历每一个可能的x
		// x = Max-n+i
		if (f[mod(i+shift, p)] == 0) {	// x
			// 当前没有与x同余的乘数
			res.push_back(i+shift);
		}
		// 去掉之前的一个因数
		f[mod(i+shift-b[i], p)]--;
		// 加上下一个因数
		f[mod(n+i+shift-b[n+i], p)]++;
	}
	
	cout << res.size() << endl;
	for (int x : res) {
		cout << x << " ";
	}
	cout << endl;
}