传送门

Sol

对于每个 \(i\) ,可以把 \(k\) 个数字分成 \((x,i-x)\) 的若干组。

那么就是求每组只能其中选择一个且可以重复的方案数。

预处理 \(f[i][j]\) 表示从 \(j\) 个组内选 \(i\) 个,每个组必须选的方案数。

\(f[i][j]=(f[i-1][j]+f[i-1][j-1]\times 2)\) 注意组有序,而选出来是无序的,所以要强制选择最后一组

那么每个组可以不选的方案数 \(g[i][j]\) 就是

\[g[i][j]=\sum_{t=0}^{j}(^{j}_{t})f[i][t]
\]

拆开组合数可以 \(NTT\) 优化。

对于询问分两种情况

  1. \(i\) 为奇数

    那么 \(x\ne i-x\) 直接枚举多少个选组内的,其它的可重组合即可
  2. \(i\) 为偶数

    存在一组 \(x=i-x\) 枚举是否选 \(x\) (因为只能选一个),然后像奇数一样做即可

    \(\Theta(n^2logn)\) 丝毫不虚
# include <bits/stdc++.h>
using namespace std;
typedef long long ll; namespace IO {
const int maxn((1 << 21) + 1); char ibuf[maxn], *iS, *iT, c;
int f; inline char Getc() {
return iS == iT ? (iT = (iS = ibuf) + fread(ibuf, 1, maxn, stdin), iS == iT ? EOF : *iS++) : *iS++;
} template <class Int> inline void In(Int &x) {
for (f = 1, c = Getc(); c < '0' || c > '9'; c = Getc()) f = c == '-' ? -1 : 1;
for (x = 0; c >= '0' && c <= '9'; c = Getc()) x = (x << 1) + (x << 3) + (c ^ 48);
x *= f;
}
} using IO :: In; const int maxn(4005);
const int mod(998244353); int f[maxn][maxn], c[maxn][maxn], fac[maxn], ifac[maxn]; inline void Inc(int &x, int y) {
x += y;
if (x >= mod) x -= mod;
} inline int Pow(ll x, int y) {
ll ret = 1;
for (; y; y >>= 1, x = x * x % mod)
if (y & 1) ret = ret * x % mod;
return ret;
} inline int Calc(int n, int num, int res) {
int ret = 0;
if (!res) ret = f[n][num];
else {
for (int j = 0; j <= n; ++j)
Inc(ret, 1LL * f[n - j][num] * c[j + res - 1][j] % mod);
}
return ret;
} int a[maxn << 2], b[maxn << 2], len, l, r[maxn << 2], w[2][maxn << 2]; inline void NTT(int *p, int opt) {
for (int i = 0; i < len; ++i) if (i < r[i]) swap(p[i], p[r[i]]);
for (int i = 1; i < len; i <<= 1)
for (int j = 0, t = i << 1; j < len; j += t)
for (int k = 0; k < i; ++k) {
int nw = w[opt == -1 ? 1 : 0][len / i * k];
int x = p[j + k], y = 1LL * p[j + k + i] * nw % mod;
p[j + k] = (x + y) % mod, p[j + k + i] = (x - y + mod) % mod;
}
if (opt == -1) {
int inv = Pow(len, mod - 2);
for (int i = 0; i < len; ++i) p[i] = 1LL * p[i] * inv % mod;
}
} inline void Init() {
f[0][0] = c[0][0] = ifac[0] = fac[0] = 1;
for (int i = 1, j; i <= 4000; ++i)
for (c[i][0] = j = 1; j <= i; ++j)
c[i][j] = (c[i - 1][j] + c[i - 1][j - 1]) % mod;
for (int i = 1; i <= 2000; ++i)
for (int j = 1; j <= 2000; ++j)
f[i][j] = (f[i - 1][j] + 1LL * f[i - 1][j - 1] * 2 % mod) % mod;
for (int i = 1; i <= 2000; ++i) fac[i] = 1LL * fac[i - 1] * i % mod;
ifac[2000] = Pow(fac[2000], mod - 2);
for (int i = 1999; i; --i) ifac[i] = 1LL * ifac[i + 1] * (i + 1) % mod;
for (len = 1; len <= 4000; len <<= 1) ++l;
for (int i = 0; i < len; ++i) r[i] = (r[i >> 1] >> 1) | ((i & 1) << (l - 1));
for (int i = 1; i < len; i <<= 1) {
int wn = Pow(3, (mod - 1) / (i << 1));
for (int nw = 1, k = 0; k < i; ++k, nw = 1LL * nw * wn % mod)
w[0][len / i * k] = nw, w[1][len / i * k] = Pow(nw, mod - 2);
}
for (int i = 0; i <= 2000; ++i) {
for (int j = 0; j < len; ++j) a[j] = b[j] = 0;
for (int j = 0; j <= 2000; ++j) a[j] = ifac[j], b[j] = 1LL * ifac[j] * f[i][j] % mod;
NTT(a, 1), NTT(b, 1);
for (int j = 0; j < len; ++j) a[j] = 1LL * a[j] * b[j] % mod;
NTT(a, -1);
for (int j = 0; j <= 2000; ++j) f[i][j] = 1LL * a[j] * fac[j] % mod;
}
} int n, k; int main() {
In(k), In(n), Init();
int r = k << 1;
for (int i = 2; i <= r; ++i) {
if (i & 1) {
int num = max(0, min(i >> 1, i - 1) - max(1, i - k) + 1);
int res = k - (num << 1);
printf("%d\n", Calc(n, num, res));
}
else {
int num = max(0, min(i >> 1, i - 1) - max(1, i - k) + 1);
int res = k - (num << 1) + 1;
if (num) --num;
printf("%d\n", (Calc(n, num, res) + Calc(n - 1, num, res)) % mod);
}
}
return 0;
}

ARG102E:Stop. Otherwise...的更多相关文章

  1. java web 开发三剑客 -------电子书

    Internet,人们通常称为因特网,是当今世界上覆盖面最大和应用最广泛的网络.根据英语构词法,Internet是Inter + net,Inter-作为前缀在英语中表示“在一起,交互”,由此可知In ...

  2. 所有selenium相关的库

    通过爬虫 获取 官方文档库 如果想获取 相应的库 修改对应配置即可 代码如下 from urllib.parse import urljoin import requests from lxml im ...

  3. In-Memory:内存数据库

    在逝去的2016后半年,由于项目需要支持数据的快速更新和多用户的高并发负载,我试水SQL Server 2016的In-Memory OLTP,创建内存数据库实现项目的负载需求,现在项目接近尾声,系统 ...

  4. 从直播编程到直播教育:LiveEdu.tv开启多元化的在线学习直播时代

    2015年9月,一个叫Livecoding.tv的网站在互联网上引起了编程界的注意.缘于Pingwest品玩的一位编辑在上网时无意中发现了这个网站,并写了一篇文章<一个比直播睡觉更奇怪的网站:直 ...

  5. 【.net 深呼吸】细说CodeDom(8):分支与循环

    有人会问,为啥 CodeDom 不会生成 switch 语句,为啥没生成 while 语句之类.要注意,CodeDom只关心代码逻辑,而不是语法,语法是给写代码的人用的.如果用.net的“反编译”工具 ...

  6. 【原】Android热更新开源项目Tinker源码解析系列之三:so热更新

    本系列将从以下三个方面对Tinker进行源码解析: Android热更新开源项目Tinker源码解析系列之一:Dex热更新 Android热更新开源项目Tinker源码解析系列之二:资源文件热更新 A ...

  7. App开发:模拟服务器数据接口 - MockApi

    为了方便app开发过程中,不受服务器接口的限制,便于客户端功能的快速测试,可以在客户端实现一个模拟服务器数据接口的MockApi模块.本篇文章就尝试为使用gradle的android项目设计实现Moc ...

  8. TODO:macOS编译PHP7.1

    TODO:macOS编译PHP7.1 本文主要介绍在macOS上编译PHP7.1,有兴趣的朋友可以去尝试一下. 1.下载PHP7.1源码,建议到PHP官网下载纯净到源码包php-7.1.0.tar.g ...

  9. In-Memory:在内存中创建临时表和表变量

    在Disk-Base数据库中,由于临时表和表变量的数据存储在tempdb中,如果系统频繁地创建和更新临时表和表变量,大量的IO操作集中在tempdb中,tempdb很可能成为系统性能的瓶颈.在SQL ...

随机推荐

  1. 从源代码解释Android事件分发机制

    在ViewRootImpl的setView方法中.用户的触摸按键消息是体如今窗体上的.而windowManagerService则是管理这些窗体,它一旦接收到用户对窗体的一些触摸按键消息,会进行对应的 ...

  2. 【css】gradient匹配ps渐变叠加效果

    CSS3 Gradient分为linear-gradient(线性渐变)和radial-gradient(径向渐变). 一.linear-gradient(线性渐变) eg:background: l ...

  3. Postman无法正常启动解决办法

    问题描述: 应用程序窗口能够打开,但就是这样一直空白,什么都不显示.接下来,主窗口以纯白色加载,不显示任何其他内容. 接下来主窗口背景米色加载和菜单栏加载和工作.应用程序将永远保持这样, 有时界面会变 ...

  4. [CF1019C]Sergey's problem[构造]

    题意 找出一个集合 \(Q\),使得其中的点两两之间没有连边,且集合中的点可以走不超过两步到达其他所有不在集合中的点.输出任意一组解. \(n\leq 10^6\) 分析 考虑构造,先从 \(1\) ...

  5. Visual Studio Package 插件开发(Visual Studio SDK)

    背景 这段时间公司新做了一个支付系统,里面有N个后台服务,每次有更新修改,拷贝打包发布包“不亦乐乎”...于是我想要不要自己定制个打包插件. 部分朋友可能会认为,有现成的可以去找一个,干嘛不用持续集成 ...

  6. git 创建标签和删除标签

    创建标签 在Git中打标签非常简单,首先,切换到需要打标签的分支上: $ git branch * dev master $ git checkout master Switched to branc ...

  7. UE4官方行为树快速入门文档解析和修改

    近学习了UE4官方文档的行为树快速入门指南,发现里面的部分逻辑稍稍有点混乱和重叠,于是加入了自己的想法,修改了部分行为树逻辑,优化了其AI寻路能力. 初始的基本操作和资源创建同官方文档一样:1个Fol ...

  8. PHP Laravel 连接并访问数据库

    第一次连接数据库 数据库配置位于config/database.php数据库用户名及密码等敏感信息位于.env文件创建一个测试表laravel_course <?php namespace Ap ...

  9. webpack4-用之初体验,一起敲它十一遍

    众所周知,webpack进入第4个大版本已经有2个月的时间了,而且webpack团队升级更新的速度也是非常的惊人 在写下如下内容的时候webpack已经出到了4.6的版本了,剑指5.0应该是指日可待了 ...

  10. PAT甲题题解-1037. Magic Coupon (25)-贪心,水

    题目说了那么多,就是给你两个序列,分别选取元素进行一对一相乘,求得到的最大乘积. 将两个序列的正和负数分开,排个序,然后分别将正1和正2前面的相乘,负1和负2前面的相乘,累加和即可. #include ...