LOJ3120

52pts

\(N - D >= 2M\)或者\(M = 0\)那么就是\(D^{N}\)

只和数字的奇偶性有关,如果有k个奇数,那么必须满足\(N - k >= 2M\)

所以设\(f[i][j]\)表示第\(i\)个数有\(j\)个奇数的方案数,\(j\cdot f[i][j] \rightarrow f[i + 1][j - 1]\)和\((D - j) \cdot f[i][j] \rightarrow f[i + 1][j + 1]\)

64pts

这个只需要把上面的矩阵快速幂优化,只不过需要一点小技巧来卡一卡常……具体看代码吧

+8pts?

如果m = 1,不合法的只有每个数都不一样的情况

m = 2除去每个数都不一样的情况,有一个数大于1的情况且这个多少为2或3

以上就有72pts了!快落!

剩下的我就不会了,然后翻翻网上的题解写了一个容斥做法

至少k个方案数

奇数的序列的指数生成函数是

\(\frac{e^{x} - e^{-x}}{2} = \frac{x^{1}}{1!} + \frac{x^{3}}{3!} + \frac{x^{5}}{5!}....\)

至少有\(i\)个数为奇数的方案数

\(f_{i} = \binom{D}{i}n![x^{n}](\frac{e^{x} - e^{-x}}{2})^{i}e^{(D - i)x}\)

把\(2^{-i}\)提出去

\(f_{i} = \frac{\binom{D}{i}}{2^{i}} n!(e^{x} - e^{-x})^{i}e^{(D - i)x}\)

方案数是

然后把那个二项式展开一下

\(f_{i} = \frac{\binom{D}{i}}{2^{i}}n!e^{(D - i)x}\sum_{j = 0}^{i} (-1)^{j}e^{-jx}e^{(i - j)x}\binom{i}{j}[x^{n}]\)

合并一下就是

\(f_{i} = \frac{\binom{D}{i}}{2^{i}}n!\sum_{j = 0}^{i} (-1)^{j}e^{(D - 2j)x}\binom{i}{j}[x^{n}]\)

由于\(e^{ax}\)的第n项生成函数是\(\frac{a^{n}}{n!}\)

所以最后就是

\(f_{i} = \frac{i!\binom{D}{i}}{2^{i}}\sum_{j = 0}^{i} (-1)^{j}\frac{(D - 2j)^n}{j!(i - j)!}[x^{n}]\)

可以卷积算出来

然后

\(g_{i} = \sum_{j = i}^{D} (-1)^{j - i}\binom{i}{j}f_{j}\)

这个卷积一下也可以算

代码

#include <bits/stdc++.h>
#define fi first
#define se second
#define pii pair<int, int>
#define mp make_pair
#define pb push_back
#define space putchar(' ')
#define enter putchar('\n')
#define eps 1e-10
#define MAXN 2005
#define ba 47
//#define ivorysi
using namespace std;
typedef long long int64;
typedef unsigned int u32;
typedef double db;
template <class T>
void read(T &res) {
res = 0;
T f = 1;
char c = getchar();
while (c < '0' || c > '9') {
if (c == '-')
f = -1;
c = getchar();
}
while (c >= '0' && c <= '9') {
res = res * 10 + c - '0';
c = getchar();
}
res *= f;
}
template <class T>
void out(T x) {
if (x < 0) {
x = -x;
putchar('-');
}
if (x >= 10) {
out(x / 10);
}
putchar('0' + x % 10);
}
const int MOD = 998244353;
int D, N, M;
int fac[1000005], invfac[1000005];
int inc(int a, int b) { return a + b >= MOD ? a + b - MOD : a + b; }
int mul(int a, int b) { return 1LL * a * b % MOD; }
int C(int n, int m) {
if (n < m)
return 0;
return mul(fac[n], mul(invfac[m], invfac[n - m]));
}
void update(int &x, int y) { x = inc(x, y); }
int fpow(int x, int c) {
int res = 1, t = x;
while (c) {
if (c & 1)
res = mul(res, t);
t = mul(t, t);
c >>= 1;
}
return res;
}
namespace task1 {
int f[2][4005];
void Main() {
int cur = 0;
f[cur][0] = 1;
for (int i = 1; i <= N; ++i) {
memset(f[cur ^ 1], 0, sizeof(f[cur ^ 1]));
for (int j = 0; j <= D; ++j) {
if (j >= 1)
update(f[cur ^ 1][j - 1], mul(f[cur][j], j));
update(f[cur ^ 1][j + 1], mul(f[cur][j], D - j));
}
cur ^= 1;
}
int ans = 0;
for (int i = 0; i <= D; ++i) {
if (N - i >= 2 * M)
update(ans, f[cur][i]);
}
out(ans);
enter;
}
} // namespace task1
namespace task2 {
vector<int> v1[305], v2[305];
struct Matrix {
int f[305][305];
Matrix() { memset(f, 0, sizeof(f)); }
friend Matrix operator*(const Matrix &a, const Matrix &b) {
Matrix c;
for (int i = 0; i <= D; ++i) {
v1[i].clear();
v2[i].clear();
}
for (int i = 0; i <= D; ++i) {
for (int j = 0; j <= D; ++j) {
if (a.f[i][j])
v1[i].pb(j);
if (b.f[i][j])
v2[i].pb(j);
}
}
for (int i = 0; i <= D; ++i) {
for (auto k : v1[i]) {
for (auto j : v2[k]) {
update(c.f[i][j], mul(a.f[i][k], b.f[k][j]));
}
}
}
return c;
}
} a, ans;
Matrix fpow(Matrix x, int c) {
Matrix res = x, t = x;
--c;
while (c) {
if (c & 1)
res = res * t;
t = t * t;
c >>= 1;
}
return res;
}
void Main() {
for (int i = 0; i <= D; ++i) {
if (i >= 1)
update(a.f[i][i - 1], i);
update(a.f[i][i + 1], D - i);
}
ans = fpow(a, N);
int res = 0;
for (int i = 0; i <= D; ++i) {
if (N - i >= 2 * M) {
update(res, ans.f[0][i]);
}
}
out(res);
enter;
}
} // namespace task2
namespace task3 {
const int MAXL = (1 << 20);
int W[MAXL + 5];
vector<int> a, b, f, g;
void NTT(vector<int> &p, int L, int on) {
p.resize(L);
for (int i = 1, j = L >> 1; i < L - 1; ++i) {
if (i < j)
swap(p[i], p[j]);
int k = L >> 1;
while (j >= k) {
j -= k;
k >>= 1;
}
j += k;
}
for (int h = 2; h <= L; h <<= 1) {
int wn = W[(MAXL + on * MAXL / h) % MAXL];
for (int k = 0; k < L; k += h) {
int w = 1;
for (int j = k; j < k + h / 2; ++j) {
int u = p[j], t = mul(p[j + h / 2], w);
p[j] = inc(u, t);
p[j + h / 2] = inc(u, MOD - t);
w = mul(w, wn);
}
}
}
if (on == -1) {
int invL = fpow(L, MOD - 2);
for (int i = 0; i < L; ++i) {
p[i] = mul(p[i], invL);
}
}
}
vector<int> operator*(vector<int> a, vector<int> b) {
vector<int> c;
int t = a.size() + b.size() - 2, l = 1;
while (l <= t) l <<= 1;
c.resize(l);
NTT(a, l, 1);
NTT(b, l, 1);
for (int i = 0; i < l; ++i) c[i] = mul(a[i], b[i]);
NTT(c, l, -1);
return c;
}
void Init() {
W[0] = 1;
W[1] = fpow(3, (MOD - 1) / MAXL);
for (int i = 2; i < MAXL; ++i) W[i] = mul(W[i - 1], W[1]);
}
void Main() {
Init();
a.resize(D + 1);
b.resize(D + 1);
int t = 1;
for (int i = 0; i <= D; ++i) {
a[i] = mul(mul(t, fpow(inc(D, MOD - 2 * i), N)), invfac[i]);
b[i] = invfac[i];
t = mul(t, MOD - 1);
}
f = a * b;
f.resize(D + 1);
t = 1;
for (int i = 0; i <= D; ++i) {
f[i] = mul(f[i], fac[D]);
f[i] = mul(f[i], invfac[D - i]);
f[i] = mul(f[i], t);
t = mul(t, (MOD + 1) / 2);
}
a.clear();
a.resize(D + 1);
t = 1;
for (int i = 0; i <= D; ++i) {
a[D - i] = mul(t, invfac[i]);
f[i] = mul(f[i], fac[i]);
t = mul(t, MOD - 1);
}
g = a * f;
int L = N - 2 * M;
int ans = 0;
for (int i = 0; i <= L; ++i) {
update(ans, mul(g[i + D], invfac[i]));
}
out(ans);
enter;
}
} // namespace task3
int main() {
#ifdef ivorysi
freopen("f1.in", "r", stdin);
#endif
fac[0] = 1;
for (int i = 1; i <= 1000000; ++i) fac[i] = mul(fac[i - 1], i);
invfac[1000000] = fpow(fac[1000000], MOD - 2);
for (int i = 999999; i >= 0; --i) invfac[i] = mul(invfac[i + 1], i + 1);
read(D);
read(N);
read(M);
if (N - D >= 2 * M || M == 0) {
out(fpow(D, N));
enter;
} else if (N <= 4000)
task1::Main();
else if (D <= 300)
task2::Main();
else if (M == 1) {
int ans = fpow(D, N);
if (D >= N) {
update(ans, MOD - mul(fac[D], invfac[D - N]));
}
out(ans);
enter;
} else if (M == 2) {
int ans = fpow(D, N);
if (D >= N) {
update(ans, MOD - mul(fac[D], invfac[D - N]));
}
for (int j = 2; j <= 3; ++j) {
if (N - j <= D - 1) {
int t = mul(D, C(D - 1, N - j));
update(ans, MOD - mul(mul(t, invfac[j]), fac[N]));
}
}
out(ans);
enter;
} else {
task3::Main();
}
return 0;
}

【LOJ】#3120. 「CTS2019 | CTSC2019」珍珠的更多相关文章

  1. @loj - 3120@ 「CTS2019 | CTSC2019」珍珠

    目录 @description@ @solution@ @accepted code@ @details@ @description@ 有 \(n\) 个在范围 \([1, D]\) 内的整数均匀随机 ...

  2. Loj #3124. 「CTS2019 | CTSC2019」氪金手游

    Loj #3124. 「CTS2019 | CTSC2019」氪金手游 题目描述 小刘同学是一个喜欢氪金手游的男孩子. 他最近迷上了一个新游戏,游戏的内容就是不断地抽卡.现在已知: - 卡池里总共有 ...

  3. LOJ 3120: 洛谷 P5401: 「CTS2019 | CTSC2019」珍珠

    题目传送门:LOJ #3120. 题意简述: 称一个长度为 \(n\),元素取值为 \([1,D]\) 的整数序列是合法的,当且仅当其中能够选出至少 \(m\) 对相同元素(不能重复选出元素). 问合 ...

  4. LOJ 3124 「CTS2019 | CTSC2019」氪金手游——概率+树形DP

    题目:https://loj.ac/problem/3124 看了题解:https://www.cnblogs.com/Itst/p/10883880.html 先考虑外向树. 考虑分母是 \( \s ...

  5. loj3120 「CTS2019 | CTSC2019」珍珠

    link .... 感觉自己太颓废了....还是来更题解吧...[话说写博客会不会涨 rp 啊 qaq ? 题意: 有 n 个物品,每个都有一个 [1,D] 中随机的颜色,相同颜色的两个物品可以配对. ...

  6. LOJ #3119「CTS2019 | CTSC2019」随机立方体 (容斥)

    博客链接 里面有个下降幂应该是上升幂 还有个bk的式子省略了k^3 CODE 蛮短的 #include <bits/stdc++.h> using namespace std; const ...

  7. LOJ #3119. 「CTS2019 | CTSC2019」随机立方体 组合计数+二项式反演

    好神的一道计数题呀. code: #include <cstdio> #include <algorithm> #include <cstring> #define ...

  8. 「CTS2019 | CTSC2019」氪金手游 解题报告

    「CTS2019 | CTSC2019」氪金手游 降 智 好 题 ... 考场上签到失败了,没想容斥就只打了20分暴力... 考虑一个事情,你抽中一个度为0的点,相当于把这个点删掉了(当然你也只能抽中 ...

  9. 「CTS2019 | CTSC2019」随机立方体 解题报告

    「CTS2019 | CTSC2019」随机立方体 据说这是签到题,但是我计数学的实在有点差,这里认真说一说. 我们先考虑一些事实 如果我们在位置\((x_0,y_0,z_0)\)钦定了一个极大数\( ...

随机推荐

  1. Cogs 746. [网络流24题] 骑士共存(最大独立集)

    [网络流24题] 骑士共存 ★★☆ 输入文件:knight.in 输出文件:knight.out 简单对比 时间限制:1 s 内存限制:128 MB 骑士共存问题 «问题描述: 在一个n*n个方格的国 ...

  2. docker部署项目: centos+python+redis+mysql+uwsgi+nginx

    一.Centos7安装docker 1.1 环境配置 先测试是否下载了docker:查看镜像:docker images没有下载,就依次执行以下环境的安装 ①curl http://mirrors.a ...

  3. T2695 桶哥的问题——吃桶 题解

    校内测试 ------T3 对于这个题,首先想到的应该就是暴力枚举了吧,看看数据范围,60就是白送的啦!(但是我也不知道怎么才20分qwq) 思路分析: 这个题要你求所有套餐的总价值,先看一眼产生套餐 ...

  4. SpringMVC——文件上传下载

    一.单文件上传 1.导入依赖 <dependency> <groupId>commons-io</groupId> <artifactId>common ...

  5. Java内置锁synchronized的实现原理及应用(三)

    简述 Java中每个对象都可以用来实现一个同步的锁,这些锁被称为内置锁(Intrinsic Lock)或监视器锁(Monitor Lock). 具体表现形式如下: 1.普通同步方法,锁的是当前实例对象 ...

  6. JavaWeb_(Hibernate框架)Hibernate论坛项目中多对多案例

    基于SSH论坛小型项目 传送门 用户和发帖进行举例 多对多关系:多个用户可以回复多个帖子 因此引入了一张回复表,用来保存用户id和帖子id CREATE TABLE `hforum`.`answer` ...

  7. reduce()之js与python

    最近在自学python基础,发现很多js中的内置函数python中都存在,甚至比js更加简洁,下面说一下reduce()在js和python中的用法,做个对比便于记忆. reduce()简介: red ...

  8. 如何简单的在linux上安装jdk并配置环境变量

    这篇文章是为了给我一会自己安装的时候方便使用的,所以内容很简单,平时在wendows系统上安装很容易,但是换到linux系统上面就蒙圈了. 一.下载jdk文件 我这提供的是官方的地址:http://w ...

  9. HashMap简单介绍

    哈希表(hash table)也叫散列表,是一种非常重要的数据结构,应用场景及其丰富,许多缓存技术(比如memcached)的核心其实就是在内存中维护一张大的哈希表. 一.什么是哈希表 在讨论哈希表之 ...

  10. oracle口令文件在windows和linux系统下的命名和位置

      分类: Oracle 1.windows系统下 oracle口令文件在:$ORACLE_HOME/database目录下: 命名规则为:PWD$SID.ora   2.linux系统下 oracl ...