题目链接
没做现场题,按着自己思路码的,有错请纠正。
第4题无思路
1、生成数列的个数(dp)
思路:定义
d
p
[
i
]
[
j
]
dp[i][j]
dp[i][j]表示取i个数,模3结果为j的组合数。那么
d
p
[
n
]
[
0
]
dp[n][0]
dp[n][0]即为所求。
#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int maxn = 100010;
const int mod = 1000000007;
int dp[maxn][3];
// dp
// dp[i][j] = sum(dp[i-1][k] * num[3 - k], 0 <= k < 3)
// num[i]表示区间[l,r]中%3为i的个数
int Mul(int a, int b) {
return 1LL * a * b % mod;
}
int Add(int a, int b) {
return (1LL * a + b) % mod;
}
int main() {
int n, l, r;
scanf("%d%d%d", &n, &l, &r);
memset(dp, 0, sizeof(dp));
int num[3];// 区间[l,r]中%3为i的个数num[i]
--l;
for (int i = 0;i < 3; ++i) {
int numr = r / 3 + ((r % 3) > i);
int numl = l / 3 + ((l % 3) > i);
num[i] = numr - numl;
}
// printf("num {%d %d %d}\n", num[0], num[1], num[2]);
int len = r - l;
dp[0][0] = 1;
for (int i = 1;i <= n; ++i) {
// dp[i][0] = dp[i-1][0] * num[0] + dp[i-1][1] * num[2] + dp[i-1][2] * num[1];
// dp[i][1] = dp[i-1][0] * num[1] + dp[i-1][1] * num[0] + dp[i-1][2] * num[2];
// dp[i][2] = dp[i-1][0] * num[2] + dp[i-1][1] * num[1] + dp[i-1][2] * num[0];
dp[i][0] = Add(Add(Mul(dp[i-1][0], num[0]), Mul(dp[i-1][1], num[2])), Mul(dp[i-1][2], num[1]));
dp[i][1] = Add(Add(Mul(dp[i-1][0], num[1]), Mul(dp[i-1][1], num[0])), Mul(dp[i-1][2], num[2]));
dp[i][2] = Add(Add(Mul(dp[i-1][0], num[2]), Mul(dp[i-1][1], num[1])), Mul(dp[i-1][2], num[0]));
// printf("dp[%d] {%d %d %d}\n", i, dp[i][0], dp[i][1], dp[i][2]);
}
printf("%d\n", dp[n][0]);
}
/*
2 1 3
*/
2、统计不能被k整除的子集个数(组合数学)
思路:两个数之和能被k整除,说明
x
+
y
=
=
0
,
m
o
d
k
x + y == 0 ,\mod k
x+y==0,modk,等价于
x
%
k
+
y
%
k
=
=
0
,
m
o
d
k
x\%k + y\%k == 0,\mod k
x%k+y%k==0,modk。那么我们直接把每个数对k取模,用数组
n
u
m
num
num统计模k为i的数量
n
u
m
[
i
]
num[i]
num[i]。对于数i,我们只要不取模k为k-i的数,就不会和相加,结果为0(模k)。
计算总数量,对一个互斥集合,我们要么取
2
n
u
m
i
2^{num_i}
2numi,要么取
2
n
u
m
k
−
i
2^{num_{k-i}}
2numk−i,要么取空集,取法有
2
n
u
m
i
+
2
n
u
m
k
−
i
−
1
2^{num_i}+2^{num_{k-i}}-1
2numi+2numk−i−1。我们把他们累乘即可。
答案为
∏
i
=
0
k
/
2
−
1
2
n
u
m
i
+
2
n
u
m
k
−
i
−
1
\prod_{i = 0}^{k/2-1}2^{num_i}+2^{num_{k-i}}-1
∏i=0k/2−12numi+2numk−i−1。如果
k
%
2
=
=
0
k\%2 == 0
k%2==0,需要特别计算。
#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int maxn = 1010;
const int mod = 1000000007;
//const int NUM = 1000;
int num[maxn];
int p[maxn];
// init 预处理 2^i
void init() {
p[0] = 1;
for (int i = 1;i < maxn; ++i) {
p[i] = p[i-1] * 2LL % mod;
}
}
int Mul(int a, int b) {
return 1LL * a * b % mod;
}
int Add(int a, int b) {
return (1LL * a + b) % mod;
}
int Sub(int a, int b) {
return (1LL * a + mod - b) % mod;
}
int main() {
init();
int n, k;
scanf("%d%d", &n, &k);
memset(num, 0 ,sizeof(num));
for (int i = 1, x;i <= n; ++i) {
scanf("%d", &x);
++num[x%k];
}
int ans = 1;
ans = Mul(ans, num[0] + 1);
for (int i = 1;i < (k + 1) / 2; ++i) {
// ans = ans * (p[num[i]] + p[num[k-i]] - 1);
ans = Mul(ans, Sub(Add(p[num[i]], p[num[k-i]]), 1));
}
if (k % 2 == 0) {
// ans = ans * p[num[k/2]];
// ans = Mul(ans, p[num[k/2]]);
ans = Mul(ans, num[k/2] + 1);
}
printf("%d\n",ans);
}
/*
4 3
1 2 3 4
*/
3、启动机器(贪心)
贪心模拟,每次尽力晚的启动机器
#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int maxn = 1010;
const int mod = 1000000007;
int n, m, t;
int main() {
scanf("%d%d%d%d", &n, &t, &m);
int pre, w;
bool flag = 1;
if (t < m) flag = 0;
scanf("%d", &pre);
int ans = m;
queue<int> q;
for (int i = pre - m;i < pre; ++i)// 贪心,尽可能放前面
q.push(i);
for (int i = 1;i < n; ++i) {
scanf("%d", &w);
if(!flag) continue;
int bound = w - t +1;// 可用边界
while(!q.empty() && q.top() < bound) q.pop();
if(q.size() < m) {
int len = m - q.size();
ans += len;
for(int j = w - len; j < w; ++j)// 贪心,尽可能放前面
q.push(j);
}
}
if(!flag) printf("-1\n");
else printf("%d\n",ans);
}
4、 线段和树
n比较大,怎么O(n)实现,暂无思路。