2019年第十届蓝桥杯A组省赛

A.平方和（5分）
答案：2658417853

试题 A: 平方和
本题总分：5 分
【问题描述】
小明对数位中含有 2、0、1、9 的数字很感兴趣，在 1 到 40 中这样的数包括 1、2、9、10 至 32、39 和 40，共 28 个，他们的和是 574，平方和是 14362。注意，平方和是指将每个数分别平方后求和。
请问，在 1 到 2019 中，所有这样的数的平方和是多少？
【答案提交】
这是一道结果填空的题，你只需要算出结果后提交即可。本题的结果为一 个整数，在提交答案时只填写这个整数，填写多余的内容将无法得分。
提示：如果你编写程序计算，发现结果是负的，请仔细检查自己的程序， 不要怀疑考场的编程软件。

思路：
简单枚举题。 从1到2019分别遍历每一个正整数，判断该正整数的每一数位是否包含”2019“中的4个数字之一，若包含，则结果ans加上它的平方和。

#include<bits/stdc++.h>
using namespace std;
#define INF 0x7fffffff
typedef long long ll;
const int maxn=50005;
bool check(int x){
	while(x){
		if(x%10==2||x%10==0||x%10==1||x%10==9) return true;
		x/=10;
	}
	return false;
}
int  main(){
	ll ans=0;
	for(int i=1;i<=2019;i++){
		if(check(i)) ans+=i*i;
	}
	cout<<ans<<endl;
	return 0;
} 

---

B.数列求值（5分）
答案：4659

试题 B: 数列求值
本题总分：5 分
【问题描述】
给定数列 1, 1, 1, 3, 5, 9, 17, …，从第 4 项开始，每项都是前 3 项的和。求第 20190324 项的最后 4 位数字。
【答案提交】
这是一道结果填空的题，你只需要算出结果后提交即可。本题的结果为一个 4 位整数（提示：答案的千位不为 0），在提交答案时只填写这个整数，填写多余的内容将无法得分。

思路：
简单模拟题。 从第4项开始，它的数值取决于它的前3项值的和，所以我们可以用a,b,c来存储第k项的前3项数值，每次求出第k项数值以后再更新一下a,b,c。最后k为20190324的时候，c存储的值即为结果。
注意：由于我们只需要输出结果的后四位，我们可以让每一项数值对10000取模，即最多只留4位整数，否则结果将会溢出。

#include<iostream>
#include<queue>
#include<vector>
using namespace std;
#define INF 0x7fffffff
typedef long long ll;
const int maxn = 50005;
const int mod = 10000;
int  main() {
	int a = 1, b = 1, c = 1,k=3;
	while (++k) {
		ll d = (a+b+c)%mod;
		a = b; b = c; c = d;
		if (k == 20190324) break;
	}
	cout << c << endl;
	return 0;
}

---

C.最大降雨量（10分）
答案：34

试题 C: 最大降雨量
本题总分：10 分
【问题描述】
由于沙之国长年干旱，法师小明准备施展自己的一个神秘法术来求雨。这个法术需要用到他手中的 49 张法术符，上面分别写着 1 至 49 这 49 个数字。法术一共持续 7 周，每天小明都要使用一张法术符，法术符不能重复使用。每周，小明施展法术产生的能量为这周 7 张法术符上数字的中位数。法术施展完 7 周后，求雨将获得成功，降雨量为 7 周能量的中位数。由于干旱太久，小明希望这次求雨的降雨量尽可能大，请大最大值是多少？
【答案提交】
这是一道结果填空的题，你只需要算出结果后提交即可。本题的结果为一 个整数，在提交答案时只填写这个整数，填写多余的内容将无法得分。

思路：
思维+手模。我们可以7行7列的方格图如下所示：(图上第i行上的 ‘o’ 代表第i周的每一天的法术符数值，’#‘代表第i周的中位数，’$‘代表7周能量的中位数，也就是最后结果)
o o o # o o o
o o o # o o o
o o o # o o o
o o o $ o o o
o o o # o o o
o o o # o o o
o o o # o o o
假设’#‘左边的元素数值都小于’#‘的数值，’#‘右边的元素数值都大于’#‘的数值，’$ ‘上边的元素数值都小于’ $ ‘的数值，’$ ‘下边的元素数值都大于’ $ ‘的数值。
为了使$的数值最大，我们可以计算一下一定会比 $ 的数值大的天数，显然就是黄色元素的个数，即3+4*3=15天，所以$ 的最大数值是49-15=34 。

---

D.迷宫（10分）
答案：DDDDRRURRRRRRDRRRRDDDLDDRDDDDDDDDDDDDRDDRRRURRUURRDDDDRDRRRRRRDRRURRDDDRRRRUURUUUUUUULULLUUUURRRRUULLLUUUULLUUULUURRURRURURRRDDRRRRRDDRRDDLLLDDRRDDRDDLDDDLLDDLLLDLDDDLDDRRRRRRRRRDDDDDDRR

试题 D: 迷宫
本题总分：10 分
【问题描述】
下图给出了一个迷宫的平面图，其中标记为 1 的为障碍，标记为 0 的为可以通行的地方。
010000
000100
001001
110000
迷宫的入口为左上角，出口为右下角，在迷宫中，只能从一个位置走到这 个它的上、下、左、右四个方向之一。
对于上面的迷宫，从入口开始，可以按DRRURRDDDR 的顺序通过迷宫， 一共 10 步。其中 D、U、L、R 分别表示向下、向上、向左、向右走。
对于下面这个更复杂的迷宫（30 行 50 列），请找出一种通过迷宫的方式，其使用的步数最少，在步数最少的前提下，请找出字典序最小的一个作为答案。请注意在字典序中D<L<R<U。（如果你把以下文字复制到文本文件中，请务 必检查复制的内容是否与文档中的一致。在试题目录下有一个文件 maze.txt， 内容与下面的文本相同）
01010101001011001001010110010110100100001000101010
00001000100000101010010000100000001001100110100101
01111011010010001000001101001011100011000000010000
01000000001010100011010000101000001010101011001011
00011111000000101000010010100010100000101100000000
11001000110101000010101100011010011010101011110111
00011011010101001001001010000001000101001110000000
10100000101000100110101010111110011000010000111010
00111000001010100001100010000001000101001100001001
11000110100001110010001001010101010101010001101000
00010000100100000101001010101110100010101010000101
11100100101001001000010000010101010100100100010100
00000010000000101011001111010001100000101010100011
10101010011100001000011000010110011110110100001000
10101010100001101010100101000010100000111011101001
10000000101100010000101100101101001011100000000100
10101001000000010100100001000100000100011110101001
00101001010101101001010100011010101101110000110101
11001010000100001100000010100101000001000111000010
00001000110000110101101000000100101001001000011101
10100101000101000000001110110010110101101010100001
00101000010000110101010000100010001001000100010101
10100001000110010001000010101001010101011111010010
00000100101000000110010100101001000001000000000010
11010000001001110111001001000011101001011011101000
00000110100010001000100000001000011101000000110011
10101000101000100010001111100010101001010000001000
10000010100101001010110000000100101010001011101000
00111100001000010000000110111000000001000000001011
10000001100111010111010001000110111010101101111000
【答案提交】
这是一道结果填空的题，你只需要算出结果后提交即可。本题的结果为一 个字符串，包含四种字母 D、U、L、R，在提交答案时只填写这个字符串，填写多余的内容将无法得分。

思路：
深搜+剪枝策略。 一看是迷宫类型题，而且和经典迷宫题基本一致，并且行数和列数不算很大，所以显然可以用dfs来解决。dfs的代码基本上就是迷宫题的模板代码，只是由于回溯次数太多了会超时，所以本题我通过记忆化搜索进行剪枝便可以大幅度地提高搜索效率，得出结果。记忆化搜索的方法是用数组d[i][j]记录从起点（1，1）到点（i，j）的最少步数，若在回溯过程发现当前步数大于d[i][j]则直接返回，避免无效访问。

#include<iostream>
#include<queue>
#include<vector>
using namespace std;
#define INF 0x7fffffff
typedef long long ll;
const int maxn = 100;
const int mod = 10000;
char str[maxn][maxn], ds[5] = {'D','L','R','U'},s[2000],ans[2000];
bool vis[maxn][maxn];
int dx[5] = { 1,0,0,-1 },d[maxn][maxn];
int dy[5] = { 0,-1,1,0 };
void dfs(int x, int y,int step) {
	if (x == 30 && y == 50) {
		if (step >= d[x][y]) return;
		d[x][y] = step;
		for (int i = 1; i <step; i++) ans[i]=s[i];
	}
	if (step > d[x][y]) return;
	for (int i = 0; i <= 3; i++) {
		int tx = x + dx[i], ty = y + dy[i];
		if (tx >= 1 && tx <= 30 && ty >= 1 && ty <= 50 && !vis[tx][ty]) {
			vis[tx][ty] = true; s[step] = ds[i]; d[x][y] = step ;
			dfs(tx, ty,step+1);  
			vis[tx][ty] = false;
		}
	}
}
int  main() {
	for (int i = 1; i <= 30; i++)
		for (int j = 1; j <= 50; j++) {
			d[i][j] = INF;
			cin >> str[i][j];
			if (str[i][j] == '0') vis[i][j] = false;
			else vis[i][j] = true;
		}
	dfs(1, 1,1);
	for (int i = 1; i < d[30][50]; i++) cout << ans[i];
	return 0;
}

---

E.RSA解密（15分）
答案：579706994112328949

试题 E: RSA 解密
本题总分：15 分
【问题描述】
RSA 是一种经典的加密算法。它的基本加密过程如下。
首先生成两个质数 p, q，令 n = p · q，设 d 与 (p − 1) · (q − 1) 互质，则可找到 e 使得 d · e 除 (p − 1) · (q − 1) 的余数为 1。
n, d, e 组成了私钥，n, d 组成了公钥。
当使用公钥加密一个整数 X 时（小于 n），计算 C = X^d mod n，则 C 是加密后的密文。
当收到密文 C 时，可使用私钥解开，计算公式为 X = C^e mod n。例如，当 p = 5, q = 11, d = 3 时，n = 55, e = 27。
若加密数字 24，得 24^3 mod 55 = 19。解密数字 19，得 19^27 mod 55 = 24。
现在你知道公钥中 n = 1001733993063167141, d = 212353，同时你截获了别人发送的密文 C = 20190324，请问，原文是多少？
【答案提交】
这是一道结果填空的题，你只需要算出结果后提交即可。本题的结果为一 个整数，在提交答案时只填写这个整数，填写多余的内容将无法得分。

思路：
熟练掌握数论中的几大常用算法。 扩展欧几里得算法+快速幂+快速乘。 由于n是由两个质因数p和q的乘积得来，显然n只有两个因子即p和q;我们可以通过遍历sqrt(n)个元素求出q和p,因为求出了一个因子p,那么另一个因子q=n/p。 求出了p和q以后我们便可以得到(p-1)*(q-1)的结果k。又因为题目已知d与k互质，且(d * e)%k== 1 。显然就是告诉我们e就是方程d*e==1(modk)的逆元。 套用扩展欧几里得算法模板便可以求出逆元e，但需要注意的是，e不能为负数，所以若e是负数的话，则需要对e进行处理。 此时数e,n,C的值我们都已知了，所以直接带入公式即可求出结果X。由于本题涉及的数值太大，求（C ^e）%n的时候肯定会溢出，所以需要用的到快速幂模板和快速乘模板。

#include<iostream>
#include<queue>
#include<vector>
using namespace std;
#define INF 0x7fffffff
typedef long long ll;
typedef unsigned long long ull;
ll n = 1001733993063167141,d= 212353,C= 20190324,p,q,e,ans;
void get_p_q() {
	for (ll i = 2; i * i <= n; i++) {
		if (n % i == 0) {
			p = i; q = n / i;
			break;
		}
	}
}
ll exgcd(ll a, ll b, ll& x, ll& y) {
	if (!b) { x = 1, y = 0; return a; }
	else {
		ll ans = exgcd(b, a % b, x, y);
		ll temp = y; y = x - (a / b) * y;
		x = temp;
		return ans;
	}
}
ll inv(ll a, ll b) {
	ll x, y;
	exgcd(a, b, x, y);
	if (x < 0) x = (x % b + b) % b;
	return x;
}
ll qmul(ll a, ll b, ll mod) {
	ll res = 0;
	while (b) {
		if (b & 1) res = (res + a) % mod;
		b >>= 1;
		a = (a + a) % mod;
	}
	return res;
}
ll qpow(ll base, ll power, ll mod) {
	ll res = 1;
	while (power) {
		if (power & 1) res = qmul(res, base, mod);
		power >>= 1;
		base = qmul(base,base,mod);
	}
	return res;
}
int  main() {
	get_p_q();
	e = inv(d, (p - 1) * (q - 1));
	ans = qpow(C, e, n);
	cout << ans << endl;
	return 0;
}

---

F.完全二叉树的权值(15分）

试题 F: 完全二叉树的权值
时间限制: 1.0s 内存限制: 256.0MB 本题总分：15 分
【问题描述】
给定一棵包含 N 个节点的完全二叉树，树上每个节点都有一个权值，按从上到下、从左到右的顺序依次是 A1, A2, · · · AN，如下图所示：
现在小明要把相同深度的节点的权值加在一起，他想知道哪个深度的节点权值之和最大？如果有多个深度的权值和同为最大，请你输出其中最小的深度。
注：根的深度是 1。
【输入格式】
第一行包含一个整数 N。
第二行包含 N 个整数 A1, A2, · · · AN 。
【输出格式】
输出一个整数代表答案。
【样例输入】
7
1 6 5 4 3 2 1
【样例输出】
2
【评测用例规模与约定】
对于所有评测用例，1 ≤ N ≤ 100000，−100000 ≤ Ai ≤ 100000。

思路：
由于是完全二叉树，所以节点个数与二叉树的深度之间是有关系的。即深度d=log[n]+1 (log[n]向下取整)。 求出了树的深度以后，我们便可以遍历每一个深度下的结点权值和。对于深度是第i层来说，前i-1层的总结点数为2^(i-1)-1;前i层的总结点数为2 ^(i)-1。所以第i层的结点有a[2 ^(i-1)]~a[2 ^i-1]。当然需要注意的是，由于完全二叉树不一定是满二叉树，所以最后一层的最后一个结点序号不一定是a[2 ^d-1]，需要加个限制条件即序号<=n。
时间复杂度：O(n)

#include<iostream>
#include<queue>
#include<vector>
using namespace std;
#define INF 0x7fffffff
const int maxn = 1e5 + 10;
typedef long long ll;
typedef unsigned long long ull;
int get_d(int x) {
	int res = 0;
	while (x/2) {
		res++; x >>= 1;
	}
	return res;
}
int qpow(int base, int power) {
	int res = 1;
	while (power) {
		if (power & 1) res *= base;
		power >>= 1;
		base *= base;
	}
	return res;
}
int  main() {
	ll n,d,ans,maxx,a[maxn]; 
	cin >> n;  
	d=get_d(n)+1;   //d为二叉树的深度
	for (int i = 1; i <= n; i++)  cin >> a[i];
	ans = 1; maxx = a[1];
	for (int i = 2; i <= d; i++) {
		int l = qpow(2, i - 1), r = qpow(2, i);
		ll res = 0;
		for (int j = l; j < r&&j<=n;j++) {
			res += a[j];
		}
		if (res > maxx) maxx = res, ans = i;
	}
	cout << ans << endl;
	return 0;
}

---

G.外卖店优先级(20分）

试题 G: 外卖店优先级
时间限制: 1.0s 内存限制: 256.0MB 本题总分：20 分
【问题描述】
“饱了么”外卖系统中维护着 N 家外卖店，编号 1 ∼ N。每家外卖店都有一个优先级，初始时 (0 时刻) 优先级都为 0。
每经过 1 个时间单位，如果外卖店没有订单，则优先级会减少 1，最低减到 0；而如果外卖店有订单，则优先级不减反加，每有一单优先级加 2。
如果某家外卖店某时刻优先级大于 5，则会被系统加入优先缓存中；如果优先级小于等于 3，则会被清除出优先缓存。
给定 T 时刻以内的 M 条订单信息，请你计算 T 时刻时有多少外卖店在优先缓存中。
【输入格式】
第一行包含 3 个整数 N、M 和 T 。
以下 M 行每行包含两个整数 ts 和 id，表示 ts 时刻编号 id 的外卖店收到一个订单。
【输出格式】
输出一个整数代表答案。
【样例输入】
2 6 6
1 1
5 2
3 1
6 2
2 1
6 2
【样例输出】
1
【样例解释】
6 时刻时，1 号店优先级降到 3，被移除出优先缓存；2 号店优先级升到 6， 加入优先缓存。所以是有 1 家店 (2 号) 在优先缓存中。
【评测用例规模与约定】
对于 80% 的评测用例，1 ≤ N, M, T ≤ 10000。
对于所有评测用例，1 ≤ N, M, T ≤ 100000，1 ≤ ts ≤ T ，1 ≤ id ≤ N。

思路：
中等模拟题。 我们用vis[maxn]数组标记第i家店在T时刻是否在优先缓存中，最后的结果就是vis[]数组中true的个数。 我们可以对每一家店i的所有订单进行分析，将i店的所有订单按照时间先后顺序存储到不定长数组vector< int >v中，然后遍历它的所有订单。 用d表示该i店的优先级，若d>5则vis[i]=true；若d<=3则vis[i]=false;对于该店的第j个订单，若它与j-1个订单的时间差dt>1，即有dt-1个时刻店i无订单，那么d=d-dt+1。为了防止d<0，所以d=max(0,d-dt+1)；然后d+=2,因为每有一个订单优先级+2；处理完店i的所有订单后，再额外判断在T时刻店i的优先级是否在优先缓存中。
时间复杂度：O(max(m,n))

#include<iostream>
#include<queue>
#include<vector>
#include<cstring>
#include<set>
#include<algorithm>
using namespace std;
#define INF 0x7fffffff
const int maxn = 1e5 + 10;
typedef long long ll;
typedef unsigned long long ull;
struct Node {
	int ts, id;
}g[maxn];
int cmp(Node a, Node b) {
	return a.ts < b.ts;
}
vector <int>v[maxn]; //q[i]存储第i家店有订单的时刻
bool vis[maxn];
int  main() {
	iostream::sync_with_stdio(false);
	int n, m, t;
	cin >> n >> m >> t;
	for (int i = 1; i <= n; i++) vis[i] = false;
	for (int i = 1; i <= m; i++) cin >> g[i].ts >> g[i].id;
	sort(g + 1, g + 1 + m,cmp);
	for (int i = 1; i <= m; i++) {
		int ts = g[i].ts, id = g[i].id;
		v[id].push_back(ts);
	}
	for (int i = 1; i <= n; i++) {
		if (v[i].size()==0) continue;
		int pre_t = v[i][0], now_t = v[i][0], d = 2,len=v[i].size();
		for (int j = 1; j < len; j++) {
			now_t = v[i][j];
			int dt = now_t - pre_t;
			if (dt > 1) d = max(0, d - dt + 1);
			if (d <= 3) vis[i] = false;
			d += 2; pre_t = now_t;
			if (d > 5) vis[i] = true;
		}
		if (v[i][len - 1] != t) d = max(0, d - t + v[i][len - 1]);
		if (d <= 3) vis[i] = false;
	}
	int cnt = 0;
	for (int i = 1; i <= n; i++) if (vis[i]) cnt++;
	cout << cnt << endl;
	return 0;
}

---

H.修改数组(20分）

试题 H: 修改数组
时间限制: 1.0s 内存限制: 256.0MB 本题总分：20 分
【问题描述】
给定一个长度为 N 的数组 A = [A1, A2, · · · AN]，数组中有可能有重复出现的整数。
现在小明要按以下方法将其修改为没有重复整数的数组。小明会依次修改A2, A3, · · · , AN。
当修改 Ai 时，小明会检查 Ai 是否在 A1 ∼ Ai−1 中出现过。如果出现过，则小明会给 Ai 加上 1 ；如果新的 Ai 仍在之前出现过，小明会持续给 Ai 加 1 ，直到 Ai 没有在 A1 ∼ Ai−1 中出现过。
当 AN 也经过上述修改之后，显然 A 数组中就没有重复的整数了。现在给定初始的 A 数组，请你计算出最终的 A 数组。
【输入格式】
第一行包含一个整数 N。
第二行包含 N 个整数 A1, A2, · · · , AN 。
【输出格式】
输出 N 个整数，依次是最终的 A1, A2, · · · , AN。
【样例输入】
5
2 1 1 3 4
【样例输出】
2 1 3 4 5
【评测用例规模与约定】
对于 80% 的评测用例，1 ≤ N ≤ 10000。
对于所有评测用例，1 ≤ N ≤ 100000，1 ≤ Ai ≤ 1000000。

思路：
第一种骗分思路： 因为对于80%的样例来说，N<=10000;所以可以直接以时间复杂度O(n^2)暴力便可通过80%的样例。当然由于样例可能太简单了，我用暴力提交的结果是90分。暴力的方法是从i=2开始遍历每一个元素，若该元素未标记，则标记当前元素并且继续遍历下一个元素，否则让当前元素一直加1直至未标记，再将当前元素标记。
第二种AC思路：并查集。 我们可以借助并查集找到从当前元素x开始第一个未出现的元素；即x=find(x)；找到以后，由于元素x要放入序列中，所以元素x不能再次使用，即fa[x]=x+1。

第一种骗分思路代码：（90分）

#include<iostream>
#include<queue>
#include<vector>
#include<cstring>
#include<set>
#include<algorithm>
using namespace std;
#define INF 0x7fffffff
const int maxn = 1e5 + 10;
typedef long long ll;
typedef unsigned long long ull;
int a[maxn], vis[maxn];
int  main() {
	iostream::sync_with_stdio(false);
	int n; cin >> n;
	for (int i = 1; i <= n; i++) {
		vis[i] = 0;
		cin >> a[i];
	}
	vis[a[1]] = 1;
	for (int i = 2; i <= n; i++) {
		int x = a[i];
		while (vis[x]) x++;
		a[i] = x; vis[x] = 1;
	}
	for (int i = 1; i <= n; i++) cout << a[i] << " ";
	return 0;
}

---

第二种AC思路代码：

#include<iostream>
#include<queue>
#include<vector>
#include<cstring>
#include<set>
#include<algorithm>
using namespace std;
#define INF 0x7fffffff
const int maxn = 1e6 + 10;
typedef long long ll;
typedef unsigned long long ull;
int fa[maxn];
inline int find(int x) {
	return x == fa[x] ? x : fa[x] = find(fa[x]);
}
int  vis[maxn];
int  main() {
	iostream::sync_with_stdio(false);
	int n; cin >> n;
	for (int i = 1; i <= maxn; i++) fa[i] = i;
	for (int i = 1; i <= n; i++) {
		int x; cin >> x;
		x = find(x);
		fa[x] = x + 1;
		cout <<x << " ";
	}
	return 0;
}

---

I.糖果(25分）

试题 I: 糖果
时间限制: 1.0s 内存限制: 256.0MB 本题总分：25 分
【问题描述】
糖果店的老板一共有 M 种口味的糖果出售。为了方便描述，我们将 M 种口味编号 1 ∼ M。
小明希望能品尝到所有口味的糖果。遗憾的是老板并不单独出售糖果，而 是 K 颗一包整包出售。
幸好糖果包装上注明了其中 K 颗糖果的口味，所以小明可以在买之前就知道每包内的糖果口味。
给定 N 包糖果，请你计算小明最少买几包，就可以品尝到所有口味的糖果。
【输入格式】
第一行包含三个整数 N、M 和 K。
接下来 N 行每行 K 这整数 T1, T2, · · · , TK，代表一包糖果的口味。
【输出格式】
一个整数表示答案。如果小明无法品尝所有口味，输出 −1。
【样例输入】
6 5 3
1 1 2
1 2 3
1 1 3
2 3 5
5 4 2
5 1 2
【样例输出】
2
【评测用例规模与约定】
对于 30% 的评测用例，1 ≤ N ≤ 20 。
对于所有评测样例，1 ≤ N ≤ 100，1 ≤ M ≤ 20，1 ≤ K ≤ 20，1 ≤ Ti ≤ M。

思路：
状压dp。 由于蒻蒻实力太差，只能看出这是一道状压dp的题目，具体如何定义状态和状态如何转移还是借鉴了大佬的思想。实现状压dp的基础就是学会二进制表示集合，二进制的一些与、或和异或的位运算。 具体内容请读者移驾它文学习一下。 我们用dp[s]表示当前状态为集合s时买的最少糖果包数。 首先初始化dp[s],并且用数组candy[i]存储第i包糖果包含的口味集合； 然后遍历每一包糖果i，枚举每一种状态j是否出现过，若出现过，则用temp表示新状态j|candy[i],该状态表示集合j加上了第i包糖果集合，再处理该状态的值；否则枚举下一个状态；

#include<iostream>
#include<queue>
#include<vector>
#include<cstring>
#include<set>
#include<algorithm>
using namespace std;
#define INF 0x7fffffff
const int maxn = 1e6 + 10;
typedef long long ll;
typedef unsigned long long ull;
int  dp[1<<21];
int  main() {
	iostream::sync_with_stdio(false);
	int n, m, k,candy[105]; cin >> n >> m >> k;
	for (int i = 1; i <= n; i++) {
		candy[i] = 0;
		for (int j = 1; j <= k; j++) {
			int x; cin >> x;
			candy[i] |= (1<<x-1);
		}
	}
	memset(dp, -1, sizeof(dp));
	dp[0] = 0;
	for (int i = 1; i <= n; i++) {
		for (int j = 0; j < (1 << m); j++) {
			if (dp[j] == -1) continue;
			int temp = candy[i] | j;
			if (dp[temp] == -1 || dp[temp] > dp[j] + 1) dp[temp] = dp[j] + 1;
		}
	}
	if (dp[(1 << m) - 1] == -1) cout << "-1" << endl;
	else cout << dp[(1 << m) - 1] << endl;
	return 0;
}

---

J.组合数问题(25分）

试题 J: 组合数问题
时间限制: 1.0s 内存限制: 256.0MB 本题总分：25 分
【问题描述】
给 n, m, k， 求 有 多 少 对 (i, j) 满 足 1 ≤ i ≤ n, 0 ≤ j ≤ min(i, m)且C≡0(mod k)，k 是质数。其中 C j 是组合数，表示从 i 个不同的数中选出 j 个组成一个集合的方案数。
【输入格式】
第一行两个数 t, k，其中 t 代表该测试点包含 t 组询问，k 的意思与上文中相同。
接下来 t 行每行两个整数 n, m，表示一组询问。
【输出格式】
输出 t 行，每行一个整数表示对应的答案。由于答案可能很大，请输出答案除以 10^9 + 7 的余数。
【样例输入】
1 2
3 3
【样例输出】
1
【样例说明】
在所有可能的情况中，只有 C1 = 2 是 2 的倍数。
【样例输入】
2 5
4 5
6 7
【样例输出】
0
7
【样例输入】
3 23
23333333 23333333
233333333 233333333
2333333333 2333333333
【样例输出】
851883128
959557926
680723120
【数据规模和约定】
对于所有评测用例，1 ≤ k ≤ 108, 1 ≤ t ≤ 105, 1 ≤ n, m ≤ 1018，且 k 是质数。评测时将使用 10 个评测用例测试你的程序，每个评测用例的限制如下：
评测用例编号 t n, m k
1, 2 ≤ 1 ≤ 2000 ≤ 100
3, 4 ≤ 105 ≤ 2000 ≤ 100
5, 6, 7 ≤ 100 ≤ 1018 ≤ 100
8, 9, 10 ≤ 105 ≤ 1018 ≤ 108

---